La segmentation en objets d'une image consiste à extraire de l'image des régions d'intérêt suivant un critère défini. Nous segmentons l'image par un algorithme de contours actifs dans le cadre d'une approche variationnelle. Partant d'un contour initial quelconque, le contour actif évolue, suivant une équation aux dérivées partielles. L'équation d'évolution du contour actif est déduite de la dérivation du critère. Au vu de la dépendance du critère à la région considérée, la dérivation du critère par rapport à la région n'est pas aisée. Nous utilisons des outils de dérivation empruntés à l'optimisation de domaine: les gradients de forme.<br />La contribution de cette thèse réside dans l'élaboration et l'étude de différents descripteurs de région. Pour chaque critère, nous calculons la dérivée du critère à l'aide des gradients de forme, et en déduisons l'équation d'évolution du contour actif.<br />Le premier descripteur définit un a priori géométrique sans contrainte paramétrique: il minimise la distance du contour actif à un contour de référence. Nous l'avons appliqué à la déformation de courbe, la segmentation et le suivi de cible.<br />Le deuxième descripteur caractérise le mouvement de l'objet par un modèle de mouvement. Le critère associé définit conjointement une région et son mouvement sur plusieurs images consécutives. Nous avons appliqué ce critère à l'estimation et la segmentation conjointe du mouvement et au suivi d'objets en mouvement.

La détection des ondes gravitationnelles par les détecteurs interférométriques terrestres, tels que VIRGO ou LIGO, et par la mission spatiale LISA sera fortement facilitée par la connaissance théorique à priori du signal. On s'intéresse dans cette thèse à l'étude d'une des sources de rayonnement gravitationnel les plus importantes, à savoir les systèmes binaires d'objets compacts. Plus précisément, on considère, dans le cadre de la relativité générale, les dernières orbites de la phase de quasi-équilibre. Elles permettent, d'une part, de fournir des données initiales aussi réalistes que possible pour la phase de coalescence et, d'autre part, d'apporter de nombreuses informations sur les objets compacts émetteurs.<br /><br />Un effort est fait pour améliorer, rendre ces données initiales les plus réalistes possible d'un point de vue astrophysique. Nous avons ainsi construits les premières séquences de binaires d'étoiles de quarks étranges, et ce pour différentes équations d'état. Contrairement au cas d'étoiles à neutrons polytropiques, la séquence se termine par une instabilité dynamique. Nous avons également calculé des configurations de binaires d'étoiles à neutrons à l'aide d'une théorie sans onde allant au delà de l'approximation communément admise de métrique spatiale conformément plate. Les solutions obtenues devraient être plus précises et de meilleures conditions initiales que celles réalisées jusqu'alors. Nous avons enfin étudié, pour des systèmes d'un seul trou noir puis des trous noirs binaires, l'influence de conditions de bords aux horizons provenant du formalisme des horizons isolés et regroupnt des ingrédients de quasi-équilibre.

