L’époque de réionisation (EoR) est une phase de grands changements qu’a subit l’Univers dans son premier milliard d’années. Suite à l’apparition des premières étoiles et à l’émission de photons énergétique par ces dernières, l’hydrogène a été réionisé. Cette transition a eu un impact sur la formation des galaxies. J’ai activement participé au développement d’EMMA, un code de simulation numérique ayant pour objectif d’étudier les processus à l’œuvre durant l’EoR. J’ai développé et implémenté un modèle de formation et d’évolution stellaire. Ces travaux ont contribué à la réalisation d’une simulation dédiée à l’étude de l’EoR parmi les plus grosses réalisées à l’heure actuelle. J’ai contribué au développement d’outils dédiés à l’exploration de simulations de ce type. J’ai étudié la façon dont le rayonnement s’échappe des galaxies en fonction des paramètres du modèle stellaire, et montré que les supernovæ peuvent augmenter la fraction de photons libérés. J’ai également étudié la propagation des fronts d’ionisation et montré qu’il était possible de réduire la vitesse de la lumière par trois (et ainsi diminuer le temps de calcul du transfert du rayonnement par 3), tout en conservant des résultats corrects ? / The epoch of reionization (EoR) is a phase of big changes in the first billion years of the Universe history. After the apparition of the first stars and the emission of energetic radiation by thoses ones, the hydrogen was reionized. This transition has an impact on the galaxies formations. I was part of the development team of EMMA, a numerical simulation code who aimed to study the processes happening during the EoR. I developed and implement a stellar formation and evolution model. These works contributed to the realisation of one of the biggest simulation dedicated to the study of the EoR yet. I contribute to the development of a tool dedicated to the exploration of this kind of simulations. I study how the radiation escaped the galaxies as a function of the parameters of the stellar model, and showed that supernovae could increase the ratio of escaping photon. I also studied the ionization fronts propagation and showed that the speed of light could be reduced by a factor 3 (and then divide the computational cost of the radiative transfer by 3), while keeping corrects results .

Ce travail s’intéresse à la phase de montée d’un aéronef civil. Les trajectoires minimisant le temps de montée ainsi que que celles minimisant la consommation de carburant sont étudiées au travers du contrôle optimal géométrique. La dynamique associée à la phase de montée possède un phénomène dit de perturbation singulière. Ce phénomène, présent dans les systèmes multi-échelle, rend difficile la résolution numérique du problème de contrôle associé. La réduction desystème hamiltonien, permettant de s’affranchir de la difficulté numérique introduite par la perturbation singulière, est étudiée d’un point de vue théorique puis numérique. Dans un second temps, le système réduit est étudié géométriquement. L’utilisation des outils du contrôle géométrique combinée à celui des synthèses à temps court permet de déterminer des familles de trajectoires localement temps-optimales pour des temps courts. Cette étude est complétée par une étude des trajectoires temps-optimales en présence de contraintes d’état. D’un point de vue plus numérique, les méthodes directes et indirectes sont utilisées pour résoudre les différents problèmes. Une synthèse locale est alors réalisée en partant des familles de trajectoires déterminées pour des temps courts. Une étude des trajectoires minimisant la consommation de carburant est également réalisée. / This work concerns the climbing phase of a civil aircraft. The trajectories which minimize the climbing time and the one which minimize the fuel consumption are studied throughout geometric optimal control. The climbing phase dynamics presents a characteristics called singular perturbation. This phenomena exists in multi-scale dynamics which makes the numerical study of the associated control problem difficult. Theoretically and numerically we study the reduction of hamiltonian system. This concept allows to remove the numerical complexity induced by the singular perturbation. Secondly, the reduced system is studied geometrically. Families of timeoptimal trajectories in small time are determined thanks to geometric control tools and small time synthesis. A study of time-optimal trajectories with active state constraints completes this work. From a more numerical point of view, direct and indirect methods are used to solve the climbing problems. A local synthesis for time-optimal trajectory is established starting from the families of trajectory determined in small time. A study of minimum fuel consumption trajectories is also realized.

