Dans cette thèse nous mettons en oeuvre une méthode vortex hybride pénalisée (HVP) afin desimuler des écoulements incompressibles autour de corps non profilés dans des milieux complexessolides-fluides-poreux. Avec cette approche particulaire hybride, le phénomène de convection estmodélisé à l’aide d’une méthode vortex afin de bénéficier du caractère peu diffusif et naturel desméthodes particulaires. Un remaillage des particules est alors réalisé systématiquement sur unegrille cartésienne sous-jacente afin d’éviter les phénomènes de distorsion. D’autre part, les effetsdiffusifs et d’étirement ainsi que le calcul de la vitesse sont traités sur la grille cartésienne, àl’aide de schémas eulériens. Le traitement des conditions de bords aux parois de l’obstacle esteffectué à l’aide d’une technique de pénalisation, particulièrement bien adaptée au traitementde milieux solides-fluides-poreux.Cette méthode HVP est appliquée au contrôle passif d’écoulement. Cette étude de contrôleest effectuée respectivement en 2D et en 3D autour d’un demi-cylindre et d’un hémisphère parl’ajout d’un revêtement poreux à la surface de l’obstacle. La présence de cette couche poreusemodifiant la nature des conditions aux interfaces, permet de régulariser l’écoulement global etde diminuer la traînée aérodynamique de l’obstacle contrôlé. A travers des études paramétriquessur la perméabilité, l’épaisseur et la position du revêtement poreux, ce travail vise à identifier desdispositifs de contrôles efficaces pour des écoulements autour d’obstacles comme des rétroviseursautomobiles. / In this work we use a hybrid vortex penalization method (HVP) to simulate incompressibleflows past bluff bodies in complex solid-fluid-porous media. In this hybrid particle approach,the advection phenomenon is modeled through a vortex method in order to benefit from thenatural description of the flow supplied by particle methods and their low numerical diffusionfeatures. A particle remeshing is performed systematically on an underlying Cartesian grid inorder to prevent distortion phenomena. On the other hand, the viscous and stretching effects aswell as the velocity calculation are discretized on the mesh through Eulerian schemes. Finally,the treatment of boundary conditions is handled with a penalization method that is well suitedfor the treatment of solid-fluid-porous media.The HVP method is applied to passive flow control. This flow control study is realized pasta 2D semi-circular cylinder and a 3D hemisphere by adding a porous layer on the surface of thebody. The presence of such porous layer modifies the characteristics of the conditions at theinterfaces and leads to a regularization of the wake and to a decrease of the aerodynamic dragof the controlled obstacle. Through parametric studies on the permeability, the thickness andthe position of the porous coating, this works aims to identify efficient control devices for flowsaround obstacles like the rear-view mirrors of a ground vehicle.

Dans cette thèse nous mettons en oeuvre une méthode vortex hybride pénalisée (HVP) afin desimuler des écoulements incompressibles autour de corps non profilés dans des milieux complexessolides-fluides-poreux. Avec cette approche particulaire hybride, le phénomène de convection estmodélisé à l’aide d’une méthode vortex afin de bénéficier du caractère peu diffusif et naturel desméthodes particulaires. Un remaillage des particules est alors réalisé systématiquement sur unegrille cartésienne sous-jacente afin d’éviter les phénomènes de distorsion. D’autre part, les effetsdiffusifs et d’étirement ainsi que le calcul de la vitesse sont traités sur la grille cartésienne, àl’aide de schémas eulériens. Le traitement des conditions de bords aux parois de l’obstacle esteffectué à l’aide d’une technique de pénalisation, particulièrement bien adaptée au traitementde milieux solides-fluides-poreux.Cette méthode HVP est appliquée au contrôle passif d’écoulement. Cette étude de contrôleest effectuée respectivement en 2D et en 3D autour d’un demi-cylindre et d’un hémisphère parl’ajout d’un revêtement poreux à la surface de l’obstacle. La présence de cette couche poreusemodifiant la nature des conditions aux interfaces, permet de régulariser l’écoulement global etde diminuer la traînée aérodynamique de l’obstacle contrôlé. A travers des études paramétriquessur la perméabilité, l’épaisseur et la position du revêtement poreux, ce travail vise à identifier desdispositifs de contrôles efficaces pour des écoulements autour d’obstacles comme des rétroviseursautomobiles. / In this work we use a hybrid vortex penalization method (HVP) to simulate incompressibleflows past bluff bodies in complex solid-fluid-porous media. In this hybrid particle approach,the advection phenomenon is modeled through a vortex method in order to benefit from thenatural description of the flow supplied by particle methods and their low numerical diffusionfeatures. A particle remeshing is performed systematically on an underlying Cartesian grid inorder to prevent distortion phenomena. On the other hand, the viscous and stretching effects aswell as the velocity calculation are discretized on the mesh through Eulerian schemes. Finally,the treatment of boundary conditions is handled with a penalization method that is well suitedfor the treatment of solid-fluid-porous media.The HVP method is applied to passive flow control. This flow control study is realized pasta 2D semi-circular cylinder and a 3D hemisphere by adding a porous layer on the surface of thebody. The presence of such porous layer modifies the characteristics of the conditions at theinterfaces and leads to a regularization of the wake and to a decrease of the aerodynamic dragof the controlled obstacle. Through parametric studies on the permeability, the thickness andthe position of the porous coating, this works aims to identify efficient control devices for flowsaround obstacles like the rear-view mirrors of a ground vehicle.

