L'objectif de cette thèse est de développer des méthodes numériques dans un code de simulation numérique direct contenant une méthode de suivi Lagrangien de phase disperse. Ces méthodes numériques ont permis d'optimiser le suivi Lagrangien parallèle des gouttes afin de pouvoir simuler de très grands domaines. Un algorithme de gestion des collisions entre les gouttes a aussi été développé et mis en place ainsi qu'une procédure permettant d'intégrer des corps immergés de formes quelconques. Ces méthodes numériques ont permis d'effectuer des simulations de référence et d'obtenir des statistiques détaillées de diverses configurations telles que le croisement de jets de gouttes qui se collisionnent ou bien la propagation d'un front de flamme dans un spray évaporant. La comparaison de ces caractéristiques avec les mêmes simulations où la prise en charge du spray était effectuée par une méthode Eulérienne multi-fluide ont été faites. Les très bonnes concordances qualitatives et quantitatives entre les prises en charge Lagrangiennes et Eulériennes ont permis de montrer les capacités du modèle Eulérien à rendre compte de phénomènes physiques complexes. Ce modèle Eulérien permettra de bien meilleurs scaling lors de son portage sur des architectures massivement parallèles.