D'Alembert est l'auteur d'un quatrième et dernier traité d'hydrodynamique méconnu des historiens des sciences : le Mémoire 57 du t. VIII de ses "Opuscules mathématiques" (1780). Il y répond, pour l'essentiel, à un écrit de Borda de 1766 contenant une sévère remise en question des principes et concepts fondateurs de son "Traité des fluides" (1744). Cette polémique entre les deux géomètres s'inscrit dans une période de crise de l'hydrodynamique touchant à la concordance entre théorie et expérience. Les équations aux dérivées partielles obtenues par D'Alembert dans son "Essai d'une nouvelle théorie de la résistance des fluides" (1752), puis par Euler dans ses célèbres mémoires de 1755, n'ont pu être résolues. L'approche unidimensionnelle initiée par Daniel Bernoulli dans l'"Hydrodynamica" (1738) puis adoptée par Jean Bernoulli dans l'"Hydraulica" (1743) et par D'Alembert dans le "Traité des fluides" (1744), repose, quant à elle, sur une approximation, l'hypothèse du parallélisme des tranches, trop éloignée des résultats expérimentaux disponibles. C'est cet écart que Borda tente de pallier en 1766 afin de faire sortir la discipline de l'impasse dans laquelle elle se trouve.<br />Dans le Mémoire 57, D'Alembert, directement visé par ce dernier, donne son point de vue sur le sujet, répond implicitement aux critiques de Borda, revient sur ses écrits antérieurs, et propose de nouvelles pistes de recherche théoriques. L'étude de ce traité, tant du point de vue physique que mathématique, engage ici un réexamen des principes et concepts fondateurs de sa théorie des écoulements, la plupart n'ayant été que partiellement abordés par l'historiographie. Il nous permet également de situer ses recherches vis-à-vis de celles des principaux artisans du développement de l'hydrodynamique à cette époque, Daniel Bernoulli, Jean Bernoulli, Euler et Lagrange, et de donner une idée de leur réception et de leur pérennité.

