À l’heure actuelle, les écoulements à deux fluides non miscibles jouent un rôle très important dans plusieurs domaines, que ça soit en science ou en ingénierie. Leur complexité est tellement élevée que les modèles actuels ne permettent de résoudre que des cas particuliers ou simplifiés avec un degré de précision qui demeurent souvent plutôt modeste. Une nouvelle approche numérique parait être une nécessité pour capturer la complexité physique du phénomène. Pour ce faire nous avons besoin d’outils robustes. Au niveau de l’interface de séparation entre les deux fluides non miscibles, les variables physiques sont discontinues, ce qui pose un défi majeur dans la description des variables et des conditions aux limites à l’interface. Le fait que les densités et les viscosités de chaque fluide soient différentes de part et d’autre de l’interface donne naissance à des défauts et des impuretés dans le champ des vitesses, ce qu’on appelle une discontinuité faible. Pour sa part, l’existence de la force de tension superficielle au niveau de l’interface crée une discontinuité sur le champ de pression, ce qu’on appelle une discontinuité forte. Un autre grand problème se pose au niveau de l’étude numérique du problème, où les méthodes numériques classiques ont une précision assez limitée dans ce genre de situation. L’objectif de ce travail est de fournir une étude complète de la dynamique de l’interface entre deux fluides non miscibles à l’aide d’outils mathématiques, physiques et numériques robustes. D’abord, une étude analytique du problème a été faite où l’équation de Navier-Stokes et les conditions de saut sur les variables physiques au niveau de l’interface de séparation entre les fluides ont été prouvées en détail. Pour traiter les discontinuités, nous avons discrétisé nos variables à l’aide de la méthode XFEM. Dû aux larges distorsions rencontrées dans ce genre d’écoulement, nous avons utilisé l’approche Eulérienne, pour corriger les oscillations des solutions dues aux choix du système de coordonnées nous avons utilisé les techniques de stabilisation SUPG/PSPG. Le traitement de la courbure des interfaces K été fait à l’aide de l’opérateur Laplace Beltrami et le suivi d’interface à l’aide de la méthode ¨Level-set¨. Pour le traitement des conditions de saut au niveau de l’interface la méthode Nitsche est développée dans différents contextes. Après avoir développé un modèle physique et mathématique dans les premières parties de notre travail, nous avons fait une étude numérique à l’aide de la plateforme de calcul FEniCS, qui est une plateforme de développement en langage C++ avec une interface Python. Un code de calcul a été développé dans le cas des écoulements de deux fluides non miscibles avec les modèles physiques et les outils mathématiques développés dans les sections précédentes. / The two-phase flow problems have an important role in the multitude of domains in science and engineering. Their complexity is so high that the actual models can solve only particular or simplified cases with a certain degree of precision. A new approach is a necessity to understand the evolution of new ideas and the physical complexity in this kind of flow, to contribute to the study of this field. A good study requires solid and robust tools to have performing results and a maximum of efficacy. At the interface of separation between the two immiscible fluids, the physical parameters are discontinuous, which gives us difficulties for the description of the physical variables at the interface and boundary conditions. The fact that the density and the viscosity are discontinuous at the interface creates kinks in the velocity, which represent a weak discontinuity. The existence of the surface tension at the interface create a discontinuity for the pressure field, it represents a strong discontinuity. The main objective of this work is to make a complete study based on strong and robust physical, mathematical and numerical tools. A strong combination, capable of capturing the physical aspect of the interface between the two fluids with a very good precision. Building such a robust, cost effective and accurate numerical model is challenging and requires lots of efforts and a multidisciplinary knowledge in mathematics, physics and computer science. First, an analytical study was made where the one fluid model of the Navier-Stokes equation was proved from Newton’s laws and jump conditions at the interface was proved and detailed analytically. To treat the problem of discontinuity, we used the XFEM method to discretize our discontinuous variables. Due to the large distortion encountered in this kind of fluid mechanic problems, we are going to use the Eulerian approach, and to correct the oscillation of solutions we will use the SUPG/PSPG stabilization technic. The treatment of the interface curvature k was done with the Laplace Beltrami operator and the interface tracking with the Level-set method. To treat the jump conditions with a very sharp precision we used the Nitsche’s method, developed in different cases. After building a strong mathematical and physical model in the first parts of our work, we did a numerical study using the FEniCS computational platform, which is a platform of computational development based on C++ with a Python interface. A numerical code was developed in this study, in the case of two-phase flow problem, based on the previous mathematical and physical models detailed in previous sections.

