Un programme d'inversion des formes d'ondes sismiques a été développé pour estimer les propriétés d'un milieu stratifié élastique dans le domaine temps-distance. Ce programme est capable de prendre en compte des données haute~ fréquences (0-50 Hz) et des milieux finement stratifiés, ce qui en fait une méthode particulièrement adaptée à l'inversion de profils de sismique marine, y compris ceux enregistrés par des sismographes fond de mer. Le problème inverse est résolu par une tchnique de moindres carrés généralisée. Cette technique nécessite de pouvoir résoudre rapidement le problème direct et calculer efficacement la matrice des dérivées de Fréchet qui représente la sensibilité au premier ordre des champs d'ondes par rapport à de faibles modifications des paramètres élastiques du modèle. Le problème direct est résolu en utilisant une technique classique de matrices de réflexion-transmission généralisée combinée avec une intégration en nombres d'ondes discrets. Les dérivées de Fréchet sont obtenues sous forme analytique à l'aide des fonctions de Green du milieu non perturbé. Nous proposons deux approches pour calculer ces dérivées. La première consiste à les calculer directement dans le domaine temps-distance en représentant les interfaces par des arrangements réguliers de points diffractants. La deuxième approche est une nouvelle méthode de calcul basée sur une formulation analytique qui permet d'évaluer avec rapidité et précision les dérivées de Fréchet pour une perturbation des paramètres élastiques d'une couche plane d'extension infinie. La méthode d'inversion proprement dite utilise une technique de gradients conjugués dans le domaine temps-distance, en intégrant toutes les techniques développées précédemment. La méthode est illustrée par des exemples d'inversion de données synthétiques réalistes, mais aussi par des enregistrements réels de type OBS de la campagne Malis en Mer Ligure. La robustesse de la méthode a été étudiée en imposant aux données ou au milieu de propagation certaines erreurs violant les hypothèses de base de l'algorithme. La méthode permet par ailleurs de considérer un grand nombre de géométries source-récepteur et peut s'appliquer en particulier à des enregistrements de profils sismiques verticaux, ainsi qu'à des modèles présentant une forte variabilité des propriétés élastiques suivant la profondeur.

Ces travaux s'articulent autour du calcul des solutions périodiques dans les systèmes dynamiques non linéaires, au moyen de méthodes numériques de continuation. La recherche de solutions périodiques se traduit par un problème avec conditions aux limites périodiques, pour lequel nous avons implémenté deux méthodes d'approximation : - Une méthode spectrale dans le domaine fréquentiel, l'équilibrage harmonique d'ordre élevé, qui repose sur une formulation quadratique des équations. Nous proposons en outre une extension de cette méthode aux cas de non-linéarités non rationnelles. - Une méthode pseudo-spectrale dans le domaine temporel, la collocation à l'aide fonctions polynômiales par morceaux. Ces méthodes transforment le problème continu en un système d'équations algébriques non linéaires, dont les solutions sont calculées par continuation à l'aide de la méthode asymptotique numérique. L'ensemble de ces outils, complétés d'une analyse linéaire de stabilité, sont intégrés au code de calcul MANLAB. Applications : Un modèle physique non-régulier de clarinette est étudié en détail : à partir de la branche de solutions statiques et ses bifurcations, on calcule les différentes branches de solutions périodiques, ainsi que leur stabilité et leurs bifurcations. Ce modèle est ensuite adapté au cas du saxophone, pour lequel on intègre une caractérisation acoustique expérimentale, afin de mieux tenir compte de la géométrie complexe de l'instrument. Enfin, nous étudions un modèle physique simplifié de violon, avec une non-régularité liée frottement de Coulomb. / Periodic solutions of nonlinear dynamical systems are the focus of this work. We compute periodic solutions through a BVP formulation, solved with two numerical methods: - a spectral method, in the frequency domain: the hogh-order Harmonic Balance Method, using a quadratic formulation of the original equations. We also propose an extension to nonrational nonlinearities. - a pseudo-spectral method, in the time domain : the arthogonal collocation at Gauss point, with piece-wise polynomial interpolation. Both methods lead to a system of nonlinear algebraic equations, and its solutions are computed by a continuation algorithm : the Asymptotic Numerical Method. These methods are embeded in the numerical package MANLAB, together with a linear stability analysis. Application We then apply these methods to physical models of several instruments : a clarinet, a saxophone, and a violin. The clarinet model contains a non-smooth contact between the reed and the mouthpiece. The study focuses on the evolution of frequency, loudness, and spectrum along the branch of periodic solutions when varying the mouth pressure. The saxophone model is very similar, but an experimental characterization of the bore is used in that case. Finally, the violin model with a non-smooth Coulomb contact law and a simplified resonator is studied, showing the variety of models that can be treated using this method.

