Um modelo numérico para simulação da resposta aeroelástica de asas fixas é proposto. A estratégia adotada no trabalho é a de tratar a aerodinâmica e a dinâmica estrutural separadamente e então acoplá-las na equação de movimento. A caracterização dinâmica de uma asa protótipo é feita pelo método dos elementos finitos e a equação de movimento é escrita em função das coordenadas modais. O carregamento aerodinâmico não-estacionário é determinado pelo método de malha de vórtices. A troca de informações entre as malhas estrutural e aerodinâmica é feita através do método de interpolação por splines de superfície e a equação de movimento é resolvida iterativamente no domínio do tempo, utilizando-se um método preditor-corretor. As teorias de aerodinâmica, dinâmica estrutural e do acoplamento entre elas são apresentadas separadamente, juntamente com os respectivos resultados obtidos. A resposta aeroelástica da asa protótipo é representada por curvas de deslocamentos modais em função do tempo para várias velocidades de vôo e a ocorrência de flutter é verificada quando estas curvas divergem (i.e. as amplitudes aumentam progressivamente). Transformadas de Fourier destas curvas mostram o acoplamento de freqüências característico do fenômeno de flutter. / A numerical model for the simulation of the aeroelastic response of fixed wings is proposed. The methodology used in the work is to treat the aerodynamic and the structural dynamics separately and then couple them in the equation of motion. The dynamic characterization of a prototype wing is done by the finite element method and the equation of motion is written in modal coordinates. The unsteady aerodynamic loads are predicted using the vortex lattice method. The exchange of information between the aerodynamic and structural meshes is done by the surface splines interpolation scheme, and the equation of motion is solved interactively in the time domain, employing a predictor-corrector method. The aerodynamic and structural dynamics theories, and the methodology to couple them, are described separately, together with the corresponding obtained results. The aeroelastic response of the prototype wing is represented by time histories of the modal coordinates for different airspeeds, and the flutter occurrence is verified when the time histories diverge (i.e. the amplitudes keep growing). Fast Fourier Transforms of these time histories show the coupling of frequencies, typical of the flutter phenomenon.

Gravel bed braided rivers are distinctive natural environments that provide a wide range of key environmental, economic and recreational services. There is, however, a growing concern that over the twentieth century, an increasing number of braided rivers have metamorphosed into wandering or single thread channels, representing a loss of key habitats, geodiversity and amenity. While in some situations, shifts in channel pattern may be unambiguously linked to abrupt changes in flow or sediment supply, the lack of a theoretical basis underpinning the development and maintenance of braiding makes identification of the cause and effect of channel metamorphosis hazardous. A growing body of research has suggested that the transition between channel patterns may depend on the poorly understood interaction between the flow regime, sediment supply and vegetation colonisation. Such interactions are governed by critical thresholds, due to changes in flow resistance and bank strength associated with the distribution, form and intensity of vegetation colonisation. Subtle changes in flow or sediment supply that promote vegetation growth or indeed remove it through inundation or attrition. This can lead to complex non-linear shifts in the balance of forces that govern sediment transport and bedform morphodynamics, ultimately resulting in one-way changes in channel morphology. There is, therefore, a critical need to develop a quantitative understanding of these feedbacks in order to design sustainable river management programmes that seek to optimize the ecological and socio-economic benefits these rivers offer. During the last three decades, significant advances in the understanding of the morphodynamics of braided rivers have been made through a combination of field and physical experimentation. More recently, the emerging field of numerical modelling has created a new avenue to investigate the processes that govern channel dynamics. While this methodology offers significant promise through the construction of virtual experiments that examine the spectrum of drivers and responses of river systems, such models require careful and critical evaluation before they can be used to guide management practice. The wider goal of this research is to explore the application of a numerical modelling to investigate the feedbacks associated with the development and maintenance of braiding. Specifically, the state-of-the-art numerical model, BASEMENT, was used to examine channel responses to steady, and unsteady flow regimes, with and without the representation of vegetation. The research focuses on four main contributions: 1. The development of a systematic framework to quantify the evolving form and processes of braided rivers that can be used as part of a comprehensive approach to model validation. 2. Simulation of braiding development and maintenance using BASEMENT under steady flow conditions. Model simulations based on the natural prototype of the braided River Feshie were used to examine the sensitivity of emergent channel morphologies to the model parameterisation, focusing in particular on the representation of bank erosion and gravity-driven sediment transport. A novel multi6metric framework for model validation is presented and the results demonstrate the critical importance of lateral bank migration processes in order to maintain braided morphologies under steady flow. 3. A critical evaluation of the simulation of braiding under different form of steady and unsteady flow regimes is presented. These experiments investigate how the morphodynamics of braiding vary under energetically-normalised flow regimes characterized by differences in hydrograph form (peak discharge and duration). This experiment provides a novel insight into the role of flow variation in the maintenance of braiding. 4. Finally, the feedback between flow regimes, sediment transport and vegetation growth are examined using a novel model of vegetation colonisation and die- back. Four scenarios are presented, a no-vegetation model, one based on low growth rate, one based on an intermediate growth rate, and finally a high growth rate parameterisation. These simulations provide a clear insight into the non-linear processes driving channel evolution and demonstrate how subtle changes in the balance between flow frequency and vegetation growth can lead to divergent channel patterns. In summary, this thesis aims to advance our understanding of the morphodynamics of braided rivers and the role numerical models may have in helping to interrogate their behaviour and governing controls. It is hoped that this work may contribute, albeit in a small way, to advancing the science that promotes the sustainability of these fascinating and valuable environments.