Soit $\mathcal{W}(d)$ un $d$-tissu non singulier du plan implicitement présenté par une équation différentielle $F(x,y,y')=0$, et de connexion associée $(E,\nabla)$. De nouveaux invariants de $\mathcal{W}(d)$ sont mis à jour ; en particulier, on montre que $(E,\nabla)$ est entièrement déterminé par la connaissance d'une $1$-forme fondamentale et du polynôme de linéarisation du tissu.\esp Nous indiquons également comment la courbure de la connexion rend compte de la linéarisation du tissu. En étudiant la trace de la courbure de la connexion, on montre que le fibré déterminant de $(E,\nabla)$ est isomorphe au produit tensoriel des fibrés en droites associés aux $3$-tissus extraits. Nous donnons ensuite une caractérisation géométrique des tissus de trace nulle, en généralisant la construction de l'hexagone de Thomsen. En outre, on présente un procédé explicite de détermination du rang de $\mathcal{W}(d)$ pour $d$ quelconque, à partir des seuls coefficients de $F$. En application, nous retrouvons des résultats connus en géométrie des tissus, et indiquons des perspectives nouvelles, notamment pour l'étude des tissus exceptionnels.

La modélisation des propulseurs ioniques de type SPT pose de nombreux <br />problèmes dans le domaine du transport des particules chargées. Nous nous <br />intéressons à deux de ces problèmes, à savoir le transport des électrons et <br />le calcul du potentiel électrique.<br /><br />Le transport des électrons résulte de l'influence conjuguée des champs <br />(électrique et magnétique) établis dans la cavité du propulseur et des <br />collisions des électrons (dans la cavité et avec la paroi limitant celle-ci). <br />Nous avons participé au développement d'un modèle SHE (Spherical Harmonics <br />Expansion) qui résulte d'une analyse asymptotique de l'équation de Boltzmann <br />munie de conditions de réflexion aux bords. Ce modèle permet d'approcher la <br />fonction de distribution en énergie des électrons en résolvant une <br />équation de diffusion dans un espace \{position, énergie\}. Plus précisément, <br />nous avons étendu une démarche existante au cas où les collisions en volume <br />(excitation, ionisation) et les collisions inélastiques à la paroi <br />(attachement et émission secondaire) sont prises en compte. Enfin, nous <br />avons écrit un code de résolution du modèle SHE, dont les résultats ont <br />été comparés avec ceux d'une méthode de Monte Carlo. <br /><br />\vspace*{1mm}<br />Dans un deuxième temps, nous avons étudié le calcul du potentiel électrique. <br />La présence du champ magnétique impose d'écrire le courant d'électrons sous <br />la forme ${\cal J}=\sigma \nabla W$<br /> où W est le potentiel électrique et le tenseur de conductivité $\sigma$<br />est fortement anisotrope compte tenu des grandeurs physiques en jeu dans <br />le SPT. Pour résoudre $\mbox{div }{\cal J}(x,y)=S(x,y)$, <br />nous avons implémenté une méthode de volumes finis <br />sur maillage cartésien permettant de résoudre ce problème elliptique <br />anisotrope, et nous avons vérifié qu'elle échouait lorsque le rapport <br />d'anisotropie devenait grand. Aussi nous avons développé une méthode de <br />paramétrisation, qui consiste à extrapoler la solution d'un problème <br />anisotrope à l'aide d'une suite de problèmes isotropes. Cette méthode a <br />donné des résultats encourageants pour de forts rapports d'anisotropie, <br />et devrait nous permettre d'atteindre des cas réels.

On s'est intéressé dans ce travail à la résolution numérique des équations d'Euler linéarisées autour d'un écoulement stationnaire, subsonique et assez régulier. Dans le but d'obtenir des matrices symétriques dans ces équations, et in fine, une équation d'équilibre énergétique, nous considérons la linéarisation d'une forme symétrique des équations d'Euler tridimensionnelles. Nous proposons un schéma non-diffusif de type Galerkin discontinu en domaine temporel (GDDT) s'appuyant sur une formulation centrée-élément avec des ux numériques totalement centrés et un schéma en temps explicite de type saute-mouton, ce qui permet d'obtenir une approximation sans dissipation et fournit une estimation précise des variations de l'énergie aéroacoustique. En effet, dans le cas général de la linéarisation autour d'un écoulement non-uniforme, il existe une équation d'équilibre énergétique au niveau continu que nous vérifions au niveau discret. Nous montrons qu'il existe un terme source discret permettant de conserver exactement l'énergie ce qui permet de prouver la stabilité de notre schéma. Ainsi, notre schéma non-diffusif de type GDDT fournit un outil précis pour contrôler des phénomènes tel que les instabilités de Kelvin-Helmholtz. Nous illustrons la capacité de notre méthode aussi bien sur plusieurs cas tests académiques que sur différentes configurations complexes grâce à une implémentation parallèle.

Dans ce travail, nous nous intéressons à des problèmes issus de la Mécanique des Fluides et plus particulièrement au cas des aérosols (ou sprays, c'est-à-dire un ensemble de particules en suspension dans un fluide environnant). Les phénomènes physiques mis en jeu sont modélisés par des équations aux dérivées partielles (EDP). La phase continue (fluide environnant) est décrite par des équations issues de la mécanique des milieux continus de type Navier-Stokes ou Euler. La phase dispersée est décrite par une équation cinétique de type Boltzmann.<br /><br />Le premier résultat que nous présentons est consacré à l'étude mathématique d'un couplage entre une équation cinétique de type Vlasov et les équations d'Euler isentropiques. Ces équations modélisent un spray fin. Nous démontrons l'existence en temps petit d'une solution régulière pour le couplage Vlasov-Euler isentropique.<br /><br />Ensuite, nous présentons les équations précises relatives à la modélisation des collisions, coalescences et fragmentations dans un spray.<br /><br />Nous décrivons par la suite la simulation numérique du couplage fluide-cinétique dans un code industriel (Commissariat à l'Énergie Atomique), en particulier l'ajout des phénomènes de collisions.<br /><br />Un deuxième modèle de fragmentation est également présenté. Ce modèle est plus pertinent dans les cas où les particules de la phase dispersée ont un grand nombre de Weber.<br /><br />Enfin, nous présentons un résultat concernant une estimation explicite de trou spectral pour l'opérateur de Boltzmann avec potentiels durs linéarisé, et pour l'opérateur de Landau avec potentiels durs linéarisé.

Dans cette thèse, nous nous intéressons à des méthodes numériques pour un modèle d'écoulements incompressibles et miscibles ayant des application dans l'hydrogéologie et l'ingénierie pétrolière. Nous étudions et analysons un schéma numérique combinant une méthode d'éléments finis mixtes (EFM) et une méthode des volumes finis (VF) pour approcher le système couplé entre une équation elliptique (pression-vitesse) et une équation de convection-diffusion-réaction (concentration). Le schéma VF considérée est de type "vertex centred" semi-implicite en temps : explicite pour la convection et implicite pour la diffusion. On utilise un schéma de Godunov pour approcher le terme convectif et une approximation élément fini P1 pour le terme de diffusion. Nous montrons des résultats de stabilité L≂ estimations BV et le principe du maximum discret sous une condition CFL appropriée. Ensuite, nous montrons la convergence de la solution approchée obtenue par le schéma combiné EFM-VF vers la solution du problème couplé. La démonstration de la convergence se fait en plusieurs étapes : premièrement, on déduit la convergence forte de la solution approchée de la concentration dans L2(Q), en utilisant la stabilité L≂, les estimations BV et des arguments de compacité. Dans l'étape suivante, on étudie le schéma découplé EFM, en donnant des résultats de convergence pour la pression et la vitesse. Enfin, le processus de convergence de la solution approchée du schéma combiné EFM-VF vers la solution exacte est obtenu par passage à la limite et par unicité de solution pour le problème continu. Des simulations numériques académiques et réalistes pour des problèmes bidimensionnels confirment la stabilité et l'efficacité du schéma combiné. Enfin, nous étudions des estimateurs d'erreur a posteriori de type résiduel pour une équation de convection-diffusion-réaction discrétisée par un schéma VF "vertex centred" semi-implicite en temps. Nous introduisons deux sortes d'indicateurs. Le premier est local en temps et en espace et constitue un outil efficace pour l'adaptation du maillage à chaque pas de temps. Le second est global en espace mais local en temps et peut être utilisé pour l'adaptation en temps. Nous montrons que l'estimateur est une borne supérieure de l'erreur. Des résultats numériques d'adaptations de maillage sont présentés et montrent l'efficacité de la méthode. La partie logiciels de ce travail porte sur deux volets. Le premier a permis de réaliser un code de calcul 2D, MFlow, écrit en C++, pour la résolution du système des écoulements miscibles considérés dans cette thèse. Le second volet concerne la collaboration avec un groupe de chercheurs pour l'élaboration de la plate-forme Homogenizer++ réalisée dans le cadre du GDR MoMaS (http://momas.univ-lyon1.fr/).

Les étoiles de la branche asymptotique (AGBS) sont surtout connues pour deux faits historiques: - le "mystère des éléments s": riches en neutrons et possédant une durée de vie inférieure à celle des AGBS, ces éléments sont néanmoins présents dans les enveloppes de ces étoiles. De plus, les abondances de ces éléments "s" sont mal reproduites par les modèles existants. - le "mystère des étoiles carbonées": ces étoiles, dont les abondances vérifient C/O>1, sont majoritairement de faible luminosité. C'est un drainage (DUP) qui permet aux AGBS de devenir carbonées: l'enveloppe convective (CE) descend là où brûle 1H et prélève du Carbone, produit plus bas par la combustion de 4He et amené là par des languettes convectives récurrentes, les pulses thermiques. Or, les profondeurs des DUP modélisés pour des AGBS de faibles masses n'expliquent la forte proportion de Carbone observée. Ces deux mystères conduisent à identifier deux paramètres majeurs de l'évolution des abondances: la perte de masse et le DUP. Le premier contrôle la fin de la phase AGB et la température à la base de la CE, à masse stellaire fixée. Le deuxième maîtrise l'apport de protons entre les couches de 1H et de 4He et la modification des abondances de surface. Et c'est de la physique non standard qui est sous-jacente à ces paramètres: couplage gaz-poussière, pulsations, ondes de chocs, dégénérescence, convection avec gradient de composition, mélange non standard (non thermiquement induit)... Plusieurs études ayant été faites, tant sur les mécanismes de production que sur le paramétrage de la perte de masse, ce travail est centré sur l'influence du traitement du DUP sur la nucléosynthèse des AGBS de masses intermédiaires. La nucléosynthèse étant connue grâce à une coopération étroite entre méthodes observationnelles et théoriques, cette thèse a été menée sur les deux fronts. Cela rend possible "l'éclairage" de chacune des méthodes par les forces et les faiblesses de l'autre et donc de mieux évaluer la fiabilité d'une mesure réalisée sur Terre lorsqu'on veut l'interpréter en terme d'abondances au centre de l'étoile. La partie observationnelle de ce travail est focalisée sur l'amélioration de la précision de rapports isotopiques délicats à interpréter dans des enveloppes particulières: rapports d'éléments lourds de sources brillantes ou rapports des CNO dans des sources présentant des anomalies d'abondances. L'interprétation des observations a requis deux modélisations avancées dont l'une, faisant partie intégrante de cette thèse, porte sur l'émission de l'enveloppe, et l'autre, portant sur le processus "s", a fait l'objet d'une collaboration. Dans le domaine théorique, le code d'évolution stellaire de Grenoble a été optimisé pour les AGBS afin d'étudier l'influence de la modélisation physique et du traitement numérique du DUP sur la nucléosynthèse. Le cœur de la modélisation d'une étoile, l'équation d'état (EOS), a été réécrite et une grande partie du traitement numérique du code a été revue. La physique du mélange non standard a été incorporée et "l'arsenal" numérique nécessaire à la mise en œuvre de cette physique a fait l'objet d'une collaboration. Ce travail a permis la précision des caractéristiques d'enveloppes d'étoiles J (riches en 13C), la détermination de 35Cl / 37Cl dans IRC+10216 et la production d'outils de calcul performants pour les étoiles (même de faibles masses) et leurs enveloppes moléculaires. Mais il a surtout montré que l'incorporation de mélange non standard permet à la fois d'augmenter la profondeur du DUP et de produire une poche de 13C propice au processus "s". Ce fait, ajouté aux difficultés pour concilier les résultats d'observations des étoiles J avec les modèles, tend à réorienter le futur des modèles d'AGBS. Il faudrait disposer de modèle d'AGBS plus généraux, notamment pour les faibles masses, cohérents avec la physique et le numérique de la phase évolutive qui les précède, elle qui est naturellement modélisée avec du mélange non standard. En prime, les changements effectués dans le code et notamment l'élaboration d'une nouvelle EOS, ont permis de modéliser le flash de l'Helium et de produire des tracés d'étoiles pré séquence principale (PMS) jusqu'à 0.1 masse solaire concurrençant les codes les plus spécialisés.

La physique de la lumière ainsi que les outils géométriques pour la Conception Assistée par Ordinateur sont à la base des logiciels de simulation des phénomènes lumineux pour la fabrication des systèmes optiques. Ce n'est pas sans difficulté que les industriels conçoivent ces logiciels dont un des principaux handicaps est que les simulations sont très coûteuses en temps. L'objectif principal de ce travail est de rechercher et développer des algorithmes de calcul plus performants. Dans un premier temps, on décrit précisément le modèle du transport des photons dans ce contexte, composé de l'équation de Boltzmann accompagné de conditions de bord, et qui, dans le cas de milieux homogènes par morceaux, se ramène à l'équation de radiosité. Ensuite, on présente les outils géométriques utilisés dans le modeleur hybride CSG (Constructive Solid Geometry) et BRep (Boundary Representation) ainsi que les algorithmes de base nécessaires à la recherche d'intersections entre des demi-droites et des objets géométriques. Puis, un tour d'horizon des méthodes d'accélération des calculs en radiométrie par localisation spatiale est présenté. En tenant compte des contraintes industrielles, une telle méthode d'accélération est alors adaptée au contexte puis développée dans un environnement logiciel existant. Des expérimentations numériques montrent l'efficacité des nouvelles bibliothèques. Enfin, une étude théorique des complexités en temps et en mémoire liées aux méthodes de localisation spatiale, faisant intervenir les sommes de Minkowski d'ensembles géométriques, débouche sur une stratégie consistant à minimiser la complexité en temps pour choisir les paramètres de localisation.

Cette thèse a pour objet le comportement asymptotique de solutions globales d'Equations aux Dérivées Partielles d'évolution paraboliques semilinéaires. A travers deux exemples distincts, on traite de la convergence en temps des solutions vers des solutions particulières (ondes progressives, solutions autosimilaires). Dans un premier temps, on étudie la stabilité asymptotique des ondes progressives à symétrie sphérique dans une équation de réaction-diffusion scalaire avec non-linéarité bistable. On obtient un résultat de stabilité pour de petites perturbations radiales et d'instabilité pour des perturbations quelconques. Dans un deuxième temps, on calcule un développement asymptotique jusqu'au second ordre des solutions, à donnée initiale petite, de Navier-Stokes et de Navier-Stokes Coriolis dans une bande tridimensionnelle. On montre notamment que leur comportement asymptotique est régi par le tourbillon d'Oseen. On généralise ensuite ce résultat à toute solution globale uniformément bornée en temps, sans aucune hypothèse de petitesse. Enfin, on met en évidence de telles solutions pour l'équation de Navier-Stokes Coriolis pour les fluides tournants dans le cas d'une rotation suffisamment rapide.