Cette thèse est consacrée au développement d’un outil de pré-dimensionnement par éléments finis pour l’estimation de l’endommagement par fatigue polycyclique de structures linéaires sous chargements multiaxiaux et stationnaires gaussiens. L’état de contraintes atteint dans ces structures étant aléatoire, il devient nécessaire de raisonner en terme de statistique et l’approche spectrale s’avère particulièrement adaptée pour cette situation. Dans ce travail, les méthodes spectrales sont améliorées par la prise en compte des largeurs de bande des spectres de réponses des structures soumises à des chargements nonproportionnels et de moyennes non nulles. Le critère d’endommagement de Sines est retenu. L’étude numérique de la distribution de l’endommagement résultant des incertitudes sur les paramètres matériaux et des applications dans le domaine de l’optimisation des structures sont également abordées. / This dissertation is devoted to the development of a finite element pre-designing tool for the estimation of high cycle fatigue damage for linear structure subject to multiaxial and Gaussian stationary loads. Given the random nature of the stress state, it becomes necessary to look at this problem from a statistical point of view and spectral methods appear to be the best approach. In this work, spectral methods are improved by taking into account the response spectral bandwidths of structures under non proportional and non zero mean loads; the Sines’ damage criterion is employed. The numerical analysis of the fatigue damage distribution resulting from the uncertainties of the material parameters and applications in the field of structural optimisation are also treated.

La quantification et la propagation des incertitudes dans les modèles physiques apparaissent comme des voies essentielles vers l'amélioration de la prédiction de leur réponse. Le développement d'outils de modélisation des incertitudes et d'estimation de leur impact sur la réponse d'un modèle a constitué un axe de recherche privilégié dans de nombreux domaines scientifiques. Cette dernière décennie, un intérêt croissant a été porté à des méthodes numériques basées sur une vision fonctionnelle des incertitudes. Ces méthodes, couramment baptisées ``méthodes spectrales stochastiques'', sont issues d'un mariage fructueux de l'analyse fonctionnelle et de la théorie des probabilités.<br /><br />Reposant sur des bases mathématiques fortes, les méthodes spectrales de type Galerkin semblent constituer une voie prometteuse pour l'obtention de prédictions numériques fiables de la réponse de modèles régis par des équations aux dérivées partielles stochastiques (EDPS). Plusieurs inconvénients freinent cependant l'utilisation de ces techniques et leur transfert vers des applications de grande taille : le temps de calcul, les capacités de stockage mémoire requises et le caractère ``intrusif'', nécessitant une bonne connaissance des équations régissant le modèle et l'élaboration de solveurs spécifiques à une classe de problèmes donnée. Un premier volet de mes travaux de recherche a consisté à proposer une stratégie de résolution alternative tentant de lever ces inconvénients. L'approche proposée, baptisée méthode de décomposition spectrale généralisée, s'apparente à une technique de réduction de modèle a priori. Elle consiste à rechercher une décomposition spectrale optimale de la solution sur une base réduite de fonctions, sans connaître la solution a priori. <br /><br />Un deuxième volet de mes activités a porté sur le développement d'une méthode de résolution d'EDPS pour le cas où l'aléa porte sur la géométrie. Dans le cadre des approches spectrales stochastiques, le traitement d'aléa sur l'opérateur et le second membre est en effet un aspect aujourd'hui bien maîtrisé. Par contre, le traitement de géométrie aléatoire reste un point encore très peu abordé mais qui peut susciter un intérêt majeur dans de nombreuses applications. Mes travaux ont consisté à proposer une extension de la méthode éléments finis étendus (X-FEM) au cadre stochastique. L'avantage principal de cette approche est qu'elle permet de traiter le cas de géométries aléatoires complexes, tout en évitant les problèmes liés au maillage et à la construction d'espaces d'approximation conformes.<br /><br />Ces deux premiers volets ne concernent que l'étape de prédiction numérique, ou de propagation des incertitudes. Mes activités de recherche apportent également quelques contributions à l'étape amont de quantification des incertitudes à partir de mesures ou d'observations. Elles s'insèrent dans le cadre de récentes techniques de représentation fonctionnelle des incertitudes. Mes contributions ont notamment porté sur le développement d'algorithmes efficaces pour le calcul de ces représentations. En particulier, ces travaux ont permis la mise au point d'une méthode d'identification de géométrie aléatoire à partir d'images, fournissant une description des aléas géométriques adaptée à la simulation numérique. Une autre contribution porte sur l'identification de lois multi-modales par une technique de représentation fonctionnelle adaptée.

