Le champ de vitesse et de pression le long des cordes vibrantes était étudié à l'aide d'un modèle numérique et par observation expérimentale. Le modèle mathématique, basé sur les équations de Navier-Stokes 2D incompressibles discrétisées avec la méthode des éléments finis en formulation d'Euler-Lagrange arbitraire, était programmé dans le langage Fortran, à l'aide de la librairie numérique Mélina. Dans les résultats des simulations numériques on peut observer le développement du jet derrière la glotte et l'évolution des tourbillons de recirculation lors du cycle d'oscillation des cordes vocales. Le modèle physique en échelle 4: 1 fournit les signaux d'accélération, pression supraglottale et le signal acoustique émis par un corps élastique, qui vibre dans le paroi du canal aérodynamique grâce au couplage avec l'écoulement. Sur le même modèle, le champ de vitesse dans le domaine supraglottale était mesuré par la PIV synchronisé avec les oscillations.

Dans la première partie, je décris mes travaux de recherche en interaction fluide-structure (résultat théorique, numérique et de contrôle), puis je présente des méthodes fréquentielles pour la contrôlabilité et la stabilisation avec applications aux problèmes de stabilisation de systèmes discrétisés. Finalement, je donne des résultats obtenus sur la dépendance par rapport au domaine de quelques EDP