Ce travail est dédié au développement d'une méthode numérique pour la simulation des écoulements liquide-bulles. La présence des bulles dans l'écoulement visqueux et incompressible est prise en compte via une méthode de pénalisation. Dans cette représentation Euler-Lagrange, les bulles supposées indéformables et parfaitement sphériques sont assimilées à des objets pénalisés interagissant avec le fluide. Une méthode VOF (Volume Of Fluid) est employée pour le suivi de la fonction de phase. Une adaptation de la discrétisation des équations de Navier-Stokes est proposée afin d'imposer la condition de glissement à l'interface entre le liquide et les bulles. Une méthode de couplage entre le mouvement des bulles et l'action du liquide est proposée. La stratégie de validation est la suivante. Dans un premier temps, une série de cas-tests est proposée; les objets pénalisés sont supposés en non-interaction avec le fluide. L'étude permet d'exhiber la convergence et la précision de la méthode numérique. Dans un second temps le couplage est testé via deux types de configurations de validation. Le couplage est d'abord testé en configuration de bulle isolée, pour une bulle en ascension dans un liquide au repos pour les Reynolds Re=17 and Re=71. Les résultats sont comparés avec la théorie établie par la corrélation de Mei pour les bulles sphériques propres décrivant intégralement la dynamique de la bulle. Enfin, des simulations en configurations de nuage de bulles sont présentées, pour des populations mono- et bidisperses dans un domaine entièrement périodique pour des taux de vide s'établissant entre 1% et 15%. Les statistiques fournies par les simulations caractérisant l'agitation induite par les bulles sont comparées à des résultats expérimentaux. Pour les simulations de nuages de bulles bidisperses, de nouveaux résultats sont présentés. / This work is devoted to the development of a numerical method for the simulation of two-phase liquid-bubble flows. We use a volume penalization method to take into account bubbles in viscous incompressible flows. The chosen Euler-Lagrange framework involves spherical and nondeformable bubbles represented as moving penalized obstacles interacting with the fluid. A VOF (Volume Of Fluid) method is used to track the phase function while a discretization of the penalized conservation equations is realized to impose slip conditions at the liquid-bubble interface. A coupling method devised from the penalized momentum equations is proposed. The validation process is set as following. First, the fluid is supposed non-acting on the bubbles; several test-cases are presented; we consider configurations with different penalized obstacles shapes (curved channel, inclined channel), the obstacles are either static or dynamic; in each configuration an analytical solution is known. The results show the compliance and the quality of our numerical closures by exposing the convergence order of the method. In order to verify the accuracy of the coupling method, numerical simulations of a 1mm diameter single bubble rising in a quiescent liquid are performed for Re=17 and Re=71. Results are compared with theory established by using Mei correlation for clean spherical bubbles describing the whole dynamics of the rising bubble. Finally, simulations of bubble swarms, in mono- and bidisperse configurations have been carried out in a fully periodic box with moderate void fractions ranging from 1% to 15%. The statistics provided by the simulations characterizing the bubble-induced agitation are compared to experimental results. For the bidisperse bubble swarm configuration, new results are presented.

Les travaux présentés dans ce manuscrit se concentrent sur l'approximation numérique des équations de la magnétohydrodynamique (MHD) et sur leur stabilisation pour des problèmes caractérisés par des nombres de Reynolds cinétique élevés ou par des écoulements multiphasiques. Nous validons numériquement un nouveau modèle de Simulation des Grandes Echelles (ou Large Eddy Simulations, LES), dit de viscosité entropique, sur des écoulements de cylindre en précession ou créés par des turbines contra-rotatives (écoulement de Von Kármán). Ces études sont réalisées avec le code MHD SFEMaNS développé par J.-L. Guermond et C. Nore depuis 2002 pour des géométries axisymétriques. Ce code est basé sur une décomposition spectrale dans la direction azimutale et des éléments finis de Lagrange dans un plan méridien. Nous adaptons une méthode de pseudo-pénalisation pour prendre en compte des turbines en mouvement, ce qui étend le code SFEMaNS à des géométries quelconques. Nous présentons aussi une méthode originale d'approximation des équations de Navier-Stokes à densité variable qui utilise la quantité de mouvement comme variable et la viscosité entropique pour stabiliser les équations de la masse et du mouvement. / The investigations presented in this manuscript focus on the numerical approximation of the magnetohydrodynamics (MHD) equations and on their stabilization for problems involving either large kinetic Reynolds numbers or multiphase flows. We validate numerically a new Large Eddy Simulation (LES) model, called entropy viscosity, on flows driven by precessing cylindrical containers or counter-rotating impellers (Von Kármán flow). These studies are performed with SFEMaNS MHD-code developed by J.-L. Guermond and C. Nore since 2002 for axisymmetric geometries. This code is based on a spectral decomposition in the azimuthal direction and a Lagrange finite element approximation in a meridian plane. We adapt a pseudo-penalization method to report the action of rotating impellers that extends the range of SFEMaNS's applications to any geometry. We also present an original approximation method of the Navier-Stokes equations with variable density. This method uses the momentum as variable and stabilizes both mass and momentum equations with the same entropy viscosity.

