Numerical weather modeling can benefit from improved wind sensing in the Earth's atmospheric boundary layer (ABL). Small, low-cost, uncrewed aircraft (drones) can be used to measure wind and a distribution of these vehicles could potentially provide measurements with much greater density and resolution, in both space and time, than current methods allow. To measure wind, a drone could be equipped with dedicated wind-measuring sensors, although these can be costly and obtrusive and must be carefully calibrated to account for interference effects. State estimation algorithms that combine a drone's operational measurements with a flight dynamic model can be used to infer wind without a dedicated wind sensor, although the sensor quality affects measurement accuracy. Previous studies have explored the effects of various sensors on wind estimate accuracy, but the effect of flight dynamic model fidelity has received less attention. This dissertation presents analysis of different aerodynamic model-free and model-based wind estimation methods, comparing six wind estimation formulations using experimental flight data from a small, fixed-wing aircraft. Each formulation is implemented using a Kalman filter, an extended Kalman filter, and an unscented Kalman filter. These filters are designed based on different assumptions related to the flight dynamic model, available sensors, and available measurements. Having identified a promising estimation approach, the dissertation next explores the value of incorporating unsteady effects into a flight dynamic model for model-based wind estimation. An unsteady aerodynamic model for a small, fixed-wing aircraft is developed, identified, and validated using experimental flight data. An extended Kalman filter is then designed and implemented for two motion models -- one that includes unsteady effects and another that does not. Analysis of the wind estimates and the estimation differences show that, while the unsteady flight dynamic model better predicts the aircraft motion, the value of incorporating this model for wind estimation is questionable. / Doctor of Philosophy / Wind velocity sensing is crucial to understanding the meteorological processes at low altitudes. The integration of low-cost drones has allowed them to be used as wind-sensing platforms. This is achieved by equipping small drones with dedicated wind-measuring sensors, often costly and infeasible, or inferring wind velocity from the drone's motion. Algorithms designed to infer wind can be used by combining onboard flight sensor measurements with a drone's flight dynamic model to infer wind. However, low-cost drones are usually equipped with low-cost flight sensors, which frequently lead to higher measurement uncertainty and degrade the accuracy of wind estimates. Previous studies have explored the effects of various sensors on wind estimates, but errors due to low-fidelity dynamic models have received less attention. This dissertation first presents a detailed analysis of different flight dynamic model-free and model-based wind estimation methods. It compares six wind estimation formulations. Each formulation is implemented in wind inferring algorithms called a Kalman filter, an extended Kalman filter, and an unscented Kalman filter. These algorithms are designed based on different assumptions related to the flight dynamic model, available flight sensors, and available measurements. Secondly, the value of incorporating a fixed-wing, unsteady flight dynamic model in a wind estimation scheme is analyzed. To this end, an unsteady flight dynamic model for a fixed-wing drone is developed, identified, and validated from data acquired from the drone's flight history. Furthermore, an extended Kalman filter is designed and implemented for two motion models -- one that includes unsteady effects and another that does not. The analysis of the time histories of the wind estimates and the wind estimate differences show that both model-based estimators perform equally well.

Yes / A problem with 1-D water hammer modelling is in the application of accurate unsteady friction. Moreover, investigating the time response of fluid dynamics and unsteady turbulence structures during the water hammer is not possible with a 1-D model. This review article provides a summary of 1-D modelling using the recent finite volume approach and the discussion extends to a quasi-2-D model and historical developments as well as recent advancements in 3-D CFD simulations of water hammer. The eddy viscosity model is excellent in capturing pressure profiles but it is computationally intensive and requires more computational time. This article reviews 3-D CFD simulations with sliding mesh, an immersed solid approach, and dynamic mesh approaches for modelling valve closures. Despite prediction accuracy, a huge computational time and high computer resources are required to execute 3-D flow simulations with advanced valve modelling techniques. Experimental validation shows that a 3-D CFD simulation with a flow rate reduction curve as a boundary condition predicted accurate pressure variation results. Finally, a brief overview of the transient flow turbulence structures for a rapidly accelerated and decelerated pipe flow using DNS (Direct numerical simulation) data sets is presented. Overall, this paper summarises past developments and future scope in the field of water hammer modelling using CFD.

