Ce travail de thèse est une contribution à l’analyse de bifurcation des écoulements fluides avec prise en compte des interactions fluide-structure. Les phénomènes d’instabilité en interaction fluide-structure apparaissent dans de nombreux domaines de la vie courante ou industriels comme, par exemple : le flottement d’un drapeau dans le vent ou bien l’écoulement au sein d’échangeurs thermiques sur les sites de production d’énergie, l’écoulement autour des câbles sous-marins pour l’extraction de matières premières ou la fixation des plateformes off-shore, l’écoulement autour des structures aéronautiques ou navales. Dans ces situations, un phénomène complexe de vibration des structures induite par vortex peut se produire. L’objectif de la thèse est de proposer un algorithme permettant l’analyse de stabilité de tels systèmes. Ainsi, le couplage original d’une méthode de perturbation d’ordre élevé (Méthode Asymptotique Numérique - MAN) à une discrétisation spatiale permettant la prise en compte des interactions fluide-structure est proposée. À cet effet, une description purement eulérienne du mouvement est retenue. L’interaction fluide- structure est décrite au moyen d’une méthode de frontières immergées (MFI) à forçage continu (méthode de pénalisation) et discret (méthode Ghost-Cell). La présence d’obstacles au sein de l’écoulement est obtenue au moyen de la méthode de Level-Set. En complément, un intégrateur temporel des équations du mouvement associant la MAN, la MFI et une technique d’homotopie est proposé. L’ensemble de ces algorithmes est appliqué à des problèmes d’écoulement incompressible, à faible nombre de Reynolds, d’un fluide visqueux newtonien en présence d’obstacles solides rigides (fixes ou mobiles). L’analyse de stabilité d’un écoulement dans une conduite avec expansion/contraction soudaine (bifurcation stationnaire), et autour d’un cylindre (bifurcation de Hopf) est traitée. L’analyse transitoire d’un écoulement autour d’un cylindre rigide et mobile est également proposée. Les résultats obtenus permettent d’évaluer la précision et la performance des algorithmes proposés. Ainsi, les résultats de cette thèse permettent de conclure sur le bien-fondé de l’approche et constituent une première étape vers l’analyse de stabilité d’écoulements en présence de structures complexes, représentatifs de situations réelles / This thesis is a first contribution to the bifurcation analysis of fluid flows by taking into account fluid-structure interactions. Instability with fluid-structure interactions appears in many areas of everyday life or industry such as, for example: flag floating in the wind, flow within heat exchangers for energy production, flow around submarine cables for the extraction of raw materials or the fixing of off-shore platforms, flow around aeronautical or naval structures. In these situations, complex vortex-induced vibrations of the structures can occur. The aim of the thesis is to propose an algorithm allowing stability analysis of such systems. Thus, an original coupling of a high order perturbation method (Asymptotic Numerical Method - ANM) to a spatial discretization which takes into account fluid-structure interactions is proposed. For this purpose, a purely Eulerian description of the motion is retained. Fluid-structure interaction is described using an immersed boundary method (IBM) with continuous forcing (penalization method) and discrete (Ghost-Cell method) forcing. The presence of bodies within the flow is obtained by means of the Level-Set method. In addition, a time integrator of the governing equations associating ANM, IBM and homotopy technique is proposed. All these algorithms are applied to analyse incompressible flows, at low Reynolds number, of a Newtonian viscous fluid in the presence of rigid solids (fixed or moving). Bifurcation analysis of flows in a channel with sudden expansion / contraction (stationary bifurcation), or around a cylinder (Hopf bifurcation) are carried out. Transient analysis of a flow around a moving rigid cylinder is also proposed. Our results make it possible to evaluate accuracy and performance of the proposed algorithms. Thus, thesis results allow to conclude on the validity of the proposed approach. Finally, this thesis work constitutes a first step towards flow stability analysis in the presence of complex structures, representative of real situations.

