Les vésicules sont des gouttes immergées dans un fluide externe visqueux, dont le rayon vautquelques dizaines de micromètres et entourées par une membrane imperméable constituée de lipides, dont l’épaisseur est approximativement 4 nm. La membrane d’une vésicule est un systèmeoriginal du point de vue mécanique : celle-ci présente à la fois des propriétés fluides (les lipidespeuvent s’écouler librement le long de la membrane, mais la surface est incompressible locale-ment) et des propriétés solides (la membrane résiste à la flexion). Les propriétés spécifiques de lamembrane rendent ce système à la fois très déformable et très contraint.Ce manuscrit s’intéresse à la modélisation d’une vésicule soumise à des efforts extérieurs d’o-rigine hydrodynamique, dans le régime de Stokes. Une attention particulière est consacrée à lasituation d’une vésicule qui sédimente. Cette situation est étudiée analytiquement dans le régimedes faibles déformations. Il est montré que plusieurs familles de solutions stationnaires non triv-iales existent, grâce aux propriétés spécifiques de la membrane. L’étude de la sédimentation d’unevésicule est poursuivie par le développement d’un code numérique capable de simuler de grandesdéformations. Pour cela, des méthodes numériques originales de calcul de prise en compte de laflexion et de l’incompressibilité surfacique sont développées. Ce code permet d’étudier la forma-tion d’un tube à l’arrière d’une vésicule en sédimentation. Ces tubes sont de fins (rapport d’aspecttypique longueur/rayon ∼ 100) cylindres connectés à la vésicule d’origine. Il est montré que cesformes tubes sont des formes stationnaires. Un modèle théorique est proposé et comparé auxsimulations numériques. Ce modèle met en lumière l’importance particulière de la tension dansces formes. Une modélisation mécanique basée sur un milieu de Cosserat surfacique courbé estégalement présentée, et permet d’identifier la contribution de la flexion au tenseur des contraintes.Cette contribution est un ingrédient indispensable pour comprendre les formes tubes. / Vesicles are drops of radius of a few tens micrometers, bounded by an impermeable lipidmembrane of approximately 4 nm thickness, and embedded in an external viscous fluid. Thevesicle membrane is an original system from the mechanical point of view : it presents bothincompressible fluid properties (the lipids can flow freely along the membrane, but membraneis incompressible locally) and solid properties (the membrane resists to bending). The specificproperties of the membrane make the system very deformable and very constrained at the sametime.This manuscript deals with the modelisation of a vesicle subjected to external stresses of hydrodynamical origin, in the Stokes regime. A particular attention is paid to the situation of asettling vesicle. This situation is studied analytically in the small deformation regime. It is foundthat several families of non-trivial stationnary shapes exist, owing to the specific properties ofthe membrane. The study of a settling vesicle is pursued by the development of a numerical codeable to deal with large deformations. Original numerical methods are developped to deal with thecomputation of the bending and with the surface incompressibility constraint. This code permitsto study the formation of tether at the rear of a settling vesicle. These tethers are thin (typicalaspect ratio : length/radius ∼ 100) cylinders of membrane connected to the original vesicle. Itis shown that these tethered shapes are stationary shapes. A theoretical model is proposed andcompared to numerical simulations. This model shows the particular importance of tension inthese shapes. A mechanical modelling based on a curved Cosserat surface is also presented, andpermits to identify the bending contribution to the stress tensor. This contribution is a salientingredient to understand tethered shapes.