La connaissance précise des champs électromagnétiques rayonnés par les dispositifs hyperfréquences nécessite des outils instrumentaux permettant la mesure directe ou indirecte de ces champs. La technique du scan champ proche fait partie de ces outils. Ce manuscrit décrit les travaux de caractérisation des sondes électromagnétiques pour une plate-forme de scan champ proche développée au laboratoire LAPLACE. L’accent a été mis sur la simulation électromagnétique des dispositifs de test utilisés dans le travail de déconvolution du facteur d’antenne des sondes de champs magnétiques ou électriques. Ces dispositifs comprennent aussi bien des structures planaires telles que des interconnexions en ligne microruban que des composants en guide d’ondes tels que des guides ouverts de formes rectangulaire ou circulaire. Pour analyser ces structures des logiciels commerciaux basés sur la méthode des différences finis ont été utilisés. Dans le cas des structures rayonnantes un programme basé sur la méthode de l’opérateur transverse a été développé, permettant la détermination de l’admittance de rayonnement et les champs rayonnés en zones proche et lointaine. Les résultats obtenus ont été validés par des simulations avec des outils commerciaux, et par des mesures réalisés au laboratoire. / The accurate knowledge of electromagnetic fields radiated by microwave devices requires instrumental tools for direct or indirect measurement of these fields. Near-field scan technique is one of those tools. This manuscript describes the work done to characterize electromagnetic probes using a near field scan platform developed in the laboratory LAPLACE. We focused our work on the electromagnetic simulation of test devices that are used in the deconvolution of antenna factor of magnetic or electric probes. These devices include both planar structures such as microstrip line and also waveguide components such as rectangular or circular open-ended waveguides. To analyze these structures, commercial software based on finite element method was used. In case of radiating structures, a program based on transverse operator method was developed. It allows the determination of the admittance of radiation and the radiated electromagnetic fields in near-field and far-field regions. The results were validated by simulations with commercial tools, and measurements made in the laboratory.

L'objectif principal de cette thèse est consacré à l'étude de la dynamique et la rhéologie d'une suspension de particules denses qui se comportent comme des fluides complexes. La premier partie de cette thèse est consacrée à l'étude de la déformation, le comportement dynamique et la rhéologie d'une suspension de vésicule (un modèle simple pour les globules rouges) sous l'action d'un écoulement externe appliqué (cisaillement simple et Poiseuille confiné) dans la limite de faible nombre de Reynolds. L'étude basée sur des simulations numériques en utilisant la méthode des intégrales de frontière. Cette étude est inspirée par le comportement des globules rouges dans le système microvasculaire. Notre étude est ensuite consacrée aux effets du confinement et du nombre capillaire sur la forme, le comportement dynamique et la viscosité effective d'une suspension de vésicules. Nous avons montré que pour des membranes rigides (nombre capillaire petit), on peut observer en plus de la forme parachute et pantoufle, les formes suivantes : (i) forme d'oscillation centrée, (ii) forme d'oscillation décentrée et (iii) la forme cacahuète. Egalement, nous avons examiné l'influence du contraste de viscosité sur la dynamique et la rhéologie d'une vésicule. Nous avons montré qu'il existe une phase de "coexistence" entre la forme pantoufle et la forme parachute. Dans la deuxième partie de cette thèse, nous avons proposé un modèle analytique et une étude numérique pour étudier les propriétés dynamiques et rhéologiques d'une suspension de particules rigides sous écoulement de Poiseuille confiné. Le débit, la dissipation et la viscosité apparente sont étudiés en fonction de la structure des plaques dans le canal. Egalement, l'étude numérique d'une suspension de particules sphériques (formes des chaînes de particules) est en accord qualitatif avec le modèle analytique qui considère les longues plaques. Cette étude numérique est basée sur une méthode de la dynamique des particules du fluide, où les particules sont représentées par un champ scalaire ayant une viscosité élevée à l'intérieur. / The principal objective of this thesis is devoted to study the dynamics and rheology of a suspension of dense particles which behave as complex fluids. The first part of this thesis is dedicated to the study of the deformation, the dynamic behavior and rheology of a single neutrally buoyant suspended vesicle (a closed phospholipid membrane), as a response to external applied flows (simple shear and confined Poiseuille flows), is studied in the limit of small Reynolds. The study based on numerical simulations using the boundary integral method. This study is inspired by the behavior of red blood cells in the microvasculature. Our study is then dedicated on the effects of confinement and capillary number on the shape, the dynamic behavior and the effective viscosity of a vesicles suspension. We have shown that for rigid membranes (small capillary number) may be exist other shapes parachute and slipper, these shapes are : (i) Centered Snaking, (ii) Off-Centered Snaking and (iii) Peanut-shaped. Equally, we examined the influence of viscosity contrast on the dynamics and rheology of a vesicle. We have shown that there is a phase of "coexistence" between the slipper and parachute shape. The coexistence of shapes seems to be supported by experiments, though a systematic experimental study is lacking. In the second part of this thesis, we proposed an analytical model and a numerical study for study the dynamics and rheology properties of a rigid particles suspended in a confined Poiseuille flow. The flow rate, the dissipation and the apparent viscosity are studied as a function of the underlying structure. Equally, the numerical study is in good qualitative agreement with the analytical theory considering long plates. The numerical study is based on a fluid dynamic particle method where the particles are represented by a scalar field having high viscosity inside.