Thesis (Ph.D.)--Massachusetts Institute of Technology, Dept. of Aeronautics and Astronautics, 1981. / MICROFICHE COPY AVAILABLE IN ARCHIVES AND AERO / Includes bibliographical references. / by Saeed Farokhi. / Ph.D.

Thesis (Ph.D.)--Massachusetts Institute of Technology, Dept. of Mechanical Engineering, 1980. / MICROFICHE COPY AVAILABLE IN ARCHIVES AND ENGINEERING. / Includes bibliographical references. / by Michael E. McClurken. / Ph.D.

Este trabalho representa a continuidade de estudos envolvendo a problemática dos escoamentos com superfície livre, contemplando a análise do fenômeno transiente em canais, a partir do modelo matemático unidimensional de Saint-Venant. Para tanto, é desenvolvido um modelo computacional em linguagem FORTRAN, capaz de avaliar o comportamento do escoamento não permanente. As equações hidrodinâmicas completas são discretizadas por um esquema completamente implícito de diferenças finitas e aplicadas no modelo computacional para a avaliação de dois casos. O modelo é previamente testado para um caso simples, cujos resultados são analisados viabilizando o modelo. No primeiro caso, o modelo é aplicado ao canal de alimentação da Usina Hidrelétrica Monjolinho em São Carlos-SP, para avaliar a necessidade de vertedouro quando se dá o fechamento brusco da turbina, e a ocorrência da entrada de ar na mesma quando da sua abertura repentina. No segundo caso, procurou-se avaliar o desenvolvimento do escoamento no Canal do Trabalhador, responsável pelo abastecimento da cidade de Fortaleza-CE. Com manobras de enchimento e esvaziamento do sistema, é possível determinar o tempo de antecedência de liga-desliga do sistema de recalque a partir das alturas dágua e velocidades de ocorrência, permitindo também a automação para as operações de controle. Em ambos os casos o modelo reproduziu resultados que ilustram com coerência os conceitos pré-estabelecidos, constituindo numa ferramenta útil para análise do fenômeno transiente nos escoamentos em condutos livres. / This work presents a computational model developed in FORTRAN language for the study of unsteady open-channel flows with the use of Saint-Venant one-dimensional equation. The discretization of hydrodynamic equations are presented in a completely implicit method of finite differences and applied in the model for the investigation of two cases, besides the one used previously to test the model. In the first case, the model is applied for a channel that supplies the Monjolinho hydroelectric plant in Sao Carlos SP, aiming to evaluate the need of a spillway when the turbine is closed and the flow abruptly stopped, as well as the occurrence of air entering the turbine when it is opened instantaneously. In the second case, the model simulates the development of the flow in the Trabalhador channel, responsible for the water supply in the city of Fortaleza - CE, in order to make possible the automation of operational control, based on data of flow velocity and water level. In both cases the model is presented as a useful tool for the analysis of unsteady open-channel flows, showing results and coherency with theory.

