Στην παρούσα διατριβή αναπτύσσεται η γραμμικοποιημένη μέθοδος των πεπερασμένων διαφορών. Η μέθοδος αυτή αποτελεί νέα αριθμητική τεχνική επιλύσεως των διαφορικών εξισώσεων δευτέρας τάξεως, κανονικών και με μερικές παραγώγους και εφαρμόζεται σε εσωτερικές και εξωτερικές ροές. Συγκεκριμένα, στο πρώτο κεφάλαιο γίνεται εισαγωγή στις μεθόδους διακριτοποιήσεως μιας διαφορικής εξισώσεως και στον τρόπο επίλυσης ενός γραμμικού συστήματος αλγεβρικών εξισώσεων. Στο δεύτερο κεφάλαιο εκτίθεται η μέθοδος των Allen και Southwell για τη λύση της εξισώσεως της στροβιλότητας προκειμένου περί δισδιάστατης ροής και η τροποποίηση και βελτίωση της μεθόδου από τον Dennis. Εν συνεχεία περιγράφεται η περαιτέρω βελτίωση της μεθόδου, διά της οποίας ο πίνακας του συστήματος των συντελεστών των αγνώστων καθίσταται διαγωνίως υπέρτερος. Στο τρίτο κεφάλαιο η βελτιωμένη μέθοδος επεκτείνεται στις τρεις διαστάσεις με μεταβλητούς συντελεστές και εφαρμόζεται στην λύση αντιστοίχου προβλήματος ροής ιδανικού ρευστού εντός σωλήνα. Στο τέταρτο κεφάλαιο η παρούσα μέθοδος γενικεύεται μεταβλητούς συντελεστές. Στο πέμπτο κεφάλαιο επιλύεται η εξίσωση της στροβιλότητας στις τρεις διαστάσεις σε σύστημα κυλινδρικών συντεταγμένων και ο αλγόριθμος εφαρμόζεται στο πρόβλημα του von Kármán. Τέλος, στο έκτο κεφάλαιο εξετάζεται η περίπτωση της αριθμητικής εξομοίωσης στροβίλων Taylor προκειμένου περί ροής ιξώδους ρευστού, το οποίο περιέχεται εντός σφαιρικού κελύφους, η δε κίνηση του ρευστού προκαλείται από την διαφορική περιστροφή των δύο σφαιρών γύρω από την κοινή διάμετρό τους. / In the present dissertation a numerical technique is developed on the solution of ordinary and/or partial differential equations of second order with variable coefficients. In particular, in the first chapter we make a general introduction to the discretisation methods and to the methods of solution of a linear system of algebraic equations. In the second chapter the method of Allen and Southwell on the solution of the vorticity equation in two dimensional flow is presented as well as its modification and improvement by Dennis. The contribution of the present work is on the further improvement of the previous method, such that the matrix of the coefficients of the unknowns becomes diagonally dominant and its extension to three dimensional problems. The consecutive steps of the method (Linearized Finite Difference Method, LFDM) are worked out on the ordinary differential equation with variable coefficients. Next, the method is applied to the solution of two problems from Fluid Mechanics. In the third chapter the improved method is extended to three dimensions with variable coefficients. Accordingly the method is applied to the solution of the equations of motion of a perfect fluid moving in a straight tube. In the fourth chapter the present method is generalized with variable coefficients. In the fifth chapter the three dimensional equation of the vorticity expressed in cylindrical coordinates is solved and the results are applied to the von Kármán problem. Finally in the sixth chapter we use this method to simulate Taylor vortices in spherical annular flow in the presence of heat, the motion of the fluid being induced by the differencial rotation of two spheres about their common diameter.