Les systèmes multi-planétaires détectés par la mission Kepler présentent un excès de paires de planètes proches de résonances de moyen mouvement d'ordre un (e.g. 2:1, 3:2) mais avec un rapport de périodes légèrement supérieur à la valeur résonante (e.g. Pext/Pint = 2.02 au lieu de 2). Différents mécanismes ont été proposés pour expliquer cette observation. Dans cette thèse, nous étudions la possibilité que ces systèmes étaient initialement en résonance et ont atteint leur configuration actuelle à cause de la dissipation par effet de marée dans les planètes. De manière plus générale, nous établissons des critères analytiques permettant de prévoir l'évolution à long terme de systèmes en résonance d'ordre quelconque en présence de dissipation de marée. Nous montrons que cette évolution dépend de l'importance relative de la dissipation dans chacune des deux planètes. L'état final d'un système n'est pas le même suivant que la planète subissant la plus forte dissipation est la planète interne ou la planète externe. Ainsi, à partir de la configuration actuellement observée et plus spécifiquement du rapport de périodes de deux planètes proches d'une résonance, on peut déduire quelle planète a subi la dissipation la plus importante. Cela donne des contraintes importantes sur la nature des planètes considérées, car il est communément admis que les planètes rocheuses dissipent plus fortement par effet de marée que les planètes gazeuses. / Multi-planetary systems discovered by the Kepler mission present an excess of planet pairs close to first-order mean-motion resonances (e.g. 2:1, 3:2) but with a period ratio slightly higher than the resonant value (e.g. Pout/Pin = 2.02 instead of 2). Several mechanisms have been proposed to explain this observation. In this thesis, we investigate the possibility that these systems were initially in resonance and reached their current configuration due to tidal dissipation in the planets. More generally, we establish analytical criteria for predicting the long term evolution of systems locked in resonances of any order undergoing tidal dissipation. We show that this evolution depends on the relative importance of the dissipation in each of both planets. The final state of a system depends on whether the planet undergoing the stronger dissipation is the inner or the outer planet. Therefore, from the currently observed configuration and more specifically the period ratio of two planets close to a resonance, one can deduce which of both planet underwent the strongest dissipation. This provides important constraints on the nature of the considered planets, since it is commonly accepted that rocky planets dissipate more strongly by tidal effect than gaseous planets.

