Mes différentes activités de recherche qui s'insèrent dans une problématique axée autour du transport multi-échelle en milieu poreux sont détaillées ici au travers de 3 volets: (1) Influence des hétérogénéités sur les processus de transfert (2) Prise en compte des effets de couplage en transport réactif (3) Impact des processus biologiques sur l'hydrodynamique et le transport Un survol des défis qui restent à surmonter et des perspectives à venir de ce travail de recherche concluent le document.

Les travaux présentés dans cette thèse portent sur l'ordonnancement d'applications sur des plate- formes hétérogènes à grande échelle. Dans la mesure où le problème général est trop complexe pour être résolu de façon exacte, nous considérons deux relaxations. Tâches divisibles : La première partie est consacrée aux tâches divisibles, qui sont des appli- cations parfaitement parallèles et pouvant être arbitrairement subdivisées pour être réparties sur de nombreux processeurs. Nous cherchons à minimiser le temps de travail total lors de l'exécution de plusieurs applications aux caractéristiques différentes sur un réseau linéaire de processeurs, sachant que les données peuvent être distribuées en plusieurs tournées. Le nombre de ces tour- nées étant fixé, nous décrivons un algorithme optimal pour déterminer précisément ces tournées, et nous montrons que toute solution optimale requiert un nombre infini de tournées, résultat restant vrai sur des plate-formes non plus linéaires mais en étoile. Nous comparons également notre méthode à des méthodes déjà existantes. Ordonnancement en régime permanent : La seconde partie s'attache à l'ordonnancement de nombreuses copies du même graphe de tâches représentant une application donnée. Au lieu de chercher à minimiser le temps de travail total, nous optimisons uniquement le cœur de l'or- donnancement. Tout d'abord, nous étudions des ordonnancements cycliques de ces applications sur des plate-formes hétérogènes, basés sur une seule allocation pour faciliter leur utilisation. Ce problème étant NP-complet, nous donnons non seulement un algorithme optimal, mais éga- lement différentes heuristiques permettant d'obtenir rapidement des ordonnancements efficaces. Nous les comparons à ces méthodes classiques d'ordonnancement, telles que HEFT. Dans un second temps, nous étudions des applications plus simples, faites de nombreuses tâches indépendantes, que l'on veut exécuter sur une plate-forme en étoile. Les caractéristiques de ces tâches variant, nous supposons qu'elles peuvent être modélisées par des variables aléatoires. Cela nous permet de proposer une ε-approximation dans un cadre clairvoyant, alors que l'ordonnan- ceur dispose de toutes les informations nécessaires. Nous exposons également des heuristiques dans un cadre non-clairvoyant. Ces différentes méthodes montrent que malgré la dynamicité des tâches, il reste intéressant d'utiliser un ordonnancement statique et non des stratégies plus dynamiques comme On-Demand. Nous nous intéressons ensuite à des applications, dont plusieurs tâches sont répliquées sur plu- sieurs processeurs de la plate-forme de calcul afin d'améliorer le débit total. Dans ce cas, même si les différentes instances sont distribuées aux processeurs tour à tour, le calcul du débit est difficile. Modélisant le problème par des réseaux de Petri temporisés, nous montrons comment le calculer, prouvant également que ce calcul peut être fait en temps polynomial avec le modèle Strict One-Port. Enfin, le dernier chapitre est consacré à l'application de ces techniques à un processeur multi- cœur hétérogène, le Cell d'IBM. Nous présentons donc un modèle théorique de ce processeur ainsi qu'un algorithme d'ordonnancement adapté. Une implémentation réelle de cet ordonnanceur a été effectuée, permettant d'obtenir des débits intéressants tout en simplifiant l'utilisation de ce processeur et validant notre modèle théorique.