Au cours de cette thèse, nous nous sommes intéressés à la modélisation des interactions fluide-structure entre un fluide visqueux, incompressible et une structure pouvant être déformable. Après avoir présenté les différentes approches possibles de modélisation, nous introduisons une nouvelle méthode de type frontière immergée : la méthode IPC ("Image Point Correction"). Combinant approches Ghost-Cell et Pénalisation, cette méthode mixte du second degré globalement et localement en vitesse, est validée sur différents cas tests (comparaisons des coefficients aérodynamiques pour des cylindres fixes ou mobiles, sédimentation 2D d'un cylindre). Nous avons ensuite appliqué la méthode IPC à la simulation de la nage. Dans un premier temps, le solveur 2D a été couplé avec un algorithme d'optimisation mathématique afin de déterminer la loi de nage optimale pour une géométrie de poisson donnée. Puis, dans un second temps, nous avons simulé la nage 3D après reconstruction approchée de la géométrie, basée sur des images du nageur. Enfin, grâce à l'outil du squelette, une reconstruction réaliste du poisson est proposée. / The aim of this thesis is to investigate the modeling of fluid-structure interactions. The fluid is viscous and incompressible, and the structure is subject to an imposed deformation. After a survey of the different existing approachs to model fluid-structure interactions, we introduce a new immersed boundary method: the IPC method (”Image Point Correction”). This ap-proach merges Ghost-Cell and Penalty concepts. It is globally and locally second order in velocity, and it is validated through several canonical simulations. Then, we apply the IPC method to fish-like swimming. First, the 2D solver is paired up with a mathematical optimization algorithm to determine the optimal swimming law for a given fish geometry. Secondly, we simulate a 3D swimmer after performing an approximated reconstruction of the geometry based on actual fish pictures. Finally, thanks to the skeleton approach, a realistic reconstruction of the fish is exposed.

This study presents an improved ghost-cell immersed boundary approach to represent a solid body in compressible flow simulations. In contrast to the commonly used approaches, in the present work ghost cells are mirrored through the boundary described using a level-set method to farther image points, incorporating a higher-order extra/interpolation scheme for the ghost cell values. In addition, a shock sensor is in- troduced to deal with image points near the discontinuities in the flow field. Adaptive mesh refinement (AMR) is used to improve the representation of the geometry efficiently. The improved ghost-cell method is validated against five test cases: (a) double Mach reflections on a ramp, (b) supersonic flows in a wind tunnel with a forward- facing step, (c) supersonic flows over a circular cylinder, (d) smooth Prandtl-Meyer expansion flows, and (e) steady shock-induced combustion over a wedge. It is demonstrated that the improved ghost-cell method can reach the accuracy of second order in L1 norm and higher than first order in L∞ norm. Direct comparisons against the cut-cell method demonstrate that the improved ghost-cell method is almost equally accurate with better efficiency for boundary representation in high-fidelity compressible flow simulations. Implementation of the improved ghost-cell method in reacting Euler flows further validates its general applicability for compressible flow simulations.

Ces travaux portent sur la modélisation et la simulation numérique des effets du déplacement d'un corps indéformable dans un écoulement multiphasique compressible. Ils se placent dans le cas où plusieurs objets sont en mouvement ou dans le cas où un objet est en mouvement dans un milieu aux géométries complexes. L'étude ne peut alors pas être placée dans le référentiel lié à l'objet en mouvement. Le modèle est basé sur une méthode multiphasique à interfaces diffuses où les différentes phases sont en équilibre mécanique. Le système régissant l'écoulement fluide est augmenté d'une équation d'advection. Cette dernière s'applique sur une fonction Level Set dont le niveau zéro permet de localiser le mobile dans l'espace. Des termes de couplage sont ajoutés au membre de droite des équations d'évolution de la quantité de mouvement et de l'énergie totale. Ces termes sont composés d'un facteur du type pénalisation et d'un facteur du type relaxation de vitesses. Cette nouvelle méthode permet de simuler des cas complexes où peuvent interagir des mobiles à hautes vitesses, des ondes de choc et des interfaces liquide/gaz. / This work deals with modelling and the numerical simulation of the effects of a moving rigid body on a multiphase flow. Here more than one object is moving, or an object is moving in a complex geometry domain. So the reference frame linked to the moving body can't be used. The model is build on a multiphase diffuse interface method with mechanical equilibrium. An advection equation is added. It applies on a Level Set function used to track the moving body. Coupling terms are added to the momentum equation and to the total energy equation. These terms are made of a penalization factor and a velocity relaxation factor. This new method allows to simulate complex cases where can interact high velocity objects, shock waves and liquid / gas interfaces.

This research is mainly focused on the development of the adaptive Cartesian grid methods for compressibl  e flow. At first, the ghost cell method and its applications for inviscid compressible flow on adaptive tree Cartesian grid are developed. The proposed method is successfully used to evaluate various inviscid compressible flows around complex bodies. The mass conservation of the method is also studied by numerical analysis. The extension to three-dimensional flow is presented. Then, an h-adaptive Runge–Kutta discontinuous Galerkin (RKDG) method is presented in detail for the development of high accuracy numerical method under Cartesian grid. This method combined with the ghost cell immersed boundary method is also validated by well documented test problems involving both steady and unsteady compressible flows over complex bodies in a wide range of Mach numbers. In addition, in order to suppress the failure of preserving positivity of density or pressure, which may cause blow-ups of the high order numerical algorithms, a positivity-preserving limiter technique coupled with h-adaptive RKDG method is developed. Such a method has been successfully implemented to study flows with the large Mach number, strong shock/obstacle interactions and shock diffraction. The extension of the method to viscous flow under the adaptive Cartesian grid with hybrid overlapping bodyfitted grid is developed. The method is validated by benchmark problems and has been successfully implemented to study airfoil with ice accretion. Finally, based on an open source code, the detached eddy simulation (DES) is developed for massive separation flow, and it is used to perform the research on aerodynamic performance analysis over the wing with ice accretion.

Despite the high cost of memory and CPU time required to resolve the boundary layer, a viscous unstructured grid solver has many advantages over a structured grid solver such as the convenience in automated grid generation and shock or vortex capturing by solution adaption. Since the geometry and flow phenomenon of a helicopter are very complex, unstructured grid-based methods are well-suited to model properly the rotor-fuselage interaction than the structured grid solver. In present study, an unstructured Cartesian grid solver is developed on the basis of the existing solver, NASCART-GT. Instead of cut-cell approach, immersed boundary approach is applied with ghost cell boundary condition, which increases the accuracy and minimizes unphysical fluctuations of the flow properties. The standard k-epsilon model by Launder and Spalding is employed for the turbulence modeling, and a new wall function approach is devised for the unstructured Cartesian grid solver. It is quite challenging and has never done before to apply wall function approach to immersed Cartesian grid. The difficulty lies in the inability to acquire smooth variation of y+ in the desired range due to the non-body-fitted cells near the solid wall. The wall function boundary condition developed in this work yields stable and reasonable solution within the accuracy of the turbulence model. The grid efficiency is also improved with respect to the conventional method. The turbulence modeling is validated and the efficiency of the developed boundary condition is tested in 2-D flow field around a flat plate, NACA0012 airfoil, axisymmetric hemispheroid, and rotorcraft applications. For rotor modeling, an actuator disk model is chosen, since it is efficient and is widely verified in the study of the rotor-fuselage interaction. This model considers the rotor as an infinitely thin disk, which carries pressure jump across the disk and allows flow to pass through it. The full three dimensional calculations of Euler and RANS equations are performed for the GT rotor model and ROBIN configuration to test implemented actuator disk model along with the developed turbulence modeling. Finally, the characteristics of the rotor-fuselage interaction are investigated by comparing the numerical solutions with the experiments.

Liquid sloshing with unrestrained free surface in a moving container is associated with various engineering problems, such as tankers on highways, liquid oscillations in large storage tanks caused by earthquakes, sloshing of liquid cargo in ocean-going vessels, and the motion of liquid fuel in aircraft and spacecraft. The purpose of this study is to develop a three-dimensional (3D) numerical wave tank with or without internal structures to investigate the mechanism of liquid sloshing and the interaction between the fluid and internal structures. The developed 3D time-independent finite difference method is applied on solving liquid sloshing in tanks with or without the influence of baffles under the ground motion of six-degrees of freedom. The 3D Navier-Stokes equations were solved and transformed to a tank-fixed coordinate system, and the fully nonlinear kinematic and dynamic free surface boundary conditions for fluid sloshing in a rectangular tank with a square base were considered. The fluid is assumed incompressible in this study. The complicated interaction in the vicinity of the fluid-structure interface was solved by implementing one dimensional ghost cell approach and the stretching grid technique near the fluid-structure boundaries were used to catch the detailed evolution of local flow field. A PC-cluster was established by linking several single computers to reduce the computational times due to the implementation of the 3D numerical model. The Message Passing Interface (MPI) parallel language and MPICH2 software were utilized to code the computer codes and to carry out the circumstance of parallel computation, respectively. The developed numerical scheme was verified by rigorous benchmark tests. Not only the reported analytical, numerical and experimental studies were compared with the present numerical results, the experimental investigation was also involved in the present work to further validate the accuracy of the numerical scheme. All the benchmark tests of this study showed excellent accuracy of the developed numerical scheme. For a tank without internal structures, the coupled motions of surge and sway are simulated with various excitation angles, excitation frequencies and water depths. The characteristics of sloshing waves are dissected in terms of the classification of sloshing wave types, sloshing amplitude, beating phenomenon, sloshing-induced forces and energy transfer of sloshing waves. Six types of sloshing waves, named single-directional, diagonal, square-like, swirling-like, swirling and irregular waves, were found and classified in the present study and the occurrence of these waves are tightly in connection with the excitation frequency of the tank. The effect of excitation angle on the characteristics of sloshing waves is explored and discussed, especially for swirling waves. The spectral analyses of sloshing displacement of various sloshing waves are examined and a clear evidence of the correlation between sloshing wave patterns and resonant modes of sloshing waves are demonstrated. The mechanism of switching direction of swirling waves is discussed by investigating the situation of circulatory flow, the instantaneous free surface, the gravitational effect and the instantaneous direction of external forcing. The coupling effects of heave, surge and sway motions were also included in this study and the result showed an unstable influence of heave motion on the kinematic and dynamic characteristics of sloshing waves when the vertical excitation frequency of the tank is twice as large as the fundamental natural frequency. Except irregular waves, the other types of sloshing waves are converted into swirling waves due to the effect of heave motion. The study related to tuned liquid damper (TLD) in 2D and 3D tanks were considered. A comprehensive investigation for a 2D tank with vertically tank bottom-mounted baffles (baffled tank) are demonstrated and discussed with respect to the influence of baffle height on the natural mode of the tank, the evolution of vortices and vortex shedding phenomenon, the relationship between the vortex shedding frequency and the excitation frequency of the tank, the vortex size generated in the vicinity of the baffle tip, the interaction of vortices inside the tank. The baffle height shows a significant influence on the shift of the first natural frequency of the baffled tank and the liquid depth also plays an important part in determining this influence. In other words, the shift of the first natural mode due to various baffle height is varied with water depths. The design of two baffles separated by 0.2 times the tank breadth is an efficient tool to not only reduce the sloshing amplitude but switch the first natural frequency of the tank. The sloshing displacement is affected distinctly by different numbers of baffles mounted vertically on the tank bottom. The more baffles mounted onto the tank bottom, the smaller the sloshing displacement is presented in both the transient and steady-state periods. The processes of the evolution of vortices near the baffle tip are categorized into four phases: the formation of separated shear layer and generation of vortices, the formation of a vertical jet and shedding of vortices, the interaction between shedding vortices and sloshing flow (the generation of snaky flow) and the interaction between snaky flow and sloshing waves. Vortex shedding phenomenon due to stronger vertical jets occurs when the excitation frequency is close to the first natural mode of the baffled tank. The size of the vortex generated near the baffle tip is discussed and the vortex size is closely correlated with the baffle height. Two types of 3D tuned liquid dampers, a vertically tank bottom-mounted baffle and a vertical plate, are discussed for a tank under coupled surge-sway motions. The wave types of diagonal and single-directional waves switch to the swirling type due to the influence of the baffle. The phenomenon of square-like waves or irregular waves coexisting with swirling waves is found in the baffled tank under diagonal excitation. The baffle and the vertical plate mounted parallel to the east (west) wall of the tank can effectively reduce the sloshing amplitude when the excitation angle is between 0 degree and 10 degree and the corresponding sloshing displacement in the sway (z) direction becomes more dominant with the increase of the excitation angle. The shift of the first natural mode of the baffled tank due to various baffle heights in the x direction is dominated in this design of baffled tank. The length of the plate can cause a significant influence on not only the variation of the natural frequencies but the type of the sloshing waves. The influence of the vertical plate on the irregular waves is insignificant and several peaks appear in the spectral analysis of the sloshing displacement for the irregular waves and the numbers of peaks are more than that of the baffled tank.

Nous nous interessons à la modélisation des interactions fluide-structure entre un fluide visqueux incompressible et une structure pouvant être déformable. Apres une approche des méthodes de type frontière immergée existantes, nous présentons une nouvelle approche : la méthode IPC (Image Point Correction) que nous validons ensuite sur différents cas tests. Puis, nous l'appliquons à la simulation 2D puis 3D de la nage d'un poisson grâce à une reconstruction utilisant l'outil du squelette.