Thesis (MScEng (Department of Mechanical and Mechatronic Engineering))--Stellenbosch University, 2006. / ENGLISH ABSTRACT: Wind conditions in South Africa determine the need for a small-scale wind turbine to produce useable power at windspeeds below 7m/s. In this project, a range of windspeeds, within which optimal performance o the wind turbine is expected, was selected. The optimal performance was assessed in terms of the Coefficient of Power(Cp), which rates the turbines blade's ability to extract energy form the avalible wind stream. The optimisation methods employed allowed a means of tackling the multi-variable problem such that the aerodynamic characteristics of the blade were ideal throughout the wind speed range. The design problem was broken down into a two-dimensional optimisaion of the airfoils used at the radial stations, and a three-dimensional optimisation of the geometric features of the wind rotor. by means of blending various standard airfoil profiles, a new profile was created at each radial station. XFOIL was used for the two-dimensional analysis of these airfoils. Three-dimensional optimisn involved representation of the rotor as a simplified model and use of the Blade Element Momentum(BEM) method for analysis. an existimg turbine blade, on which the design specifications were modelled, was further used for comparative purposes throughout the project. The resulting blade design offers substantial improvements on the reference design. The application of optimisation methods has successfully aided the creation of a wind turbine blade with consistent peak performance over a range of design prints. / Sponsored by the Centre for Renewable and Sustainable Energy Studies, Stellenbosch University

