Les tourbillons sont des structures importantes pour une large gamme d'écoulements de fluides, notamment les sillages, l'interaction fluide-structure, les décollements de couche limite et la turbulence. Cependant, les méthodes numériques classiques n'arrivent généralement pas à donner une représentation précise des tourbillons. Ceci est principalement lié à la dissipation numérique des schémas qui, si elle n'est pas spécifiquement calibrée pour le calcul des écoulements tourbillonnaires, conduit à une diffusion artificielle très rapide des tourbillons dans les calculs. Parmi d'autres approches, la méthode "Vorticity Confinement" (VC) de J. Steinhoff permet de compenser la dissipation des schémas au sein des tourbillons en introduisant une anti-dissipation non-linéaire, mais elle n’est précise qu’au premier ordre. D’autre part, des progrès significatifs ont récemment été accomplis dans le développement de méthodes numériques d’ordre élevé. Celles-ci permettent de réduire ce problème de dissipation excessive, mais la diffusion des tourbillons reste importante pour de nombreuses applications. La présente étude vise à développer des extensions d’ordre élevé de la méthode VC pour réduire cette dissipation excessive des tourbillons, tout en préservant la précision d'ordre élevé des schémas. Tout d'abord, les schémas de confinement sont analysés dans le cas de l'équation de transport linéaire, à partir de discrétisations couplées et découplées en espace et en temps. Une analyse spectrale de ces schémas est effectuée analytiquement et numériquement en raison de leur caractère non linéaire. Elle montre des propriétés dispersives et dissipatives améliorées par rapport aux schémas linéaires de base à tous les ordres de précision. Dans un second temps, des schémas VC précis au troisième et cinquième ordre sont développés pour les équations de Navier-Stokes compressibles. Les termes correctifs restent conservatifs, invariants par rotation et indépendants du schéma de base, comme la formulation originale VC2. Les tests numériques valident l'ordre de précision et la capacité des extensions VC d’ordre élevé à réduire la dissipation dans les tourbillons. Enfin, les schémas avec VC sont appliqués au calcul des écoulements turbulents, dans une approche de simulation de grandes échelles implicite (ILES). Les schémas numériques avec VC présentent une résolvabilité améliorée par rapport à leur version linéaire de base, et montrent leur capacité à décrire de façon cohérente ces écoulements tourbillonnaires complexes. / Vortices are flow structures of primary interest in a wide range of fluid dynamics applications including wakes, fluid-structure interaction, flow separation and turbulence. Albeit their importance, standard Computational Fluid Dynamics (CFD) methods very often fail to provide an accurate representation of vortices. This is primarily related to the schemes’ numerical dissipation which, if inadequately tuned for the calculation of vortical flows, results in the artificial spreading and diffusion of vortices in numerical simulations. Among other approaches, the Vorticity Confinement (VC) method of J. Steinhoff allows balancing the baseline dissipation within vortices by introducing non-linear anti-dissipation in the discretization of the flow equations, but remains at most first-order accurate. At the same time, remarkable progress has recently been made on the development of high-order numerical methods. These allow reducing the problem of excess dissipation, but the diffusion of vortices remains important for many applications. The present study aims at developing high-order extensions of the VC method to reduce the excess dissipation of vortices, while preserving the accuracy of high-order methods. First, the schemes are analyzed in the case of the linear transport equation, based on time-space coupled and uncoupled formulations. A spectral analysis of nonlinear schemes with VC is performed analytically and numerically, due to their nonlinear character. These schemes exhibit improved dispersive and dissipative properties compared to their linear counterparts at all orders of accuracy. In a second step, third- and fifth-order accurate VC schemes are developed for the compressible Navier-Stokes equations. These remain conservative, rotationally invariant and independent of the baseline scheme, as the original VC2 formulation. Numerical tests validate the increased order of accuracy and the capability of high-order VC extensions to balance dissipation within vortices. Finally, schemes with VC are applied to the calculation of turbulent flows, in an implicit Large Eddy Simulation (ILES) approach. In these applications, numerical schemes with VC exhibit improved resolvability compared to their baseline linear version, while they are capable of producing consistent results even in complex vortical flows.