The present study is an integral part of a broader study focused on the design and implementation of self-cleaning culverts, i.e., configurations that prevent the formation of sediment deposits after culvert construction or cleaning. Sediment deposition at culverts is influenced by many factors, including the size and characteristics of material of which the channel is composed, the hydraulic characteristics generated under different hydrologic events, the culvert geometry design, channel transition design, and the vegetation around the channel. The multitude of combinations produced by this set of variables makes the investigation of practical situations challenging. In addition to the above considerations, the field observations, and the laboratory and numerical experiments have revealed additional complexities of the flow and sediment transport through culverts that further increase the dimensions of the investigation. The flow complexities investigated in this study entail: flow non-uniformity in the areas of transition to and from the culvert, flow unsteadiness due to the flood wave propagation, and the complex correlation between the flow and sediment hydrographs produced during storm events. To date, the literature contains no systematic studies on sediment transport through multi-box culverts. Similarly, there is limited knowledge about the non-uniform, unsteady sediment transport in channels of variable geometry. Furthermore, there are few readily useable numerical models that can reliably simulate flow and sediment transport in such complex situations. Given the current state of knowledge, the main goal of the present study is to investigate the above flow complexities in order to provide the needed insights for optimizing the culvert design. The research was phased so that field observations were conducted first to understand the culvert behavior in Iowa landscape. Modeling through complementary hydraulic model and numerical experiments was subsequently carried out to gain the practical knowledge for the development of the self-cleaning culvert designs.

Unsteady flow separation represents a highly complex and important area of study within fluid mechanics. The extent of separation and specific time scales over which it occurs are not fully understood and has significant consequences in numerous industrial applications such as helicopters, jet engines, hydroelectric turbines and wind turbines. A direct consequence of unsteady separation is the erratic movement of the separation point which causes highly dynamic and unpredictable loads on an airfoil. Current computational models underestimate the aerodynamic loads due to the inaccurate prediction of the emergence and severity of unsteady flow separation especially in response to a sudden change in the effective angle of attack. To capture the complex flow phenomena over wind turbine blades during stall development, a scaled three-dimensional non-rotating blade model is designed to be dynamically similar to a rotating full-scale NREL 5MW wind turbine blade. A time-resolved particle image velocimetry (PIV) investigation of flow behavior during the stall cycle examines the processes of stall development and flow reattachment. The flow fields are examined through the application of Eulerian techniques such as proper orthogonal decomposition and empirical mode decomposition to capture unsteady separation characteristics within the flow field. Then, for a higher order description, coherent structures such as vortices and material lines are resolved to fully characterize the flow during a full pitching cycle in a Lagrangian framework. The Eulerian and Lagrangian methods described in the present analysis is extended to investigate the spanwise characteristics within the root section of a three dimensional airfoil. Furthermore, statistical information of the separation point is pursued along four spanwise positions during two cases of unsteady separation. The results of the study describe a critical role of surface vorticity accumulation in unsteady separation and reattachment. Evaluation of the unsteady characteristics of the shear layer reveal evidence that the build-up and shedding of surface vorticity directly influence the dynamic changes in separation point. The quantitative characterization of surface vorticity and shear layer stability enables improved aerodynamic design, but also has broader implications on the larger discipline of unsteady fluid dynamics.