Dans les propulseurs de fusée, des instabilités aéroacoustiques et des interactions de type fluide-structure sont à l'origine de fortes oscillations de poussées pouvant déranger la poussée du moteur mais également causer des dommages non négligeables. On trouve dans les moteurs de fusée des protections thermiques de face (PTF) coincées entre les pains de propergol. Leurs déplacements se trouvent être la principale cause des interactions fluide structure (IFS) présentes dans les booster. Dans ce contexte, nous avons développé une approche numérique visant à simuler les problèmes d'IFS. Notre méthode se base sur le couplage de deux codes dissociés : l'écoulement est simulé avec CARBUR tandis que la dynamique des structures déformables est traitée par MARCUS. Une loi de comportement hyperélastique a été implémentée dans CARBUR afin de simuler le mouvement des PTF. Une campagne expérimentale a été menée dans notre laboratoire sur le tube à chocs T80 et en guide de validation du couplage des codes, les résultats numériques et expérimentaux ont été confrontés. / In a solid rocket motor, high pressure oscillations induced by aeroacoustic instabilities and fluid structure interaction (FSI) may lead to disturb rocket thrust and cause damages. In the rocket motors, flexible inhibitors made of insulating material are initially bonded to the propellant, and FSI is mainly induced by their displacement. In this context, a numérical approach to simulate FSI problems has been developped. Our method is based on the coupling of two dissociated codes : fluid flow is computed with CARBUR, while the dynamics of deformable structures is simulated by MARCUS. A hyperelastic behaviour law has been implemented in MARCUS in order to simulate the movement of flexible inhibitors. An experimental approach has been leaded in the shock waves tubes (T80) in our laboratory and as a validation of FSI coupling codes, numerical and experimental results have been compared.

Cette thèse contribue à la vérification et à la validation du modèle bifluide de Baer-Nunziato, pour modéliser les phénomènes de transitoires hydrauliques dans les réseaux de tuyauteries industrielles. Il s’agit d’abord de modéliser les écoulements de transitoires hydrauliques avec le modèle bifluide en représentation eulérienne, puis d’étendre ce modèle en formalisme ALE (Arbitrary Lagrangian Eulerian) pour prendre en compte l’interaction fluide-structure (IFS). Pour modéliser les écoulements, des lois de fermetures du modèle bifluide concernant les termes interfaciaux, les termes sources et les lois thermodynamiques (EOS) ont d’abord été étudiées. Ensuite, le système complet a été simulé avec une méthode à pas fractionnaires qui admet deux étapes, l’une pour la résolution de la partie convective, l’autre pour les termes sources. L’ensemble de schémas a été vérifié et étendu aux EOS ‘Stiffened Gas généralisées’ afin de représenter le changement de phase eau-vapeur. Après avoir retrouvé certains phénomènes typiques associés aux transitoires hydrauliques, le modèle bifluide a été validé avec l’expérience de Simpson, l’expérience Canon, et comparé avec deux modèles homogènes sur ces deux expériences. Enfin, une version ALE du modèle bifluide a été mise en œuvre et vérifiée sur un cas de propagation d’ondes de pression dans une conduite flexible. La variation de la célérité des ondes dans le fluide liée au couplage fluide/structure a été bien retrouvée. La validation a été effectuée sur un cas expérimental d'explosion dans une tuyauterie en eau. Les simulations sont en bon accord avec les données expérimentales. / This thesis contributes to the verification and the validation of the Baer-Nunziato (BN) model, to modelize water hammer phenomena in industrial piping systems. It consists of two parts, the first is to modelize water hammer flows with the BN model in Eulerian representation and the second is to extend this model to the ALE (Arbitrary Lagrangian Eulerian) formalism so as to take into account fluid-structure interaction (FSI). To modelize water hammer flows, closure laws of the BN model concerning the interfacial/source terms and the equations of state (EOS) were first studied. Then the whole system was simulated with a fractional step method including two steps, one for the resolution of the convective part, the other for the source terms. All schemes have been extended to ‘generalized Stiffened Gas’ EOS in order to represent phase-change. After regaining some typical phenomena associated with water hammer flows, the BN model was validated with the Simpson experiment, a classical water hammer test case, and the Canon experience, a rapid decompression of fluid in a high pressure duct. Moreover, the model was compared with two homogeneous models on both experiments. Finally, an ALE version of the BN model was implemented, and verified on a case of wave propagation in a ‘single’ phase flow and a two-phase flow in a flexible pipe. The variation of wave propagation speed due to the coupling between the fluid and the structure has been well retrieved. The validation was performed on an experiment which examines the response of a pipe filled with water, subjected to a violent pressure peak (140 bar). The simulations are in good agreement with the experimental data.