Ce travail de thèse porte principalement sur l'utilisation du calcul moulien et de la technique d'arborification introduits par Jean Ecalle dans les années 70 et leurs applications à l'étude des systèmes dynamiques discrets ou continus.L'une des contributions est une étude systématique des conditions sous lesquelles l'arborification permet de restaurer la convergence de séries formelles via l'introduction d'une notion d'invariance d'un moule sous arborication. Ces résultats permettent de donner une preuve détaillée du théorème de Brjuno de linéarisation analytique des champs de vecteurs telle qu'elle est proposée par Jean Ecalle dans son article "Singularités non abordables par la géométrie". Ces résultats ont été obtenus en collaboration avec Dominique Manchon (Université de Clermont Ferrand) et Jacky Cresson.La puissance du calcul moulien est ensuite illustrée par la résolution presque complète de la conjecture de Jarque-Villadelprat sur les centres isochrones Hamiltoniens. Cette conjecture stipule qu'il n'existe pas de champs de vecteurs polynomiaux du plan de degré pair qui soit hamiltonien. L'examen de la structure algébrique de la correction, introduite dans les années 90 par G. Gallavotti et généralisée ensuite par Jean Ecalle et Bruno Vallet, et son calcul explicite via le calcul moulien, nous ont permis d'obtenir des conditions explicites d'obstructions à l'isochronisme. L'aspect algébrique et combinatoire de ces objets et méthodes conduisent naturellement à une classication des conditions de centre via une notion de complexité. L'arborication quand à elle permet l'unification de nombreuses approches et une simplication de divers travaux, notamment ceux de J.C.Butcher autour de la structure algébrique des méthodes de Runge-Kutta qui a induit ce que les numériciens appellent des B-séries. En étudiant la structure algébrique de l'opérateur de substitution associé à un difféomorphisme, en particulier celui relié à une méthode de Runge-Kutta et celui associé à la solution de l'équation diérentielle sous-jacente, on présente le codage de Butcher comme une traduction particulière de l'arborification directe de l'opérateur de substitution. Notons que ce phénomène est large et permet d'inclure les travaux plus récents sur l'approche par trajectoires rugueuses des solutions d'équations différentielles stochastiques.Une seconde partie de la thèse concerne la recherche des groupes de symétries de Lie des tissus du plan en suivant une approche d'Alain Hénaut (Université de Bordeaux). Ce travail nous a permis de préciser la relation entre la dimension de ces groupes de symétries et le caractère linéarisable ou hexagonale des tissus du plan. Dans le cas des arrangements de droites, on obtient ainsi une relation profonde entre le module de dérivations de Saito associé à l'arrangement et le groupe de symétrie du tissu associé. / This thesis work mainly focuses on the use of the mould calculus and the technic of arborification which had been introduced both by J.Ecalle in the seventies and theirs applications to the study of continuous or discrete systems.One of the contributions is the systematic study of conditions under which the arborification allows to reestablish the convergence of formal series via introduction of a notion of invariance of mould under arborification. These results allow to give a detailed proof of Brjuno Theorem of analytic linearizability of vector fields as it is proposed by J.Ecalle in his article "Singularité non abordable par la géométrie". These results were obtained jointly with Dominique Manchon (University of Clermont Ferrand) and Jacky Cresson.The power of the mould calculus is then illustrated by an almost complete resolution of the Jarque-Villadelprat's conjecture about Hamiltonian Isochronous centers. This conjecture states that there is not existing polynomial vector fields in the plane of odd degree which are Hamiltonian. The study of the algebraic structure of the correction, introduced in the nineties by G.Gallavotti and then generalized by J.Ecalle and B.Vallet and its explicit computation via mould calculus, enables us to obtain explicit conditions of obstruction to isochronicity. The algebraic and combinatoric aspect of these objects and methods brings naturally to the classification of center conditions through a notion of complexity. The arborification allows to the unification of different approaches and a simplicification of different works, especially those of J.C.Butcher about algebraic structures of Runge-kutta methods, who had introduced that is called B-series by numerical mathematicians. Studying the algebraic structure of the substitution operator associated to a diffeomorphism, especially the one related to a Runge-Kutta method and the one which is associated to the solution of the underlying differential equations, we present the Butcher's encoding as a special translation of a direct arborification of the substitution automorphism. We can conclude that this phenomenon is wide and allows to include more recent studies on the approach by rough path of stochastic differential equations.A second part of this thesis involves the research of Lie group of symmetries of planar webs following Hénaut's approach (University of Bordeaux).This work allows to precise the relation between the dimension of the groups of symmetries and the linearizability or hexagonal character of planar webs. In the the case of line arrangement, we obtain a depthful relation between the modulus of derivations of Saito associated to the line arrangement and the group of symmetries of the associated web.

