Dans ce projet, nous étudions certaines classes d équations non-linéaires aux différences pour toutes les valeurs non-négatives admissibles de paramètres et de conditions initiales. Dans la première partie, nous présentons quelques définitions initiales et nécessaires de la stabilité dans la littérature des équations aux différences. Nous exprimons également plusieurs résultats connus et quelques théorèmes qui seront utiles dans notre recherche à la suite. Dans le deuxièllle chapitre, on étudie l'intervalle invariant, le caractère des sellli-cycles, la stabilité globale, et le "boundedness" de l'équation aux différences [Formule mathématique]. Dans la deuxième partie de ce travail, nous étudions la famille de modèles de la population, d'ordre supérieur, de l'équation logistique non-autonorne [Formule mathématique]. En particulier, le comporternent périodique, l'attractivité des solutions et la stabilité de solutions sont examinés en détail. Ce modèle représente les différentes croissances de population telles que la plupart des saumons et "spruce budworms" ayant assez de nourriture (feuillage). L'impact du caractère saisonnier , qui est indubitablenlent présent dans la croissance de la population, sur le comporternent des solutions (croissance de la population) est également envisagé. Enfin, nous présentons un modèle déterministe de l'infection par le VIH (virus d'immunodéficience hurnaine) en présence de la trithérapie. Puis la stabilité globale asymptotique de l'équilibre "disease-free" et l'équilibre endémique sont étudiés pour le modèle continu. En outre, le modèle est implicitement estimé par la méthode des différences finies, un système d'équations aux différences, afin de résoudre le système d 'IVP (problème aux valeurs initiales) . La dynamique du modèle estimé est complètement déterminée par une quantité de seuil R appelée le "basic reproduction number". Lorsque le nombre associé à la reproduction est inférieur à l'unité, le résultat indique que l'infection par le VIH peut être éliminée de la personne infectée. Un équilibre stable endémique existe lorsque le nombre associé à la reproduction est supérieur à l'unité (conduisant à la persistance et l'existence du VIH au sein de la personne infectée, puis de la communauté) .

L'objectif des travaux présentés dans cette thèse est d'explorer de nouvelles méthodes numériques efficaces et robustes pour la résolution par la méthode des éléments finis des problèmes d'écoulements des fluides viscoélastiques pour des nombres de Weissenberg élevés. La motivation de cette étude provient en grande partie au fait que les fluides viscoélastiques sont à la base de nombreuses applications industrielles, notamment dans l'industrie des polymères, l'industrie des pâtes de papiers, l'industrie alimentaire, etc... Ces écoulements sont régis par un système d'équations aux dérivées partielles fortement non linéaire composé des équations de Navier Stokes (ou Stokes pour les écoulements lents) et d'une équation constitutive modélisant la contribution du polymère. Les rhéologues ont mis à la disposition de la communauté scientifique plusieurs modèles. Dans ce travail, afin de pouvoir effectuer des comparaisons avec les travaux effectués par d'autres chercheurs déjà disponibles dans la vaste littérature du domaine, nos simulations numériques ont été faites pour les modèles de Oldroyd-B, PTT et Giesekus. La méthodologie utilisée est toutefois assez générale. Une des difficultés majeures rencontrées dans la simulation numérique des fluides viscoélastiques est la perte de convergence des algorithmes pour des problèmes à faibles nombres de Weissenberg. Ainsi, dans l'optique d'améliorer la convergence des schémas numériques, nous nous sommes proposés de développer et d'implémenter de nouveaux algorithmes en nous basant sur la nouvelle formulation en logarithmique du tenseur de conformation. Cette formulation a été utilisée pour la première fois dans le contexte des fluides viscoélastiques par Fattal et Kupferman et al. (21). Par la suite, en nous inspirant de Coronado et al. (35), nous avons introduit deux autres formulations logarithmiques. La première est une variante de celle proposée par Coronado et al. (35) avec un traitement spécial pour le terme convectif. La deuxième formulation permet une linéarisation complète des équations pour la résolution par la méthode de Newton. Les deux méthodes que nous proposons ont l'avantage d'être relativement simple à mettre en oeuvre. Pour tester la robustesse de ces nouvelles formulations, quelques calculs ont été effectués pour le problème d'écoulement autour d'un cylindre et l'inévitable problème de la contraction 4:1 sur des maillages non structurés anisotropes obtenus par la stratégie générale de remaillage adaptatif développée au GIREF. Nous avons ensuite jugé utile d'appliquer cette nouvelle formulation logarithmique pour la résolution de quelques problèmes de surfaces libres. Nous avons étudié la présence des écoulements secondaires et leur influence sur la morphologie des écoulements complexes des problèmes multi-couches dans une filière. En effet, nous avons montré qu'il possible de résoudre des problèmes multi-couches dans une filière tridimensionnelle en résolvant un problème bidimensionnel de surface libre instationnaire seulement sur une section. Nous avons également traité un problème de cisaillement d'une goutte dans une matrice fluide. Il s'agit ici d'un écoulement régi par la tension superficielle. Le calcul des sufaces libres est rendu possible grâce à la méthode des surfaces de niveau (level set method). Cette méthode requiert la résolution d'une équation de transport pour trouver la nouvelle position de l'interface. L'interface constitue une zone de transition de largeur 2 ? des paramètres rhéologiques des deux fluides. Ces paramètres sont régularisés à l'interface pour éviter de gérer des discontinuités. Une technique de réinitialisation a été mise en oeuvre pour remédier à d'éventuels problèmes de conservation de masse. Ces techniques sont falicitées par l'introduction de la méthode de remaillage adaptatif. Celle-ci permet de concentrer les éléments dans l'interface permettant ainsi de faire les calculs avec une meilleure précision. Pour la résolution des équations hyperboliques, notamment l'équation constitutive, de transport et de régularisation de l'interface, nous avons utilisé la méthode de stabilisation de type SUPG. Il résulte de ces développements une méthode de résolution nettement plus performante que les méthodes classiques.

La thèse présentée ici contribue à l'étude des rotations et des déplacements d'un corps solide dans l'espace. Elle se compose de trois articles qui traitent de la modélisation statistique de données directionnelles, rotationnelles et de déplacements. Le premier article propose un nouveau modèle pour les données directionnelles qui sont des vecteurs unitaires dans R p [matrice de rotation R, p x p]. On modélise la dispersion et l'orientation de vecteurs unitaires aléatoires par une nouvelle distribution indexée par des paramètres d'échelle et de localisation et on suggère d'utiliser la méthode des moments pour estimer ses paramètres. Le deuxième article propose une nouvelle distribution pour les déplacements dans R p [matrice de rotation R, p x p]. On traite de modèles stochastiques pour les variations d'une séquence de déplacements autour d'une valeur moyenne. Sous quelques spécifications pour les erreurs expérimentales, on montre que l'estimateur de la rotation moyenne associée au déplacement moyen est meilleur comparé à celui obtenu par des modèles marginaux qui utilisent la rotation seule. Le troisième article étudie le mouvement relatif d'un corps solide par rapport à un autre qui s'obtient comme une rotation selon un axe fixe et une translation le long de cet axe. On suggère un modèle pour ce type de mouvement avec une valeur prédite qui dépend de la position et de l'orientation de l'axe de rotation. On suppose que le déplacement observé est égal à la valeur modale perturbée par des erreurs qui suivent une distribution suggérée dans le deuxième article. On montre, par une étude de Monte-Carlo, que le modèle proposé est plus performant que le modèle marginal de rotation ou de translation pour l'estimation de l'axe de rotation. Ainsi, le traitement combiné de données de rotation et de translation permet d'obtenir des estimation plus précises.