Dans cette thèse, nous avons étudié la simulation numérique des phénomènes d’interaction fluide-structure (IFS) par la méthode des éléments finis pour un fluide incompressible et non visqueux en interaction avec une structure flexible. Les modèles numériques développés sont basé sur une approche d’IFS partitionnée. Une amélioration basée sur une compensation des effets de massé ajoutée est proposée au cours de la thèse afin d’assurer la convergence et la stabilité du schéma de couplage partitionné indépendamment de la densité du fluide impliqué. L’approche corrective nécessite une estimation de la matrice de masse ajoutée et demande une légère modification de l’algorithme itératif. Les méthodes proposées ont été validées sur les cas académiques en comparaison avec les solutions analytiques et sont appliqués au cas d’une nouvelle conception de pompe pour tout type de fluides (gaz, liquides, fluide chargé…), en vue d’affiner la compréhension de son fonctionnement et ainsi mieux la caractériser. Les méthodes ainsi que les validations sont publiées sur un article qui a été accepté par le revue scientifique « Computers & Fluids ». Une présentation orale a effectuée pendant la conférence internationale ACE-X2012 à Istanbul et une autre a été accepté par la conférence nationale CSMA-2013 à Giens. / The numerical simulation of fluid-structure interaction (FSI) by the finite element method has been studied in the context of an incompressible and inviscid flow interacting with a very flexible structure.The numerical models developed in this work are based on a partitioned FSI approach. An improvement based on a compensation of the added-mass effect is proposed during the PhD research in order to ensure the convergence and the stability of the partitioned coupling scheme for all fluids regardless of its density. This simple correction requires to estimate the added-mass matrix and to modify slightly the iterative algorithm.The proposed methods were validated by comparing with analytical solutions for several academic cases and are applied to a novel pumping technology, which is applicable to all kinds of fluid (gas, liquid, slurry...). The main objective is to provide a better understanding about its operations and to improve the designing of pump. The methods and the validation cases are published in an article which has been accepted by the scientific review Computers & Fluids. They were also presented during the international conference ACE-X2012 in Istanbul and have been accepted and scheduled for oral presentation during the national conference CSMA-2013 in Giens.

Ce travail de recherche est dédié à la compréhension des mécanismes physiques de l’érosion de cavitation dans un fluide compressible à l’échelle fondamentale de l’implosion d’une bulle de cavitation. Suite à l’implosion d’une bulle de vapeur à proximité d’une surface solide, des très hautes pressions sont générées. Ces pressions sont considérées responsables de l’endommagement (érosion) des surfaces solides observé dans la plupart des applications. Notre approche numérique démarre avec le développement d’un solveur compressible capable de résoudre les bulles de cavitation au sein du code volumes finis YALES2 en utilisant un simple modèle de mélange homogène des phases fluides. Le solveur est étendu à une approche ALE (Arbitraire Lagrangien Eulérien) dans le but de mener des simulations d’interaction fluide-structure sur un maillage mobile. La réponse du matériau solide est calculée avec le code de calcul éléments finis Cast3M, et nous a permis de mener des simulation avec un couplage d’abord monodirectionnel, ensuite bidirectionnel, entre le fluide et le solide. On compare des résultats obtenus à deux dimensions, puis à trois, avec des observations expérimentales. On discute les chargements de pression estimés, et les réponses de différents matériaux pour des implosions de bulle à des différentes distances de la surface. Enfin, à travers l’utilisation de simulations avec couplage bidirectionnel entre fluide et solide, on identifie l’amortissement des chargements de pression pour les différents matériaux. / This research is devoted to understanding the physical mechanism of cavitation erosion in compressible liquid flows on the fundamental scale of cavitation bubble collapse. As a consequence of collapsing bubbles near solid wall, high pressure impact loads are generated. These pressure loads are believed to be responsible for the erosive damages on solid surface observed in most applications. Our numerical approach begins with the development of a compressible solver capable of resolving the cavitation bubbles in the finite-volume solver YALES2 employing a simplified homogenous mixture model. The solver is extended to Arbitrary Lagrangian-Eulerian formulation to perform fluid structure interaction simulation with moving mesh capabilities. The material response is resolved with the finite element solver Cast3M, which allowed us to perform one-way and two-way coupled simulations between the fluid and solid domains. In the end, we draw comparisons between 2D and 3D vapor bubble collapse dynamics and compare them with experimental observations. The estimated pressure loads on the solid wall and different responses of materials for attached and detached bubble collapses are discussed. Finally, the damping of pressure loads by different materials is identified with two-way coupled fluid-structure interaction.

