Nous considèrerons ici les modèles linéaires de multiflots, en mettant l'accent sur leurs multiples applications, notamment dans les domaines de l'ordonnancement et de la gestion de production. Il est bien connu que ces problèmes, présentés sous forme de programmes linéaires, sont difficiles à résoudre, contrairement à leurs homologues en flot simple. Les méthodes de résolution classiques proposent, déjà dans le cas continu, des solutions approchées. On distingue: les méthodes de décomposition par les prix, par les ressources, ainsi que les techniques de partitionnement. Si l'on rajoute la contrainte d'intégralité sur les flots, ces problèmes deviennent extrêmement difficiles. Nous nous sommes intéressés à un cas particulier des problèmes de multiflots, à savoir: le biflot entier de coût minimum. Nous avons développé une approche de résolution heuristique basée sur un principe de décomposition mixte, opérant itérativement, à la fois par une allocation de ressources et par un ajustement des coûts. L'implémentation de cette approche met en évidence des résultats prometteurs, obtenus sur des problèmes de biflot purs, générés aléatoirement. Nous avons donc envisagé une deuxième application sur des problèmes de biflot plus structurés. Ces problèmes de biflot ont été proposés pour la modélisation du problème de voyageur de commerce. Cette application débouche d'une part, sur l'utilisation d'un algorithme de recherche d'un circuit hamiltonien dans un graphe, et d'autre part, sur le développement de techniques heuristiques pour la construction de tournées intéressantes