Cette thèse est dediée à l'alanyse numérique de méthodes numériques dans le domaine du contrôle optimal, et est composée de deux parties. La première partie est consacrée à des nouveaux résultats concernant des méthodes numériques pour le contrôle optimal de systèmes hybrides, qui peuvent être contrôlés simultanément par des fonctions mesurables et des sauts discontinus dans la variable d'état. La deuxième partie est dédiée è l'étude d'une application spécifique surl'optimisation de trajectoires pour des lanceurs spatiaux avec contraintes en probabilité. Ici, on utilise des méthodes d'optimisation nonlineaires couplées avec des techniques de statistique non parametrique. Le problème traité dans cette partie appartient à la famille des problèmes d'optimisation stochastique et il comporte la minimisation d'une fonction de coût en présence d'une contrainte qui doit être satisfaite dans les limites d'un seuil de probabilité souhaité. / This thesis is devoted to the analysis of numerical methods in the field of optimal control, and it is composed of two parts. The first part is dedicated to new results on the subject of numerical methods for the optimal control of hybrid systems, controlled by measurable functions and discontinuous jumps in the state variable simultaneously. The second part focuses on a particular application of trajectory optimization problems for space launchers. Here we use some nonlinear optimization methods combined with non-parametric statistics techniques. This kind of problems belongs to the family of stochastic optimization problems and it features the minimization of a cost function in the presence of a constraint which needs to be satisfied within a desired probability threshold.

Les architectures de calcul haute performance les plus récentes intègrent de plus en plus de nœuds de calcul qui contiennent eux-mêmes plus de cœurs. Les bus mémoires et les réseaux de communication sont soumis à un niveau d'utilisation critique. La programmation parallèle sur ces nouvelles machines nécessite de porter une attention particulière à ces problématiques pour l'écriture de nouveaux algorithmes. Nous analysons dans cette thèse un code de simulation de turbulences de plasma et proposons une refonte de la parallélisation de l'opérateur de gyromoyenne plus adapté en termes de distribution de données et bénéficiant d'un schéma de recouvrement calcul -- communication efficace. Les optimisations permettent un gain vis-à-vis des coûts de communication et de l’empreinte mémoire. Nous étudions également les possibilités d'évolution de ce code à travers la conception d'un prototype utilisant un modèle programmation par tâche et un schéma de communication asynchrone adapté. Cela permet d'atteindre un meilleur équilibrage de charge afin de maximiser le temps de calcul et de minimiser les communications entre processus. Un maillage réduit adaptatif en espace est proposé, diminuant le nombre de points sans pour autant perdre en précision, mais ajoutant de fait une couche supplémentaire de complexité. Ce prototype explore également une distribution de données différente ainsi qu'un maillage en géométrie complexe adapté aux nouvelles configurations des tokamaks. Les performances de différentes optimisations sont étudiées et comparées avec le code préexistant et un cas dimensionnant sur un grand nombre de cœurs est présenté. / Recent high performance computing architectures come with more and more cores on a greater number of computational nodes. Memory buses and communication networks are facing critical levels of use. Programming parallel codes for those architectures requires to put the emphasize on those matters while writing tailored algorithms. In this thesis, a plasma turbulence simulation code is analyzed and its parallelization is overhauled. The gyroaverage operator benefits from a new algorithm that is better suited with regard to its data distribution and that uses a computation -- communication overlapping scheme. Those optimizations lead to an improvement by reducing both execution times and memory footprint. We also study new designs for the code by developing a prototype based on task programming model and an asynchronous communication scheme. It allows us to reach a better load balancing and thus to achieve better execution times by minimizing communication overheads. A new reduced mesh is introduced, shrinking the overall mesh size while keeping the same numerical accuracy but at the expense of more complex operators. This prototype also uses a new data distribution and twists the mesh to adapt to the complex geometries of modern tokamak reactors. Performance of the different optimizations is studied and compared to that of the current code. A case scaling on a large number of cores is given.