On s'intéresse dans cette thèse à l'étude de méthodes numériques pour les systèmes de réaction-diffusion. Tout d'abord, on étudie le schéma par régularisation du résidu et ses extrapolations; ce schéma introduit un préconditionneur en espace lors de la discrétisation en temps. On prouve la stabilité en norme usuelle et la convergence en norme d'énergie de cette méthode et on l'applique au préconditionnement de méthodes spectrales par des méthodes d'éléments finis. Cette application nécessite le calcul d'asymptotiques précises des polynômes de Legendre et de leurs extrema. On prouve aussi la convergence et l'ordre deux d'une méthode de splitting semi-discrétisée en temps pour les systèmes de réaction-diffusion, l'approximation de Peaceman-Rachford. Enfin, on applique ces méthodes à la simulation d'une équation parabolique non linéaire pour modéliser la croissance de grains et à une équation parabolique non locale venant de la mécanique statistique et modélisant les systèmes autogravitants de fermions.

Les supernovas gravitationnelles sont un sujet d'étude passionnant, aussi bien par la diversité et la richesse de la physique mise en jeu que par l'émerveillement que peut exercer l'apparition d'un nouvel objet brillant dans le ciel. Les supernovas incitent les chercheurs aussi bien que les profanes au rêve et à l'imagination. Cette thèse traite des simulations numériques des supernovas gravitationnelles, et en particulier du problème du transport des neutrinos. Les neutrinos jouent un rôle crucial dans le chauffage et le refroidissement de la matière, et donc dans l'explosion d'une supernova. Leur implémentation numérique est difficile, à cause des connaissances très détaillées nécessaires pour le traitement des interactions, mais surtout à cause de la très grande exigence en temps de calcul des équations à résoudre. Après une première partie sur le contexte physique dans lequel se trouvent les supernovas gravitationnelles actuellement, je fais le point sur leur traitement numérique, avec un accent sur les méthodes numériques présentes dans le code que j'ai utilisé tout au long de ma thèse, le code CoCoNuT. Pour le traitement des neutrinos, deux pistes sont proposées. La troisième partie de cette thèse propose de fortement simplifier le traitement des neutrinos. C'est ce qui est fait dans l'implémentation du schéma de fuite, que j'ai pu par la suite utiliser pour étudier la formation de trous noirs. Dans cette optique, je montre l'influence de l'ajout de particules exotiques (pions et hypérons) sur la formation du trou noir. La présence de pions ou d'hypérons déclenche l'effondrement plus vite que dans les simulations où leur présence n'est pas prise en compte. De plus, les hypérons présentent une transition de phase, intéressante à étudier, par exemple pour la contrepartie en ondes gravitationnelles. La deuxième piste proposée pour le traitement des neutrinos est la résolution de l'équation de Boltzmann en relativité générale. La dérivation de cette équation dans le cadre de la relativité générale, et plus précisément en formalisme 3+1, fait l'objet de la quatrième partie de cette thèse. Enfin, la cinquième partie de cette thèse introduit un nouveau code de résolution numérique de l'équation de Boltzmann. Je montre d'une part les méthodes utilisées, et d'autre part des tests numériques, qui valident le bon fonctionnement du code.