La microscopie holographique utilise des techniques d’interférométrie pour mesurer les changements de volume des neurones et des astrocytes avec une précision sans précédent. Un défi important serait de relier les changements de phase mesurés aux changements de volume du neurone et plus encore de relier l’étendue de ces changements de volumes à certaines propriétés des neurones comme le niveau d’activité des cotransporteurs cation-chlorure (CCC) et certaines propriétés biomécaniques des membranes. L’objectif à plus long terme est d’utiliser des changements de phase pour détecter des modifications dans la réponse volumique des neurones à un choc osmotique par exemple, modifications qui pourraient éventuellement permettre de détecter des pathologies. Pour comprendre l’information que l’on peut tirer des mesures expérimentales, il est important de comprendre le lien entre différentes variables : force de la pompe Na⁺ – K⁺ATPase, la perméabilité de la membrane à l’eau, les propriétés biomécaniques de la membrane et les changements de phase observés par l’expérimentateur. Pour y arriver, nous aborderons quelques notions sur les systèmes dynamiques, plus précisément nous utiliserons les Equations Différentielles Ordinaires (E.D.O.) afin d’éffectuer la modélisation mathématique du phénomène illustrant la variation du volume de la membrane cellulaire, ainsi que les variations des quantités de K⁺, Na⁺ et Cl⁻, qui constituent la principale composition ionique des astrocytes, qui sont les cellules étudiées dans ce projet. Dans ce même régistre de rappel mathématique sur les systèmes dynamiques, nous parlerons des bifurcations, pour lesquelles nous décrirons quand et comment est ce qu’elles apparaîssent tout en les illustrant par des exemples, ceci dans l’optique de se préparer à une meilleure compréhension des résultats à venir après l’étude de notre modèle, puisqu’on espère y observer des bifurcations. Nous serons ainsi amenés à étudier profondémént le système d’E.D.O obtenu, notamment la recherche des points d’équilibre et leurs comportements dans l’espace des phases, voir s’il ya lieu des points de bifurcation et leurs interprétations pour la cellule concernée. Le but visé étant d’obtenir des bifurcations, ce qui expliquerait le dysfonctionnement des astrocytes, et expliquerait certainement l’origine de certaines maladies maladies neurodégénératives ; nous verrons finalement après étude du modèle qu’il n’existe pas de bifurcation, néanmoins la simplicité du modèle utilisé ouvre des portes à de futurs projets plus complexes qui permettront peut-être d’atteindre les objectifs visés. / The holographic microscopy uses interferometry techniques for measuring changes in volume of neurons with an unprecedented accuracy. A major challenge is to relate the measured phase changes with the neuron volume changes and more to relate the extent of these changes volumes to certain properties of neurons such as the activity level of Cation-Chloride Cotransporter (CCC) and some biomechanical properties membranes. The longer term objective is the use of phase changes for detecting changes in the density response of neurons to an osmotic shock which could possibly allow the detection of many kind of pathologies. To understand the information that can be derived from experimental measurements, it is important to understand the relationship between different variables: force pump Na⁺ – K⁺ ATPase, membrane permeability of water, biomechanical properties of the membranes and the phase changes observed by the experimenter. To achieve this, we need some dynamical system skills, we will use the Ordinary Differential Equations (E.D.O) in order to perform the mathematical modeling of the phenomenon illustrating the variation of the membrane volume, as well as the variations in quantities of K⁺, Na⁺ and Cl⁻, which constitute the main ionic composition of astrocytes, which are the cells studied in this project. In this mathematical recall on dynamical systems, we will talk about the bifurcations for a better understanding of the incoming results since we are expecting bifurcations for our model. We will study deeply the E.D.O. system obtained including the search of equilibrium points and their behavior in the phase space, and we will see if there are bifurcations and what is their meaning. The aim being to obtain bifurcations, which would explain the dysfunction of the astrocytes, and would certainly explain the origin of certain neurodegenerative diseases; we will finally see, after studying the model, that there is no bifurcation, nevertheless the simplicity of the model used opens doors to more complex future projects that will perhaps achieve the desired objectives.