Dans divers domaines de la physique, certains phénomènes sont modélisés par des systèmes hyperboliques non-conservatifs. En particulier, dans le domaine de la physique des plasmas, dont l'un des champs d'application majeur est la Fusion par Confinement Inertiel, le système d'Euler bi-température, modélisant les phénomènes de transport de particules chargées, en est un exemple. La difficulté de l'étude de ces systèmes réside dans la présence de termes non-conservatifs, qui empêchent la définition classique des solutions faibles. Pour parvenir à une définition de ce type de solutions, on a recours à l'emploi de systèmes cinétiques sous-jacents. Dans ce manuscrit, on s'intéresse à l'étude numérique de ces systèmes cinétiques pour la résolution du système d'Euler bi-température.Ce manuscrit se divise en deux parties. La première partie contient l'étude numérique du système d'Euler bi-température. Dans un premier chapitre, on résout numériquement les équations en dimension 1 d'espace par le biais d'un système sous-jacent issu de la physique des plasmas: le système de Vlasov-BGK-Ampère. On présente une méthode numérique préservant l'asymptotique pour ce système sous-jacent et on montre, par des simulations numériques, que le schéma limite obtenu donne des résultats consistants avec Euler bi-température. Dans un second chapitre, on résout le même modèle en dimension 2 d'espace par un système sous-jacent de type BGK discret. On démontre une inégalité d'entropie pour les solutions issues du modèle sous-jacent, ainsi qu'une inégalité discrète de dissipation d'entropie pour le schéma.Dans la deuxième partie de ce manuscrit, on s'intéresse au développement de méthodes numériques pour quelques modèles cinétiques. On considère ici le cas des écoulements raréfiés de mélanges de gaz, dans l'optique d'une application aux cas des plasmas. Premièrement, on présente un schéma cinétique adaptatif et dynamique en vitesse pour les gaz inertes. Par l'emploi de lois de conservation discrètes, la solution est approchée sur un ensemble de vitesses discrètes local et dynamique. Dans un second temps, on propose une extension de cette méthode visant à améliorer les performances de celle-ci. Puis, ces deux versions de la méthode sont comparées à la méthode classique sur grille fixe uniforme sur une série de cas tests.Enfin, dans le dernier chapitre, on propose une méthode numérique pour la résolution d'une extension de ces équations, prenant en compte la présence de réactions chimiques au sein du mélange. Le contexte considéré est celui des réactions chimiques bi-moléculaires réversibles lentes. La méthode proposée, de type implicite-explicite, est linéaire, stable et conservative. / In various domains of physics, several phenomena can be modeled via the use of nonconservative hyperbolic systems. In particular, in plasma physics, in the process of developping and understanding the phenomena leading to Inertial Confinement Fusion, the bi-temperature Euler sytem can be used to model particle transport phenomena in a plasma. The difficulty of the mathematical study of such systems dwells in the presence of so-called non-conservative products, which prevent the classical definition of weak solutions via distribution theory. To attempt to define these quantities, it is useful to supplement the hyperbolic system with an underlying kinetic model. In this work, the objective is the numerical study of such kinetic systems in order to solve the bi-temperature Euler system.This manuscript is split in two parts. The first one contains the study of the bi-temperature Euler system. In the first chapter, this system in dimension 1 is solved by the use of an underlying kinetic model sprung from plasma physics: the Vlasov-BGK-Ampère system. An asymptotic-preserving numerical method is introduced, and it is shown that the scheme obtained in the limit is consistant with a scheme for teh bi-temperature Euler system. In the following chapter, the same hyperbolic model in dimension 2 is studied, this time via a discrete-BGK type underlying model. An entropy inequality is proved for solutions coming from the kinetic model, as well as a discrete entropy dissipation inequality.In the second part of the manuscript, we are interested in the development of numerical schemes for gas mixture rarefied flows. Firstly, an adaptive kinetic scheme is introduced for inert gas mixtures. By the use of discrete conservation laws, the solution is approximated on a set of discrete velocities that depends on space, time and species. Secondly, an extension of the method is proposed in order to improve the efficiency of the first method. Finally, the two methods are compared to the classical fixed grid method on a series of test cases.In the last chapter, a numerical method is proposed for rarefied flows of reacting mixtures. The setting considered is the case of slow bimolecular reversible chemical reactions. The method introduced is an explicit-implicit treatment of the relaxation operator, which is shown to be stable, linear and conservative.

Ce mémoire est consacré à l'analyse des effets du confinement de l'air et du caractère multidimensionnel de l'écoulement lors des essais de caractérisation hydrodynamique d'un sol par la méthode du drainage interne. L'étude expérimentale et numérique de la compression de l'air au cours de l'infiltration montre que les caractéristiques du sol ne sont pas affectées par l'écoulement de l'air, dès lors que la pression capillaire est mesurée et non celle de l'eau uniquement. L'étude numérique du caractère axisymétrique de l'écoulement lors d'un essai à l'infiltromètre double anneaux montre que malgré l'importance du transfert radial, son influence sur la caractérisation d'un milieu homogène reste faible. En revanche, la caractérisation de milieux stratifiés horizontalement est très délicate.