Thesis (PhD) -- Stellenbosch University, 2014. / ENGLISH ABSTRACT: With increasing environmental concerns and the drive for a greener economy comes an increased desire to improve turbine engine fuel efficiency and reduce emissions. Unfortunately weight reduction techniques used increase the blade loading, which in turn increases the losses. Non-axisymmetric endwall contouring is one of several techniques being investigated to reduce loss in a turbine. An investigation at Durham University produced a non-axisymmetric endwall design for a linear cascade. An adaption of the most promising endwall was investigated in an annular rotating test rig at the CSIR using steady state instrumentation. The current investigation extends those investigations into the unsteady time domain. Previous investigations found that a generic rotor endwall contour improved efficiency by controlling the endwall secondary flow vortex system in both a linear cascade and an annular 1½ stage rotating test turbine. The current research was aimed at determining if there were any unsteady effects introduced by the contoured endwall. The approach was unique in that it investigated the unsteady effects of an endwall contour originally designed for a linear cascade both experimentally and numerically at three incidence angles (positive, zero and negative to represent increased load, design load and decreased load respectively), the results of which are openly available. Unsteady experimental hotfilm results showed that the endwall contour made the velocity profile more radially uniform by reducing the strength of the endwall secondary flow vortex system. The fluctuations in the velocity were also reduced producing a more temporally uniform velocity profile. The FFT magnitude of the velocity at the blade passing frequency below midspan was also reduced. It was found that the reduction in the endwall secondary flow vortex system due to the contour increased with increasing loading. Numerical results showed that the oscillations in the flow were small and did not penetrate the boundary layer. The contoured rotor was forward and aft loaded when compared to the annular rotor, resulting in a weaker cross passage pressure gradient which allowed the endwall secondary flow vortex system to be less tightly wrapped. Numerical results did not show a significant difference in the oscillations observed in the annular and contoured rotor. A new objective function for use in the endwall optimisation process was proposed that acts as a proxy for efficiency, but is less prone to uncertainty in the results. When used on the current results it shows the same trend as efficiency. It remains to be used to design an endwall for full validation. / AFRIKAANSE OPSOMMING: Met ’n toenemende omgewingsbesorgdheid en die strewe na ’n groener ekonomie kom ’n toenemende behoefte om turbine enjin brandstofdoeltreffendheid te verbeter en vrystellings te verlaag. Ongelukkig het gewigsbesparingstegnieke wat gebruik is die lemlading verhoog, wat op sy beurt die verliese verhoog. Nie-assimmetriese endwandprofilering is een van verskeie tegnieke wat ondersoek word om verliese in ’n turbine te verminder. ’n Ondersoek by die Universiteit van Durham het ’n nie-assimmetriese endwandontwerp vir ’n lineêre kaskade gelewer. ’n Aanpassing van die mees belowende endwand is in ’n annulêre roterende toetsopstelling by die WNNR getoets, deur gebruik te maak van bestendige toestand instrumentasie. Die huidige ondersoek brei daardie ondersoeke uit na die nie-bestendige verwysingsraamwerk . Vorige ondersoeke het bevind dat die generiese rotor endwandprofiel doeltreffendheid verbeter as gevolg van die beheer van die endwand sekondêre vloei draaikolkstelsel in beide ’n lineêre kaskade sowel as ’n annulêre 1½ stadium roterende toetsturbine. Die huidige navorsing was daarop gemik om vas te stel of die endwandprofiel enige onbestendige effekte tot gevolg gehad het. Die benadering was uniek in die sin dat dit die onbestendige effekte ondersoek het van ’n endwandprofiel wat oorspronklik ontwerp is vir ’n lineêre kaskade beide eksperimenteel en numeries op drie invalsshoeke (positief, nul en negatief om onderskeidelik verhoogde lading, ontwerplading en verlaagde lading te verteenwoordig), waarvan die resultate algemeen beskikbaar is. Onbestendige eksperimentele warmfilm resultate het getoon dat die endwandprofiel die snelheidsprofiel meer radiaal uniform gemaak het deur die vermindering van die sterkte van die endwand sekondêre vloei werwelstelsel. Die skommelinge in die snelheid is ook verminder wat ’n meer tyduniforme snelheidsprofiel gelewer het. Die FFT (Fast Fourier Transform) grootte van die snelheid van die lem verbygaan frekwensie onder lem midbestek het ook verminder. Daar was bevind dat die vermindering in die endwand sekondêre vloei draaikolkstelsel as gevolg van die endwandprofiel toeneem met toenemende lading. Numeriese resultate het getoon dat die ossilasie in die vloei klein was en nie die grenslaag binnegedring het nie. Die rotor met gevormde wand het ’n voor- en agterlading gehad in vergelyking met die rotor met annulêre wand, wat tot ’n laer drukgradient dwarsop die vloeirigting gelei het, die endwand sekondêre vloei draaikolkstelsel minder beperk het. Numeriese resultate het nie ’n beduidende verskil in die ossilasies tussen die annulêre en gevormde rotorwand getoon nie. ’n Nuwe doelwitfunksie vir gebruik in die endwand optimersproses is voorgestel wat dien as ’n plaasvervanger vir doeltreffendheid, maar minder geneig is tot onsekerheid in die resultate. Wanneer dit gebruik word op die huidige resultate toon dit dieselfde tendens as doeltreffendheid. Dit moet nog gebruik word in die ontwerp van ’n endwand vir volledige bevestiging.