La dissipation d'énergie turbulente est un processus clef dans le milieu interstellaire (MIS) froid, non seulement pour comprendre les voies de formation des étoiles, mais aussi en tant que source d'énergie supra-thermique et donc capable d'ouvrir de nouveaux chemins chimiques pour expliquer les abondances d'espèces soumises à des barrières endoénergétiques, telles que $\mathrm{CH}^+$ et $\mathrm{SH}^+$ qui sont observées dans le MIS. Dans ce contexte, l'intermittence spatio-temporelle du taux de dissipation d'énergie joue un rôle crucial car elle conduit à une injection d'énergie supra-thermique qui peut localement être très intense par rapport au taux moyen. Les caractéristiques détaillées de la distribution spatiale et les propriétés géométriques des lieux de dissipation intense peut fournir des indications précieuses pour les modèles chimiques.Nous étudions ici ces structures à l'aide de simulations numériques directes, avec un soin sans précédent pour résoudre les processus dissipatifs numériquement. Le nombre de Mach dans le MIS diffus peut prendre des valeurs aussi bien grandes que petites et nous encadrons les effets de la compression entre deux catégories de turbulence en déclin. Dans le cas extrême incompressible, nous réalisons des simulations pseudo-spectrales de magnétohydrodynamique visqueuse et résistive, avec un éclairage particulier sur la diffusion ambipôlaire due à la vitesse relative qui existe entre les ions et les neutres. Du côté de la compressibilité extrême, nous considérons des simulations isothermes basées sur des schémas sur grille (type Godunov) qui incluent dissipation visqueuse et résistive: nous nous intéressons ici particulièrement à la dissipation numérique.Nos simulations incompressibles montrent que la diffusion ambipôlaire conduit le champ magnétique à petite échelle dans une configuration libre de force de Lorentz. En conséquence, l'échelle caractéristique du chauffage par friction ion-neutre se déplace à plus grande échelle, aux échelles inertielles de la turbulence, bien plus grandes que la taille estimée par le raisonnement dimensionnel classique. Les structures dissipatives dans nos simulations sont des feuilles cohérentes spatiallement, chacune portant une nature dissipative bien distincte (visqueuse, ohmique ou bien ambipôlaire). Nous révélons les lois statistiques qui gouvernent leurs caractéristiques et nous calculons les exposants des fonctions de structure qui quantifient l'intermittence.Nous montrons que les simulations compressibles sont sujettes à une forte dissipation numérique: dans nos expériences, à peu près la moitié de la dissipation totale est attribuée aux termes de dissipation physique, le reste est produit par le schéma numérique. Nous avons mis au point une méthode pour estimer localement l'énergie perdue dans le schéma et nous l'utilisons pour examiner la structure en feuilles du champ de dissipation comme dans les simulations incompressibles. Par exemple, bien que nous confirmons que la dissipation physique visqueuse est dominée par les feuilles de cisaillement plutôt que par les chocs, nous ne pouvons exclure que la dissipation numérique ne renverse cet équilibre si celle-ci était dominée par les chocs.Pour finir, nous examinons l'efficacité de diagnostiques observationnels variés pour tracer les structures de forte dissipation. En particulier, nous trouvons que les incréments de centroïdes de vitesse ou des paramètres de Stokes corrèle très bien sur le plan du ciel avec certaines structures de forte dissipation. Nous calculons aussi les exposants de l'intermittence mesurés pour ces mêmes traceurs et nous trouvons qu'ils s'étendent sur une large plage de valeurs. Enfin, nous mélangeons les phases dans l'espace de Fourier associé à la boîte de simulation périodique, et nous démontrons le rôle crucial de la cohérence de phase pour produire la structure filamentaire observée dans les cartes d'incréments commecelles obtenues récemment par la collaboration Planck / Turbulent energy dissipation is a key process in the cold interstellar medium (ISM), not only on the road to star formation but also as a source of suprathermal energy able to open new chemical routes, otherwise inactive at the low gas temperature. Such routes are required, though, to explain the high abundance of species such as $\mathrm{CH}^+$ and $\mathrm{SH}^+$ observed in the ISM. In this context, the space-time intermittency of energy dissipation is particularly relevant because it drives injection of suprathermal energy in the ISM locally far above the average level. The detailed characteristics of the spatial distribution and the geometrical properties of the energy dissipation rate can provide valuable inputs to chemical models. We study them here with the aid of direct numerical simulations with unprecedented dedication to resolve the dissipation processes numerically. As the sonic Mach number in the diffuse interstellar medium can take values in a wide range, we bracket the possible physics by considering two categories of decaying turbulence models. On the incompressible extreme, we perform pseudo-spectral simulations of viscous and resistive magnetohydrodynamics, with a particular emphasis on ambipolar diffusion due the ion-neutral drift. On the compressible extreme we consider grid-based (Godunov) simulations of isothermal resisitive and viscous magnetohydrodynamics where our focus is on numerical dissipation. Our incompressible simulations show that ambipolar diffusion leads to force-free magnetic fields at small scales. As a result, the typical scale of ion-neutral friction heating is displaced to large scales in the inertial range, much greater than dimensional analysis would predict. The structures of high dissipation are spatially coherent sheets, each with a single nature of dissipation (viscous, ohmic or ambipolar). We reveal their statistical scaling laws and compute their intermittency exponents. We show that compressible simulations are subject to a lot of numerical dissipation: in our set up, less than half of the total dissipation is accounted for by the physical terms, the rest is produced by the numerical scheme. We design a method to recover locally the energy lost in the scheme and we use it to examine the sheet-like structure of the dissipation field as in our incompressible simulations. We show that numerical dissipation prevents us to assess the nature of dissipative structures. For instance, although we confirm previous results that physical dissipation in shearing sheets rapidly dominates over shocks, the balance could be reversed if numerical dissipation were shock dominated. Finally, we examine the efficiency of various observational tracers to characterize the structures of high dissipation. In particular, we find that increments of molecular line centroid velocity or of polarization Stokes parameters correlate very well on the plane of the sky with specific structures of high dissipation. We also compute the intermittency exponents measured for these tracers and find they span a broad range of possible values. At last, we mix the Fourier phases to demonstrate the crucial role of coherence in producing the filamentary structure of observable maps of increments such as recently produced by the Planck collaboration.

