O presente trabalho consiste no desenvolvimento de um código computacional baseado no método dos elementos finitos (MEF), para análise bidimensional de interação fluido-estrutura. Desenvolve-se um código bidimensional para dinâmica de fluidos compressíveis, viscosos ou não, em formulação Euleriana, com base no algoritmo CBS – characteristic based split. Então o código desenvolvido é adaptado para poder ser acoplado a um programa de formulação Lagrangeana para análise dinâmica de estruturas, o que é feito através do emprego da descrição Lagrangeana - Euleriana arbitrária (ALE). Por fim procede-se o acoplamento com um código para análise de estruturas, de formulação posicional e não linear geométrica, baseado no método dos elementos finitos. / The present work consists of the development of a computational code based on the element finite method for fluid-structure interaction analysis. A two-dimensional fluid dynamic Eulerian code is developed based on the CBS algorithm – characteristic based split. Then, the computational code is modified to be coupled with a Lagrangean structures dynamical code by using the arbitrary Lagrangean – Eulerian description (ALE). At the end, the coupling is made with a positional nonlinear geometrical structural dynamics code based on the finite element method.

Apresentam-se as bases teóricas do Método de Elementos Finitos (MEF) piezoelétrico, e a sua aplicação na modelagem de transdutores de ultra-som piezoelétricos, que consiste na determinação das características vibracionais (frequências de ressonância e anti-ressonância, modos de vibrar e coeficiente de acoplamento eletromecânico), obtenção da curva de admitancia, análise transiente da estrutura piezoelétrica sujeita a uma excitação pulsada e análise da influência da variação das constantes piezoelétricas do transdutor com o raio. Utilizando-se o MEF aplicado a acústica obteve-se o campo acústico gerado pelo transdutor operando em onda contínua, bem como iniciou-se o estudo da propagação de ondas num líquido, analisando-se as ondas geradas pela excitação pulsada de um pistão plano em contato com o fluido. Os modos de vibrar e os valores de frequências de ressonância obtidos para um transdutor, foram comparados com os resultados experimentais. / The theoretical basis of piezoelectric finite element method (FEM), and its application in piezoelectric ultrasonic transducer modelling is presented. Among these applications we have the calculation of resonance and antiresonance frequencies, vibration modes, piezoelectric coupling coefficient, admittance curve and transient analysis of piezoelectric structure excited by a short pulse. By means of piezoelectric FEM the influence of variation of piezoelectric constant with radius is analysed. It is discussed three kind of functions (linear, cosinoidal and Gaussian). This technique is called apodization. The acoustic filed generated by the transducer operating in continuous wave (CW) was calculated by using FEM applied to acoustic, considering the fluid-structure coupling. The study of wave propagation in liquids is started by using FEM, analyzing the waves generated by a plane piston in contact with the fluid, excited by a short pulse. For each case discussed above, all boundary conditions and hypothesis assumed in the construction of finite element models are discussed. Although the models considered are circular transducers, the concepts acquired can be expanded to other geometries. The vibrational modes were visualized by means of a laser interferometry technique (ESPI), and the admittance curves were measured by using an impedometer. These results were compared with the FEM results, and the models precision was discussed.

