Dans le but d'établir un modèle prédictif pour la croissance tumorale in-vivo et la thérapie, le modèle multi-échelle doit être élaboré et calibré par étape et de façon individuelle pour chaque type de cellule ciblé et pour di érents environnements (in-vitro et in-vivo). Nous présenterons, en tant que preuve de concept et à partir de di érentes sources de données, les étapes de la construction et de la paramétrisation du modèle de la croissance avasculaire des lignées de cellules EMT6/Ro et SK-MES-1. Dans une première étape, un modèle multi-échelle à base d'agents a été construit et validé avec des données provenant de la littérature sur les sphéroïdes multicellulaires de carcinomes mammaires de souris EMT6/Ro. Pour cette lignée de cellules, il a pu prédire que la cinétique de croissance est contrôlée par une combinaison de contraintes spatiales et de limitation des nutriments. Il a été trouvé que l'ATP est la ressource critique que les cellules essayent de garder constante en permutant d'un métabolisme aérobique à anaérobique et ce pour de larges plages de concentrations d'oxygène et de glucose. La saturation de la croissance a été observé uniquement dans le cas de faibles concentrations d'oxygène et de glucose ce que le modèle a pu expliqué par une migration guidée par l'adhésion de cellule à cellule. Dans une seconde étape, le modèle a été adapté à la lignée cellulaire SK-MES-1. Nous avons calibré la cinétique de croissance qualitativement en analysant des images de cryosections de sphéroïdes marquées pour l'apoptose et la prolifération et quantitativement en la comparant des courbes de croissance. Au delà de l'ATP, le lactate a été identi é comme contrôlant la taille du noyau nécrotique. Pour rendre compte de la situation in-vivo, nous proposons une extension du modèle qui prend en compte un réseau de vaisseaux sanguins et le phénomène de l'angiogenèse associé. A n de paramétrer les propriétés des vaisseaux fonctionnels et dans le but de valider les lois de l'angiogenèse, nous menons à partir d'images de perfusion d'agents de contraste une étude de sensibilité aux paramètres. Dans un premier temps, nous résolvons le problème direct de la perfusion des agents de contraste dans un réseau de vaisseaux perméables ou non. Ensuite, nous résolvons le problème inverse rigoureusement et, grâce à des comparaisons directes entre les paramètres originaux et ceux récupérés, nous étudions la capacité de prédiction du modèle dans di érents cas.

Dans les fours de raffinerie, l'efficacité du transfert énergétique vers le pétrole brut avant sa distillation est altérée par la formation d'un composé carboné dans les tubes, appelé coke. Cela conduit à l'augmentation des coûts de production et de maintenance, et exige une compréhension accrue ainsi qu'un meilleur contrôle de ce phénomène. Cet encrassement est de type chimique et induit par les fortes températures. Dans les fours de cette dimension, le transfert de chaleur s'effectue principalement par rayonnement des produits de combustion. Le flux radiatif net sur les surfaces d'échange des tubes dépend de la température de toutes les surfaces solides et a donc besoin d'être prédit avec une précision suffisante. La température sur les tubes est le résultat d'un équilibre entre le rayonnement thermique et la conduction. Le comportement thermique de l'ensemble du système est un problème de couplage entre le rayonnement et la conduction. Une méthodologie complète de couplage est exposée et validée de la manière suivante. Dans ce problème, la flamme est décrite par un modèle analytique axisymétrique avec chimie complexe. Le couplage avec la conduction dans les tubes est réalisé par l'utilisation d'une méthode aux ordonnées discrètes (DOM) avec un modèle spectral de type bandes étroites pour le rayonnement des gaz de combustion. Un bilan énergétique confirme que les transferts de chaleur sont dominés par le rayonnement thermique. Un bon accord avec les mesures disponibles sur un four réel montre que l'approche proposée est capable de prédire le rayonnement thermique. L'étape suivante consiste à coupler le calcul de la température du tube à une loi d'encrassement. Un modèle chimique simple est utilisé. Il est validé à l'aide d'une expérience de laboratoire. La comparaison entre les températures obtenues avec la simulation et celles mesurées par des sondes thermiques montre que la simulation est capable de capturer l'évolution de la température dans le tube avec précision. Enfin, un modèle d'encrassement pour la configuration réelle est trouvé puis appliqué dans une simulation couplée complète. Cette simulation montre un bon accord entre l'évolution de la température sur site et dans la simulation. Une analyse plus poussée est réalisée sur les profils de température, de flux radiatif et de dépôt de coke et montre l'impact de ce dépôt sur l'installation.