Ces travaux s'articulent autour du calcul des solutions périodiques dans les systèmes dynamiques non linéaires, au moyen de méthodes numériques de continuation. La recherche de solutions périodiques se traduit par un problème avec conditions aux limites périodiques, pour lequel nous avons implémenté deux méthodes d'approximation : - Une méthode spectrale dans le domaine fréquentiel : l'équilibrage harmonique d'ordre élevé, qui repose sur une formulation quadratique des équations. Nous proposons en outre une formulation originale permettant d'étendre cette méthode aux cas de non-linéarités non rationnelles. - Une méthode pseudo-spectrale par éléments dans le domaine temporel : la collocation à l'aide fonctions polynômiales par morceaux. Ces méthodes transforment le problème continu en un système d'équations algébriques non linéaires, dont les solutions sont calculées par continuation à l'aide de la méthode asymptotique numérique. L'ensemble de ces outils, intégrés au code de calcul MANLAB et complétés d'une analyse linéaire de stabilité, sont alors utilisés pour l'étude des régimes périodiques d'une classe particulière de systèmes dynamiques non linéaires : les instruments de musique auto-oscillants. Un modèle physique non-régulier de clarinette est étudié en détail : à partir de la branche de solutions statiques et ses bifurcations, on calcule les différentes branches de solutions périodiques, ainsi que leur stabilité et leurs bifurcations. Ce modèle est ensuite adapté au cas du saxophone, pour lequel on intègre une caractérisation acoustique expérimentale, afin de mieux tenir compte de la géométrie complexe de l'instrument. Enfin, nous étudions un modèle physique simplifié de violon, avec une non-régularité liée frottement de Coulomb. Cette dernière application illustre ainsi la polyvalence des outils développés face aux différents types de non-régularité.