Les réseaux de tuyauterie industriels sont le siège de niveaux importants de vibrations induites par l’écoulement qui peuvent mener à la rupture par fatigue des installations. La présente étude se concentre sur l’analyse et modélisation simplifiée de la source vibratoire associée au passage d’un écoulement liquide turbulent par un coude à 90°. Une approche combinant expériences et simulation a été conduite. Une boucle de circulation d’eau munie d’un coude transparent a été conçue afin de permettre des me-sures de vitesse à l’intérieur du coude. Pour cette finalité, les techniques de Vélocimétrie par Imagerie de Particules (PIV) plane et stéréoscopique ont été utilisées. La pression pariétale et les vibrations du banc d’essais ont été mesurées simultanément. Plusieurs configurations d’écoulement ont été testées afin d’obtenir une riche base de données couplées reliant l’écoulement fluide à l’excitation dynamique des parois et, finalement, à la réponse vibratoire de la structure. En parallèle, l’écoulement instationnaire d’eau dans le coude a été simulé au moyen d’une approche du type Simulation des Grandes Echelles (LES). La simulation fluide a permis d’étudier en détails la topologie de l’écoulement turbulent au passage du coude ainsi que le champ instationnaire de pression fluctuante induit sur la paroi. Finale-ment, un ensemble d’outils statistiques a été appliqué aux données expérimentales et numériques afin de proposer un modèle simplifié des transferts qui relient l’écoulement turbulent à la sollicitation dynamique de la structure contenant le coude. / The flow of dense fluids within thin-walled piping systems may lead to significant levels of Flow-Induced Vibration, mainly in the vicinities of singularities such as obstacles inserted into the flow, sudden changes of cross-sectional area or flow direction. This study focuses on the analysis and reduced-order modelling of the vibrational source associated with the turbulent flow of liquids through a 90° elbow.A mixed experimental-computational approach is undertaken. A closed water loop containing a transparent elbow was designed in order to allow for fluid velocity measurements inside the singularity. To this purpose, planar and stereoscopic Particle Image Velocimetry (PIV) were employed. Wall pressure fluctuations and structural vibrations were measured simultaneously. Several flow configurations were tested in order to obtain a large coupled database linking the flow to the dynamic excitation, and then to the vibration response of the structure.In parallel, the unsteady water flow through the elbow was computed using Large-Eddy Simulation (LES). The fluid simulation allowed for a detailed study of the turbulent flow through the singularity and the unsteady pressure field induced on the piping walls. Finally, a set of statistical tools was applied to both experimental and computational data in order to propose a reduced-order model of the transfer function that links the tur-bulent flow to the dynamic excitation of the elbowed piping structure.