<p>The continuous trend in gas turbine design towards lighter, more powerful and more reliable engines on one side and use of alternative fuels on the other side renders flutter problems as one of the paramount challenges in engine design. Flutter denotes a self-excited and self-sustained aeroelastic instability phenomenon that can lead to material fatigue and eventually damage of structure in a short period of time unless properly damped. The design for flutter safety involves the prediction of unsteady aerodynamics as well as structural dynamics that is mostly based on in-house developed numerical tools. While high confidence has been gained on the structural side unanticipated flutter occurrences during engine design, testing and operation evidence a need for enhanced validation of aerodynamic models despite the degree of sophistication attained. The continuous development of these models can only be based on the deepened understanding of underlying physical mechanisms from test data.</p><p>As a matter of fact most flutter test cases treat the turbomachine flow in two-dimensional manner indicating that the problem is solved as plane representation at a certain radius rather than representing the complex annular geometry of a real engine. Such considerations do consequently not capture effects that are due to variations in the third dimension, i.e. in radial direction. In this light the present thesis has been formulated to study three-dimensional effects during flutter in the annular environment of a low-pressure turbine blade row and to describe the importance on prediction of flutter stability. The work has been conceived as compound experimental and computational work employing a new annular sector cascade test facility. The aeroelastic response phenomenon is studied in the influence coefficient domain having one blade oscillating in various three-dimensional rigid-body modes and measuring the unsteady response on several blades and at various radial positions. On the computational side a state-of-the-art industrial numerical prediction tool has been used that allowed for two-dimensional and three-dimensional linearized unsteady Euler analyses.</p><p>The results suggest that considerable three-dimensional effects are present, which are harming prediction accuracy for flutter stability when employing a two-dimensional plane model. These effects are mainly apparent as radial gradient in unsteady response magnitude from tip to hub indicating that the sections closer to the hub experience higher aeroelastic response than their equivalent plane representatives. Other effects are due to turbomachinery-typical three-dimensional flow features such as hub endwall and tip leakage vortices, which considerably affect aeroelastic prediction accuracy. Both effects are of the same order of magnitude as effects of design parameters such as reduced frequency, flow velocity level and incidence. Although the overall behavior is captured fairly well when using two-dimensional simulations notable improvement has been demonstrated when modeling fully three-dimensional and including tip clearance.</p>