L’objectif de ma thèse est d'étudier la propagation des ondes acoustiques dans les milieux poreux. La théorie de l'homogénéisation (Boutin et Auriault, 1990; Malinouskaya, 2007; Li, 2010) est utilisée avec les modèles de réseaux tels que LBM, LSM, LSM2S, LBM-LSM, LBM-LSM2S. Ces modèles nous permettent de déterminer les propriétés macroscopiques, les vitesses acoustiques et les effets d'atténuation dans les échantillons Fontainebleau avec deux composants (pores et quartz) et les échantillons STATOIL avec trois composantes (pores, quartz et d'argile). La modélisation numérique vise à résoudre 3 problèmes. Le premier problème est la caractérisation des échantillons; par la détermination de la porosité et des fonctions de corrélation avec les composantes de Fourier correspondantes (Adler, 1992; Nguyen, 2013). Le second porte sur la propagation des ondes dans les échantillons secs; les vitesses sont dérivées du tenseur de rigidité efficace C(eff) qui peut être calculé par LSM (Pazdniakou, 2012) ou LSM2S. Le troisième problème correspond aux échantillons saturés par un fluide incompressible ou compressible; les vitesses sont obtenues par résolution de l'équation de Christoffel après les déterminations de C(eff), fr la perméabilité dynamique K et ses réactions à la pression de fluide α et β. Pour les échantillons de Fontainebleau, les calculs sont effectués avec des modèles pré-existents tels que LSM, LBM, LBM-LSM. Ces modèles de bases sont étendus pour milieux avec plusieurs composants solides; ils sont validés via des comparaisons avec d’autres méthodes (Nemat-Nasser et Iwakuma, 1982; Torquato, 1998, 2000; Cohen, 2004). Les vitesses, le module d’élasticité et cisaillement efficace de tous les échantillons secs ainsi que les vitesses et l'atténuation des ondes dans les échantillons saturés sont déterminés. Les séries de résultats obtenues concordent bien avec des corrélations empiriques et théoriques, tels que le modèle d’IOS d’ Arns (1996), les modèles empiriques de Nur et al. (1991), Krief (1990) et avec le modèle de Gassmann. Les résultats numériques sont un peu plus grands que les données expérimentales d’ Han (1986) et de Gomez et al. (2010); les raisons de cette petite différence ont été provisoirement analysées, mais sa cause n'a pas été identifiée sans ambiguïté. / The purpose of this Ph.D. thesis is to study acoustic waves in porous media. The homogenization theory (Boutin and Auriault, 1990; Malinouskaya, 2007; Li, 2010) is used together with the lattice models such as LBM, LSM, LSM2S, LBM-LSM, LBM-LSM2S in order to determine the macroscopic properties, the acoustic velocities, the attenuation effects in Fontainebleau samples with two components (pore and quartz) and in STATOIL samples with three components (pore, quartz and clay). Three problems are studied numerically in this work. The first problem is devoted to characterizations of samples; this is done with the determination of the porosity and of the correlation functions with the corresponding Fourier components (Adler, 1992; Nguyen, 2013). The second one addresses wave propagation in dry samples; the velocities are derived from the effective stiffness tensor C(eff) which can be calculated by LSM (Pazdniakou, 2012) or LSM2S. The third one corresponds to samples saturated by incompressible or compressible fluids; the velocities can be obtained from the Christoffel equation after determining C(eff) , the dynamic permeability K and the reactions to fluid pressure α and β. For Fontainebleau samples, the calculations are performed with basic existing models such as LSM, LBM, LBM-LSM. These basic models are extended to solids with multiple components; they are validated by comparisons with others (Nemat-Nasser and Iwakuma, 1982; Torquato, 1998, 2000; Cohen, 2004). The velocities, the effective bulk and shear modulus of all the dry samples as well as the velocities and the attenuation effected in saturated samples are determined. These results are in good agreement compared with existing models and results, such as the IOS model of Arns (1996), the empirical models of Nur et al. (1991), Krief (1990) and with Gassmann’ s model. The numerical results are slightly larger than the experimental data of Han (1986) and Gomez et al. (2010); the origin of this small discrepancy has been tentatively analyzed, but its cause has not been unambiguously identified.

L'objet de ce travail de recherche est d'estimer le rôle des structures tourbillonnaires pour améliorer l'efficacité des refroidisseurs microfluidiques. Cette problématique a été étudiée en deux parties. L'objectif de la première partie est de quantifier les transferts thermiques engendrés par des tourbillons lors de leur interaction avec une paroi chaude. Ceci a été réalisé numériquement en deux étapes. La première étape a pour vocation d'améliorer la compréhension des processus de transferts thermiques d'un tourbillon transporté dans un écoulement laminaire. Deux processus majeurs ont été mis en avant : l'advection des particules froides vers la paroi chaude et le mélange par advection.Suite à ces observations, les transferts thermiques induis par une structure tourbillonnaire ont été optimisés en fonction de ses caractéristiques initiales. Cette optimisation a permis d'augmenter l'efficacité d'un système de refroidissement microfluidique de plus de 50% dans le cas d'une allée tourbillonnaire.La seconde partie couple une étude expérimentale à un modèle numérique. Cette étude se focalise sur le formation des structures tourbillonnaires à l'échelle micrométrique. Les structures tourbillonnaires sont générées par le couplage d'un micro jet synthétique à un écoulement transverse. Dans cette configuration, plusieurs topologies de structures tourbillonnaires ont été observées en fonction du nombre de Reynolds du jet et de l'écoulement. Ces différentes topologies ont été cartographiées en fonction des nombres de Reynolds du système. Finalement, pour une certaine topologie de tourbillon il a été possible de lier leurs caractéristiques aux paramètres de contrôle du système. Grâce à cette relation, il est possible de contrôler le type de structures formées par le système de refroidissement / The purpose of this research paper is to determine the role of vortical structures in the improvement of microfluidic cooling systems’ efficiency.This study is presented in two parts. The aim of the first part is to measure heat transfers produced by vortices while interacting with a hot wall. This part was carried out numerically following two steps. The aim of the first step was to improve the understanding of the processes involved during heat transfers in a vortex convected in a laminar flow. As a result, two major processes can be highlighted : the advection of cold particles to the hot wall and an advective mixing. Following these observations, heat transfers produced by a vortical structure were optimized in accordance with its initial characteristics. Consequently, the optimization of heat transfers increased microfluidic cooling systems’ efficiency by more than 50% in the case of vortex streets.The second part combines an experimental study with a numerical model. This study focuses on the creation of vortical structures on a microscopic scale. For this purpose, vortical structures are produced by combination of a synthetic micro jet with a crossflow. The outcome is that several topologies of vortical structures can be observed depending on the Reynolds number of the jet and the crossflow. These topologies can be mapped according to the Reynolds numbers of the system. As a result, a connection can be made between characteristics of some topologies of vortices and control parameters of the system. As a conclusion, the type of structures produced by the cooling system can be controlled