Pour améliorer la conception des turbomachines, les industriels doivent appréhender des phénomènes aéroélastiques complexes présents dans les compresseurs comme les cycles limites d’interaction fluide-structure des fans. La compréhension et la modélisation de ces phénomènes impliquent de développer des modèles numériques complexes intégrant des phénomènes multi-physique et de valider ces modèles à l’aide de bancs d’essais. Le banc d’essai du compresseur CREATE est instrumenté pour étudier des instabilités aérodynamiques couplées à des vibrations, notamment sur le rotor du premier étage, et permet de valider des modèles numériques. La modélisation de l’écoulement en amont du premier étage du compresseur à l’aide du logiciel Turb’Flow, développé pour l’étude des écoulements dans les compresseurs aéronautiques, a permis de mettre en évidence l’importance des conditions limites d’entrée pour l’obtention de résultats précis. En particulier, il a été possible de modéliser correctement l’ingestion d’une alimentation non-homogène en entrée de la roue directrice d’entrée. Ce phénomène peut se produire en amont des fans et interagir avec un mode de la structure. Une stratégie de couplage partitionné fort explicite dans le domaine temporel a été introduite dans le logiciel Turb’Flow. Comme cette méthode présente un risque de décalage temporel à l’interface fluide-structure, une attention particulière a été portée à la modélisation de la conservation de l’énergie à cette interface. La conservation de l’énergie à l’interface est cruciale quand les déplacements sont importants et quand un comportement non-linéaire fort apparaît entre le fluide et la structure (onde de choc et amortissement structurel nonlinéaire). Parallèlement au développement du module aéroélastique, le schéma implicite de Runge- Kutta d’ordre 3 en temps (RKI-3) a été développé et évalué sur un cas de dynamique (vibration d’une aube de turbine transsonique) et sur un cas de propagation d’onde de choc. L’utilisation du schéma RKI-3 permet d’augmenter, à iso-précision, d’un ordre le pas de temps par rapport aux schémas de Gear et de Newmark. S’il apporte un gain en temps CPU pour l’étude de la dynamique des structures, il est pénalisant dans le cadre de simulation URANS. Cependant, le schéma RKI-3 est utilisable dans le cadre de simulations couplées fluide-structure. / To increase turbomachinery design, manufacturers have to comprehend complex aeroelastic phenomena involving compressors like fluid-structure interaction limit cycles of fans. The understanding and the modeling of these phenomena involve developing complex solvers coupling techniques and validating these techniques with bench tests. The bench test of the CREATE compressor is instrumented to study the coupling between aerodynamic instabilities and structure vibration, in particular on the first stage rotor, and allows to validate numerical techniques. The flow modeling upstream to the first stage with the Turb’Flow flow solver (targeting turbomachinery applications) shows that, to have accurate results, inlet limit conditions must take into account. The ingestion of non-homogeneous flow upstream to the inlet guide vane is accurately modeled. This phenomenon can appear upstream to fans and interact with structure Eigen-modes. Explicit partitioned strong coupling considered in time domain was implemented in a Turb’Flow flow solver. As there is a risk of time shift at the fluid-structure interface, careful attention should be paid to energy conservation at the interface. This conservation is crucial when displacements are large and when strong non-linear behaviors occur in both fluid and structure domains, namely shock waves, flow separations and non-linear structural damping. In parallel with coupling technique development, the three-order implicit Runge-Kutta scheme (RKI-3) was implemented and validated on a structure dynamic case (transonic turbine blade vibration) and on a case of shock waves propagation. The RKI-3 scheme allows increasing the time step of one order of magnitude with the same accuracy. There is a CPU time gain for structure dynamics simulations, but no for URANS simulations. However, the RKI-3 scheme can be to use for fluid-structure coupling simulations. The coupling technique was validated on a test case involving tube in which the shock wave impinges on a cross flow flexible panel, initially at rest. This case allows modeling an interaction between sonic flow and a panel movement with a tip clearance. Some numerical simulations were carried out with different temporal schemes. The RKI-3 scheme has no influence on results (compared with Gear and/or Newmark scheme) on the energy conservation at the fluid-structure interface. Compared to experimental results, pressure is in fairly good ix Liste des publications agreement. The analysis of numerical results highlighted that a vertical shock tube with up and down waves creates pressure fluctuation. Frequency is under predicted and amplitude is not in fairly good agreement. The panel root modeling might be questionable.