Ce travail est consacré à l'étude de structures cellulaires bidimensionnelles et notamment à leur évolution au cours du temps. Après une phase transitoire (dont la durée dépend de l'ordre initial de la structure) l'évolution atteint généralement un régime stationnaire (où l'aire moyenne des cellules varie linéairement avec le temps et où les propriétés sans dimension, telles que le désordre topologique et la distribution des nombres de côtés des cellules, sont invariantes). Une méthodologie d'analyse d'images a été mise au point pour caractériser les structures à un instant donné. Leurs propriétés métriques et topologiques sont déterminées en tenant compte d'une correction de biais statistique. Nous avons réalisé des expériences portant sur des mousses de savon bidimensionnelles (réalisées entre deux plaques). Un système de drainage a également été mis en place dans le but de conserver l'épaisseur des arêtes constante au cours de l'évolution. Lorsque l'état initial de la structure est très ordonné, le stade transitoire est très long, ce qui limite notre étude à cette seule phase de l'évolution. En revanche, nous avons pû étudier le régime stationnaire en partant de structures initiales désordonnées et en particulier déterminer la valeur du désordre topologique dans cette phase. Cette valeur semble dépendre de la composition du liquide moussant, de la taille de la boîte contenant la mousse, mais apparemment pas du drainage. Des simulations ont été effectuées à l'aide d'un programme développé par H. Telley à l'EPFL. Ce programme est fondé sur l'utilisation des complexes de Laguerre bidimensionnels et périodiques. Ces simulations ont fourni des résultats comparables à ceux observés pour les mousses, mais également pour les polycristaux, grâce à l'ajustement d'un paramètre distributif. Celui-ci est relié de façon simple aux transformations topologiques élémentaires intervenant au cours de l'évolution. La validation du programme a été effectuée non seulement pour le régime stationnaire mais aussi pour le régime transitoire (pour les mousses de savon)

Un polymère thermoplastique fondu glisse sur une surface solide lisse et non-adsorbante : la vitesse est non nulle à l'interface. Le glissement est réduit si l'on greffe des macromolécules à la surface. Dans ce travail théorique, nous modélisons cette réduction du glissement à l'aide de mécanismes moléculaires, donc à l'échelle nanométrique. Le comportement microscopique des molécules et la réponse rhéologique macroscopique de l'interface sont décrits. Les prédictions du modèle présenté et celles d'autres modèles moléculaires sont confrontées aux résultats des expériences menées par une équipe du lavoratoire. Un certain nombre des lois de comportement sont expliquées, mais la gamme explorée des paramètres ne permet pas de départager tous les modèles. Un test différent est proposé, utilisant des polymères en étoile. D'autres problèmes ont été abordés : conformation d'une macromolécule unique dans un fondu chimiquement différent, statique d'une brosse polymère et d'une étoile dans les mêmes conditions, pénétration partielle d'un fondu dans une brosse chimiquement identique très dense, généralisation aux polymères branchés statistiques d'une méthode de séparation de polymères en étoile en solution diluée, détermination du point de gel pour une structure constituée de polymères en anneau (gel "olympique"), dynamique d'étalement d'une goutte d'hélium superfluide analogue à celle d'une goutte de polymère liquide. Présentation (PDF p. 6, English p. 8) Partie I. Statique et dynamique de polymères liquides (PDF p. 10) Chapitre 1. Généralités sur les polymères (PDF p. 11, English p. 12, français p. 16) Chapitre 2. Polymères fondus : propriétés et mise en oeuvre (PDF p. 39, English p. 40, français p. 42) Chapitre 3. Polymères fondus : éléments de théorie (PDF p. 58, English p. 59, français p. 66) Partie II. Statique et dynamique de chaînes greffées (PDF p. 153) Chapitre 4. Statique de chaînes greffées (PDF p. 154, English p. 155, français p. 158) Chapitre 5. Glissement d'un polymère fondu sur une surface greffée (PDF p. 200, English p. 201, français p. 223) Partie III. Autres problèmes (PDF p. 313) Chapitre 6. Caractérisation de polymères branchés par perméation (PDF p. 314, français p. 315, article in English p. 318) Chapitre 7. Gel olympique (PDF p. 332, English p. 333, français p. 338) Chapitre 8. Étalement d'une goutte par effet Josephson (PDF p. 358, français p. 359, article in English p. 361) Conclusion (PDF français p. 370, English p. 372) Bibliographie (PDF p. 374) Table des matières / Table of contents (PDF p. 380)