Les travaux de cette thèse ont pour but de développer un prototype d’éco-pilote, nommé EcoNav, permettant d’optimiser la vitesse d’un bateau afin de réduire sa consommation de carburant. EcoNav est composé de plusieurs modules dont : un modèle hydraulique 2D simulant l’écoulement hydrodynamique (vitesse du courant et hauteur d’eau) le long du trajet du bateau; - un modèle de résistance à l’avancement servant à alimenter un modèle de prédiction de la consommation de carburant; - un algorithme d’optimisation permettant de trouver le profil optimal de vitesse. Afin de pouvoir estimer la consommation de carburant, un modèle numérique de la résistance à l’avancement en milieu confiné a été développé durant la première partie de cette thèse. Ce modèle numérique 3D simule l’écoulement du fluide autour du bateau et permet de calculer les forces agissant sur sa coque. La résolution des équations RANS est couplée avec un algorithme de quasi-Newton afin de trouver la position d’équilibre du bateau et calculer son enfoncement. Cette méthode est validée en comparant les résultats numériques avec des résultats expérimentaux issus d’essais en bassin de traction. L’influence de l’enfoncement sur la résistance à l’avancement et la précision de la méthode est étudiée en comparant les résultats numériques obtenus avec et sans enfoncement. La précision des modèles empiriques de prédiction de la résistance à l’avancement est également comparée à celle du modèle numérique. Enfin, le modèle numérique est utilisé afin de déterminer si le confinement en largeur ou en profondeur ont une influence identique sur l’augmentation de résistance à l’avancement. Les résultats de cette étude permettent d’établir si le confinement de la voie d’eau peut être caractérisé à l’aide d’un paramètre unique (coefficient de blocage par exemple) ou bien deux paramètres permettant de distinguer le confinement latéral et vertical. Dans la seconde partie de cette thèse, les méthodes numériques utilisées pour le modèle d’éco-pilote sont décrites et comparées afin de sélectionner celles qui sont le plus adaptées à chaque module. EcoNav est ensuite utilisé afin de modéliser un cas réel : celui du bateau automoteur Oural navigant sur la Seine entre Chatou et Poses (153 km). La consommation optimisée est comparée à la consommation non optimisée, calculée à partir des vitesses AIS observées sur le tronçon étudié. L’influence de la trajectoire du bateau et de son temps de parcours sur sa consommation sont également étudiés. Les résultats de ces investigations ont montré qu’optimiser la vitesse du bateau permet d’obtenir une réduction de la consommation de carburant de l’ordre de 8 % et qu’optimiser la trajectoire du bateau ainsi que prendre en compte des informations en temps réel (disponibilité des écluses, trafic sur le fleuve) peuvent permettre de réaliser des économies de carburant supplémentaires. / An eco-driving prototype, named EcoNav, is developed with the aim of optimizing a vessel speed in order to reduce fuel consumption for a given itinerary. EcoNav is organized in several modules : - a 2D hydraulic model simulating the flow conditions (current speed and water depth) along the itinerary; - a ship resistance model calculating the thrust necessary to counteract the hydrodynamic forces ; - a fuel consumption model calculating the fuel consumption corresponding to the thrust input; - a non linear optimization algorithm calculating the optimal speed profile. In order to evaluate the fuel consumption of an inland vessel, a ship resistance numerical model is developed in the first part of this PhD. This 3D numerical model simulates the flow around an inland self-propelled vessel and evaluates the hydrodynamic forces acting on the hull. A RANS solver is coupled with a quasi-Newton approach to find the equilibrium position and calculate ship sinkage. This method is validated by comparing the results of numerical simulations to towing tank tests. The numerical results with and without sinkage are also compared to study the influence of sinkage on ship resistance and on the accuracy of the method. Additionally, some empirical models are investigated and compared with the accuracy of the numerical method. Finally, the numerical model is used to determine if channel with and water depth restriction contribute to the same amount of ship resistance increase for the same level of restriction. The results of that investigation give insight to whether channel restriction can be characterized by a unique parameter (for instance the blockage ratio) or two parameters to distinguish water depth and channel with effects. In the second part of this PhD, the numerical methods used in the speed optimization model are described and validated. The speed optimization model is then used to simulate a real case: the itinerary of the self-propelled ship Oural on river Seine, between Chatou and Poses (153 km). The optimized fuel consumption is compared with the non-optimized fuel consumption, based on AIS speed profile retrieved on this itinerary. The effects of the ship trajectory and travel duration on fuel consumption are also investigated. The results of those investigations showed that optimizing the ship speed lead to an average fuel saving of 8 % and that using an optimal track and including real time information such as lock availability and river traffic can lead to additional fuel savings.