Le palier hydrodynamique est une solution de guidage en rotation particulièrement appréciée pour ses caractéristiques d’amortissement à hautes vitesses de rotation. Cependant les performances des machines tournantes lubrifiées par un film fluide sont impactées par des effets non linéaires difficiles à analyser. La prédiction du comportement du système par la simulation nécessite une modélisation avancée de l’écoulement de lubrifiant dans le palier hydrodynamique. Enfin, l’état de surface semble avoir un impact important sur l’écoulement du fluide lubrifiant, lui-même agissant sur les caractéristiques statiques et dynamiques des parties tournantes. Cette étude vise à améliorer les modèles numériques liés à l’impact de l’état de surface des paliers hydrodynamiques sur la dynamique de ligne d’arbre. La méthode d’homogénéisation multi-échelles a été utilisée à cet effet dans un algorithme multi-physiques pour décrire l’interaction entre la structure flexible en rotation et les films fluides des supports de lubrification. Différents modèles ont été utilisés pour prendre en compte la présence de zone de rupture de film lubrifiant. Des méthodologies non-linéaires fréquentielles ont été mises en place afin de permettre l’étude paramétrique des solutions périodiques d’un tel système et de leur stabilité. Afin de confronter ce modèle complexe à la réalité, un banc d’essai miniature a également été conçu. Différents échantillons présentant des états de surface modifiés par ablation à l’aide de LASER femto-seconde ont été testés. L’étude expérimentale a permis de vérifier certaines tendances prévues par la simulation. Des améliorations des performances des paliers hydrodynamiques par rapport aux vibrations auto-entretenues du système ont été démontrées pour certaines textures. En revanche toutes les améliorations ne sont pas prédites par les algorithmes d’homogénéisation multi-échelles. La présence de recirculation dans les aspérités du motif a été mise en évidence à partir de la résolution locale des équations de Navier-Stokes. Ce résultat participe à la remise en question des hypothèses classiques utilisées en texturation, et peut justifier les améliorations obtenues expérimentalement avec les paliers texturés. / The hydrodynamic bearing provides good damping properties in rotating machineries. However, the performances of rotor-bearings systems are highly impacted by nonlinear effects that are difficult to analyze. The rotordynamics prediction requires advanced models for the flow in the bearings. The surface of the bearings seems to have a strong impact on the lubricant flow, acting on the static and dynamic properties of the rotating parts. This study aims to enhance the simulation of the bearings’ surface state effect on the motion of the rotating shaft. The flexible shaft interacts with textured hydrodynamic bearings. Multi-scales homogenization is used in a multi-physics algorithm in order to describe the fluid-structure interaction. Different models are used to account for the cavitation phenomenon in the bearings. Nonlinear harmonic methods allow efficient parametric studies of periodic solutions as well as their stability. Moreover, a test rig has been designed to compare predictions to real measurements. Several textured shaft samples modified with femto-seconds LASER surface texturing are tested. In most cases the experimental study showed similar results than the simulation. Enhancements of the vibration behaviors of the rotor-bearing system have been revealed for certain texturing patterns. The self-excited vibration, also known as "oil whirl" phenomenon, is stabilized on a wide rotating frequency range. However, the simulation tool does not predict well the enhancements that are observed. Vortices in surface texturing patterns have been revealed numerically with Navier-Stokes equation resolution. These results are opposed to the classical lubrication hypothesis. It is also a possible explanation of the enhancements that are experimentally measured with textured bearings.

Esta investigación busca evaluar la Interacción entre el Suelo y la Estructura (ISE) en reservorios elevados tipo INTZE de tres diferentes capacidades, cimentados en diferentes tipos de suelos, ubicados en la zona con más peligro sísmico en el Perú. Para esto se planteó realizar una muestra de 27 configuraciones estructurales, con variaciones de altura y diámetro del fuste, y diámetro de losa de cimentación, las cuales están cimentadas en cuatro tipos de suelos de acuerdo con la norma Sismorresistente E.030-16. Se evaluará la influencia de la carga sísmica, mediante un análisis dinámico espectral definido por los lineamientos de las normas nacionales E.030-16 y complementariamente por códigos y estándares internacionales como el ACI 350.3-06. Para el dimensionamiento de los reservorios se aplicó el criterio de Otto Intze, análisis estático de esfuerzos y criterios del código ACI 350.3-06 y ACI 307-48. Estos reservorios fueron modelados en el software SAP2000, contando con 216 simulaciones en total. Se aplicó el modelo de Housner (1963) para la Interacción Fluido-Estructura y para la ISE se aplicó el modelo recomendado por el código FEMA P-750. Para analizar los resultados de interés, como las fuerzas internas en el fuste, la cortante basal, momento de volteo, desplazamientos y periodos, se aplicó la prueba de normalidad de Shapiro Wilk para la variación porcentual que se genera al evaluar la ISE, con el fin de determinar intervalos de variación con una probabilidad del 95%. Además, se aplicó la prueba de independencia Chi-Cuadrado para determinar cuantitativamente si la ISE influye en la reducción de las respuestas de interés. Se concluye de la prueba de independencia que, para un nivel de significancia de 15% la ISE influye en la reducción de la fuerza cortante basal y momento de volteo. Así mismo, para un nivel de significancia de 5% la ISE influye en la reducción de las fuerzas internas locales de fuste. / This research seeks to evaluate the soil structure interaction (SSI) in elevated water tanks of three different capacities, founded on different types of soil, all located in the most seismic hazard zone in Peru. Twenty-seven models were evaluated, all with different structural adjustments, including variations in the capacity, height, and diameter of the foundation slab, which are founded on four types of soil, according to the Peruvian Code for Earthquake Resistant Design E.030-16. We are going to analyze the seismic load influence through a spectral dynamic analysis, according to the Peruvian Code E.030-16, and the ACI 350.3-06 code for elevated tanks. For the sizing of the elevated tanks, the Otto Intze criterion, the static stress analysis, the ACI 350.3-06 and the ACI 307-48 code were applied. These elevated water tanks were modeled in the SAP2000 software, with 216 simulations. The Housner model (1963) was applied for the Fluid Structure Interaction, and for the Soil-Structure Interaction (SSI), we applied the model recommended by the FEMA P-750 code. To analyze the results of interest, such as internal forces in the shaft, base shear, rocking moment, displacements and periods, we applied the Shapiro Wilk normality test for the percentage variation that is generated when we evaluate the SSI, in order to determinate variations intervals with a probability of 95%. Also, we applied the Chi-Square test of independence to determinate, quantitatively, if the results are influenced by the SSI. We concluded, from the independence test that, for a significance level of 15%, the SSI influences the reduction of the basal shear force and the rocking moment. Likewise, for a significance of 5%, the SSI influences the reduction of the internal shaft forces. / Tesis