Dans ce projet de recherche, on s'intéresse au développement et à l'évaluation de nouvelles méthodes numériques pour les écoulements peu profonds. De nouvelles techniques de discrétisation spatiales et temporelles des équations sont proposées. La première partie de la thèse est dédiée au développement d'une méthode des volumes finis explicite d'ordre élevé et d'une famille de schémas semi-implicites qui sont efficaces pour la modélisation des processus lents et rapides dans les écoulements océaniques et atmosphériques. La deuxième partie du projet de recherche concerne la construction d'un schéma numérique efficace sans solveur de Riemann pour les écoulements peu profonds avec une topographie variable sur un maillage non structuré. Dans cette partie de la thèse, une nouvelle approche est proposée pour l'analyse de stabilité des schémas numériques non structurés pour les équations en eaux peu profondes. Dans la troisième partie de la thèse, deux schémas de volumes finis sont développés pour les lois de conservation sur des surfaces courbes qui ont un large potentiel d'être appliqués aux écoulements peu profonds sur la sphère. Dans ces cas, les schémas numériques sont développés en adoptant la démarche suivie par Stanley Osher. Cette démarche consiste à utiliser des systèmes hyperboliques simples qui génèrent des phénomènes d'ondes complexes et des solutions qui ont différentes structures. Ces solutions sont très efficaces pour tester les méthodes numériques. Dans notre cas, nous avons utilisé les équations de Burgers qui ont joué un rôle très important dans le développement des schémas numériques à capture de chocs en mécanique des fluides. / This research project focuses on the development and evaluation of numerical methods for shallow flows by proposing new spatial and temporal discretization techniques. First, a new high-order explicit finite volume method and a class of semi-implicit schemes are introduced which are effective for modelling fast and slow waves in oceanic and atmospheric flows. In the second part of the research project, a central-upwind scheme is proposed for shallow water flows on variable topography using unstructured grids. In this part of the project, a new approach is proposed for the stability analysis of unstructured numerical schemes for shallow water equations. In the third part of the thesis, two finite volume methods are developed for the conservation laws on curved geometries which are potentially applicable to shallow flows on a sphere. For such cases, numerical schemes are developed by using the approach followed by Stanley Osher. This approach employs simple hyperbolic systems which generate complex wave phenomena, and solutions that are effective for assessing numerical methods. In our case, Burgers’ equations are used since they have played an important role in the development of shock-capturing schemes in fluid mechanics.