Le polyèdre associé à un problème d'optimisation combinatoire est l'enveloppe convexe des (vecteurs d'incidence des) solutions réalisables de ce problème. De nombreux problèmes d'optimisation combinatoire se formulent comme une maximisation de fonctions linéaires sur les polyèdres qui leurs sont associés. La description du polyèdre par un système d'inéquations linéaires est intimement liée à la résolution du problème correspondant, par le biais de la programmation linéaire. Afin de déterminer un tel système, une approche classique consiste à décomposer le problème en sous-problèmes tels que les polyèdres associés soient connus ; une composition ultérieure de ces derniers conduit à une description du polyèdre associé au problème considéré. L'objet principal de cette thèse est l'étude de la composition des polyèdres. Dans un premier temps, une approche de composition, basée sur la programmation dynamique et les méthodes de projection polyédrale, est étudiée et des résultats généraux sont proposés, permettant ainsi d'unifier des recherches existantes dans ce domaine. Cette approche est, ensuite, appliquée à la composition de polyèdres associés au problème du voyageur de commerce. En seconde partie, considérant le problème du stable, des opérations sur les graphes (composition par identification de sous-graphes de deux graphes donnés, adjonction d'une nouvelle arête) sont traitées. Des résultats polyédraux sont donc donnés, et des conséquences concernant la perfection et la h-perfection des graphes sont montrés

Cette thèse concerne le développement et la mise en œuvre d'un environnement pour le calcul scientifique et la modélisation. L'approche retenue est celle d'une décomposition stratifiée des problèmes, ceci dans un double but: marquer le cheminement progressif des étapes de description, allant de l'énoncé informel vers un langage cible en passant par des langages intermédiaires plus ou moins formalisés ; et, d'obtenir une décomposition structurée, modulaire, pour aller du problème initial vers le programme. Dans le but de vérifier la cohérence des descriptions, des schémas de résolutions, des décompositions, nous associons à tout énoncé des conditions logiques dépendant du langage de description. Pour cela, il nous a paru nécessaire d'étudier les formulations logiques décrites par des systèmes polynômiaux sur les corps finis de la forme Z/pZ. L'étude de ces systèmes nous conduisent à traiter le problème de l'élimination des quantificateurs sur un corps fini, le problème du calcul du résultat sur Z/pZ: des algorithmes sont proposés, ainsi qu'une généralisation de la méthode de Dixon-Biard. Le problème de la déduction est aussi abordé. Ces algorithmes nous permettent de vérifier localement la cohérence d'un énoncé mais aussi d'une décomposition de problème. Ceci rend envisageable une vérification globale. Un éditeur de strates sous Grif est présenté

Ce manuscrit comporte trois parties distinctes, la première concerne la simulation numérique d'écoulements de fluides. La deuxième partie porte sur les équations de l'électromagnétisme et leur couplage avec les équations cinétiques de Vlasov dans le cadre de simulations numériques en physique des plasmas. La dernière partie évoque rapidement des travaux qui ont donné lieu à des publications mais qui ne rentrent pas complètement dans un des deux cadres abordés précédemment. Dans la première partie, l'objectif est d'étendre aux maillages triangulaires non structurés une méthode numérique éprouvée pour résoudre les équations de Stokes bidimensionnelles : la méthode Marker And Cell qui a été développée sur des maillages en quadrilatères quasi-réguliers dans les années 60. L'idée proposée pour cela est de résoudre le problème de Stokes avec pour variables le tourbillon, la vitesse et la pression. Alors que les résultats numériques obtenus sur des maillages réguliers sont satisfaisants, ceux sur des maillages non structurés ne le sont pas. Il s'est avéré lors de l'étude théorique que ce problème est un problème de stabilité. On montre théoriquement et numériquement que la formulation tourbillon-vitesse-pression est une généralisation de la formulation fonction courant-tourbillon permettant la prise en compte de conditions limites plus générales. L'instabilité, due à des fonctions harmoniques discrètes, peut être levée en utilisant de véritables fonctions harmoniques, calculées à l'aide de leur représentation intégrale, dans le schéma numérique. On résout ainsi l'instabilité de la formulation fonction courant-tourbillon et on améliore les précédents résultats de convergence connus de cette formulation. En particulier, on démontre alors une convergence en moyenne quadratique du tourbillon de l'ordre de 3/2 à 2 dans les cas les plus réguliers (contre un-demi avant). Puis on utilise le fait que la formulation tourbillon-vitesse-pression est équivalente à la formulation fonction courant-tourbillon pour redéfinir une nouvelle formulation tourbillon-vitesse-pression. En effet, il est bien connu que la formulation classique en fonction courant-tourbillon (n'utilisée qu'en deux dimensions d'espace) est mal posée lorsque l'on cherche le tourbillon dans l'espace de Sobolev H1 car son gradient n'est alors pas contrôlé, mais bien posée dans un autre espace de fonctions moins régulières. On étend alors ce résultat au cas tri-dimensionnel et l'on obtient une nouvelle formulation en tourbillon-vitesse-pression bien posée dans un nouvel espace. On démontre aussi théoriquement que ce nouvel espace est bien celui introduit en 2D, ce qu'on confirme par des résultats numériques. La deuxième partie concerne la résolution d'équations cinétiques intervenant dans la simulation directe des plasmas et des faisceaux de particules chargées (modèles de Vlasov-Poisson ou Vlasov-Maxwell). Grâce à l'augmentation de la puissance de calculs des ordinateurs, la simulation de l'évolution des plasmas et des faisceaux de particules basée sur une résolution directe de l'équation de Vlasov sur un maillage de l'espace des phases devient une alternative aux méthodes particulaires (Particle In Cell) habituellement employées. La force de ces simulations directes réside dans le fait qu'elles ne sont pas bruitées (contrairement aux méthodes PIC) et que l'approximation est de même résolution sur tout l'espace des phases, en particulier dans les régions à faible densité de particules où des phénomènes physiques importants ont lieu. L'inconvénient principal est que beaucoup de points sont inutiles car la fonction de distribution des particules y est nulle, ce qui rend ces méthodes directes coûteuses en temps de calcul. On introduit alors une méthode de résolution directe de l'équation de Vlasov sur un maillage {\bf{mobile}} de l'espace des phases. Ce qui permet de ne mailler que la partie de l'espace des phases sur laquelle la fonction de distribution des particules est {\emph{a priori}} non nulle. Nous avons utilisé avec succès cette méthode de maillage mobile en 1D afin de simuler un problème d'interaction laser-plasma. Nous introduisons donc un maillage mobile de l'espace des phases qui suit parfaitement le développement des instabilités et permet de réduire drastiquement le temps de calcul. Nous avons aussi obtenu les premiers résultats d'une méthode de maillage mobile en 4D en couplant le maillage mobile et une méthode de décomposition de domaines. \\ En ce qui concerne la résolution des équations de Vlasov-Maxwell : nous travaillons sur le développement d'un solveur Maxwell éléments finis d'arêtes d'ordre élevé couplé à un code PIC en trois dimensions d'espaces (6D de l'espace des phases). Un point important afin que le couplage fonctionne est que l'équation de conservation de la charge doit être vérifiée au niveau numérique à chaque pas de temps. Ce qui implique que le courant obtenu à partir de l'évolution de l'équation de Vlasov doit être calculé d'une façon bien particulière. Les premiers résultats obtenus en 2D confirment que la méthode de calcul du courant proposée conserve bien la charge comme attendu.

Cette thèse se propose d'étudier la thermodynamique de la matière nucléaire qui constitue l'écorce des proto-étoiles et étoiles à neutrons. Une connaissance détaillée des propriétés thermodynamiques de la matière d'étoile est nécessaire afin d'améliorer la compréhension des phénomènes physiques impliqués dans le refroidissement des proto-étoiles, et dans la formation de supernovae de type II. L'un des objectifs fondamentaux est d'extraire le diagramme des phase de la matière d'étoile afin de déterminer si celui-ci présente des instabilités et/ou des points critiques. Le travail présenté ici se divise en deux parties, une première portant sur des approches classiques, et une seconde sur une approche quantique. Les approches classiques sont basées sur le modèle d'Ising et le groupe de renormalisation. Elles vont nous permettre d'obtenir des informations qualitatives sur la phénoménologie des transitions de phase de la matière d'étoile, et de discuter l'influence de la frustration Coulombienne sur le diagramme des phases. L'approche quantique est basée sur un modèle de dynamique moléculaire fermionique que nous avons dérivé à partir du formalisme de la fonctionnelle de la densité, et implémenté numériquement avec des forces de Skyrme optimisée pour la matière riche en neutrons. Le travail de cette thèse montre des premières applications à l'étude thermodynamique de systèmes finis, et à des calculs de structure nucléaire pour des noyaux légers. Nous proposerons également une ébauche du formalisme qui permettra à terme de traiter numériquement le problème quantique de la matière infinie d'étoile à l'aide du modèle de dynamique moléculaire.

.

On étudie la distance de Prokhorov, la distance de Geffroy et la distance de Fortet-Mourier-Wasserstein. On résout en particulier le problème du calcul des distances. On traite quelques problèmes relatifs à l'estimation. Puis on donne une condition nécessaire et suffisante de non séparation asymptotique de deux chaines de Markov lorsque l'espace des états est de cardinal M

.

.

L'objet de cette thèse est l'étude de quelques propriétés de l'équation d'évolution de Schrödinger. Dans un premier temps, on s'intéresse à un problème inverse concernant cette équation posée en domaine borné, avec potentiel, lequel dépend uniquement de la variable d'espace, et donnée de Dirichlet sur le bord. On démontre, à l'aide d'une inégalité de Carleman, que le problème inverse de la détermination du potentiel à partir de la mesure du flux de la solution à travers une partie du bord est un problème bien posé. Dans un deuxième temps, il est question de l'équation de Schrödinger considérée dans $\mathbb R^3$ avec un potentiel coulombien, localement singulier, et un potentiel électrique non borné, tous deux dépendant des variables d'espace et de temps. On montre successivement l'existence d'une unique solution régulière pour l'équation linéaire et pour l'équation avec non-linéarité de Hartree. Ce sont des étapes préliminaires à l'étude d'un système couplant à travers le potentiel coulombien, cette équation de Hartree-Fock et une équation issue de la dynamique newtonienne. Les résultats obtenus ici sont indispensables à l'étude finale des problèmes de contrôle optimal bilinéaire posés à partir de ces différentes équation, le contrôle de la solution étant effectué par le potentiel électrique. On démontre l'existence d'un contrôle optimal et on donne la condition d'optimalité correspondante dans les cas appropriés\vspace(0,5cm)