Cette thèse a pour thème central l'influence stabilisatrice de politiques de taxation appropriées sur la dynamique économique. La stabilisation est ici entendue comme la possibilité d'atténuer voire d'éliminer, par de telles politiques, la probabilité d'apparition de fluctuations endogènes dans l'économie. Le premier chapitre examine, dans le cadre d'une économie réelle, les propriétés stabilisatrices de taxes linéairement progressives, proches des codes d'imposition par tranches prévalant dans de nombreux pays développés. Le deuxième chapitre met en évidence que l'influence stabilisatrice (contracyclique) de la taxation progressive diminue lorsqu'on introduit une disposition fiscale réaliste permettant aux entreprises de déduire de leur profit avant impôt les dépenses variables de maintenance. Le troisième chapitre évalue l'impact de la taxation des revenus dans une économie monétaire avec agents hétérogènes. Cette étude confirme que la progressivité fiscale réduit la probabilité d'apparition de fluctuations endogènes dans l'économie. Cependant on montre dans ce cadre monétaire, à la différence des économies réelles étudiées dans les deux premiers chapitres, que même un niveau élevé de progressivité ne peut éliminer complètement la possibilité de fluctuations endogènes. Enfin, le dernier chapitre analyse comment la politique fiscale peut affecter la volatilité agrégée et la croissance dans des économies sujettes à des contraintes de crédit. Dans un modèle de croissance endogène combinant des frictions sur le marché du capital et une inégalité d'accès aux opportunités d'investissement entre les individus, il est montré qu'une politique fiscale bien choisie, même linéaire, taxant les revenus du travail et organisant des transferts vers l'investissement innovant, peut remédier au défaut de financement de ces investissements lors des crises, stabiliser la dynamique et placer l'économie sur un sentier de croissance permanente et soutenue.

Un défi courant de la Mécanique Linéaire Elastique de la Rupture (MLER) est de prendre en considération les non-linéarités induites par les déformations du front de ssure. Pour ceci, une approche appropriée est la méthode de perturbation du front initiée par Rice (1985). Elle permet d'actualiser les Facteurs d'Intensité des Contraintes (FIC) quand le front d'une ssure plane est perturbé dans son plan. Cette approche, et ses prolongements ultérieurs à des cas plus complexes, sont rappelés dans cette thèse. Les applications concernent la déformation du front de ssure quand il se propage dans des milieux homogènes ou hétérogènes en rupture fragile, en fatigue ou propagation sous-critique. Les formes de ssures correspondant à des FIC uniformes le long du front ont été obtenues : fissure à front droit ou circulaire, mais également quand les bifurcations existent, à front ondulé. Pour un front initialement droit, on a montré que lors de la propagation dans un milieu homogène, d'éventuelles perturbations du front disparaîssent progressivement quelles que soient leurs longueurs d'onde, tandis que dans un milieu à propriétés de rupture désordonnées, le front devient rugueux de type Family-Vicsek (1985). La prolongation de ces approches de perturbation à des géomètries plus réalistes et à la coalescence de fissures est également envisagée. La comparaison quantitative avec des expériences est à l'étude. En outre, l'application de la MLER et de son principe énergétique de minimisation à l'initiation et aux morphologies de fissures apparaissant dans différents problèmes physiques (séchage de suspensions, impact de plaque, indentation de film et phénomènes géologiques) est considérée.

Cette étude concerne le delta du Nestos (Grèce du Nord) et les deltas dont la formation est plutôt influencée par l'action fluviatile et l'action de la houle. Les forages et les sondages électriques ont permis de constater, pendant le dernier stade de l'évolution du delta, une translation de l'écoulement du Nestos d'Ouest en Est. Ceci est aussi vérifié par l'analyse des conditions géomorphologiques et fut déterminant pour le régime hydrogéologique des deux secteurs du delta. Après l'aménagement de la rivière, commence la phase de destruction du delta. Les barrages naturels des îles de Thassopoula et Thassos n'ont pas permis la formation des profils sous-marins typiques, sauf dans le cas des côtes orientales. L'effet de la houle a donné trois types de côtes. Le transport des sédiments et la formation de la plate-forme sous-marine sont influencés par la présence des îles ci-dessus mentionnées. L'analyse morphométrique du delta exprime la distribution du volume des sédiments et indique que la convexité ou la concavité des côtes et des profils sous-marins résulte de cette distribution. Plus encore, cette analyse peut être utilisée pour la connaissance des conditions hydrogéologiques, avant la réalisation de recherches précises. Les analyses minéralogiques des sédiments et la détermination des fossiles ont complété la connaissance sur la formation du delta, et les conditions paléoécologiques. L'étude des conditions hydrogéologiques a été faite, pour les marbres, des bordures du delta, à l'aide de la géologie, des mesures du d~bit des sources et l'analyse des courbes, des analyses chimiques et isotopiques et de l'analyse tectonique des affleurements des marbres. Cette étude pour les systèmes aquifères de la plaine deltaïque, qui résultent par la formation du delta, est basée sur les essais de pompage, ainsi que sur les analyses chimiques et isotopiques. Dans ce domaine, on a essayé de critiquer les méthodes d'essais de pompage dans les conditions hétérogènes du delta et dans les conditions d'utilisation des forages qui exploitent toutes les couches aquifères traversées. En tenant compte des pertes de charge dans les forages, en essayant d'établir une corrélation entre le Débit Spécifiques des forages et la Transmissivité on est arrivé à établir un cadre d'application des méthodes des essais de pompage sous les conditions ci-dessus mentionnées. L'étude inclut encore la description des conditions hydrogéologiques et hydrochimiques de la nappe phréatique et du système des nappes profondes, leurs variations et leurs conditions d'alimentation.