L'étude de la propagation de la lumière dans les fibres optiques monomodes requiert la prise en compte de plusieurs phénomènes compliqués tels que la dispersion modale de polarisation et l'effet Kerr. Il s'est avéré que l'évolution de l'enveloppe lentement variable du champ électrique est bien décrite par un système couplé d'équations de Schrödinger non linéaires à coefficients aléatoires : l'équation de Manakov PMD. Cette équation fait intervenir différentes échelles dont le ratio est donné par un petit paramètre. La première partie de ce travail consiste à étudier le comportement asymptotique de la solution de l'équation de Manakov PMD lorsque ce petit paramètre tend vers zéro. En généralisant la théorie de l'Approximation-Diffusion au cadre de la dimension infinie, on a montré que la dynamique asymptotique est donnée par une équation aux dérivées partielles stochastiques dirigée par un mouvement brownien de dimension trois. Dans une seconde partie, nous proposons un schéma de différences finies de type Crank Nicolson pour cette équation pour lequel nous obtenons un ordre de convergence en probabilité d'ordre 1/2. La discrétisation du bruit doit être implicite afin d'obtenir un schéma conservatif et stable. Enfin la dernière partie est relative à la simulation numérique de la dispersion modale de polarisation et à ses effets sur la propagation et la collision de solitons de Manakov. Dans ce cadre, on propose une méthode de réduction de variance valable pour les équations aux dérivées partielles stochastiques.

Cette thèse est consacrée à l'étude de certains aspects numériques et mathématiques liés à la simulation d'écoulements incompressibles triphasiques à l'aide d'un modèle à interfaces diffuses de type Cahn-Hilliard/Navier-Stokes. La discrétisation spatiale est effectuée par éléments finis. La présence d'échelles très différentes dans le système suggère l'utilisation d'une méthode de raffinement local adaptatif. La procédure mise en place permet de tenir compte implicitement des non conformités des maillages générés, pour produire in fine des espaces d'approximation conformes. Nous montrons, en outre, qu'il est possible d'exploiter cette méthode pour construire des préconditionneurs multigrilles. Concernant la discrétisation en temps, notre étude a commencé par celle du système de Cahn-Hilliard. Pour remédier aux problèmes de convergence de la méthode de Newton utilisée pour résoudre ce système (non linéaire), nous proposons un schéma semi-implicite permettant de garantir la décroissance de l'énergie. Nous montrons l'existence et la convergence des solutions discrètes. Nous poursuivons ensuite cette étude en donnant une discrétisation en temps inconditionnellement stable du modèle complet Cahn-Hilliard/Navier-Stokes ne couplant pas fortement les deux systèmes. Nous montrons l'existence des solutions discrètes et, dans le cas où les trois fluides ont la même densité, nous montrons leur convergence. Nous étudions, pour terminer cette partie, diverses problématiques liées à l'utilisation de la méthode de projection incrémentale. Enfin, la dernière partie présente plusieurs exemples de simulations numériques, diphasiques et triphasiques, en deux et trois dimensions.

La Simulation aux Grandes Echelles (SGE) est de plus en plus utilisée dans les processus de développement et la conception des réacteurs aéronautiques industriels. L'une des raisons pour ce besoin résulte dans la capacité de la SGE à fournir des informations instantanées d'un écoulement turbulent augmentant la quantité des prédictions de la composition des gaz d'échappement. Ce manuscrit de thèse aborde deux sujets récurrents de la SGE. D'une part, les schémas numériques pour la SGE nécessitent certaines propriétés, notamment une précision élevée avec une diffusivité faible pour ne pas nuire aux modèles de turbulence. Afin de répondre à ce pré requis, une famille de schémas d'intégration temporelle d'ordre élevée est proposée, permettant de modifier la diffusion numérique du schéma. D'autre part, la SGE étant intrinsèquement instationnaire, elle est très consommatrice en temps CPU. De plus, une géométrie complexe prend beaucoup de temps de simulation même avec les super calculateurs d'aujourd'hui. Dans le cas particulier d'intérêt et souvent rencontré dans les applications industrielles, l'approche bas-Mach est constitue une alternative intéressante permettant de réduire le coût et le temps de retour d'une simulation LES. L'impact et la comparaison des formalismes compressible et incompressible sont toutefois rarement quantifiés, ce qui est proposé dans ce travail pour une configuration représentative d'un brûleur swirlé industriel mesuré au CORIA.