L'adaptation de maillage peut être très utile dans l'application de la méthode des éléments finis à certains problèmes physico-mathématiques. L'adaptation est une procédure qui nécessite, entre autres, des algorithmes efficaces de calcul des dérivées de fonctions de type Pk ou Qk- Des méthodes d'estimation des dérivées ont donc été élaborées pour résoudre ce problème. Ce travail consiste à expliquer au lecteur ce qu'est exactement la récupération des dérivées et à dresser une liste des méthodes classiques qui sont utilisées de nos jours. Nous développons aussi une nouvelle méthode de récupération pour les fonctions de type Pi qui est, à plusieurs égards, beaucoup plus performante que les méthodes classiques. Un aspect intéressant de cette nouvelle méthode est qu'elle est globale et qu'elle n'a ainsi pas besoin de traitement spécial sur le bord des maillages.

Ce travail est consacré à l'étude du couplage des équations de Navier-Stokes et de la chaleur posées dans un domaine polygonal non convexe. Après avoir analysé le comportement singulier des solutions de ces équations près des coins du domaine considéré, nous présentons une formulation mixte duale de ces équations basée sur l'introduction de deux nouvelles inconnues : a qui est le tenseur gradient de la vitesse et le champ vectoriel ^ qui désigne le gradient de la température. Ensuite, en considérant une famille de maillage Th du domaine fi, nous analysons une méthode d'éléments finis mixte duale basée sur cette dernière formulation en utilisant l'élément fini de Raviart-Thomas de plus bas degré pour approximer les nouvelles inconnues cr et <^ sur chaque triangle K de la triangulation Th. Tandis que les variables n, p et T seront approximées par des polynômes de degré zéro sur chaque triangle K. En particulier, nous montrons que l'on peut retrouver l'ordre de convergence quasi-optimal si le maillage est raffiné suivant certaines règles qui sont essentiellement celles introduites par G. Raugel et basées sur le fait que les solutions sont régulières dans des espaces de Sobolev à poids. Nous discutons les aspects d'implémentation de la méthode d'éléments finis mixte duale raffinée pour ces équations en utilisant un algorithme de point fixe combiné à une formulation hybride de deux systèmes issus du découplage du problème discret : l'un correspondant aux équations de Navier-Stokes et l'autre à l'équation de la chaleur. Nos résultats numériques obtenus, en plus de confirmer l'ordre de convergence optimal pour un problème test posé dans un domaine polygonal non convexe, sont tout-à-fait comparables avec ceux existants dans la littérature pour la convection naturelle dans une cavité carrée.