Dans le domaine du calcul de haute performance, les architectures d'exécution sont en continuelle évolution. L'augmentation du nombre de nœuds de calcul, ou le choix d'une topologie réseau plus rapide représentent un investissement important tant en temps qu'en moyen financier. Les méthodes de prédiction de performances permettent de guider ce choix. En parallèle à ce développement, les systèmes HPC pair-à-pair (P2P) se sont également développés ces dernières années. Ce type d'architecture hétérogène permettrait la résolution des problèmes scientifiques pour un coût très faible par rapport au coût d'une architecture dédiée.Ce manuscrit présente une méthode nouvelle de prédiction de performances pour les applications réelles de calcul distribué, exécutées dans des conditions réelles. La prédiction prend en compte l'optimisation du compilateur. Les résultats sont extrapolables et ils sont obtenus pour un ralentissement réduit. Ce travail de recherche est implémenté dans un logiciel nouveau nommé dPerf. dPerf est capable de prédire les performances des applications C, C++ ou Fortran qui communiquent en utilisant les normes MPI ou P2P-SAP et qui s'exécutent sur une architecture cible pair à pair, hétérogène et décentralisée. La précision de cette contribution a été étudiée sur (i) la transformée Laplace, pour l'aspect séquentiel, (ii) le benchmark IS de NAS, pour l'aspect MPI, (iii) et le code de l'obstacle pour l'aspect calcul P2P décentralisé et l'extrapolation du nombre de nœuds.

L'objectif de la thèse effectuée dans le cadre du projet européen « Intelligent system for NET shape FORming of Sheet MEtal Product » est la recherche de procédures optimisées pour la mise au point de simulations numériques d'opérations industrielles d'emboutissage. Ce travail a pour but d'apporter des réponses efficaces, applicables au sein d'une entreprise en regard de la mauvaise perception de ce type de procédé dans le milieu industriel.<br />La procédure met en avant le rôle clé de l'identification des lois de comportement élastoplastiques définissant le comportement de la tôle pendant la simulation.<br />Dans un premier temps on s'intéresse à l'identification classique à partir d'essais de traction utilisés habituellement à cet effet.<br />Différentes techniques et protocoles sont alors retenus pour tirer le meilleur parti de ce type d'essais en utilisant notamment les techniques expérimentales d'analyse d'images.<br />Ces choix conduisent à une amélioration notable des résultats concernant l'effet du retour élastique, l'état de contraintes résiduelles et la prédiction de la localisation de la déformation.<br />Cependant les limites de ce type d'identification réalisé sur des essais homogènes apparaissent clairement.<br />C'est pourquoi, nous proposons une méthode originale d'identification des paramètres de la loi de comportement à partir d'essais hétérogènes voisins du cas réel de mise en forme. L'optimisation des coefficients de la loi de comportement est alors fondée sur l'utilisation de réseaux de neurones. Les avantages de cette méthode sont démontrés sur un exemple d'emboutissage de pièces axisymétriques à flan bloqué.