Les modèles de champ moyen en QED apparaissent naturellement dans la modélisation du nuage électronique des atomes lourds. Cette modélisation joue un rôle croissant en physique et chimie quantique, les effets relativistes ne pouvant pas être négligés pour ces atomes. En physique quantique relativiste, le vide est un milieu polarisable, susceptible de réagir à la présence de champ électromagnétique.On se place dans le cadre du modèle variationnel de Bogoliubov-Dirac-Fock (BDF) qui est une approximation de champ moyen de la QED sans photon (en particulier, les interactions considérées sont purement électrostatiques).Il est à noter que pour donner un sens au modèle BDF, il est nécessaire d'introduire une régularisation ultra-violette. Il se produit un phénomène de renormalisation de charge due à la polarisation du vide : la charge de l'électron observée dépend de la charge « nue » de l'électron et du paramètre de régularisation. On étudie rigoureusement ce phénomène ainsi que le problème de la renormalisation de la masse. Cette dernière est en lien avec l'existence d'un état fondamental pour le système d'un électron dans le vide, en l'absence de tout champ extérieur. En revanche, on montre l'absence de minimiseurs dans le cas de deux électrons.Enfin, on exhibe des points critiques de l'énergie BDF, interprétés comme des états excités du vide. On met en évidence le positronium, système métastable d'un électron et de son antiparticule le positron, ainsi que le dipositronium, molécule métastable constituée de deux électrons et de deux positrons.Les méthodes utilisées sont variationnelles (concentration-compacité, lemme de Borwein et Preiss). / In QED, mean-field models appear in the modelling of the electron clouds of heavy atoms. This modelling plays a increasing role in physics and in quantum chemistry: relativistic effects cannot be neglected in these atoms. In relativistic quantum physics the vacuum is a polarizable medium that can react to the presence of an electromagnetic field.We consider the so-called Bogoliubov-Dirac-Fock (BDF) model, a variational model which is a mean-field approximation of no-photon QED (in particular the interactions are purely electrostatic).We point out that an ultraviolet regularisation is necessary to properly define the BDF model. The vacuum polarisation leads to a \emph{renormalisation} phenomenon, the "observed" charge of the electron depends on its "bare" charge and the regularisation parameter. We rigorously study both the problem of charge renormalisation and mass renormalisation. This last one is linked to the existence of ground state in the case of an electron in the vacuum, without any external field. In contrast, we show there is no ground state in the case of two electrons.Finally we exhibit some critical points of the BDF energy which are interpreted as vacuum excited states. In particular, there are the positronium (a metastable system constituted by an electron and its antiparticle called the positron) and the dipositronium (a metastable molecule constituted by two electrons and two positrons).The methods that we use are variational (concentration-compactness, Borwein and Preiss's Lemma).

Dans le cadre du contrôle quantique bilinéaire, cette thèse étudie la possibilité de retrouver l'Hamiltonien et/ou le moment dipolaire à l'aide de mesures d'observables pour un ensemble grand de contrôles. Si l'implémentation du contrôle fait intervenir des bruits alors les mesures prennent la forme de distributions de probabilité. Nous montrons qu'il y a toujours unicité (à des phases près) des Hamiltoniens de du moment dipolaire retrouvés. Plusieurs modèles de bruit sont étudiés: bruit discrète constant additif et multiplicatif ainsi qu'un modèle de bruit dans les phases sous forme de processus Gaussien. Les résultats théoriques sont illustrés par des implémentations numériques. / The problem of recovering the Hamiltonian and dipole moment, termed inversion, is considered in a bilinear quantum control framework. The process uses as inputs some measurable quantities (observables) for each admissible control. If the implementation of the control is noisy the data available is only in the form of probability laws of the measured observable. Nevertheless it is proved that the inversion process still has unique solutions (up to phase factors). Several models of noise are considered including the discrete noise model, the multiplicative amplitude noise model and a Gaussian process phase model. Both theoretical and numerical results are established.

On considère des germes de champs de vecteurs holomorphes singuliers trimimensionnels, appelés noeud-cols doublement résonants. Ces champs de vecteurs correspondent à des systèmes différentiels bidimensionnels à singularité irrégulière, et dont la partie linéaire possède deux valeurs propres non-nulles opposées. Ce type de singularité apparait par exemple à l'infini dans les équations de Painlevé PI,...,PV après compactification à poids de l'espace, pour des valeurs génériques des paramètres. Depuis Boutroux, l'étude de ces singularités a générè de nombreux travaux de recherche. Récemment, plusieurs auteurs ont fournis des informations nouvelles, en étudiant notamment les phénomènes de Stokes non-linéaires et quasi-linéaires associés, en donnant des formules de connexion. Les coefficients de Stokes quasi-linéaires sont invariants sous l'action de changement de coordonnées analytiques locaux, mais ne forment pas un système complet d'invariants analytiques. L'objectif de ce travail de thèse est de fournir une classification analytique générale et complète des noeud-cols doublement résonants. L'idée pour cela est d'adapter les travaux de Martinet et Ramis, généralisés ensuite par Stolovitch. Dans une première partie on fournit une classification formelle, i.e. sous l'action de changements de coordonnées formels, en exhibant des formes normales formelles. Dans un second temps, on étudiera l'existence de normalisations sectorielles (analytiques sur des secteurs), généralisant ainsi un théorème de Hukuhara-Kimura-Matuda. Enfin, on étudiera les recollements entre ces applications normalisantes dans les domaines d'intersections: c'est ce que l'on appellera les difféomorphismes de Stokes. Il s'agira là d'étudier des isotropies sectorielles de la forme normale. On verra que la donnée d'une forme normale formelle et d'un couple de difféomorphismes de Stokes fournira un système complet d'invariants analytiques. Enfin, dans une quatrième et dernière partie, nous calculerons certains de ces invariants pour la singularité irrégulière à l'infini de la première équation de Painlevé. / We consider germs of analytic singular vector fields in dimension three, called doubly-resonant saddle-nodes. These vector fields correspond to irregular two-dimensional systems with a pair of two opposite non-zero eigenvalues. This king of singularity appears for instance at infinity in Painlevé equations PI,...,PV, after a weighted compactifcation, for generic values of the parameters. Since Boutroux, the study of these singularities has generated many researches. Recently, several authors provided new informations, by studying for instance the associated non-linear and quasi-lineair Stokes phenomenas and by giving connection formulas. Quasi-linéaire Stokes coefficients are invariant under local analytic change of coordinates, but do not form a complete set of invariants for analytic classification. The goal of this work is to provide a complete analytic classification of doubly-resonant saddle-nodes. The idea for this is to adapt the works of Martinet and Ramis, generalized then by Stolovitch. In the first part, we give a formal classification, based on the existence on unique formal normal forms. In the second part, we prove the existence of sectorial nomalizing maps (analytic over sectors), generalizing a theorem by Hukuhara-Kimura-Matuda. In the third part, we study the Stokes diffeomorphisms, and more generaly the sectorials isotropies of the normal form. We obtain a complet set of analytic invariants. Finally, in the fourth part, we compute some of these invariants in the case of the first Painlevé equation.

Cette thèse est composée de deux parties qui sont relativement indépendantes l’un de l’autre. Dans la première partie,une autre théorie de la gravitation que l’on appelle la gravité de f(R), est étudiée. Une première analyse mathématique est discutée sur cette théorie, y compris la formulation mathématique du problème de Cauchy, la discussion sur le choix du couplage, et la formulation mathématique des équations différentielles. Ce système des équations différentielles est de quatrième ordre et très impliqué. Pour pouvoir établir l’existence locale, une série de transformations et reformulation et introduites. Elles nous amènent à une formulation que l’on l’appelle la formulation conforme augmenté. Avec cette formulation, l’existence locale est établie. La deuxième partie est consacrée à l’analyse d’un type de système non-linéaire composé des équations d’onde et équations de Klein-Gordon. Ce type de système apparaît naturellement dans de nombreux modèles physiques: le plus important, l’équation d’Einstein couplé avec un champ scalaire réel du massif et le système de la formulation conforme augmentée de la théorie de f(R). La difficulté principale est le manque de la symétrie: un des champs de vecteur de Killing conforme de l’opérateur d’onde, le champ de vecteur de scaling S := t∂ t +r∂ r, n’est pas un champ de vecteur de Killing conforme de l’opérateur de Klein Gordon. Pour franchir cette difficulté, un nouveau cadre, appelé la méthode de feuilletage hyperboloïdal, est introduit. Avec ce cadre, nous pouvons encadrer les équations d’onde et les équations de Klein-Gordon dans le même cadre. Cela nous permet d’établir un résultat d’existence globale pour les données initiales petites et localisées dans un compact. / This these is composed by two parts which are relatively independent to each other. In the first part an alternative theory of the gravitation, the so-called f(R) gravity, is studied. A first mathematical analysis is discussed on this theory, including the mathematical formulation of the Cauchy problem, the discussion on the choice of coupling, the mathematical formulation of the differential system. This system is four-order and highly involved. To establish the local well-posedness result, a series of transformations ans re-formulations is introduced and we finally arrived at a formulation, called the augmented conformal formulation with which we have managed to establish the local well-poseness theory.The second part is devoted to the analysis of a type of coupled wave and Klein-Gordon system. This kind of system arises naturally in many physical model, especially in the Einstein equation coupled with a real massive scalar field and the augmented conformal formulation of the f(R) gravity. The main difficulty to treat this type of system is the lack of symmetry: one of the conformal Killing vector filed of the linear wave operator, the scaling vector field S := t∂t+r∂r is not a conformal Killing vector field of the linear Klein-Gordon operator. To overpass this difficult, a new framework, called the hyperboloidal foliation method is introduced. With this framework we can encompass the wave equations and the Klein-Gordon equations in the same framework. This allowed us to establish a global well-posedness result for compactly supported, small amplitude initial data.

Ce travail est consacré à l’étude de la dynamique d’un système proie-prédateur de type Leslie-Gower défini par un système d’équations différentielles ordinaires (EDO) ou d’équations différentielles stochastiques (EDS), ou par des systèmes couplés d’EDO ou d’EDS. L’objectif principal est de faire l’analyse mathématique et la simulation numérique des modèles construits. Cette thèse est divisée en deux parties : La première partie est consacrée à un système proie-prédateur où les proies utilisent un refuge, le modèle est donné par un système d’équations différentielles ordinaires ou d’équations différentielles stochastiques. Le but de cette partie est d’étudier l’impact du refuge ainsi que la perturbation stochastique sur le comportement des solutions du système. Dans la deuxième partie, nous considérons un système proie-prédateur couplé en réseau. Il s’agit d’étudier comment des couplages plus ou moins forts entre plusieurs systèmes affectent l’existence et la position des points d’équilibre, et la stabilité de ces systèmes. / This work is devoted to the study of the dynamics of a predator-prey system of Leslie-Gower type defined by a system of ordinary differential equations (EDO) or stochastic differential equations (EDS), or by coupled systems of EDO or EDS. The main objective is to do mathematical analysis and numerical simulation of the models built. This thesis is divided into two parts : The first part is dedicated to a predator-prey system where the prey uses a refuge, the model is given by a system of ordinary differential equations or stochastic differential equations. The purpose of this part is to study the impact of the refuge as well as the stochastic perturbation on the behavior of the solutions of the system. In the second part, we consider a networked predator-prey system. We show that symmetric couplings speed up the convergence to a stationary distribution.

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

Dans cette thèse nous considérons l’équation des ondes amorties vectorielle sur une variété riemannienne compacte, lisse et sans bord. L’amortisseur est ici une fonction lisse allant de la variété dans l’espace des matrices hermitiennes de taille n. Les solutions de cette équation sont donc à valeurs vectorielles. Nous commençons dans un premier temps par calculer le meilleur taux de décroissance exponentiel de l’énergie en fonction du terme d’amortissement. Ceci nous permet d’obtenir une condition nécessaire et suffisante la stabilisation forte de l’équation des ondes amorties vectorielle. Nous mettons aussi en évidence l’apparition d’un phénomène de sur-amortissement haute fréquence qui n’existait pas dans le cas scalaire. Dans un second temps nous nous intéressons à la répartition asymptotique des fréquences propres de l’équation des ondes amorties vectorielle. Nous démontrons que, à un sous ensemble de densité nulle près, l’ensemble des fréquences propres est contenu dans une bande parallèle à l’axe imaginaire. La largeur de cette bande est déterminée par les exposants de Lyapunov d’un système dynamique défini à partir du coefficient d’amortissement. / In this thesis we are considering the vectorial damped wave equation on a compact and smooth Riemannian manifold without boundary. The damping term is a smooth function from the manifold to the space of Hermitian matrices of size n. The solutions of this équation are thus vectorial. We start by computing the best exponential energy decay rate of the solutions in terms of the damping term. This allows us to deduce a sufficient and necessary condition for strong stabilization of the vectorial damped wave equation. We also show the appearance of a new phenomenon of high-frequency overdamping that did not exists in the scalar case. In the second half of the thesis we look at the asymptotic distribution of eigenfrequencies of the vectorial damped wave equation. Were show that, up to a null density subset, all the eigenfrequencies are in a strip parallel to the imaginary axis. The width of this strip is determined by the Lyapunov exponents of a dynamical system defined from the damping term.

Un système est défini par les entités et leurs interrelations dans un environnement qui est déterminé par une limite arbitraire. Les systèmes complexes présentent un comportement émergent sans un contrôleur central. Les systèmes concrets désignent ceux qui sont observables dans la réalité. Un modèle nous permet de comprendre, de contrôler et de prédire le comportement du système. Un modèle différentiel à partir d'un système pourrait être compris comme une sorte de loi physique sous-jacent représenté par l'un ou d'un ensemble d'équations différentielles. Ce travail vise à étudier et mettre en œuvre des méthodes pour effectuer la modélisation des systèmes automatisée par l'ordinateur. Cette thèse pourrait être divisée en trois étapes principales, ainsi: (1) le développement d'un solveur numérique automatisé par l'ordinateur pour les équations différentielles linéaires, partielles ou ordinaires, sur la base de la formulation de matrice pour une personnalisation propre de la méthode Ritz-Galerkin; (2) la proposition d'un schème de score d'adaptation qui bénéficie du solveur numérique développé pour guider l'évolution des modèles différentiels pour les systèmes complexes concrets; (3) une implémentation préliminaire d'une application de programmation génétique pour effectuer la modélisation des systèmes automatisée par l'ordinateur. Dans la première étape, il est montré comment le solveur proposé utilise les polynômes de Jacobi orthogonaux comme base complète pour la méthode de Galerkin et comment le solveur traite des conditions auxiliaires de plusieurs types. Solutions à approximations polynomiales sont ensuite réalisés pour plusieurs types des équations différentielles partielles linéaires, y compris les problèmes hyperboliques, paraboliques et elliptiques. Dans la deuxième étape, le schème de score d'adaptation proposé est conçu pour exploiter certaines caractéristiques du solveur proposé et d'effectuer l'approximation polynômiale par morceaux afin d'évaluer les individus différentiels à partir d'une population fournie par l'algorithme évolutionnaire. Enfin, une mise en œuvre préliminaire d'une application GP est présentée et certaines questions sont discutées afin de permettre une meilleure compréhension de la modélisation des systèmes automatisée par l'ordinateur. Indications pour certains sujets prometteurs pour la continuation de futures recherches sont également abordées dans ce travail, y compris la façon d'étendre ce travail à certaines classes d'équations différentielles partielles non-linéaires. / A system is defined by its entities and their interrelations in an environment which is determined by an arbitrary boundary. Complex systems exhibit emergent behaviour without a central controller. Concrete systems designate the ones observable in reality. A model allows us to understand, to control and to predict behaviour of the system. A differential model from a system could be understood as some sort of underlying physical law depicted by either one or a set of differential equations. This work aims to investigate and implement methods to perform computer-automated system modelling. This thesis could be divided into three main stages: (1) developments of a computer-automated numerical solver for linear differential equations, partial or ordinary, based on the matrix formulation for an own customization of the Ritz-Galerkin method; (2) proposition of a fitness evaluation scheme which benefits from the developed numerical solver to guide evolution of differential models for concrete complex systems; (3) preliminary implementations of a genetic programming application to perform computer-automated system modelling. In the first stage, it is shown how the proposed solver uses Jacobi orthogonal polynomials as a complete basis for the Galerkin method and how the solver deals with auxiliary conditions of several types. Polynomial approximate solutions are achieved for several types of linear partial differential equations, including hyperbolic, parabolic and elliptic problems. In the second stage, the proposed fitness evaluation scheme is developed to exploit some characteristics from the proposed solver and to perform piecewise polynomial approximations in order to evaluate differential individuals from a given evolutionary algorithm population. Finally, a preliminary implementation of a genetic programming application is presented and some issues are discussed to enable a better understanding of computer-automated system modelling. Indications for some promising subjects for future continuation researches are also addressed here, as how to expand this work to some classes of non-linear partial differential equations.

Dans cette thèse, nous nous sommes intéressés à l'analyse locale de systèmes différentiels linéaires singulièrement perturbés et de systèmes de Pfaff complètement intégrables et multivariés à croisements normaux. De tels systèmes ont une vaste littérature et se retrouvent dans de nombreuses applications. Cependant, leur résolution symbolique est toujours à l'étude. Nos approches reposent sur l'état de l'art de la réduction formelle des systèmes linéaires singuliers d'équations différentielles ordinaires univariées (ODS). Dans le cas des systèmes différentiels linéaires singulièrement perturbés, les complications surviennent essentiellement à cause du phénomène des points tournants. Nous généralisons les notions et les algorithmes introduits pour le traitement des ODS afin de construire des solutions formelles. Les algorithmes sous-jacents sont également autonomes (par exemple la réduction de rang, la classification de la singularité, le calcul de l'indice de restriction). Dans le cas des systèmes de Pfaff, les complications proviennent de l'interdépendance des multiples sous-systèmes et de leur nature multivariée. Néanmoins, nous montrons que les invariants formels de ces systèmes peuvent être récupérés à partir d'un ODS associé, ce qui limite donc le calcul à des corps univariés. De plus, nous donnons un algorithme de réduction de rang et nous discutons des obstacles rencontrés. Outre ces deux systèmes, nous parlons des singularités apparentes des systèmes différentiels univariés dont les coefficients sont des fonctions rationnelles et du problème des valeurs propres perturbées. Les techniques développées au sein de cette thèse facilitent les généralisations d'autres algorithmes disponibles pour les systèmes différentiels univariés aux cas des systèmes bivariés ou multivariés, et aussi aux systèmes d''equations fonctionnelles. / In this thesis, we are interested in the local analysis of singularly-perturbed linear differential systems and completely integrable Pfaffian systems in several variables. Such systems have a vast literature and arise profoundly in applications. However, their symbolic resolution is still open to investigation. Our approaches rely on the state of art of formal reduction of singular linear systems of ordinary differential equations (ODS) over univariate fields. In the case of singularly-perturbed linear differential systems, the complications arise mainly from the phenomenon of turning points. We extend notions introduced for the treatment of ODS to such systems and generalize corresponding algorithms to construct formal solutions in a neighborhood of a singularity. The underlying components of the formal reduction proposed are stand-alone algorithms as well and serve different purposes (e.g. rank reduction, classification of singularities, computing restraining index). In the case of Pfaffian systems, the complications arise from the interdependence of the multiple components which constitute the former and the multivariate nature of the field within which reduction occurs. However, we show that the formal invariants of such systems can be retrieved from an associated ODS, which limits computations to univariate fields. Furthermore, we complement our work with a rank reduction algorithm and discuss the obstacles encountered. The techniques developed herein paves the way for further generalizations of algorithms available for univariate differential systems to bivariate and multivariate ones, for different types of systems of functional equations.