La description magnétohydrodynamique est utilisée pour étudier les plasmas de fusion par confinement magnétique dans deux configurations: tokamak et reversed field pinch. Une méthode de Fourier pseudo-spectrale et une technique de pénalisation en volume sont employées pour résoudre les équations. La méthode de pénalisation permet d’introduire des conditions aux limites de Dirichlet et donc de faire varier facilement la géométrie considérée. Les simulations dans des géométries toroïdales de type tokamak montrent l’apparition spontanée de vitesses. Une importante composante toroïdale se développe si le système est peu dissipatif. Il est aussi montré que la brisure de symétrie dans la forme de la section du tore fait apparaitre un moment angulaire toroïdal. Pour le Reversed Field Pinch on montre l’émergence de structures hélicoïdales. La forme de ces structures varie en fonction des coefficients de transport ainsi que du paramètre de pincement du champ magnétique imposé. Pour compléter l’étude on compare les résultats du tore aux calculs dans un cylindre périodique. Les différences dans la dynamique des deux cas sont mises en avant. Finalement les simulations sont confrontées à des expériences et un meilleur accord est observé entre simulation et expérience pour la géométrie toroïdale que pour la géométrie cylindrique. / A magnetohydrodynamic description is used to study magnetic fusion plasmas in two different configurations: tokamak and reversed field pinch. A Fourier pseudo-spectral method with a volume penalization technique are used to solve the system of equations. The penalization method is used to introduce Dirichlet boundary conditions and it al- lows to easily modify the consider geometry. The simulations of a tokamak configuration in a toroidal geometry show the spontaneous appearance of velocities. These velocities are dominated by their toroidal component if the system is little dissipative. It is also shown that the symmetry breaking of the cross section of the torus causes a toroidal angular momentum to develop. For the Reversed Field Pinch configuration we show the appearance of helical structures. The shape of these structures varies with the value of the transport coefficients and with the pinch ratio parameter of the imposed magnetic field. To complete the study, we compare the results of simulations obtained in toroidal and in periodic cylindrical geometries. The differences in the dynamics of these two cases are highlighted. Finally, simulations are compared to experimental data and a significant better agreement is observed between the simulation and the experiment for the toroidal geometry than for the cylindrical case.

Les polymères jouent un rôle important sur la réduction de la traînée et la modification de la structure des écoulements. Nous avons utilisé le modèle Oldroyd-B pour étudier l’effet de la viscoélasticité de solutions de polymères dilués sur des écoulements en deux dimensions autour d’obstacles dans un canal. Les obstacles sont pris en compte par la méthode de pénalisation volumique et des condition aux limites artificielles sans réflexion sont imposées à la sortie du canal. La discrétisation est effectuée par des schémas performants en différences finies et la résolution par une méthode multigrille. Les simulations numériques sont effectuées pour une large gamme de nombres de Reynolds et de nombres de Weissenberg. Les caractéristiques détaillées des écoulements viscoélastiques sont analysées et comparées entre elles et celles de l’écoulement du fluide sans polymère. En particulier les jeux de paramètres conduisant à une augmentation ou à une baisse de la traînée. Enfin les effets du polymère sur des écoulements turbulents sont aussi analysés. / Polymer plays an important role on the drag reduction and modification of the structure of flow. In this thesis, Oldroyd B model is employed to study the effectof viscoelasticity for polymer solutions diluted by two-dimensional direct numerical simulation. The obstacles are taken into account by penalization method. The artificial boundary condition is imposed without any reflection on the channel outlet. Flow pasta cylinder is investigated in detailed by present numerical methods. The numerical codes are valid for predicting the drag force and capturing the important character of viscoelastic flow by comparing with experimental and other numerical results. The drag map of the cylinder is obtained at a wide range of Reynolds number and Weissenberg number space. The detailed characteristics of viscoelastic flow are reported in the thesis. The effects of polymer on two-dimensional turbulent flow are also discussed by grid turbulent flow.