The marine current turbine (MCT) is an exciting proposition for the extraction of renewable tidal and marine current power. However, the numerical prediction of the performance of the MCT is difficult due to its complex geometry, the surrounding turbulent flow and the free surface. The main purpose of this research is to develop a computational tool for the simulation of a MCT in turbulent flow and in this thesis, the author has modified a 3D Large Eddy Simulation (LES) numerical code to simulate a three blade MCT under a variety of operating conditions based on the Immersed Boundary Method (IBM) and the Conservative Level Set Method (CLS). The interaction between the solid structure and surrounding fluid is modelled by the immersed boundary method, which the author modified to handle the complex geometrical conditions. The conservative free surface (CLS) scheme was implemented in the original Cgles code to capture the free surface effect. A series of simulations of turbulent flow in an open channel with different slope conditions were conducted using the modified free surface code. Supercritical flow with Froude number up to 1.94 was simulated and a decrease of the integral constant in the law of the wall has been noticed which matches well with the experimental data. Further simulations of the marine current turbine in turbulent flow have been carried out for different operating conditions and good match with experimental data was observed for all flow conditions. The effect of waves on the performance of the turbine was also investigated and it has been noticed that this existence will increase the power performance of the turbine due to the increase of free stream velocity.

The velocity field in the wake of a small scale flexible parachute canopy was measured using two-dimensional particle image velocimetry. The experiments were performed in a water tunnel with the Reynolds number ranging from 3.0-6.0 x 104. Both a fully inflated canopy and the inflation phase were investigated in a constant freestream (i.e. an infinite mass condition). The fully inflated canopy experienced a cyclic“breathing" which corresponded to the shedding of a vortex ring from the canopy. The normalized breathing frequency had a value of 0.56 +/- 0.03. The investigation of the canopy inflation showed that during the early stages of the inflation, the boundary layer on the canopy surface remains attached to the canopy while the canopy diameter increases substantially. The boundary layer begins to separate near the apex region when the diameter is ~68% of the fully inflated diameter. The separation point then progresses upstream from the canopy apex region toward the canopy skirt. During this time period, the force rapidly increases to its maximum value while the separation point of the boundary layer moves upstream towards the skirt. The force then declines rapidly and the separated boundary layer rolls-up into a large vortex ring near the canopy skirt. At the same time, the canopy is drawn into an over-expanded state after which the cyclic breathing initiates. The unsteady potential force was estimated from the rate of change of the canopy volume. It contributed no more than 10% of the peak opening force and was only significant during the early stages of inflation. The majority of the opening force was the result of the time rate of change of the fluid impulse. It accounts for approximately 60% of the peak opening force. This result shows that the formation of the viscous wake is the primary factor in the peak drag force of the canopy.

La maîtrise du bruit et des vibrations est un objectif fréquemment rencontré dans le domaine industriel. Qu’il s’agisse de questions de confort ou de sécurité, les domaines d’applications sont nombreux et variés : transport, BTP, ingénierie civile et militaire… Dans cette thèse, un problème de vibroacoustique interne avec couplage fluide-structure est étudié. Il s’agit d’une cavité remplie de fluide dont les parois sont constituées d’une structure sandwich viscoélastique. Les difficultés numériques associées à ce modèle portent sur la non linéarité du matériau et sur les propriétés des opérateurs matriciels manipulés (conditionnement, non symétrie). Le calcul des vibrations du système dissipatif couplé nécessite une valeur initiale, choisie comme la solution du problème conservatif. Cette solution n’étant pas aisée à déterminer, deux solveurs aux valeurs propres basés sur la Méthode Asymptotique Numérique (MAN) sont proposés pour résoudre le problème des vibrations libres du système conservatif. Associant des techniques de perturbation d'ordre élevé et de continuation, la MAN permet de transformer le problème non linéaire de départ en une suite de problèmes linéaires, plus simples à résoudre. Les solutions obtenues sont ensuite utilisées comme point initial pour déterminer la réponse libre du système dissipatif. Un solveur de Newton d’ordre élevé, basé sur les techniques d’homotopie et de perturbation est développé pour résoudre ce problème. Enfin, le régime forcé est étudié. Pour toutes les configurations envisagées, les résultats obtenus mettent en évidence des performances numériques améliorées par rapport aux méthodes classiquement utilisées (Arpack, Newton…). / Noises and vibrations control is a common objective in the industrial field. Whether it is a question of comfort or safety, the fields of application are numerous and varied: transport, building, civil and military engineering… In this thesis, a vibroacoustics interior problem with fluid-structure coupling is studied. A cavity filled of fluid whose walls are made of a sandwich viscoelastic structure is considered. The numerical difficulties associated with this model relate to the non-linearity of the viscoelastic material and the properties of the matrix operators used (conditioning, non-symmetry). The calculation of the vibrations of the coupled dissipative system requires an initial value, chosen as the solution to the conservative problem. Since this solution is difficult to determine, two eigenvalue algorithms based on the Asymptotic Numerical Method (ANM) are proposed to solve the problem of free vibrations of the conservative system. Combining high order perturbation and continuation techniques, ANM transforms the initial non-linear problem into a set of linear problems that are easier to solve. The solutions obtained are then used as the initial point to determine the free vibrations of the dissipative problem. A high order Newton solver, based on homotopy and perturbation techniques, is developed to solve this problem. Finally, the forced harmonic response of the damped system is computed. For all the configurations tested, the results obtained show improved numerical performance compared to the methods conventionally used (Arpack solver, Newton algorithm…).