La dynamique sédimentaire estuarienne joue un rôle très important pour la qualité de l'eau, les écosystèmes et la navigation. Les estuaires macrotidaux comme la Gironde se caractérisent par la formation de régions très chargées en matière en suspension (MES), appelées zones de turbidité maximale (ZTM), qui influencent le transport et le dépôt des sédiments fins, l'envasement des chenaux, la consommation d'oxygène dissous et le devenir des polluants. L'objectif de ce travail est de comprendre la dynamique hydro-sédimentaire, particulièrement de la ZTM, dans la section fluviale, encore peu étudié, de l'estuaire de la Gironde en lien avec les facteurs de forçage environnementaux et les perturbations du système (changements hydrologiques et morphologiques naturels et anthropiques). La méthodologie de ce travail est basée sur l'analyse de 10 années de données continues de turbidité enregistrées par le réseau de surveillance MAGEST. L'exploitation de telles séries de données, assez novatrice dans les estuaires, a notamment impliqué le développement d'une méthode d'analyse basée sur la combinaison de plusieurs méthodes spectrales. Cette approche est complétée par l'analyse des profils de turbidité et de vitesse de courant lors de cycles de marée, l'analyse de séries temporelles historiques de marée et l'exploitation d'un modèle semi-analytique 2DV. La dynamique sédimentaire de l'estuaire fluvial est d'abord détaillée à toutes les échelles de temps représentatives. A l'échelle de temps intratidale, la distribution verticale des MES et des courants, en deux points d'une même section transversale, a permis de détailler les mécanismes de transport sédimentaire. Les flux particulaires résiduels (totaux, advection, pompage tidal) ainsi estimés pour plusieurs conditions hydrologiques, démontrent le contrôle du pompage tidal sur les flux de MES lors de l'étiage. A l'échelle de temps subtidale, la réponse de la ZTM aux fluctuations hydrologiques (crues, périodes d'augmentation ou diminution continue du débit, variabilité inter-anuelle) est analysée. Ceci a permis de définir plusieurs indicateurs hydrologiques des caractéristiques de la ZTM, qui suggèrent l'intensification de la ZTM au cours des dernières décennies en lien avec la diminution des débit. La contribution relative des facteurs de forçage à la variabilité de la turbidité a été quantifiée pour différentes régions estuariennes et échelles de temps (saisonnière et plurianuelle). L'application de la méthodologique développée à l'estuaire de la Loire, qui dispose de séries de données similaires (réseau SYVEL), a permis de généraliser ces résultats. Enfin, l'effet des changements pluri-décennaux hydrologiques et morphologiques sur la propagation de la marée et la dynamique sédimentaire est détaillé dans la Garonne tidale. Il ressort une amplification du marnage et de la asymétrie de la marée au cours des six dernières décennies, principalement liée aux changements morphologiques naturels dans la Gironde en aval, les extractions de granulat et le changement de régime hydrologique. L'implémentation d'un modèle semi-analytique a permis de vérifier ces résultats et d'analyser leurs implications sur les concentrations de MES et la limite amont de la ZTM. / Estuarine suspended sediment dynamics play an important role in water quality, ecosystems and navigation. The formation of regions of high suspended sediments (SS) concentrations, called turbidity maximum zones (TMZ), is a characteristic feature of macrotidal estuaries, such as the Gironde. The TMZ influences the transport and deposition of fine sediments, channel siltation, oxygen conditions and the particulate transport of pollutants. This work aims to understand the hydro-sedimentary dynamics of the fluvial Gironde estuary, still poorly studied, in relation with environmental forcings and system perturbations (natural and anthropic hydrological and morphological changes). The methodology of this work is based on the analysis of the 10-years continuous time series of turbidity recorded by the MAGEST monitoring network. The exploitation of such time series, quite innovative in estuaries, required the development of an analysis method based on the combination of spectral techniques. This approach is completed by the analysis of turbidity and current velocity profiles over tidal cycles, the analysis of historical tide time series, and the exploitation of a 2DV semi-analytical model. First, SS dynamics of the fluvial Gironde is detailed at all representative time scales. At the intratidal time scale, the mechanisms of SS transport were described from the vertical depth of SSC and current velocities at two points of the same section. Residual fluxes (total, advective and tidal pumping), estimated for different hydrological conditions, demonstrated the control of tidal pumping on SS fluxes during periods of low river flow. At the subtidal time scale, the TMZ response to hydrological fluctuations (floods, periods of continuous river flow increase and decrease, interannual changes) was analyzed. Hydrological indicators of the TMZ features were thus defined, which suggest the TMZ intensification over the last decades in relation to the river flow decrease. The relative contributions of environmental forcings to the turbidity variability were quantified for different estuarine regions and time scales (seasonal and multiannual). The application of the same methodology to the Loire estuary, which counts on equivalent time series (SYVEL network), allowed the generalization of these results. Finally, the impact of pluri-decades hydrological and morphological changes on tidal propagation and suspended sediments dynamics is detailed in the tidal Garonne. Both tidal range and asymmetry appear to be amplified over the last six decades, mainly due to natural changes of the down Gironde, gravel extraction in the tidal Garonne and hydrological regime shifts. The implementation of an idealized model allowed verifying such results and analyzing their implications for SS concentrations and the upper TMZ limit.

