Birds and bats employ a variety of advanced wing motions in the efficient production of thrust. The purpose of this thesis is to quantify the benefit of these advanced wing motions, determine the optimal theoretical wing kinematics for a given flight condition, and to develop a methodology for applying the results in the optimal design of flapping-wing aircraft (ornithopters). To this end, a medium-fidelity, combined aero-structural model has been developed that is capable of simulating the advanced kinematics seen in bird flight, as well as the highly non-linear structural deformations typical of high-aspect ratio wings. Five unique methods of thrust production observed in natural species have been isolated, quantified and thoroughly investigated for their dependence on Reynolds number, airfoil selection, frequency, amplitude and relative phasing. A gradient-based optimization algorithm has been employed to determined the wing kinematics that result in the minimum required power for a generalized aircraft or species in any given flight condition. In addition to the theoretical work, with the help of an extended team, the methodology was applied to the design and construction of the world's first successful human-powered ornithopter. The Snowbird Human-Powered Ornithopter, is used as an example aircraft to show how additional design constraints can pose limits on the optimal kinematics. The results show significant trends that give insight into the kinematic operation of natural species. The general result is that additional complexity, whether it be larger twisting deformations or advanced wing-folding mechanisms, allows for the possibility of more efficient flight. At its theoretical optimum, the efficiency of flapping-wings exceeds that of current rotors and propellers, although these efficiencies are quite difficult to achieve in practice.

Birds and bats employ a variety of advanced wing motions in the efficient production of thrust. The purpose of this thesis is to quantify the benefit of these advanced wing motions, determine the optimal theoretical wing kinematics for a given flight condition, and to develop a methodology for applying the results in the optimal design of flapping-wing aircraft (ornithopters). To this end, a medium-fidelity, combined aero-structural model has been developed that is capable of simulating the advanced kinematics seen in bird flight, as well as the highly non-linear structural deformations typical of high-aspect ratio wings. Five unique methods of thrust production observed in natural species have been isolated, quantified and thoroughly investigated for their dependence on Reynolds number, airfoil selection, frequency, amplitude and relative phasing. A gradient-based optimization algorithm has been employed to determined the wing kinematics that result in the minimum required power for a generalized aircraft or species in any given flight condition. In addition to the theoretical work, with the help of an extended team, the methodology was applied to the design and construction of the world's first successful human-powered ornithopter. The Snowbird Human-Powered Ornithopter, is used as an example aircraft to show how additional design constraints can pose limits on the optimal kinematics. The results show significant trends that give insight into the kinematic operation of natural species. The general result is that additional complexity, whether it be larger twisting deformations or advanced wing-folding mechanisms, allows for the possibility of more efficient flight. At its theoretical optimum, the efficiency of flapping-wings exceeds that of current rotors and propellers, although these efficiencies are quite difficult to achieve in practice.

L’objectif de ce travail porte sur le développement d’un modèle d’interaction fluide-structure adapté à la dynamique des éoliennes de grandes tailles avec des pales flexibles qui se déforment de manière significative sous l’effet de la pression exercée par le vent. Le modèle développé est basé sur une approche efficace d’IFS partitionnée pour un fluide incompressible et non visqueux en interaction avec une structure flexible soumise a des grandes transformations. Il permet de fournir une meilleure estimation de la charge aérodynamique et de la réponse dynamique associée du système (pales, mat, attachements, câbles) avec un temps de calcul raisonnable et pour des simulations sur des longues périodes. Pour la modélisation structurale, un élément fini de type solide 3D est développé pour l’étude dynamique des pales d’éolienne soumises à des grands déplacements et des grandes rotations. Une amélioration du comportement en flexion est proposée par l’introduction des degrés de liberté en rotation et l’enrichissement du champ de déplacements afin de décrire plus précisément la flexibilité des pales. Cet élément solide est apte de capter des modes de hautes fréquences qui peuvent s’avérer néfastes pour la stabilité du calcul. Deux techniques sont donc proposées pour les contrôler : la régularisation de la matrice masse et le développement des schémas d’intégration robustes de conservation et de dissipation d’énergie. Les chargements aérodynamiques sont modélisés en utilisant la Panel Method. Il s’agit d’une méthode aux frontières, relativement rapide par rapport à la CFD mais suffisamment précise pour calculer la distribution de la pression exercée sur la pale. Les modèles fluide et structure interagissent via un algorithme de couplage partitionné itératif dans lequel des considérations particulières sont prises en compte dans le contexte des grandes transformations. Dans un effort visant à instaurer un indicateur de fatigue dans la méthodologie proposée, des câbles précontraints sont introduits reliant le mat de l’éolienne au support. Une nouvelle formulation complémentaire en termes de contraintes est ainsi développée pour l’analyse dynamique des câbles 3D en comportement élasto-visco-plastique. Chaque méthode proposée a été d’abord validée sur des cas tests pertinents. Par la suite, des simulations numériques d’éoliennes avec des pales flexibles sont effectuées en vue d’affiner la compréhension de leur comportement dynamique et l’intérêt que la flexibilité des pales peut apporter à leur fonctionnement. / In this work, a numerical model of fluid-structure interaction is developed for dynamic analysis of giant wind turbines with flexible blades that can deflect significantly under wind loading. The model is based on an efficient partitioned FSI approach for incompressible and inviscid flow interacting with a flexible structure undergoing large transformations. It seeks to provide the best estimate of true design aerodynamic load and the associated dynamic response of such system (blades, tower, attachments, cables). To model the structure, we developed a 3D solid element to analyze geometrically nonlinear statics and dynamics of wind turbine blades undergoing large displacements and rotations. The 3D solid bending behavior is improved by introducing rotational degrees of freedom and enriching the approximation of displacement field in order to describe the flexibility of the blades more accurately. This solid iscapable of representing high frequencies modes which should be taken under control. Thus, we proposed a regularized form of the mass matrix and robust time-stepping schemes based on energy conservation and dissipation. Aerodynamic loads are modeled by using the 3D Vortex Panel Method. Such boundary method is relatively fast to calculate pressure distribution compared to CFD and provides enough precision. The aerodynamic and structural parts interact with each other via a partitioned coupling scheme with iterative procedure where special considerations are taken into account for large overall motion. In an effort to introduce a fatigue indicator within the proposed framework, pre-stressed cables are added to the wind turbine, connecting the tower to the support and providing more stability. Therefore, a novel complementary force-based finite element formulation is constructed for dynamic analysis of elasto-viscoplastic cables. Each of theproposed methods is first validated with differents estexamples.Then,several numerical simulations of full-scale wind turbines are performed in order to better understand its dynamic behavior and to eventually optimize its operation.

This thesis deals with computing the pressure distribution around a swimming slender fish and the thrust generated by its flapping motion. The body of the fish is modeled as a missile like slender body to which a tail is attached that is modeled as a swept wing. The tail is attached to the tip of the slender body and maintains its slope with it. The motion for the swimming fish is prescribed. The fluid flow is modeled as an unsteady potential flow problem with the flow around the slender body modeled as flow over an array of cylinders of varying radii and the flow over the swept wing modeled using the vortex panel method. The pressure distribution is computed using the unsteady Bernoulli equation. The overall thrust & drag for different parameters are studied and compared