L’optimisation de la performance des voiliers est un problème difficile en raison de la complexité du systèmemécanique (couplage aéro-élastique et hydrodynamique) et du nombre important de paramètres à optimiser (voiles, gréement,etc.). Malgré le fait que l’optimisation des voiliers est empirique dans la plupart des cas aujourd’hui, les approchesnumériques peuvent maintenant devenir envisageables grâce aux dernières améliorations des modèles physiques et despuissances de calcul. Les calculs aéro-hydrodynamiques restent cependant très coûteux car chaque évaluation demandegénéralement la résolution d’un problème non linéaire d’interaction fluide-structure. Ainsi, l’objectif central de cette thèseest de proposer et développer des méthodes originales dans le but de minimiser le coût numérique de l’optimisation dela performance des voiliers. L’optimisation globale par méta-modèles Gaussiens est utilisée pour résoudre différents problèmesd’optimisation. La méthode d’optimisation par méta-modèles est étendue aux cas d’optimisations sous contraintes,incluant de possibles points non évaluables, par une approche de type classification. L’utilisation de méta-modèles à fidélitésmultiples est également adaptée à la méthode d’optimisation globale. Les applications concernent des problèmesd’optimisation originaux où la performance est modélisée expérimentalement et/ou numériquement. Ces différentes applicationspermettent de valider les développements des méthodes d’optimisation sur des cas concrets et complexes, incluantdes phénomènes d’interaction fluide-structure. / Sailing yacht performance optimization is a difficult problem due to the high complexity of the mechanicalsystem (aero-elastic and hydrodynamic coupling) and the large number of parameters to optimize (sails, rigs, etc.).Despite the fact that sailboats optimization is empirical in most cases today, the numerical optimization approach is nowconsidered as possible because of the latest advances in physical models and computing power. However, these numericaloptimizations remain very expensive as each simulation usually requires solving a non-linear fluid-structure interactionproblem. Thus, the central objective of this thesis is to propose and to develop original methods aiming at minimizing thenumerical cost of sailing yacht performance optimization. The Efficient Global Optimization (EGO) is therefore appliedto solve various optimization problems. The original EGO method is extended to cases of optimization under constraints,including possible non computable points, using a classification-based approach. The use of multi-fidelity surrogates isalso adapted to the EGO method. The applications treated in this thesis concern the original optimization problems inwhich the performance is modeled experimentally and/or numerically. These various applications allow for the validationof the developments in optimization methods on real and complex problems, including fluid-structure interactionphenomena.