Cette thèse s'intéresse au recalage de modèles dynamiques, processus d'optimisation au cours duquel on corrige des paramètres du modèle de manière à minimiser une distance calcul-essai. A chaque itération, de nombreuses quantités doivent être déterminées (modes de la structure, sensibilités, localisation d'erreur, etc.). Ainsi le recalage fait intervenir des calculs particulièrement coûteux. C'est pourquoi on développe dans un premier temps une méthode de résolution approchée, applicable à divers problèmes rencontrés en analyse modale et plus particulièrement dans le recalage. Cette méthode, dénommée "méthode d'itération sur les résidus", est ainsi appliquée au calcul de modes (de structures élastiques ou amorties, couplées ou non avec un fluide), à la résolution d'une méthode d'expansion ou encore au calcul de sensibilités. Des comparaisons avec des méthodes existantes montrent, sur des modèles industriels, des performances extrêmement intéressantes tout en gardant une méthodologie unifiée, simple de mise en oeuvre. Dans un deuxième temps, on expose les divers outils utilisés pour le recalage. Sur ce sujet, les contributions principales de cette thèse sont l'introduction d'une démarche systématique pour la détermination de la visibilité d'erreurs, ainsi qu'une méthode de placement de capteurs pour localiser un défaut donné. Finalement, ces outils sont appliqués au recalage des modèles de l'étage principal cryotechnique et de l'étage supérieur cryotechnique d'Ariane 5.

L'athérosclérose est une maladie inflammatoire qui commence quand les lipoprotéines de faible densité (LDL) entrent dans l'intima du vaisseau sanguin où elles sont oxydées (ox-LDL). Le ox-LDL est considéré comme un agent dangereux par le système immunitaire provoquant ainsi une réponse immunitaire. Cette réponse immunitaire déclenche le recrutement des monocytes dans l'intima où elles se transforment en macrophages et ensuite en cellules spumeuses. Ce dernier amplifie la production des cytokines inflammatoires et davantage de recrutement des monocytes. Ce processus auto-amplifié est compensé par la sécrétion de cytokines anti-inflammatoires (anti-inflammation biochimique) et la migration des cellules musculaires lisses pour former une chape fibreuse qui couvre le noyau lipidique. Cette chape fibreuse avec le noyau lipidique s'appellent la plaque d'athérosclérose. Celle-ci change la géométrie du vaisseau sanguin en le rétrécissant et interagit avec du flux sanguin. Cette interaction peut avoir des conséquences dangereuses liées à la rupture de plaque ou à la formation du caillot de sang. La thèse est consacrée à la modélisation mathématique de ces phénomènes. Elle est composée de deux parties : Nous développons des modèles mathématiques basés sur des équations de réaction diffusion afin de décrire le processus inflammatoire. Le premier modèle est unidimensionnel. Il nous permet d'expliquer comment le développement de l'athérosclérose dépend de la concentration en cholestérol (ox-LDL). Si cette concentration dans l'intima est basse, alors la maladie ne se développera pas. Les concentrations intermédiaires de ox-LDL peuvent mener au développement de la maladie dans certaines conditions. Nous montrons que l'inflammation se propage en front d'ondes de réaction-diffusion. Les concentrations élevées de ox-LDL engendre le développement de la maladie. Même une petite perturbation du cas non inflammatoire mène à une propagation d'ondes qui correspond à l'inflammation. Ensuite nous étudions un modèle bidimensionnel qui représente un système d'équations type réaction-diffusion sur une bande. La deuxième dimension correspond à la section transversale de l'intima et une condition aux limites non-linéaire décrit le recrutement des monocytes. Cette condition aux limites est une fonction des concentrations des cytokines. Nous démontrons l'existence des fronts de propagation d'onde et confirmons les résultats précédents qui montrent que l'athérosclérose se développe en tant qu'onde de réaction-diffusion. Les résultats théoriques des deux modèles sont confirmés par des simulations numériques qui montrent que le cas bidimensionnel converge vers le cas unidimensionnel quand l'épaisseur de l'intima tend vers zéro. Une fois la plaque se forme, elle interagit avec le flux sanguin engendrant de différentes conséquences mécaniques et biochimiques. Nous développons un modèle d'interaction fluide-structure. La plaque d'athérome composée d'un dépôt lipidique couvert par une chape fibreuse, les deux étant modélisés en tant que matériaux hyper-élastiques. Le sang est considéré comme un fluide non-Newtonien avec une viscosité variable modélisée selon la loi de Carreau. Les paramètres utilisés dans nos simulations sont tirés de données expérimentales mentionnées dans la littérature. Nous étudions les effets non-Newtoniens sur les recirculations du sang en aval de la plaque d'athérome et aussi sur les contraintes sur celle-ci. Les simulations montrent que le modèle Newtonien surestime les recirculations de manière significative par rapport au modèle non-Newtonien. Elles montrent aussi que le modèle Newtonien sous-estime légèrement les contraintes sur la plaque pour des taux de cisaillement usuels, mais cette sous-estimation devient importante pour des taux de cisaillement bas.

Habilitation à diriger des recherches : Sciences physiques : Lille 1 : 2001. / Bibliogr. f. 133-136.

Reproduction de : Thèse doctorat : Transferts thermiques : Reims : 2006. / Titre provenant de l'écran titre. Bibliogr. f. 179-182.

Thèse de doctorat : sciences. Mécanique : Ecully, Ecole centrale de Lyon : 2006. / 82 réf.

Thèse de doctorat : sciences. Mécanique : Ecully, Ecole centrale de Lyon : 2006. / Titre provenant de l'écran-titre. 82 réf.

L'objectif industriel de cette thèse est d'étudier l'impact de l'environnement moyennes et hautes fréquences sur les lanceurs de missiles et de fusées. Ce champ acoustique est aléatoire et large bande et génère des vibrations à l'intérieur des lanceurs, qui peuvent détériorer les structures et les équipements fragiles, notamment les composants électroniques des équipements embarqués (premiers modes de vibration à 200-500 Hz, situés dans la zone des moyennes fréquences des lanceurs). Il est nécessaire de réaliser un modèle numérique du lanceur, afin de déterminer les niveaux de vibrations se propageant dans les structures internes des lanceurs sur l'ensemble du spectre fréquentiel. Dans le domaine des moyennes fréquences, les méthodes par éléments-finis ou énergétiques -- comme la SEA (Analyse Energtique Statistique) -- ne sont ni efficaces ni précises. Une nouvelle méthode est alors introduite : la méthode SEA-like. Cette méthode est déterministe et énergétique - pour prendre en compte le caractère aléatoire de l'excitation et les incertitudes des paramètres de la structure - mais s'appuie également sur des analyses basses fréquences, comme l'analyse modale ou par élément fini. Cette méthode est donc valable sur des domaines fréquentielles o'u la méthode SEA ne peut plus être utilisée. Dans une première partie, un système couplé plaque-cavité soumis à une excitation rain-on-theroof est étudié. Ce système a été modélisé avec trois méthodes : la méthode SEA-like basée sur une méthode par couplage modal - qui permet l'étude des échanges énergétiques en interaction fluide-structure dans le domaine des moyennes fréquences - la méthode SEA - qui permet de valider le modèle SEA-like dans le domaine des hautes fréquences, - et enfin une méthode SEA-like basée sur une modélisation par éléments-finis du système couplé et qui permet la validation dans le domaine des moyennes fréquences. Dans une deuxième partie, la méthode SEA-like est étendue à des cas de charges aéroacoustiques. En effet, deux types d'excitations sont fréquemment rencontrées dans le domaine de l'aérospatial et notamment dans l'étude des lanceurs : l'excitation par bruit diffus, qui intervient lors de la première phase de décollage du lanceur et une excitation par couche limite turbulente, qui intervient lors du vol atmosphérique. Ces excitations sont appliqués au système couplé plaque-cavité et modélisés avec la méthode SEA-like. Pour ce faire, les différentes excitations aéroacoustiques sont assimilées à des excitations delta-corrélées équivalentes. Les caractéristiques de l'excitation apparaissent alors dans le vecteur de puissance injectée et la modélisation SEA-like du système proprement dit est indépendante du type d'excitation. La méthode SEA-like est ainsi directement généralisable à différentes excitations. Ces modèles sont validés par comparaison avec des modèles connus d'excitations aéroacoustiques.

Cette maîtrise s'inscrit dans le projet Tr-Francis du laboratoire Heki de l'Université Laval. Ce mémoire présente l'exploration de stratégies de simulation fluide-structure bidirectionnelle à l'aide du code Star-CCM+ dans le but d'étudier les résonances de turbines hydrauliques. Concrètement, le travail vise à déterminer si la stratégie disponible dans Star-CCM+ permet de reproduire la physique d'une résonance liée à l'interaction rotor-stator d'une géométrie tournante à symétrie cyclique, analogue à celle d'une turbine Francis. Pour ce faire, une géométrie plus simple, mais semblable à une aube de turbine est d'abord étudiée, un hydrofoil. Une simulation FSI bidirectionnelle est produite sur l'hydrofoil F99 en condition de synchronisation de la réponse structurelle avec le détachement tourbillonnaire. Cette première simulation permet de déterminer les paramètres généraux à utiliser. Cette simulation ne parvient pas à reproduire les résultats expérimentaux obtenus à l'Université de science et technologie de Norvège, principalement en ce qui concerne la fréquence d'émission des tourbillons. Étant donné que l'accroche n'est pas atteinte, la résonance n'est pas observée. Afin d'assurer que la stratégie est apte à reproduire la physique de la résonance, le module d'élasticité est modifié afin de faire concorder la fréquence naturelle de la structure à la fréquence des tourbillons. L'accroche des tourbillons est alors reproduite avec cohérence en ce qui concerne la déformation de la structure et de l'écoulement. La stratégie est ensuite adaptée afin qu'elle soit utilisée sur le cas test RSI conçu par Dussault au laboratoire Heki. Les modifications sont majoritairement liées au domaine tournant nécessaire pour simuler la rotation de la roue du cas test. Les conditions de simulation sont choisies afin d'obtenir une résonance du mode à deux diamètres nodaux de la roue. Ici, l'apparition d'éléments à volume négatif dans le maillage empêche la simulation d'atteindre plus de 1.3 tour de roue simulé. À ce point, les fluctuations de pression liées à l'interaction rotor-stator sont clairement observables et une déformation liée au mode 2ND est obtenue. Toutefois, la déformation de la roue est dominée par la torsion. Il est difficile d'expliquer l'apparition de la torsion, la meilleure hypothèse portant sur les conditions initiales de la simulation. Également, certaines variations de pression incohérentes sont détectées, possiblement liées à des instabilités du solveur fluide. / This master project is part of Tr-Francis project at Heki Laboratory of Université Laval. This thesis presents the exploration of two-way fluid-structure simulation strategies with Star-CCM+ to study the resonances of hydraulic turbines. Concretely, the objective is to establish if the approach of Star- CCM+ can reproduce a resonance induced by rotor-stator interaction in the case of a simplified turbine model. To do so, a simpler case is studied, a hydrofoil. This geometry is similar to a runner blade of a hydraulic turbine without the complexity of the rotating domain. A two-way simulation is performed on the hydrofoil F99 in lock-in conditions. This first simulation is used to explore the general parameters of the software. This simulation fails to reproduce experimental data from the Norwegian University of Science and Technology, mainly regarding the vortex shedding frequency. The lock-in is not achieved and so is the resonance. To verify if the strategy can reproduce the resonance physics, the young modulus of the material is modified to match the natural frequency to the vortex shedding frequency simulated. The lock-in is now reproduced with consistent structural deformation and flow. Afterward, the strategy is adapted for the RSI test case developed by Dussault at Heki laboratory. Main modifications are related to the rotating domain used to simulate the runner rotation. Operating conditions are selected to obtain resonance of the two nodal diameters mode of the runner. The simulation fails to run more than 1.3 runner turns before the appearance of negative volume cells in the runner mesh. In the simulation, the RSI pressure fluctuations are identified, and a two nodal diameter mode is observed. However, the main runner deformation is in torsion. The main hypothesis to explain this behavior is related to the fluid initial conditions. Also, inconsistent pressure fluctuations are detected in the runner, possibly related to some solver instabilities.