Um dos grandes desafios da engenharia atualmente é viabilizar soluções computacionais que reduzam o tempo de processamento e forneçam respostas ainda mais precisas. Frequentemente surgem propostas com as mais diversas abordagens que exploram novas formas de resolver tais problemas ou tentam, ainda, melhorar as soluções existentes. Uma das áreas que se dedica a propor tais melhorias é a computação paralela e de alto desempenho - HPC (High Performance Computing). Técnicas que otimizem o tempo de processamento, algoritmos mais eficientes e computadores mais rápidos abrem novos horizontes possibilitando realizar tarefas que antes eram inviáveis ou levariam muito tempo para serem concluídas. Neste projeto propõe-se a implementação computacional de uma arquitetura de computação paralela com o intuito de resolver, de forma mais eficiente, em comparação com a arquitetura sequencial, problemas de Dinâmica dos Fluidos e Interação Fluido-Estrutura e que também seja possível estender esta arquitetura para a resolução de outros problemas relacionados com o Método dos Elementos Finitos. O objetivo deste trabalho é desenvolver um algoritmo computacional eficiente em linguagem de programação científica C++ e CUDA - de propriedade da NVIDIAr - tendo como base trabalhos anteriores desenvolvidos no LMC (Laboratório de Mecânica Computacional) e, posteriormente, com a arquitetura desenvolvida, executar e investigar problemas de Dinâmica dos Fluidos e Interação Fluido-Estrutura (aplicando o método dos Elementos Finitos com Fronteiras Imersas e a solução direta do sistema de equações lineares com PARDISO) com o auxílio dos computadores do LMC. Uma análise de sensibilidade para cada problema é realizada de forma a encontrar a melhor combinação entre o número de elementos da malha de elementos finitos e o speedup, e posteriormente é feita uma análise comparativa de desempenho entre a arquitetura paralela a sequencial. Com uma única GPU conseguiu-se uma considerável redução no tempo para o assembly das matrizes globais e no tempo total da simulação. / One of the biggest challenges of engineering is enable computational solutions that reduce processing time and provide more accurate numerical solutions. Proposals with several approaches that explore new ways of solving such problems or improve existing solutions emerge. One of the biggest areas dedicated to propose such improvements is the parallel and high performance computing. Techniques that improve the processing time, more efficient algorithms and faster computers open up new horizons allowing to perform tasks that were previously unfeasible or would take too long to complete. We can point out, among several areas of interest, Fluid Dynamics and Interaction Fluid-Structure. In this work it is developed a parallel computing architecture in order to solve numerical problems more efficiently, compared to sequential architecture (e.g. Fluid Dynamics and Fluid-Structure Interaction problems) and it is also possible to extend this architecture to solve different problems (e.g. Structural problems). The objective is to develop an efficient computational algorithm in scientific programming language C ++, based on previous work carried out in Computational Mechanics Laboratory (CML) at Polytechnic School at University of São Paulo, and later with the developed architecture, execute and investigate Fluid Dynamics and Fluid-Structure Interaction problems with the aid of CML computers. A sensitivity analysis is executed for different problems in order to assess the best combination of elements quantity and speedup, and then a perfomance comparison. Using only one GPU, we could get a 10 times speedup compared to a sequential software, using the Finite Element with Immersed Boundary Method and a direct solver (PARDISO).