O presente trabalho tem por objetivo desenvolver uma formulação numérica baseada em Análise Isogeométrica para o estudo de problemas de interação fluido-estrutura (IFE) em aplicações envolvendo corpos rígidos submersos, onde escoamentos incompressíveis de fluidos Newtonianos com superfície livre são considerados. Propõe-se o emprego da Análise Isogeométrica por permitir a unificação entre os procedimentos de pré-processamento e análise, melhorando assim as condições de continuidade das funções de base empregadas tanto na discretização espacial do problema como na aproximação das variáveis do sistema de equações. O sistema de equações fundamentais do escoamento é formado pelas equações de Navier-Stokes e pela equação da conservação de massa, descrita segundo a hipótese de pseudo-compressibilidade, em uma formulação cinemática ALE (Arbitrary Lagrangean- Eulerian). A consideração da superfície livre no escoamento se dá tratando o fluido como um meio bifásico, através do método Level Set. O corpo rígido apresenta não linearidade na rotação e restrições representadas por vínculos elásticos e amortecedores viscosos, sendo a equação de equilíbrio dinâmico resolvida através do método de Newmark. O esquema de acoplamento sólido-fluido adotado é o particionado convencional, que impõe condições de compatibilidade cinemáticas e de equilíbrio sobre a interface sólido-fluido, analisando ambos os meios de maneira sequencial. A discretização das equações governantes é realizada através do esquema explícito de dois passos de Taylor-Galerkin, aplicado no contexto da Análise Isogeométrica. Por fim, são analisados alguns problemas da Dinâmica de Fluidos Computacional, de onde se concluiu que os resultados obtidos são bastante consistentes com os fenômenos envolvidos, com as ferramentas exclusivas da Análise Isogeométrica, como o refinamento k, melhorando a convergência dos resultados. Para escoamentos bifásicos, verificou-se que o método Level Set obteve resultados bastante promissores apresentando, entretanto, uma dissipação numérica excessiva. Propõe-se, para estudos futuros, a elaboração de esquemas numéricos que conservem melhor o volume da fase líquida do escoamento. / The present work aims to development of a numerical formulation based on Isogeometric Analysis for the study of Fluid-Structure Interaction problems in applications involving rigid bodies submerged, considering incompressible Newtonian flows with free surface. The use of the Isogeometric Analysis allows unification between the preprocessing and analysis steps, improving then the continuity of the base functions employed, both in the spatial discretization and approximation of the variables in the system of equations. The fundamental flow equations are formed by the Navier-Stokes and the mass conservation, described by de pseudo-compressibility hypothesis, in an ALE (Arbitrary Lagrangean-Eulerian) kinematic formulation. The free surface consideration of the flow is handled treating the fluid like a two- phase medium, using the Level Set method. The rigid body considers nonlinearity in rotation, and restrictions represented by elastic springs and viscous dampers, with the dynamic equilibrium equation being resolved using the Newmark’s method. The solid-fluid coupling scheme is the conventional partitioned, which imposes kinematics and equilibrium compatibility conditions on the solid-fluid interface, analyzing both mediums in a sequential manner. The governing equations are discretized using the explicit two step Taylor-Galerkin method, applied in an Isogeometric Analisys context. Finally, some Computational Fluid Dinamics problems are analysed, from which it was concluded that the results obtained are quite consistent with phenomena involved, with the unique tools of Isogeometric Analysis, such as k-refinement, improving the convergence of the results. For biphasic flows, it was verified that the Level Set method obtained very promising results, presenting, however, an excessive numerical dissipation. For future studies, it is proposed the elaboration of numerical schemes that better preserve the volume of the liquid phase of the flow.