Le modèle multi-fluide permet de décrire par une approche Eulérienne les sprays polydispersés et apparaît donc comme une méthode indiquée pour les applications de combustion diphasique. Sa pertinence pour la simulation à l’échelle d’applications industrielles est évaluée dans ce travail, par sa mise en oeuvre dans des configurations bi-dimensionnelle et tri-dimensionnelle plus représentatives de ce type de simulations. Cette évaluation couple une étude de faisabilité en terme de coût de calcul avec une analyse de la précision obtenue, par des comparaisons avec les résultats de méthodes de références pour la description des sprays. Afin de définir une telle référence, une hiérarchisation des modèles de spray est proposée dans ce travail, soulignant les niveaux de modélisation associée aux diverses méthodes. Une première configuration d’écoulements tourbillonnaires est utilisée pour caractériser la méthode multi-fluide. L’étude de la structure mathématique du système de lois de conservation permet d’analyser la formation de singularités et de fournir les outils permettant d’évaluer leur impact sur la modélisation. Cette étude permet également de dériver un schéma numérique robuste et efficace pour des configurations bi- et tri-dimensionnelle. La description des dynamiques de gouttes conditionnées par la taille est évaluée dans ces configurations tourbillonnaires au moyen de comparaisons quantitatives, sur des champs instantanés, où le multi-fluide est confronté à une méthode Lagrangienne, ainsi qu’à des résultats expérimentaux. Afin d’évaluer le comportement de la méthode multi-fluide dans des configurations plus représentatives des problématiques industrielles, le solveur MUSES3D est développé, permettant, entre autres, une évaluation fine des méthodes de résolution des sprays. Une implémentation originale de la méthode multi-fluide, conciliant généricité et efficacité pour le calcul parallèle, est réalisée. Le couplage de ce solveur avec le code ASPHODELE, développé au CORIA, permet d’effectuer une évaluation opérationnelle des approches Euler/Lagrange et Euler/Euler pour la description des écoulements diphasiques à inclusions dispersées. Finalement, le comportement de la méthode multi-fluide dans des jets bi-dimensionnels et dans une turbulence homogène isotrope tri-dimensionnelle permet de montrer sa précision pour la description de la dynamique de sprays évaporant dans des configurations plus complexes. La résolution de la polydispersion du spray permet de décrire précisément la fraction massique de combustible en phase vapeur, un élément clé pour les applications de combustion. De plus, l’efficacité du calcul parallèle par décomposition de domaine avec la méthode multi-fluide permet d’envisager son utilisation à l’échelle d’applications industrielles. / The multi-fluid model, providing a Eulerian description of polydisperse sprays, appears as an interesting method for two-phase combustion applications. Its relevance as a numerical tool for industrial device simulations is evaluated in this work. This evaluation assesses the feasibility of multi-fluid simulations in terms of computational cost and analyzes their precision through comparisons with reference methods for spray resolution. In order to define such a reference, the link between the available methods for spray resolution is provided, highlighting their corresponding level of modeling. A first framework of 2-D vortical flows is used to assess the mathematical structure of the multi-fluid model governing system of equations. The link between the mathematical peculiarities and the physical modeling is provided, and a robust numerical scheme efficient for 2-D/3-D configurations is designed. This framework is also used to evaluate the multi-fluid description of evaporating spray sizeconditioned dynamics through quantitative, time-resolved, comparisons with a Lagrangian reference and with experimental data. In order to assess the multi-fluid efficiency in configurations more representative of industrial devices, a numerical solver is designed, providing a framework devoted to spray method evaluation. An original implementation of the multifluid method, combining genericity and efficiency in a parallel framework, is achieved. The coupling with a Eulerian/Lagrangian solver for dispersed two-phase flows, developed at CORIA, is conducted. It allows a precise evaluation of Euler/Lagrange versus Euler/Euler approaches, in terms of precision and computational cost. Finally, the behavior of the multi-fluid model is assessed in 2D-jets and 3-D Homogeneous Isotropic Turbulence. It illustrates the ability of the method to capture evaporating spray dynamics in more complex configurations. The method is shown to describe accurately the fuel vapor mass fraction, a key issue for combustion applications. Furthermore, the method is shown to be efficient in domain decomposition parallel computing framework, a key issue for simulations at the scale of industrial devices.