Les ingénieurs sont constamment mis au défi de concevoir des avions plus légers et moins polluants. En même temps, ils se doivent de respecter un certain nombre de critères de dimensionnement établis pour assurer l'intégrité des structures aéronautiques. Pour répondre à ce double défi, les chercheurs du milieu industriel et académique travaillent à l'unisson pour constamment repousser les limites du savoir scientifique. Ces collaborations aboutissent notamment au développement de nouveaux outils mathématiques, numériques, expérimentaux et de production.On trouve de nombreuses pièces métalliques parmi les composants aéronautiques les plus critiques, comme les trains d’atterrissages, les arbres moteurs ou les chapes. Ces structures subissent des chargements extrêmes de par leur environnement thermique ou leur intense sollicitation mécanique ponctuelle ou cyclique. Pour ces structures élastoplastiques, la ruine peut généralement être causée par l'apparition et la propagation de fissures, ou l'émergence d'instabilités locales, comme la localisation de la déformation plastique, ou d'instabilités globales, comme le flambement.Dans les travaux présentés dans ce manuscrit, nous nous sommes intéressés à la détection d’instabilités locales et globales dans des structures élastoplatiques en grandes déformations. Alors qu'elles sont classiquement considérées indépendamment, nous avons fait le choix d'étudier les deux types d'instabilités simultanément pour mieux comprendre la compétition entre ces deux phénomènes. Pour ce faire, plusieurs développements analytiques sont présentés et fondés sur le critère de stabilité de Hill (1958) et le critère de localisation de la déformation de Rice (1976), ainsi que la proposition d'une nouvelle méthode d'"analyses de stabilité affaiblie" permettant d'investiguer la sensibilité aux types de conditions limites imposées. Tous ces critères ont été implémentés dans le code de calcul par éléments finis Zset, ce qui nous a permis, entre autres, d'analyser l’apparition d'instabilités élastoplastiques dans différentes éprouvettes de traction/torsion, des tubes en torsion et une chape en traction. / Aeronautical engineers are constantly challenged to provide lighter structures in order to reduce fuel consumption, and thus the environmental impact and flight costs. At the same time, the design of aeronautical structures is subjected to strict regulation aimed at ensuring the integrity of the aircraft and the safety of the passengers. To tackle this challenge, the limits of structural and material mechanics are consistently explored which in turn leads to the development of new, mathematical, numerical, experimental and manufacturing tools.There are numerous metallic parts in the most critical aeronautical structures, like landing gears, engine shafts, or mechanical lugs. These parts are subjected to extreme loading conditions due to the thermal environment or to the intense mechanical ultimate or cyclic loading. The failure of these elastoplastic structures is generally caused by the initiation and propagation of cracks or by the emergence of local instabilities, such as plastic strain localization, or global instabilities, such as buckling.In the present work, we focus on the detection of local and global elastoplastic instabilities in a finite deformation framework. While they are generally studied separately, it was chosen to study both phenomena together in order to analyze and better understand the competition between localization and buckling in elastoplastic structures. For this purpose, multiple analytical developments are presented founded on Hill's global stability criterion (1958) and Rice's strain localization criterion (1976). The new "weakened stability analysis" has been introduced in order to analyze the sensitivity to the type of prescribed boundary conditions. All these criteria have been implemented in the finite element software Zset, which allowed us to analyze the emergence of elastoplastic instabilities in various experimental samples, tubes loaded in torsion, and a lug loaded in tension.

Les graphes, ou réseaux, sont des structures mathématiques représentant des relations entre des éléments. Ces systèmes peuvent être analysés dans le but d’extraire des informations sur la structure globale ou sur des composants individuels. L'analyse de graphe conduit souvent à des problèmes hautement complexes à résoudre. À grande échelle, le coût de calcul de la solution exacte est prohibitif. Heureusement, il est possible d’utiliser des méthodes d’approximations itératives pour parvenir à des estimations précises. Lesméthodes historiques adaptées à un petit nombre de variables ne conviennent pas aux matrices creuses de grande taille provenant des graphes. Par conséquent, la conception de solveurs fiables, évolutifs, et efficaces demeure un problème essentiel. L’émergence d'architectures parallèles telles que le GPU ouvre également de nouvelles perspectives avec des progrès concernant à la fois la puissance de calcul et l'efficacité énergétique. Nos travaux ciblent la résolution de problèmes de valeurs propres de grande taille provenant des méthodes d’analyse de graphe dans le but d'utiliser efficacement les architectures parallèles. Nous présentons le domaine de l'analyse spectrale de grands réseaux puis proposons de nouveaux algorithmes et implémentations parallèles. Les résultats expérimentaux indiquent des améliorations conséquentes dans des applications réelles comme la détection de communautés et les indicateurs de popularité / Graphs, or networks, are mathematical structures to represent relations between elements. These systems can be analyzed to extract information upon the comprehensive structure or the nature of individual components. The analysis of networks often results in problems of high complexity. At large scale, the exact solution is prohibitively expensive to compute. Fortunately, this is an area where iterative approximation methods can be employed to find accurate estimations. Historical methods suitable for a small number of variables could not scale to large and sparse matrices arising in graph applications. Therefore, the design of scalable and efficient solvers remains an essential problem. Simultaneously, the emergence of parallel architecture such as GPU revealed remarkable ameliorations regarding performances and power efficiency. In this dissertation, we focus on solving large eigenvalue problems a rising in network analytics with the goal of efficiently utilizing parallel architectures. We revisit the spectral graph analysis theory and propose novel parallel algorithms and implementations. Experimental results indicate improvements on real and large applications in the context of ranking and clustering problems