Le bruit généré par l'impact du sillage de la soufflante avec les aubes du redresseur est une source prédominante des turboréacteurs en conditions d'approche, et la composante à large-bande contribue fortement au niveau sonore global. Une méthode numérique basée sur un code CAA résolvant les équations d'Euler est développée dans cette thèse afin d'estimer le bruit d'interaction entre un sillage turbulent et un aubage de géométrie quelconque. Le sillage amont est modélisé à l'aide d'une méthode stochastique supposant un spectre de turbulence homogène isotrope et une représentation spatiale simplifiée du champ de vitesse. Ces fluctuations de vitesse sont injectées dans le code CAA via des conditions aux limites adéquates. La méthode ainsi mise en place est validée dans un premier temps sur des cas d'interaction avec une plaque plane en comparant les résultats numériques aux solutions du modèle d'Amiet. Un chaînage avec une méthode intégrale est aussi réalisé pour estimer le rayonnement acoustique. La méthode numérique est ensuite mise en œuvre pour simuler la réponse acoustique d'un profil isolé avec un bord d'attaque ondulé conçu et testé dans la soufflerie de l'ISVR dans le cadre du projet européen FLOCON. Enfin, la méthodologie est étendue à des configurations de grilles d'aubes annulaires en conduit. Après une étape de validation sur des cas monofréquentiels issus d'un benchmark CAA de la NASA, des simulations large bande sont réalisées, tout d'abord sur une configuration de grille d'aubes sans incidence (écoulement axial uniforme) testée en soufflerie au LMFA, puis sur une configuration plus complexe de grille d'aubes non chargée en écoulement tournant proposée par Atassi.

Cette thèse est consacrée à la conception d'un contrôle en boucle fermée d'un écoulement de cavité subsonique. L'objectif est de réaliser un contrôleur qui dépend seulement de grandeurs observables expérimentalement et qui gère des situations où les écoulements sont excités par des perturbations aléatoires extérieures. Pour faire face à ces deux aspects essentiels, deux stratégies ont été définies: l'identification d'un modèle non-linéaire reproduisant la dynamique de l'écoulement à partir seulement d'informations mesurables et la conception d'un compensateur linéaire robuste, basée sur la théorie du contrôle H∞, qui incorpore des propriété de robustesse dans la définition de la fonction objectif. La première partie de la thèse est consacrée à l'identification d'un modèle non-linéaire grâce à des données obtenues à partir d'une expérience menée dans la soufflerie subsonique (M = 0.1) S19 sur le site Chalais-Meudon de l'ONERA. Afin de décrire la dynamique de cet écoulement, et en particulier son contenu fréquentiel, l'écoulement sans contrôle a été caractérisé par des mesures par fil chaud et de pression instationnaire et par des clichés de vélocimétrie par images des particules (PIV) résolue en temps. Un filtrage temporel a été appliqué avec succès aux clichés PIV afin d'extraire la dynamique basse fréquence de l'écoulement. Cette étape est indispensable pour pouvoir gérer des écoulements turbulents caractérisés par un spectre fréquentiel très étendu. Les modes POD obtenus ont été utilisés comme base de projection pour le champ de vitesse et les trajectoires associées ont été interpolées (apprentissage statistique) sur une structure de modèle non-linéaire autorégressif exogène (NLARX). Il s'avère que les modèles obtenus ne sont pas robustes, dans le sens où ils ne parviennent pas à reproduire la dynamique d'un ensemble de données de validation, une fois adaptés à un ensemble de données d'apprentissage. Il a été démontré que cet échec est dû aux fortes non-linéarités observées dans l'écoulement de cavité, qui rendent impraticables les méthodes d'identification. La deuxième partie de la thèse est consacrée à la conception d'un contrôleur robuste à partir de simulations numériques d'un écoulement de cavité carrée, incompressible et en régime transitionnel, pour différents nombres de Reynolds. Diverses méthodes de synthèse de contrôleur ont été testées et évaluées en utilisant plusieurs mesures de robustesse. On a constaté que la technique traditionnelle de contrôle linéaire quadratique gaussien (LQG) présente une faible robustesse aux perturbations extérieures, tandis que d'autres, comme la technique LTR (Loop Transfer Recovery) et les contrôleurs basés sur les perturbations "les pires" (worst-case), améliorent la robustesse, mais pas suffisamment pour faire face à la forte non-linéarité de l'écoulement. Dans ce but, on met en place un contrôleur qui optimise les propriétés de robustesse par rapport à des incertitudes de type "entrée-multiplicative" et de type "entrée vers sortie". Celui-ci présente des marges de robustesse fortement augmentées par rapport à l'introduction de perturbations de la partie stable de la dynamique entré-sortie, même si le prix à payer en terme de performance est significatif. Une stratégie pour prendre en compte également des perturbations de la partie instable de la dynamique entrée-sortie, comme celles obtenues par un changement du nombre de Reynolds, a été présentée.