Les résultats présentés dans cette thèse concernent l'étude de la<br />théorie de la diffusion pour des champs de Dirac dans plusieurs<br />espaces-temps de la relativité générale. Les méthodes complètement<br />dépendantes du temps développées par Enss, Sigal, Soffer, Graf,<br />Derezi\'nski et Gérard constituent le fil conducteur de ce<br />travail. Ces méthodes sont basées sur des estimations de propagation<br />comme les estimations de vitesse minimale (obtenues par une théorie de<br />Mourre) qui correspondent à une version faible du principe de Huygens <br />et sur l'étude d'observables asymptotiques naturelles comme les<br />opérateurs de vitesse asymptotiques. Dans un premier temps, on teste<br />ces méthodes en étudiant la propagation de champs de Dirac, massifs ou<br />non, perturbés par des potentiels à longue portée, en espace-temps<br />plat. On montre ainsi <br />l'existence et la complétude asymptotique des opérateurs d'onde<br />modifiés. Dans un deuxième temps, on s'intéresse à des situations<br />géométriques plus compliquées en étudiant la propagation de ces champs<br />à l'extérieur de trous noirs de Reissner-Nordström (à symétrie<br />sphérique) et de Kerr-Newman (en rotation) du point de vue<br />d'observateurs lointains. L'originalité de ce type d'étude réside dans<br />le fait que les observateurs distinguent deux régions asymptotiques<br />(l'horizon du trou noir et l'infini spatial) aux structures<br />géométriques bien différentes ce qui entraîne l'existence de deux<br />canaux de diffusion. Dans le cas de trous noirs à symétrie<br />sphérique, une décomposition sur une base d'harmoniques sphériques<br />permet de se ramener à un problème à une dimension d'espace, du type<br />espace-temps plat. La difficulté essentielle provient alors de<br />l'absence de symétrie sphérique des trous noirs de Kerr-Newman qui<br />rend impossible une telle simplification. Dans les deux cas, on montre <br />l'existence et la complétude asymptotique des opérateurs d'onde<br />(modifiés à l'infini) à l'aide des méthodes dépendantes du temps.

Ce travail porte sur l'étude des méthodes de décomposition de domaine et leur application pour résoudre des problèmes de contrôle optimal régis par des équations aux dérivées partielles. Le principe de ces méthodes consiste à ramener des problèmes de grande taille sur des géométries complexes en une suite de sous-problèmes de taille plus petite sur des géométries plus simples. En considérant une décomposition sans recouvrement, l'intérêt de ces méthodes pour les problèmes de contrôle optimal réside au niveau de l'intégration de l'équation d'état, puisqu'il est possible de partitionner le problème en une suite de problèmes plus petits, quitte à contraindre les interfaces entre les sous-domaines à obéir à des conditions de raccordement afin de déduire la solution globale à partir des solutions locales. Dans une première partie, nous étudions le cas elliptique. Nous considérons simultanément la minimisation de la fonction coût et des raccordements sur les frontières entre les sous-domaines. Cette combinaison de problèmes de minimisation et de méthodes de décomposition de domaine est traitée par des techniques de Lagrangien augmenté. Nous montrons que, sur le domaine décomposé, le problème initial se réduit à la recherche d'un point-selle. Une étude des méthodes de Lagrangien nous a permis de choisir une variante d'algorithmes existants dans la littérature et de les combiner avec un algorithme de décomposition de domaine. Dans la seconde partie, nous développons l'extension de cette approche aux problèmes de contrôle optimal régis par des systèmes paraboliques en considérant uniquement une décomposition en espace du domaine de calcul. Dans une dernière partie, nous considérons une décomposition de domaine avec recouvrement à chaque pas de la minimisation. D'une part, nous construisons un algorithme parallèle en utilisant la méthode de Schwarz multiplicative en tant que solveur. Ceci permet de déduire naturellement l'état adjoint par transposition des systèmes directs locaux. L'algorithme global défini par la méthode de minimisation de type quasi-Newton et ce solveur de Schwarz constitue une méthode robuste de résolution du problème de contrôle optimal, mais coûteuse. D'autre part, et plus particulièrement, pour des problèmes de grande taille, l'algorithme de type quasi-Newton, combiné avec le solveur de Krylov BiCGSTAB préconditionné par une méthode de Schwarz additive, est plus compétitif dans la mesure oû l'on obtient de bonnes performances parallèles. De nombreux résultats sont présentés pour préciser le comportement des algorithmes d'optimisation quand ils sont utilisés avec des méthodes de Schwarz.

Nous nous intéressons à l'étude d'un problème d'analyse des transferts de chaleur qui modélise une opération de soudage. L'approche que nous considérons ne s'occupe que de la partie solide de la plaque. Elle consiste à résoudre un problème à frontière libre. Pour cela, nous proposons une formulation en optimisation de forme. Le problème d'état est gouverné par un opérateur qui, pour certaines données, n'est pas coercif. Cela complique l'étude de la continuité du problème d'état. Nous surmontons cette difficulté en utilisant le degré topologique de Leray-Shauder, ainsi nous montrons l'existence d'un domaine optimal. Ensuite, nous considérons une discrétisation de ce problème basée sur les éléments finis linéaires. Nous prouvons alors que le problème discret admet une solution et nous montrons qu'une sous-suite des solutions de ce problème convergence vers la solution du problème continu. Enfin, nous présentons des résultats numériques réalisés par deux méthodes : la méthode déterministe basée sur le calcul du gradient de forme, et les algorithmes génétiques combinés avec la logique floue et le calcul parallèle. Ainsi une étude comparative de ces deux méthodes aux niveaux qualitatif et quantitatif a été présentée.

La motivation de cette thèse de Doctorat est de modéliser quelques problèmes biologiques avec des systèmes et des équations de réaction-diffusion. La thèse est divisée en sept chapitres: 1. On modélise des ions de calcium et des protéines dans une épine dendritique mobile (une microstructure dans les neurones). On propose deux modèles, un avec des protéines qui diffusent et un autre avec des protéines fixées au cytoplasme. On démontre que le premier problème est bien posé, que le deuxième problème est presque bien posé et qu'il y a un lien continu entre les deux modèles. 2. On applique les techniques du Chapitre 1 pour un modèle d'infection virale et réponse immunitaire dans des cellules cultivées. On propose comme avant deux modèles, un avec des cellules qui diffusent et un autre avec des cellules fixées. On démontre que les deux problèmes sont bien posés et qu'il y a un lien continu entre les deux modèles. On Žtudie aussi le comportement asymptotique et la stabilité des solutions pour des temps larges, et on fait des simulations dans Matlab. 3. Dans le Chapitre 3 on montre que la croissance a deux effets positives dans la formation de motifs ou patterns. Le premier est un effet anti-explosion (anti-blow-up) car les solutions sur un domaine croissant explosent plus tard que celles sur un domaine fixé, et si la croissance est suffisamment rapide alors elle peut même empêcher l'explosion. Le deuxième est un effet stabilisant car les valeur propres sur un domaine croissant ont des parties réelles plus petites que celles sur un domaine fixé. 4. On étend la définition de front progressif à des variétés et on en étudie quelques propriétés. 5. On étudie des front progressifs sur la droite réelle. On démontre qu'il y a deux fronts progressifs qui se déplacent dans des directions opposées et qu'ils se bloquent mutuellement, générant ainsi une solution stationnaire non-triviale. Cet exemple montre que pour des modèles à diffusion non-homogène les fronts progressifs ne sont pas nécessairement des invasions. 6. On étudie des fronts progressifs sur la sphère. On démontre que pour des sous-domaines de la sphère avec des conditions aux limites de Dirichlet le front progressif est toujours bloqué, tandis que pour la sphère complète le front peut ou bien invahir ou bien être bloqué, tout en fonction des conditions initiales. 7. On étudie un problème elliptique aux valeurs propres nonlinéaires. Sur la sphère de dimension 1 on démontre l'existence de multiples solutions non-triviales avec des techniques de bifurcation. Sur la sphère de dimension n on utilise les mêmes arguments pour dŽmontrer l'existence de multiples solutions non-triviales à symétrie axiale, i.e. qui ne dépendent que de l'angle vertical.

Cette thèse s'intéresse à la modélisation et à la commande optimale de réseaux de distribution d'eau chaude avec pertes de chaleur.<br/><br/>Développement de logiciels d'aide à la compréhension du comportement dynamique du réseau, à sa gestion optimale face à des variations de données extérieures (consommations, tarifs énergétiques,...). Modélisation et optimisation du fonctionnement hydraulique et thermique du réseau, par rapport à la somme du coût de production de l'eau chaude, et du coût de fonctionnement des pompes. Exemples de simulation et d'optimisation réalisées sur des petits réseaux

Cette thèse traite des inéquations variationnelles stochastiques et de leurs applications aux vibrations de structures mécaniques. On considère d'abord un algorithme numérique déterministe pour obtenir le régime stationnaire d'une inéquation variationnelle stochastique modélisant un oscillateur élasto-plastique excité par un bruit blanc. Une famille de solutions d'équations aux dérivées partielles définissant la mesure invariante par dualité est étudiée comme alternative à la simulation probabiliste. Puis, nous présentons une nouvelle caractérisation de l'unique mesure invariante. Dans ce contexte, nous montrons une relation liant des problèmes non-locaux et des problèmes locaux en introduisant la définition des cycles courts. Dans un cadre orienté vers les applications, nous démontrons que la variance de la déformation plastique cro^it linéairement avec le temps et nous caractérisons rigoureusement le coefficient de dérive en introduisant la définition des cycles longs. Dans la suite, nous étudions un processus approché de la solution de l'inéquation comportant des sauts aux instants de transition de l' état plastique vers l' état élastique. Nous prouvons que la solution approchée converge sur tout intervalle de temps ni vers la solution de l'inéquation, lorsque la taille du saut tend vers 0. Ensuite, nous défi nissons une inéquation variationnelle stochastique pour modéliser un oscillateur élasto-plastique excité par un bruit blanc filtré. Nous prouvons la propriété ergodique du processus sous-jacent et nous caractérisons sa mesure invariante. Nous étendons la méthode de A.Bensoussan et J.Turi avec une difficulté supplémentaire due à l'accroissement de la dimension. Finalement, dans un chapitre orienté vers l'expérimentation numérique, nous mettons en évidence par les simulations probabilistes le phénomène de phases micro-élastiques. Leur impact concerne des grandeurs utiles a l'ingénieur comme la fréquence des déformations plastiques. Un critère empirique qui peut ^etre utile à l'ingénieur est fourni afin de ne pas prendre en compte les phases micro-élastiques et ainsi évaluer d'une façon réaliste, à partir de la mesure invariante, les statistiques de la déformation plastique d'un oscillateur élasto-plastique excité par un bruit blanc.

Nous nous intéressons d'abord à l'équation aux dérivées partielles stochastique de Cahn-Hilliard en dimension 1 avec une seule singularité. C'est une équation d'ordre 4 dont la non linéarité est de type logarithmique ou en puissance négative $x^{-\alpha}$, à laquelle on ajoute la dérivée d'un bruit blanc en espace et en temps. On montre l'existence et l'unicité des solutions en utilisant les solutions d'équations approchées aux non linéarités Lipschitz. La présence d'une mesure de réflexion permet d'assurer l'existence de solutions. On étudie ces mesures à l'aide des mesures de Revuz associées et, grâce à une formule d'intégration par parties, on montre qu'elles sont identiquement nulles lorsque alpha est plus grand ou égal à 3. Dans un deuxième temps, on considère la même équation mais avec deux singularités logarithmiques en +1 et -1. Il s'agit du modèle complet de l'équation de Cahn-Hilliard. Cette fois-ci on utilise des équations approchées aux non linéarités polynomiales pour montrer l'existence et l'unicité de solutions. Deux mesures de réflexion doivent ici être ajoutées pour assurer l'existence. De plus, on montrera que la mesure invariante est ergodique. Enfin, on étudie l'équation déterministe : des simulations numériques basées sur une méthode d'élements finis de hauts degrés permettent d'illustrer plusieurs résultats théoriques. La capture des interfaces et des états stationnaires requiert une attention particulière. On s'intéressera également aux bifurcations autour de la première valeur propre du Laplacien sur des domaines généraux. Par ailleurs, quelques simulations stochastiques permettent de mettre en évidence les instants de contact avec les singularités, les évolutions stochastiques en temps long et les changements d'états stationnaires.

La thèse est consacrée au développement de nouvelles méthodes de modélisations mathématiques en biologie et en médecine, du type "off-lattice" modèles hybrides discret-continus, et de leurs applications à l'hématopoïèse et aux maladies sanguines telles la leucémie et l'anémie. Dans cette approche, les cellules biologiques sont considérées comme des objets discrets alors que les réseaux intracellulaire et extracellulaire sont décrits avec des modèles continus régis par des équations aux dérivées partielles et des équations différentielles ordinaires. Les cellules interagissent mécaniquement et biochimiquement entre elles et avec le milieu environnant. Elles peuvent se diviser, mourir par apoptose ou se différencier. Le comportement des cellules est déterminé par le réseau de régulation intracellulaire et influencé par le contrôle local des cellules voisines ou par la régulation globale d'autres organes. Dans la première partie de la thèse, les modèles hybrides du type "off-lattice" dynamiques sont introduits. Des exemples de modèles, spécifiques aux processus biologiques, qui décrivent au sein de chaque cellule la concurrence entre la prolifération et l'apoptose, la prolifération et la différenciation et entre le cycle cellulaire et de l'état de repos sont étudiés. L'émergence des structures biologiques est étudiée avec les modèles hybrides. L'application à la modélisation des filamente de bactéries est illustrée. Dans le chapitre suivant, les modèle hybrides sont appliqués afin de modéliser l'érythropoïèse ou production de globules rouges dans la moelle osseuse. Le modèle inclut des cellules sanguines immatures appelées progéniteurs érythroïdes, qui peuvent s'auto-renouveler, se différencier ou mourir par apoptose, des cellules plus matures appelées les réticulocytes, qui influent les progéniteurs érythroïdes par le facteur de croissance Fas-ligand, et des macrophages, qui sont présents dans les îlots érythroblastiques in vivo. Les régulations intracellulaire et extracellulaire par les protéines et les facteurs de croissance sont précisées et les rétrocontrôles par les hormones érythropoïétine et glucocorticoïdes sont pris en compte. Le rôle des macrophages pour stabiliser les îlots érythroblastiques est montré. La comparaison des résultats de modélisation avec les expériences sur l'anémie chez les souris est effectuée. Le quatrième chapitre est consacré à la modélisation et au traitement de la leucémie. L'érythroleucémie, un sous-type de leucémie myéloblastique aigüe (LAM), se développe à cause de la différenciation insuffisante des progéniteurs érythroïdes et de leur auto-renouvellement excessif. Un modèle de type "Physiologically Based Pharmacokinetics-Pharmacodynamic" du traitement de la leucémie par AraC et un modèle de traitement chronothérapeutique de la leucémie sont examinés. La comparaison avec les données cliniques sur le nombre de blast dans le sang est effectuée. Le dernier chapitre traite du passage d'un modèle hybride à un modèle continu dans le cas 1D. Un théorème de convergence est prouvé. Les simulations numériques confirment un bon accord entre ces deux approches.

L'étude de la propagation de la lumière dans les fibres optiques monomodes requiert la prise en compte de plusieurs phénomènes compliqués tels que la dispersion modale de polarisation et l'effet Kerr. Il s'est avéré que l'évolution de l'enveloppe lentement variable du champ électrique est bien décrite par un système couplé d'équations de Schrödinger non linéaires à coefficients aléatoires : l'équation de Manakov PMD. Cette équation fait intervenir différentes échelles dont le ratio est donné par un petit paramètre. La première partie de ce travail consiste à étudier le comportement asymptotique de la solution de l'équation de Manakov PMD lorsque ce petit paramètre tend vers zéro. En généralisant la théorie de l'Approximation-Diffusion au cadre de la dimension infinie, on a montré que la dynamique asymptotique est donnée par une équation aux dérivées partielles stochastiques dirigée par un mouvement brownien de dimension trois. Dans une seconde partie, nous proposons un schéma de différences finies de type Crank Nicolson pour cette équation pour lequel nous obtenons un ordre de convergence en probabilité d'ordre 1/2. La discrétisation du bruit doit être implicite afin d'obtenir un schéma conservatif et stable. Enfin la dernière partie est relative à la simulation numérique de la dispersion modale de polarisation et à ses effets sur la propagation et la collision de solitons de Manakov. Dans ce cadre, on propose une méthode de réduction de variance valable pour les équations aux dérivées partielles stochastiques.