La recherche sur les écoulements à surface libre turbulents, et sur le cisaillement turbulent spécifiquement progresse depuis les années 1970. Les équations décrivant le mouvement d'agitation turbulente sont bien connues. Elles sont reliées directement aux équations non moyennées de continuité et de conservation de la quantité de mouvement qui forment le système d'équations de Navier-Stokes ainsi que les différentes équations de corrélations du mouvement turbulent. Pour un écoulement turbulent, ce système d'équations présente un problème de fermeture relié à la présence des termes associés aux contraintes de Reynolds. L'échelle du phénomène et plusieurs processus importants du mouvement turbulent ne peuvent ainsi être résolus de façon exacte, et un degré d'approximation devient donc nécessaire. La recherche sur le phénomène de la turbulence a permis de développer plusieurs concepts et approches reliés, notamment, à la viscosité turbulente, à l'isotropie locale de la turbulence et aux équations de transport des termes fluctuants afin de résoudre le problème de fermeture du système d'équations de Navier-Stokes. Ces approches présentent évidemment certaines limitations. De fait, la turbulence est fortement dépendante des conditions aux frontières par rapport à l'écoulement moyen et ne permet pas de prédétermination de la forme des distributions des vitesses moyennes de l'écoulement dans le cadre de la recherche d'une solution approximative. Ce travail de recherche présente un modèle hydrodynamique turbulent 2D vertical, en vue de résoudre une problématique d'écoulement à surface libre à partir de la méthode des éléments finis de frontière ainsi qu'une démarche en laboratoire, afin de caractériser les fluctuations de vitesses et de pression d'un écoulement turbulent à l'aide d'une méthode de vélocimétrie par imagerie de particules. Pour notre modèle numérique, notre approche procède à partir d'une formulation intégrale des équations moyennées de conservation de la quantité de mouvement et de continuité et l'équation associée des corrélations de vitesses turbulentes. En spécifiant les conditions aux limites pour les vitesses, la corrélation et les tractions à la frontière, le calcul du champ de vitesses peut être déterminé en termes des variables primaires à la frontière. Dans ces conditions, il n'est pas requis de discrétiser le domaine à l'aide d'un maillage ou de cellules ce qui pennet une utilisation plus efficiente des ressources et du temps de calcul. Des applications à des cas d'écoulements simples permettent de démontrer la validité du modèle bien que l'imposition des contraintes à la frontière peut être difficile à interpréter et que dans certains cas, l'intégration numérique pour des fonctions singulières peut être hasardeuse. À tenne, notre modèle prédit correctement les solutions analytiques de l'équation de Navier-Stokes pour un écoulement de Couette et ce, sans imposer de distribution à priori des variables primaires de l'écoulement sur le domaine, dans la recherche d'une solution numérique. Le système d'équations de Navier-Stokes, pour un écoulement turbulent, peut donc être envisagé à partir d'une formulation intégrale sur la frontière uniquement. Cette situation étant en accord et cohérente avec le comportement physique attendu du phénomène où les propriétés d'un écoulement turbulent sont directement reliées aux conditions aux limites et que l'action des contraintes fluctuantes est inter-reliée avec l'écoulement moyen. D'autre part, la démarche expérimentale à l'aide de vélocimétrie par imagerie de particules a permis une meilleure compréhension des spécificités des écoulements turbulents à surface libre, à Reynolds modéré. Notre analyse sur les conditions turbulentes à la frontière de l'écoulement tend à corroborer que les corrélations de vitesses turbulentes sont en équilibre avec le gradient de pression et bien que l'écoulement moyen puisse être stationnaire, il apparaît que les corrélations varient continuellement dans le temps selon une action dissipative, ce qui signifie que le gradient de pression est en constante évolution pour équilibrer le système. De plus, il apparaît que les conditions de pression moyenne dans un modèle, ne peuvent représenter adéquatement la réalité physique de l'action de la pression instantanée selon l'agitation turbulente. Enfin, remarquons que l'étude de la turbulence en laboratoire peut être entreprise de manière adéquate, aujourd'hui, avec des moyens limités malgré l'échelle du phénomène. Mots clés: équations de Navier-Stokes, écoulement turbulent, corrélations de vitesses, éléments finis de frontière, vélocimétrie par imagerie de particules.

Dans cette thèse, nous nous intéressons à la construction d'éléments finis d'ordre élevé adaptés aux maillages hybrides, pour la résolution de systèmes hyperboliques linéaires en régimes harmonique et temporel. L'accent est plus particulièrement porté sur la construction d'éléments pyramidaux. On étudie trois formulations pour lesquelles on cherche des éléments finis "optimaux" au sens de la convergence dans la norme de l'espace considéré pour la formulation. Pour les formulations H^1 et H(rot), on construit des éléments finis "optimaux" nodaux et hp. Les matrices élémentaires sont évaluées grâce à des formules de quadrature adaptées et des estimations d'erreur sont effectuées pour vérifier la convergence des éléments optimaux construits. Pour la formulation discontinue LDG (Local Discontinuous Galerkin), on présente des éléments utilisant des fonctions de base orthogonales permettant de mettre au point une construction de la matrice de masse et un produit matrice-vecteur rapides. Dans le cas des trois formulations, on étudie les propriétés numériques des éléments construits, on vérifie que l'on retrouve bien numériquement la convergence théorique et on compare nos éléments avec d'autres éléments trouvés dans la littérature. Finalement, on présente des expériences numériques en 3D avec l'équation des ondes ou de Helmholtz, et les équations de Maxwell dans le cas des régimes temporels et harmoniques. On montre ainsi l'efficacité des maillages hybrides par rapport aux maillages purement tétraédriques ou aux maillages hexaédriques obtenus en découpant chaque tétraèdre d'un maillage purement tétraédrique en quatre hexaèdres.