Il est bien rapporté dans la littérature que les caractéristiques du joint soudé ainsi que le transfert de chaleur à l'interface liquide/solide sont fortement dépendants de l'écoulement du métal liquide. Ce travail porte sur l'étude numérique de l'influence des paramétres opératoires de soudage sur les bains de fusion TIG. Un modéle 2D incluant la cathode, le plasma d'arc et l'anode est développé pour l'étude du soudage TIG spot. Il permet également de quantifier les effets de la composition chimique du gaz, et de simuler des techniques innovantes tel que le soudage avec gaz pulsé. Le soudage TIG avec torche mobile et apport de matiére est modélisé à l'aide d'une approche hybride 2D-3D. Les simulations numériques pour des positions de soudages horizontales révèlent une dissymétrie de la soudure dans la direction opposée à la gravité, ce qui est dû à l'action de la flottabilité et à la déformation assymétrique de la surface libre. Ce défaut peut être considérablement réduit lors du soudage de matériaux à taux de soufre différents, en plaçant la pièce bas soufre en bas.

L'objet de ce travail concerne l'application de la méthode de simulation des grandes échelles au mélange de deux fluides à propriétés thermodynamiques différentes.<br />L'origine des erreurs numériques liées à la discrétisation des équations de Navier-Stokes ainsi que leur interaction avec un modèle sous-maille sont étudiées pour une turbulence homogène et isotrope en auto-amortissement. Un code de calcul à haut pouvoir de solution est alors développé pour la simulation de couches de mélange bi-espèces. La réduction de l'amplitude des ondes acoustiques initiales est assurée par l'utilisation d'un champ en similitude temporelle.<br />L'importance relative des termes sous-maille issus des équations filtrées est mesurée à partir du filtrage explicite des champs de simulations numériques directes des couches de mélange temporelles $N_2/O_2$ et $H_2/O_2$.<br />L'utilisation d'une fermeture implicite autour d'un schéma numérique dissipatif est ensuite évoquée.

Ce travail a pour but le développement de techniques de séparation des contributions acoustique et aérodynamique de champs de pression pariétale sous des écoulements turbulents. Ces techniques de séparation, dites en nombre d'onde, se basent sur les propriétés spatiales de ces fluctuations. La compréhension des phénomènes de chacun des deux types fluctuations de pression pariétale s'exerçant sur la vitre latérale d'un véhicule vient servir un enjeu industriel de réduction du niveau sonore intérieur du bruit de forme des automobiles. Pour ce faire, des expérimentations en soufflerie anéchoïque sont menées afin d'obtenir des mesures de pression pariétale fluctuante pour divers types d'écoulements. Ainsi un écoulement de couche limite turbulente pleinement développée sur plaque plane, des écoulements décollés en aval de trois géométries de marches montantes, ainsi qu'en aval d'un montant de baie sur véhicule réel sont considérés. Une représentation des spectres en nombre d'onde-fréquence (obtenue par la méthode du corrélogramme spatial) est choisie pour séparer les deux contributions. Cette technique ne s'appliquant que sur des écoulements stationnaires et homogènes, de nouvelles techniques basées sur la Décomposition Modale Empirique (EMD) sont utilisées, d'une part pour améliorer la séparation des deux contributions (grâce à l'EMD spatiale), d'autre part pour obtenir une décomposition multi-échelle au cours du temps (grâce à l'Ensemble-EMD). Grâce à ces méthodes, une observation plus précise des phénomènes pariétaux convectifs et acoustiques est proposée.

In this work we review 12 years of developments in the field of residual based discretizations for hyperbolic problems and their application to the solution of the shallow water equations. Fundamental concepts related to the topic are recalled and he construction of second and higher order schemes for steady problems is presented. The generalization to time dependent problems by means of multi-step implicit time integration, space-time, and genuinely explicit techniques is thoroughly discussed. Finally, the issues of C-property, super consistency, and wetting/drying are analyzed in this framework showing the power of the residual based approach.