Le contexte général de cette thèse s'inscrit dans le cadre de la modélisation eulérienne instationnaire des écoulements turbulents anisothermes gaz - particules. La modélisation de ces écoulements est cruciale pour de nombreuses applications industrielles et pour la prédiction de certains phénomènes naturels. Par exemple, la combustion diphasique dans les moteurs automobiles et aéronautiques est précédée par l'injection et la dispersion de carburant liquide dans la chambre de combustion. Les phénomènes mis en jeu exigent alors une prédiction locale tenant compte du caractère instationnaire de l'écoulement turbulent et de la présence de géométries complexes. De plus, de nombreuses études expérimentales et numériques récentes ont mis en évidence le rôle prépondérant de l'inertie des particules sur les mécanismes de dispersion et de concentration préférentielle en écoulement turbulent. Ceci rend donc indispensable la prise en compte de ces mécanismes dans la modélisation diphasique. Au cours de ce travail de thèse, une approche eulérienne locale et instantanée a été développée pour prédire les écoulements gaz-particules anisothermes et turbulents. Elle est basée sur l'approche statistique du Formalisme Eulérien Mésoscopique (MEF) introduite par Février et al. (JFM, 2005). Cette approche a été ici étendue aux variables thermiques pour la prise en compte du caractère anisotherme de l'écoulement. Cette approche a été ensuite utilisée dans le cadre de la méthode des moments (Kaufmann et al., JCP, 2008), et un système d'équations locales et instantanées pour la phase dispersée a été proposé. La modélisation au premier ordre exige la fermeture des moments de second ordre apparaissant dans les équations de la quantité de mouvement et de l'énergie. La proposition de telles relations constitutives fait l'objet d'une partie de la thèse. Afin de fournir une méthode capable de prédire le comportement local, instantané et anisotherme de la phase dispersée dans des configurations `a une échelle réaliste, les équations pour la phase dispersée ont été filtrées et une modélisation aux grandes échelles (LES) est effectuée. Cette modélisation étends, par la prise en compte des variables thermiques, le travail de Moreau et al. (FtaC, 2010) sur l'approche LES Euler-Euler en conditions isothermes. L'approche complète est enfin appliquée aux résultats de simulation numérique d'un jet plan turbulent gazeux froid, chargé en particules, dans une turbulence homogène isotrope chaude monophasique. / The aim of this thesis is to provide an Eulerian modeling for the dispersed phase interacting with unsteady non-isothermal turbulent flows. The modeling of these flows is crucial for several industrial applications and for predictions of natural events. Examples are the combustion chambers of areo engines where the combustion is preceded by the injection and dispersion of liquid fuel. The prediction of such phenomena involves a local modeling of the mixture for taking into account the unsteady behavior of the turbulent flow and the presence of complex geometries. Moreover, many experimental and numerical studies have recently highlighted the significant role of the particle inertia on the mechanisms of dispersion and preferential concentration. Accounting for such mechanisms is therefore essential for modeling the particle-laden turbulent flows. In this thesis, a local and instantaneous Eulerian approach able to describe and to predict the local behavior of inertial particles interacting with non-isothermal turbulent flows has been developed. It is based on the statistical approach known as Mesoscopic Eulerian formalism (MEF) introduced by Février et al. (JFM, 2005). The statistical approach has been extended to the thermal quantities in order to account for the non-isothermal conditions into the modeling. This formalism is then used in the framework of the moment approach (Kaufmann et al., JCP, 2008) and a system of local and instantaneous equations for the non-isothermal dispersed phase has been suggested. The first order modeling requires to close second-order moments appearing in momentum and energy equations. The proposal of such constitutive relations makes the object of a part of this study. In order to provide an Eulerian approach usable in real configurations at industrial scale, the equations of the dispersed phase are filtered and the approach developed in the framework of the Large-Eddy Simulations. From the work of Moreau et al. (FTaC, 2010), the Eulerian-Eulerian LES approach is then extended to non-isothermal conditions. The whole modeling is then a priori tested against numerical simulations of a cold planar turbulent particle-laden jet crossing a homogeneous isotropic decaying hot turbulence.

The Mississippi River (MR) has been engineered with the development of the levee system, dams for flood control measures, jetties, revetments and dredging of the navigation channel. These alterations have reduced the replenishment of the sediment to the Louisiana Coastal area. To aid in the restoration planning, 1-D numerical models have been calibrated and validated to predict the river response to various changes such as channel modifications, varied flow conditions and hurricane situations. This study utilized the HEC-RAS 4.1 and the CHARIMA (Dr. Forrest Holly, University of Iowa). The models were calibrated for hydrodynamics and sediment using Tarbert Landing discharges (HEC-RAS), Belle Chasse sand concentrations (CHARIMA), and Gulf of Mexico (GOM) stages. The models showed that a large percentage of the river flow is lost over the East Bank downstream of Bohemia which reduces the sand transport capacity of the river. This reach is subject to flow reversals during hurricanes.

This project is intended to demonstrate the current state of knowledge in the prediction of the tonal and broadband noise radiation from a Sevik rotor. The rotor measurements were made at the Virginia Tech Stability Wind Tunnel. Details of the rotor noise and flow measurements were presented by Wisda et al(2014) and Murray et al(2015) respectively. This study presents predictions based on an approach detailed by Glegg et al(2015) for the broadband noise generated by a rotor in an inhomogeneous flow, and compares them to measured noise radiated from the rotor at prescribed observer locations. Discrepancies between the measurements and predictions led to comprehensive study of the flow in the wind tunnel and the discovery of a vortex upstream of the rotor at low advance ratios. The study presents results of RANS simulations. The static pressure and velocity profile in the domain near the rotor's tip gap region were compared to measurements obtained from a pressure port array and a PIV visualization of the rotor in the wind tunnel. / Includes bibliography. / Thesis (M.S.)--Florida Atlantic University, 2015. / FAU Electronic Theses and Dissertations Collection

As estruturas aeronáuticas estão sujeitas a diversas solicitações, devido principalmente às interações com o escoamento aerodinâmico, que podem causar distúrbios e vibrações, comprometendo seu desempenho. As medidas aerodinâmicas aplicadas em uma aeronave podem ser obtidas por simulações computacionais ou testes experimentais. No entanto, podem existir imperfeições na simulação computacional, como por exemplo, se conseguir reproduzir algumas condições de vôo real. Sendo assim, diversas pesquisas vêm sendo realizadas para solucionar estes problemas. Dentre elas estão os testes experimentais feitos em túnel de vento com modelos de escala real em diversas condições de vôo. Desta forma, a construção de um modelo físico de um aerofólio em escala reduzida e a implementação de sensores a este modelo torna-se uma ferramenta bastante importante para validar resultados teóricos e experimentais. Assim, nesse trabalho realizou-se a construção de um modelo de aerofólio NACA0012, o desenvolvimento de um mecanismo de fixação do modelo ao túnel de vento e a implementação de um controlador de oscilação forçada. O modelo físico realiza oscilações harmônicas, em regime não-estacionário. O objetivo do trabalho foi mapear as medidas de pressão atuantes sobre modelo ensaiado em regime estacionário e não-estacionário e fazer a comparação entre os dois casos. / Aeronautical structures are affected by many loads, most of them given by the aerodynamic flow interactions. These flow interactions may cause vibration leading to structural failure, such as cracks and fatigue. The aerodynamic flow interactions can be measured by experiment or predicted by computational simulation. Otherwise, computational simulations on its own are not reliable and can not reproduce a real flight condition, such as the mean atmospheric turbulence dynamic. Many researches has been done to solve these problems for computational simulations. One of them are the wind tunnel experiments with a full scale models in many flight conditions for posterior comparison. For a smaller wind tunnel, a small scale physical prototype well instrumented becomes an important solution to validate theoretical and experimental results. In the present work the construction of a NACA 0012 airfoil model, the development of a constraint mechanism and the implementation of a forced oscilation control system were done. The physical model oscilates with a given frequency. The aim of present work is to map the pressure measurements actuating on the model, testing it under a steady state condition and a transient condition for posterior comparison of both conditions.

Os recentes desenvolvimentos na tecnologia de sistema aeronáutico de geometria variável têm sido motivados principalmente pela necessidade de melhorar o desempenho de aeronaves. O conceito de Morphing Aircraft, por meio da variação da linha de arqueamento, representa uma alternativa para sistemas aeronáuticos mais eficientes. No entanto, para aeronaves de alto desempenho, projetos com estes novos conceitos podem gerar reações aeroelásticas adversas, o que representa uma questão importante e pode vir a limitar esses novos projetos. A compreensão adequada do comportamento aeroelástico devido à variação da linha de arqueamento, particularmente em regimes transônico, compreende uma questão importante. Este trabalho consiste num estudo preliminar das consequências aeroelásticas de um sistema aeronáutico de geometria variável. O objetivo desse trabalho é explorar as repostas aeroelásticas transônicas de um aerofólio com arqueamento variável no tempo. A metodologia para análise aeroelástica é baseada num modelo de seção típica. A integração no tempo do sistema aeroelástico é obtida pelo método de Runge-Kutta de quarta ordem. A representação do escoamento transônico não estacionário foi computada por um código CFD em um contexto de malhas não estruturadas com uma formulação dada pelas equações de Euler-2D. Esses resultados preliminares podem fornecer aos projetistas informações importantes sobre as respostas aeroelásticas de um sistema aeronáutico com variação da linha de arqueamento, permitindo estabelecer um quadro adequado para futuras investigações de controle aeroelástico de sistema aeronáutico de geometria variável. / Recent developments on aircraft variable geometry technologies have been mainly motivated by the need for improving the flight performance. The morphing wing concept, by means of variable camber, represents an alternative towards more efficient lifting surfaces. However, for higher performance aircraft, this technology may lead to designs that create unsteady loads, which may result in adverse aeroelastic responses, which represents an important and limiting issue. Proper understanding of the aeroelastic behavior, particularly in transonic flight regimes, due to variations in camber comprises an important matter. This work is a primary study of aeroelastic consequences of an real-time adaptive aircraft. The objective of this work is to investigate prescribed variations to airfoil camberline and their influence to the aeroelastic response in transonic flight regime. The methodology is based on computational simulations of typical section with unsteady transonic aerodynamics solved with a Computational Fluid Dynamics (CFD) code. The time integration of the aeroelastic system is obtained by Runge-Kutta fourth order. The unsteady transonic flow was computed by a CFD code based on the 2D-Euler equations with unstructured mesh. Prescribed camber variation of a symmetrical airfoil is transferred to the CFD mesh, and aeroelastic responses and loading is assessed. These preliminary results may provide the designers valuable information on the interaction between changes in camber during airfoil aeroelastic reactions, allowing establishing an adequate framework for further aeroelastic control investigations of morphing wings.