Nous nous intéressons, dans cette thèse, à l'étude des méthodes d'extrapolation polynômiales et à l'application de ces méthodes dans l'accélération de méthodes de points fixes pour des problèmes donnés. L'avantage de ces méthodes d'extrapolation est qu'elles utilisent uniquement une suite de vecteurs qui n'est pas forcément convergente, ou qui converge très lentement pour créer une nouvelle suite pouvant admettreune convergence quadratique. Le développement de méthodes cycliques permet, deplus, de limiter le coût de calculs et de stockage. Nous appliquons ces méthodes à la résolution des équations de Navier-Stokes stationnaires et incompressibles, à la résolution de la formulation Kohn-Sham de l'équation de Schrödinger et à la résolution d'équations elliptiques utilisant des méthodes multigrilles. Dans tous les cas, l'efficacité des méthodes d'extrapolation a été montrée.Nous montrons que lorsqu'elles sont appliquées à la résolution de systèmes linéaires, les méthodes d'extrapolation sont comparables aux méthodes de sous espaces de Krylov. En particulier, nous montrons l'équivalence entre la méthode MMPE et CMRH. Nous nous intéressons enfin, à la parallélisation de la méthode CMRH sur des processeurs à mémoire distribuée et à la recherche de préconditionneurs efficaces pour cette même méthode.

Dans cette thèse, nous nous intéressons à la résolution des <br />équations de Maxwell en régime fréquentiel, afin de calculer<br />précisément la signature radar de cibles diverses. Pour avoir<br />une grande précision nécessaire pour des expérience de grande taille,<br /> nous utilisons des méthodes d'ordre élevé.<br /><br />Dans le cas scalaire, les éléments finis spectraux hexaédriques<br />avec condensation de masse, permettent d'obtenir un produit matrice vecteur <br />rapide et peux coûteux en stockage. Dans le cas vectoriel, les hexaèdres<br />de la première famille ne réalisent pas la condensation de masse, mais on peut<br />écrire un algorithme rapide de produit matrice-vecteur. Des résultats<br />numériques 3-D montrent la performance de l'algorithme proposé.<br /><br />Nous traitons également le cas où la géométrie présente<br />une symétrie de révolution. On est alors ramenés à une succession<br />de problèmes 2-D indépendants.<br />Nous proposons une méthode éléments finis d'ordre élevé <br />couplée à des équations intégrales d'ordre élevé.

L'objet de cette thèse concerne la proposition de modèles paramétriques des processus de branchement uni et multi-types. Nous mettons en valeur l’intérêt de la théorie des processus de branchement et du développement nécessaire des différents outils et de concepts propres à plusieurs domaines. Pour cela, nous commençons par rappeler quelques définitions et résultats de la théorie des processus de branchement uni et multi-types, et ce en temps discret comme en temps continu. On se consacre par la suite au développement méthodologique de ces modèles.Dans la deuxième partie de ce mémoire, nous étudions seulement l'évolution d’une seule population en temps continu, et présentons quelques familles de lois paramétriques, associées à des processus de branchement homogènes particuliers. Des méthodes récursives de calcul, ainsi que des propriétés pertinentes, concernant ces distributions de probabilité, sont dérivées des fonctions génératrices satisfaisant certaines équations aux dérivées partielles linéaires précisés. Les familles proposées seront utiles à la modélisation de systèmes plus cohérents en dynamique de populations, puisqu'on y montre que les hypothèses usuelles de distributions de Poisson ne peuvent être argumentées.Dans la troisième partie, nous étudions le comportement de l'évolution de plusieurs populations en interactions. Nous y présentons aussi des modèles paramétriques de lois, associés à des processus de branchement multi-types en temps continu et homogènes en temps. Nous considérons ensuite un modèle particulier, où une population ``mère donneuse" autonome alimente en individus K populations filles, qui sont, elles, en interaction. Ce modèle est bien adapté à l'étude des systèmes dynamiques des populations en interaction qui reste à la fois simple, mais riche en variétés de comportement. L'étude du système multi-types se fait via l'évolution des fonctions génératrices de la loi multidimensionnelles des effectifs. Pour cela, utilisant les équations différentielles ordinaires et aux dérivées partielles, nous établissons les équations implicites des distributions temporelles et multidimensionnelles, et discutons des méthodes analytiques ou numériques de leur résolution. Nous développons ensuite des exemples de modèles et en particulier celui concernant 3 et 4 populations.En conclusion, nous argumentons la pertinence de cette approche, et l’interprétation des paramètres, qui sont d'un grand intérêt pour le développement de méthodes d'inférence statistique, pour de nombreux domaines d'applications. / This thesis aims to propose parametric models for single and multi-type branching processes. The importance of the theory of branching processes is pointed out. Hence, developing various tools and specific concepts in several domains is important for applications. For those purpose, we recall some definitions and results of the single-and-multi-type branching processes theory in discrete and continuous case. Afterward, we focus on the methodological development of those models.In the second part, the evolution of a single population in the continuous case has been studied. Then, some parametric distribution families associated to particular branching mechanisms are explored. Recursive computational procedure and relevant properties concerning the associted probability distributions are derived from generating functions that satisfy specified linear partial differential equations. The suggested families are useful for the modeling of systems that are more coherent with population dynamics, contrarily to the usual hypothesis of Poisson distributions, that cannot be argued.In the third part, the evolution of different populations with interaction is explored. Similarly, some parametric models of homogeneous multi-type branching processes in continuous time are proposed. Afterwards, we consider a particular model where an autonomous donor parent population feeds in individuals, K types progeny populations that interacts. This model is well adapted to the study of dynamical systems of populations in interaction. This simple model, but has a rich variety of behaviors.The study of such systems is also done regarding the evolution of generating functions of multidimensional ndividual countrings. To achievea such study, ordinary and partial differential equations are used to establish the implicit equations of temporal and multidimensional distributions. Analytical and numerical methods for equation resolution are then discussed, and examples of particular models are developed.In conclusion, the relevancy of this approach is argumed, censidering parameters interpretation in the development of inference methods for the various applied domains.