Depuis la publication de l'ouvrage de Choquet (1955), la théorie d'espérance non linéaire a attiré avec grand intérêt des chercheurs pour ses applications potentielles dans les problèmes d'incertitude, les mesures de risque et le super-hedging en finance. Shige Peng a construit une sorte d'espérance entièrement non linéaire dynamiquement cohérente par l'approche des EDP. Un cas important d'espérance non linéaire cohérente en temps est la G-espérance, dans laquelle le processus canonique correspondant (B_{t})_{t≥0} est appelé G-mouvement brownien et joue un rôle analogue au processus de Wiener classique. L'objectif de cette thèse est d'étudier, dans le cadre de la G-espérance, certaines équations différentielles stochastiques rétrogrades (G-EDSR) à croissance quadratique avec applications aux problèmes de maximisation d'utilité robuste avec incertitude sur les modèles, certaines équations différentielles stochastiques (G-EDS) réfléchies et équations différentielles stochastiques rétrogrades réfléchies avec générateurs lipschitziens. On considère d'abord des G-EDSRs à croissance quadratique. Dans le Chapitre 2 nous fournissons un resultat d'existence et unicité pour des G-EDSRs à croissance quadratique. D'une part, nous établissons des estimations a priori en appliquant le théorème de type Girsanov, d'où l'on en déduit l'unicité. D'autre part, pour prouver l'existence de solutions, nous avons d'abord construit des solutions pour des G-EDSRs discretes en résolvant des EDPs non-linéaires correspondantes, puis des solutions pour les G-EDSRs quadratiques générales dans les espaces de Banach. Dans le Chapitre 3 nous appliquons les G-EDSRs quadratiques aux problèmes de maximisation d'utilité robuste. Nous donnons une caratérisation de la fonction valeur et une stratégie optimale pour les fonctions d'utilité exponentielle, puissance et logarithmique. Dans le Chapitre 4, nous traitons des G-EDSs réfléchies multidimensionnelles. Nous examinons d'abord la méthode de pénalisation pour résoudre des problèmes de Skorokhod déterministes dans des domaines non convexes et établissons des estimations pour des fonctions α-Hölder continues. A l'aide de ces résultats obtenus pour des problèmes déterministes, nous définissons le G-mouvement Brownien réfléchi et prouvons son existence et son unicité dans un espace de Banach. Ensuite, nous prouvons l'existence et l'unicité de solution pour les G-EDSRs multidimensionnelles réfléchies via un argument de point fixe. Dans le Chapitre 5, nous étudions l'existence et l'unicité pour les équations différentielles stochastiques rétrogrades réfléchies dirigées par un G-mouvement brownien lorsque la barrière S est un processus de G-Itô. / Since the publication of Choquet's (1955) book, the theory of nonlinear expectation has attracted great interest from researchers for its potential applications in uncertainty problems, risk measures and super-hedging in finance. Shige Peng has constructed a kind of fully nonlinear expectation dynamically coherent by the PDE approach. An important case of time-consistent nonlinear expectation is G-expectation, in which the corresponding canonical process (B_{t})_{t≥0} is called G-Brownian motion and plays a similar role to the classical Wiener process. The objective of this thesis is to study, in the framework of the G-expectation, some backward stochastic differential equations (G-BSDE) under a quadratic growth condition on their coefficients with applications to robust utility maximization problems with uncertainty on models, Reflected stochastic differential equations (reflected G-SDE) and reflected backward stochastic differential equations with Lipschitz coefficients (reflected G-BSDE). We first consider G-BSDE with quadratic growth. In Chapter 2 we provide a result of existence and uniqueness for quadratic G-BSDEs. On the one hand, we establish a priori estimates by applying the Girsanov-type theorem, from which we deduce the uniqueness. On the other hand, to prove the existence of solutions, we first constructed solutions for discrete G-BSDEs by solving corresponding nonlinear PDEs, then solutions for the general quadratic G-BSDEs in the spaces of Banach. In Chapter 3 we apply quadratic G-BSDE to robust utility maximization problems. We give a characterization of the value function and an optimal strategy for exponential, power and logarithmic utility functions. In Chapter 4, we discuss multidimensional reflected G-SDE. We first examine the penalization method to solve deterministic Skorokhod problems in non-convex domains and establish estimates for continuous α-Hölder functions. Using these results for deterministic problems, we define the reflected G-Brownian motion and prove its existence and its uniqueness in a Banach space. Then we prove the existence and uniqueness of the solution for the multidimensional reflected G-SDE via a fixed point argument. In Chapter 5, we study the existence and uniqueness of the reflected backward stochastic differential equations driven by a G-Brownian motion when the obstacle S is a G-Itô process.