La conception d’une coque de sous-marin est cruciale pour son opérabilité et la sécurité de l’équipage, mais également complexe. En effet, les ingénieurs doivent prendre en compte à la fois la légèreté, la discrétion acoustique et la résistance de la coque à la pression d’immersion et aux attaques extérieures. Les explosions sous-marines représentent une menace de premier ordre pour l’intégrité de la coque, dont le comportement doit être correctement analysé. Les travaux présentés dans cette thèse portent sur le développement d’une méthode simplifiée, basée sur des formulations analytiques, pour étudier la réponse mécanique d’un cylindre profondément immergé à une explosion sous-marine. Le but de cette méthode est de fournir aux ingénieurs une estimation rapide des dommages subis par la coque cylindrique, leur permettant de simuler un grand nombre de scénarios d’explosion. Dans ce travail de thèse, le modèle de la corde plastique sur fondation plastique est repris et adapté à l’étude d’un cylindre immergé soumis à un chargement explosif, pour lequel les effets de pression d’immersion et d’interaction fluide-structure sont à prendre en compte. Une modélisation simplifiée de l’interaction fluide-structure est couplée avec le modèle de corde plastique sur fondation plastique, en considérant d’une part les effets de la pression d’immersion et d’autre part le raidissage circonférentiel du cylindre. Des expressions analytiques sont développées pour calculer l’enfoncement de la coque ainsi que son énergie de déformation plastique. Les résultats obtenus sont comparés à des résultats d’essais et de simulations numériques, ce qui permet de valider progressivement la méthode simplifiée, mais aussi de mettre en évidence ses limitations. / The design of a hull of submarine is crucial for its operability and the safety of the crew, but also complex. Indeed, the engineers have to take into account at the same time lightness, acoustic discretion and resistance of the hull to immersion pressure and to environmental attacks. Underwater explosions represent a first-rate threat to the integrity of the hull, whose behavior needs to be properly analyzed. The works presented in this thesis concern the development of a simplified method, based on analytical formulations, to study the mechanical behavior of a cylinder deeply immersed subjected to an underwater explosion. The purpose of this method is to give engineers a fast estimation of the damage undergone by the cylindrical shell, allowing them to compute a large number of scenarios of explosion.In the scope of this thesis, the model of the plastic string on plastic foundation is picked up and adapted to the study of an immersed cylinder subjected to an explosive load, for which the effects of immersion pressure and fluid-structure interaction are to be taken into account. This simplified model of the fluid structure interaction is coupled with the model of a plastic string resting on plastic foundation, by considering on one hand the effects of the immersion pressure and on the other hand the circumferential stiffening of the cylinder. Analytical expressions are derived in order to calculate the final deflection of the shell as well as its energy of plastic deformation. The obtained results are compared with results obtained with numerical simulations, which allows to validate gradually the simplified method, but also to highlight its limitations.

Les Moteurs à Propergol Solide sont le siège d'instabilités de combustion qui se manifestent par des fluctuations de la pression interne et, en conséquence, de la poussée induite. Les phénomènes oscillatoires observés résultent d'un couplage entre l'acoustique de la chambre de combustion et une ou plusieurs instabilités hydrodynamiques de l'écoulement interne.Dans cette thèse, on s'intéresse à l'étude de l'écoulement induit par injection pariétale en géométrie axisymétrique. Plus particulièrement, on se concentre sur l'interaction qui peut avoir lieu entre cet écoulement et une structure : les Protections Thermiques de Face.Une étude expérimentale reposant sur le montage gaz froid VALDO de l'ONERA, modélisation représentative des MPS P230, a permis d'examiner l'influence d'obstacles, rigides et en élastomère, sur les instationnarités de l'écoulement. Le caractère instable de l'écoulement induit par injection pariétale est alors confirmé et le rôle primordial joué par les PTF sur les fréquences émergeant au sein de l'écoulement a été identifié.L'analyse de stabilité linéaire de l'écoulement, conduite en approche globale, permet l'étude de sa dynamique en tant qu'amplificateur de bruit. Les m¶mécanismes mis en jeu dans l'instabilité de l'écoulement induit par injection pariétale ainsi que son interaction complexe avec la couche de cisaillement issue de l'obstacle sont analysés. La réponse de l'écoulement à un forçage harmonique est alors examinée à la lumière des résultats expérimentaux.Cette étude a été complétée par une approche numérique de l'Interaction Fluide-Structure sur une configuration jugée critique. Une étape préliminaire à la simulation multi-physique est l'élaboration d'un modèle réaliste du comportement des PTF en élastomère. L'identification et la prise en considération des mécanismes propres aux élastomères dans la loi de comportement est ainsi détaillée. Les résultats issus des simulations sont confrontés aux résultats expérimentaux et rendent possible la validation d'un scénario de couplage dédié à l'Interaction Fluide-Structure au sein des MPS. / Fluid Structure Interaction of an inhibitor with the internal flow induced by wall injection was studied in an axisymmetric cold flow apparatus. Experiments were carried out over a wide range of injection velocities in order to underline how the obstacle not only modifies the mean flowfield but also its entire dynamic behaviour.The resulting instability (from the interaction of the unstable shear layer with the Taylor-Culick flow) exhibits a significant shift with respect to the Taylor-Culick instability and therefore emphasizes the strong impact of the inhibitor on hydrodynamics.The mecanisms responsible of such behaviour are studied in this thesis with the help of global linear stability analysis and multi-physics numerical computations.