The current standard of simulating flood flow in natural river reaches is based on solving the 1-D or 2-D St. Venant equations or using hybrid 1-D/2-D models based on the same equations. These models are not always able to accurately predict floodwave propagation, especially around and downstream of regions where 3-D effects become important, or at times when the main assumptions associated with these models are violated (e.g. flow becomes pressurized due to presence of a hydraulic structure like a bridge or a culvert). A 3-D modeling approach, though computationally much more expensive, is not subject to such limitations and should be able to predict accurately predict floodwave propagation even in regions where 3-D effects are expected to be significant. This dissertation describes the development and validation of a 3-D time-accurate RANS-based model to study flood-related problems in natural environments. It also discusses how results from these 3-D simulations can be used to better calibrate lower order models. Applications are included where the flow becomes pressurized during high flow conditions and the sediment entrainment potential of the flow during the flooding event is estimated. Another important category of applications discussed in the present study are floodwave propagation induced by a sudden dam break failure. Results show that 2-D models show fairly large differences with 3-D model predictions especially in regions where 3-D effects are expected to be significant (e.g. near channel-floodplain transition, in highly curved channels, near hydraulic structures). The study also discusses the use of the validated 3-D model as an engineering design tool to identify the optimum solution for flood protection measures intended to reduce flooding in the Iowa River near Iowa City. 3-D simulation results are also used to discuss hysteresis effects in the relationship between bed shear stress and the stage/discharge. Such effects need to be taken into consideration to accurately estimate erosion associated with the passage of a floodwave.

In this thesis, the seminal work of Stokes concerning the ow of a Navier-Stokesuid due to a suddenly accelerated or oscillating plate and the ow due to torsionaloscillations of an innitely long rod in a Navier-Stokes uid is extended to a uid withpressure dependent viscosity. The viscosity of many uids varies signicantly withpressure, a fact recognized by Stokes; and Barus, in fact, conducted experiments thatshowed that the variation of the viscosity with pressure was exponential. Given sucha tremendous variation, we study how this change in viscosity aects the nature of thesolution to these problems. We nd that the velocity eld, and hence the structureof the vorticity and the shear stress at the walls for uids with pressure dependentviscosity, are markedly dierent from those for the classical Navier-Stokes uid.

A prerequisite for aeroelastic stability prediction inturbomachines is the understanding of the fluctuatingaerodynamic forces acting on the blades. Unsteady transonicflows are complex because of mutual interactions betweentravelling pressure waves, outlet disturbances, shock motion,and fluctuating turbulent boundary layers. Complex phenomenaappear in the shock/boundary layer region and produce phaselags and high time harmonics, which can give a significantcontribution to the overall unsteady lift and torque, andtherefore affect flutter boundaries, cause large localstresses, or even severely damage the turbomachine. The present research work is concerned with theunderstanding of phenomena associated with travelling waves innon-uniform transonic flows and how they affect the unsteadypressure distribution on the surface as well as the far fieldradiated sound. In similitude with turbomachines potentialinteraction, the emphasis was put on the unsteady interactionof upstream propagating acoustic waves with an oscillatingshock in 2D and 3D nozzle flows. Both numerical andexperimental studies are carried out and compared with eachother. Results shows that the unsteady pressure distribution, bothon the bump surface and within the channel, results from thesuperposition of upstream and downstream propagating waves. Itis believed that outlet pressure perturbations propagateupstream in the nozzle, interact in the high subsonic flowregion according to the acoustic blockage theory, and arepartly reflected or absorbed by the oscillating shock,depending on the frequency of the perturbations and theintensity of the SBLI. Furthermore the shock motion amplitudeis found to be related to the mean flow gradients and localwave length of the perturbations rather than to the shockboundary layer interaction. The phase angle between incomingpressure perturbations and the shock motion increases with theperturbation frequency but also depends on the intensity of theSBLI. Additionally the phase angle "shift" observed underneaththe shock location linearly increases with the perturbationfrequency and the shock strength. Such phase shift is criticalregarding aeroelastic stability and might have a significantimpact on the phase angle of the overall aerodynamic forceacting on the blade and shift the aerodynamic damping fromstable to exciting. <b>Keywords:</b>Shock Boundary Layer Interaction, ShockMotion, Unsteady Flows, Nozzle Flows, Potential Interaction,Back Pressure Perturbations.