Face au déluge actuel de données principalement non structurées, les graphes ont démontré, dans une variété de domaines scientifiques, leur importance croissante comme language abstrait pour décrire des interactions complexes entre des objets complexes. L’un des principaux défis posés par l’étude de ces réseaux est l’inférence de propriétés macroscopiques à grande échelle, affectant un grand nombre d’objets ou d’agents, sur la seule base des interactions microscopiquesqu’entretiennent leurs constituants élémentaires. La physique statistique, créée précisément dans le but d’obtenir les lois macroscopiques de la thermodynamique à partir d’un modèle idéal de particules en interaction, fournit une intuition décisive dans l’étude des réseaux complexes.Dans cette thèse, nous utilisons des méthodes issues de la physique statistique des systèmes désordonnés pour mettre au point et analyser de nouveaux algorithmes d’inférence sur les graphes. Nous nous concentrons sur les méthodes spectrales, utilisant certains vecteurs propres de matrices bien choisies, et sur les graphes parcimonieux, qui contiennent une faible quantité d’information. Nous développons une théorie originale de l’inférence spectrale, fondée sur une relaxation de l’optimisation de certaines énergies libres en champ moyen. Notre approche est donc entièrement probabiliste, et diffère considérablement des motivations plus classiques fondées sur l’optimisation d’une fonction de coût. Nous illustrons l’efficacité de notre approchesur différents problèmes, dont la détection de communautés, la classification non supervisée à partir de similarités mesurées aléatoirement, et la complétion de matrices. / In an era of unprecedented deluge of (mostly unstructured) data, graphs are proving more and more useful, across the sciences, as a flexible abstraction to capture complex relationships between complex objects. One of the main challenges arising in the study of such networks is the inference of macroscopic, large-scale properties affecting a large number of objects, based solely on he microscopic interactions between their elementary constituents. Statistical physics, precisely created to recover the macroscopic laws of thermodynamics from an idealized model of interacting particles, provides significant insight to tackle such complex networks.In this dissertation, we use methods derived from the statistical physics of disordered systems to design and study new algorithms for inference on graphs. Our focus is on spectral methods, based on certain eigenvectors of carefully chosen matrices, and sparse graphs, containing only a small amount of information. We develop an original theory of spectral inference based on a relaxation of various meanfield free energy optimizations. Our approach is therefore fully probabilistic, and contrasts with more traditional motivations based on the optimization of a cost function. We illustrate the efficiency of our approach on various problems, including community detection, randomized similarity-based clustering, and matrix completion.