Dans cette thèse, nous étudions différents aspects mathématiques et numériques des équations de Gross-Pitaevskii et de Schrödinger non linéaire. Nous commençons (chapitre 1) par introduire différents modèles à partir des systèmes physiques que sont les condensats de Bose-Einstein et les impulsions lumineuses dans les fibres optiques. Cette modélisation conduit aux équations aux dérivées partielles stochastiques suivantes : l'équation de Gross-Pitaevskii stochastique et l'équation de Schrödinger non linéaire avec dispersion aléatoire. Ensuite, dans le second chapitre, nous nous intéressons au problème de l'existence et l'unicité d'une solution de ces équations. On montre notamment que le problème de Cauchy a une solution pour l'équation de Gross-Pitaevskii stochastique avec rotation grâce à la construction de la solution associée au problème. Nous abordons ensuite dans le troisième chapitre le problème du calcul des états stationnaires pour l'équation de Gross-Pitaevskii. Nous développons une méthode pseudo-spectrale de discrétisation du Continuous Normalized Gradient Flow, associée à une résolution itérative préconditionnée des sous-espaces de Krylov. Le quatrième chapitre concerne l'étude de schémas pseudo-spectraux pour la dynamique de l'équation de Gross-Pitaevskii et de Schrödinger non linéaire. On procède à une étude numérique de ces schémas (schéma de splitting de Lie et de Strang, ainsi qu'un schéma de relaxation). De plus, on analyse le schéma de Lie dans le cadre de l'équation de Schrödinger non linéaire avec dispersion aléatoire. Finalement, nous présentons, dans le cinquième chapitre, une boîte à outils Matlab (GPELab) développée dans le but de fournir les méthodes numériques que nous avons étudiées / The aim of this Thesis is to study various mathematical and numerical aspects related to the Gross-Pitaevskii and nonlinear Schrödinger equations. We begin (chapter 1) by introducing a few models starting from the physics of Bose-Einstein condensates and optical fibers. This naturally leads to introducing a stochastic Gross-Pitaevskii equation and a nonlinear Schrödinger equation with random dispersion. Next, in the second chapter, we analyze the existence and uniqueness problem for these two equations. We prove that the Cauchy problem admits a solution for the stochastic Gross-Pitaevskii equation with a rotational term by constructing the solution associated with the linear. The third chapter is concerned with the computation of stationary states for the Gross-Pitaevskii equation. We develop a pseudo-spectral approximation scheme for the Continuous Normalized Gradient Flow formulation, combined with preconditioned Krylov subspace methods. This original approach leads to the robust and efficient computation of ground states for fast rotations and strong nonlinearities. In the fourth chapter, we consider some pseudo-spectral schemes for computing the dynamics of the Gross-Pitaevskii and nonlinear Schrödinger equations. These schemes (the Lie's and Strang's splitting schemes and the relaxation scheme) are numerically studied. Moreover, we proceed to a rigorous numerical analysis of the Lie scheme for the associated stochastic PDEs. Finally, we present in the fifth chapter a Matlab toolbox (called GPELab) that provides computational solutions based on the schemes previously introduced in the Thesis