Cette thèse présente une contribution aux développements de méthodes numériques pour la simulation d'écoulements en milieu poreux, en particulier par des méthodes de décomposition de domaine espace--temps qui permettent l'utilisation de pas de temps différents dans les différents sous--domaines. Nous étudions deux types de méthodes: la première est basée sur une généralisation de l'opérateur de Steklov-Poincaré au cas de problèmes dépendants du temps, et la seconde est basée sur la méthode de Relaxation d'Onde Optimisée de Schwarz (OSWR) dans laquelle des conditions de transmission plus générales (Robin ou Ventcell) sont utilisées pour accélérer la convergence de l'algorithme. Ces deux méthodes sont étudiées sur une formulation mixte qui est bien adaptée à la modélisation de l'écoulement et du transport en milieu poreux. Nous considérons tout d'abord un problème de diffusion et formulons, pour chaque méthode, un problème sur l'interface espace -temps entre les sous-domaines. Le caractère bien posé de ces problèmes, avec des conditions aux limites de Dirichlet ou de Robin, est démontré. Les preuves de convergence de l'algorithme OSWR et de sa version semi-discrète sous forme mixte sont également données. Des expériences numériques sont menées en 2D pour comparer les performances des deux méthodes sur des problèmes fortement hétérogènes, et un préconditionneur de Neumann--Neumann dépendant du temps permet d'accélérer la première méthode. Les deux méthodes sont ensuite étendues au cas d'une équation d'advection-diffusion, l'advection et la diffusion étant traitées séparément grâce une technique de séparation d'opérateurs, ce qui permet d'utiliser des pas de temps différents pour les deux phénomènes dans chaque sous-domaine. Des conditions de transmission sont proposées séparément pour l'advection et pour la diffusion. La convergence des méthodes est étudiée sur des exemples numériques, pour des problèmes en régime d'advection dominante ou de diffusion dominante, et leur précision en temps est étudiée dans le cas de grilles non-conformes en temps. Deux exemples inspirés de la simulation du stockage de déchets nucléaires sont étudiés, et la simulation sur des temps longs est réalisée par l'intermédiaire de fenêtres en temps. Nous considérons également la méthode OSWR avec des conditions de transmission de Ventcell, étendues à la formulation mixte. Nous démontrons que les problèmes de sous--domaine avec des conditions aux limites de Ventcell sont bien posés. Nous comparons les performances des paramètres optimisés pour Ventcell et Robin dans le cas de problèmes hétérogènes pour une décomposition en deux sous-domaines. Enfin, nous étudions l'extension des deux méthodes au cas où l'interface représente une fracture pour un modèle réduit d'écoulement dans un milieu poreux fracturé.

Les applications auto-adaptatives, dont le comportement évolue en fonction de l'environnement, sont un élément clé des systèmes de demain. L'utilisation de matériel reconfigurable, combiné à la parallélisation des unités de calcul, permettent d'envisager de nouveaux niveaux de performances pour ces mêmes applications. L'objectif de cette thèse est de mettre en place un ensemble d'outils permettant la description et le déploiement d'applications parallèles auto-adaptatives. Nous proposons à la fois un modèle d'application parallèle et une architecture de plateforme reconfigurable destinée au déploiement des applications conçues en utilisant ce modèle. Notre modèle d'applications sépare le contrôle du calcul en isolant ce dernier en différents noyaux virtuels. Associée à une représentation du contrôle indépendante de la plateforme, la structure de l'application est donc totalement portable. Indépendamment de celle-ci, les noyaux de calcul peuvent être distribués selon plusieurs implémentations, selon la plateforme, mais également pour proposer différents niveaux de performances. Une couche de virtualisation permet de faire le lien entre la partie contrôle et les noyaux, en traduisant les ordres génériques en actions adaptées à l'implémentation. Concernant la plateforme, nous proposons une architecture permettant l'intégration de ressources de calcul logicielles et matérielles, pouvant être implémentées tant statiquement qu'en utilisant du matériel reconfigurable. Cette architecture parallèle, inspirée du modèle des supercalculateurs, doit permettre d'utiliser tout type d'unités d'exécution et de matériel reconfigurable comme base matérielle pour la plateforme.

L'objectif de cette thèse était de proposer, puis de valider expérimentalement, un ensemble de méthodes algorithmiques permettant le remaillage parallèle de maillages distribués, en s'appuyant sur une méthode séquentielle de remaillage préexistante. Cet objectif a été atteint par étapes : définition de structures de données et de schémas de communication adaptés aux maillages distribués, permettant le déplacement à moindre coût des interfaces entre sous-domaines sur les processeurs d'une architecture à mémoire distribuée ; utilisation d'algorithmes de répartition dynamique de la charge adaptés aux techniques parallèles de remaillage ; conception d'algorithmes parallèles permettant de scinder le problème global de remaillage parallèle en plusieurs sous-tâches séquentielles, susceptibles de s'exécuter concurremment sur les processeurs de la machine parallèle. Ces contributions ont été mises en oeuvre au sein de la bibliothèque parallèle PaMPA, en s'appuyant sur les briques logicielles MMG3D (remaillage séquentiel de maillages tétraédriques) et PT-Scotch (repartitionnement parallèle de graphes). La bibliothèque PaMPA offre ainsi les fonctionnalités suivantes : communication transparente entre processeurs voisins des valeurs portées par les noeuds, les éléments, etc. ;remaillage, selon des critères fournis par l'utilisateur, de portions du maillage distribué, en offrant une qualité constante, que les éléments à remailler soient portés par un unique processeur ou bien répartis sur plusieurs d'entre eux ; répartition et redistribution de la charge des maillages pour préserver l'efficacité des simulations après remaillage.