Le but de cette recherche consiste à la modélisation numérique et expérimentale d’écoulement sanguin dans les artères, en utilisant les nouvelles méthodes de couplage fluide structure. Plusieurs approches numériques peuvent être utilisées pour ce type de couplage, la méthode des éléments finis pour la modélisation de la structure. Pour la modélisation du domaine fluide la méthode de type éléments fnis ou méthode particulaire, Smooth Particle Hydrodynamic method, (méthode SPH) peuvent utilisées. Afin de valider ces deux méthodes numériques, Eléments Finis et SPH, une première application consiste à la modélisation numérique du gonflement d’une membrane en caoutchouc. Pour cette application des données expérimentales sont disponibles pour la validation des résultats numériques. Pour le matériau caoutchouc, une loi de comportement de type Mooney Rivlin est utilisée. Pour la structure la méthode des éléments finis, formulation coque est utilisée, pour le fluide nous avons opté pour la méthode particulaire SPH. Concernant la sensibilité par rapport au maillage, plusieurs maillages pour la modélisation des domaines fluide et structure sont testés. Une bonne corrélation en terme de déplacement et de vitesse du centre de la membrane, entre les résultats numériques et données expérimentales a été observée. Cette application est publiée dans un journal international de rang A.La seconde formulation développée dans le manuscrit consiste à la modélisation et simulation numérique du problème en vue. La méthode de couplage de type Euler Lagrange, une formulation Eulerienne à maillage fixe pour le fluide et une formulation Lagrangienne pour la structure, est utilisée pour la modélisation de l’écoulement sanguin dans les artères. Différentes méthodes de pénalisation pour le couplage fluide structure ont été testées. Pour une meilleure consistance, plusieurs types de maillages fluide et structure ont été analysés, and des pas de temps vérifiant la condition de stabilité. Comme très peu de résultats expérimentant dans le domaine biomécanique sont disponibles. Les résultats numériques sont comparés aux résultats théoriques publiés dans la littérature. Une bonne corrélation entre les résultats numériques et les données théorique a été signalée dans le manuscrit. Ces résultats sont publiés dans un Journal International de rang A. / The aim of this thesis is to investigate blood flow in arteries using Fluid-Structure Interaction (FSI) numerical approach. There are different approaches which could be used, Finite Element (FE) method for modelling artery wall and either FE or Smoothed Particles Hydrodynamic (SPH) method for blood modeling. In order to investigate the appropriate numerical method to simulate the biomedical problem. both SPH and Finite Element methods were applied. Both methods were first validated for applications where experimental data are available. In order to validate the SPH and FEM method used to simulate the fluid domain and arteries, we investigate a specific problem concerning membrane inflation. Finite Element method with shell formulation was used for the membrane made of rubber material, and SPH Particle method for the fluid. This application has been selected since experimental data are available. Mooney Rivlin constitutive material law for hyper elastic incompressible material is used for the membrane and compressible air for the fluid. For mesh sensitivity and consistency, different meshes for the membrane and different particle numbers for SPH have been investigated, Good agreement with experimental data were obtained in term of displacement and velocity of the membrane center. The application is published in an international Journal with Index Citation. The second formulation developed in the manuscript concerns the fluid structure coupling problem we are interested in. Coupled Eulerian Lagrangian (CEL) approach for the modelisation of Pulse Wave Velocity (PWV) in large arteries. Different penalty methods for the fluid structure coupling have been investigated for good accuracy and consistency of the results using different fluid and structure meshes and time step satisfying stability condition. Since experimental data are not available for these approaches, numerical results were compared to the theoretical theory and previous work published in the literature. This work is published in an International Journal with Index Citation.

Cette thèse est dédiée à l’étude analytique et numérique des instabilités d’origine thermique et thermodiffusive de fluides visco-élastiques. L’objectif recherché est de contribuer à la compréhension de la dynamique qui résulte de la compétition entre plusieurs moteurs d’instabilités. En plus du caractère visco-élastique du fluide et de la présence d’un gradient de température vertical déstabilisant, d’autres sources d’instabilités viennent s’y ajouter : Le couplage "convection/écoulement de Poiseuille" d’une part, et le couplage "convection/effet Soret" inhérent aux mélanges binaires d’autre part. Deux configurations physiques sont alors considérées. La première partie de cette thèse est consacrée aux écoulements de convection mixte de type Rayleigh-Bénard-Poiseuille de fluides visco-élastiques, alors que la deuxième partie concerne l’effet de la thermodiffusion sur les instabilités de ces fluides saturant un milieu poreux. Le choix d’un milieu poreux est essentiellement motivé par la suggestion d’un protocole industriel de séparation possible des constituants d’une solution de polymères. / This thesis is dedicated to analytical and numerical study of thermal and thermodiffusive instabilities of viscoelastic fluids. The objective is to contribute to the understanding of the dynamics that results from the competition between different origins of instabilities. In addition to the viscoelastic nature of the fluid and the presence of a destabilizing vertical temperature gradient, other sources of instabilities have to be added: the coupling "convection/Poiseuille flow" on the one hand, and the coupling "convection/Soret effect" inherent to binary mixtures on the other hand. Two physical configurations are then considered. The first part of this thesis will be devoted to the Rayleigh-Bénard-Poiseuille mixed viscoelastic fluid convection, while the second part aims to identify the effect of thermodiffusion and viscoelasticity on convective instabilities in a porous medium. The choice of a porous medium in the second part is primarily motivated by the suggestion of an industrial protocol for separating the constituents of a polymer solution.

Ce travail présente les résultats numériques et une analyse des écoulements des fluides newtoniens et viscoélastiques du type Phan-Thien-Tanner (PTT), dans une conduite courbée de 180◦ avec section carrée et une région de sortie linéaire. Les équations de Navier-Stokes sont écrites en coordonnées orthogonales généralisées, et résolues en utilisant la méthode des volumes finis. On montre les résultats numériques des lignes de courant, du développement de tourbillons et des profils de vitesses dans la région de sortie, pour trois nombres de Dean (125 ; 137 et 150) et quatre nombres de Deborah différents (0,1 ; 0 2 ; 0 3 et 0, 4). Nous avons observé que la formation de paires supplémentaires de tourbillons contre rotatifs et les transitions des modèles de deux tourbillons au modèle de quatre tourbillons et au modèle de six tourbillons, ainsi que les dimensions des tourbillons et le moment de leur apparition dépendent du nombre de Dean et du nombre de Deborah. / This work presents the numerical results and the analysis of the flow behavior of Newtonian and viscoelastic fluids, type Phan-Thien-Tanner (PTT), in an 180◦ curved duct of square cross section and straight output zone. The Navier-Stokes equations were written in generalized orthogonal coordinates, and solved using the finite volume method. The numerical results of the streamlines, development of vortex and velocity profiles in the output region are shown, for Dean numbers (125 ; 137 ; and 150) and four different Deborah numbers (0,1 ; 0,2 ; 0,3 ; and 0,4). We observed that the formation of supplementary twin counter-rotating vortex and the transitions of two-vortex patterns to four-vortex patterns or six-vortex patterns, as well as the vortex dimensions and the apparition time depend of number Dean and number Deborah.

Ce travail de thèse sur les milieux poreux est axé sur deux parties. La première concerne l'étude numérique de l'écoulement d'un fluide Newtonien ou non-Newtonien au sein d'un système fluide/poreux. L'approche à un seul domaine qui consiste à écrire l'équation de Navier-Stokes incluant le terme de Darcy-Brinkman-Forchheimer est adoptée dans cette étude. La relation entre le gradient de pression et la vitesse débitante linéaire dans le cas de Darcy où le fluide est Newtonien, est obtenue. Cette relation est étendue dans le cas non-Darcy où le fluide est non Newtonien. L'influence des nombres de Darcy et de Forchheimer sur la structure de l'écoulement est montrée. Dans la seconde partie, une étude de stabilité linéaire et numérique de la convection naturelle de fluides viscoélastiques saturant une couche poreuse horizontale chauffée par un flux constant est réalisée. Une étude d'instabilité primaire et secondaire nous a permis de montrer que pour un fluide Newtonien, la convection monocellulaire perd sa stabilité au profit des rouleaux longitudinaux. Dans le cas des fluides viscoélastiques, on trouve que l'élasticité du fluide induit la sélection des rouleaux transversaux propagatifs. Une solution numérique basée sur un schéma aux différences finies est venue conforter ces résultats analytiques. / The present thesis on porous media concentrates in two parts. The first concerns the numerical study of the flow of a Newtonian or a non-Newtonian fluid within a fluid/porous system. The approach of a single domain, which consists of/in writing the Navier-Stokes equation including the Darcy-Brinkham-Forchheimer term, is chosen in this study. The linear relation between the pressure gradient and the bulk velocity in the Darcy case, for which the fluid is Newtonian, is obtained. This relation is extended to the non-Darcy case, for which the fluid is non-Newtonian. The influence of Darcy and Forchheimer numbers on the structure of the flow is presented. In the second part, linear stability and numerical analysis of the natural convection of viscoelastic fluids saturating a horizontal porous layer heated by a constant flux is performed. A primary and secondary instability study allowed to show that, for a Newtonian fluid, the unicellular convection loses its stability to the benefit of longitudinal rolls. In the case of viscoelastic fluids, the elasticity of the fluid leads to the selection of propagation transverse rolls. A numerical solution based on a finite difference scheme has reinforced these analytical results.

Dans ce travail de thèse tant théorique que numérique, on étudie les différentes instabilités pouvant se développer dans un milieu poreux saturé par un fluide viscoélastique et chauffé par le bas. La formulation mathématique des équations de ce problème repose sur la loi phénoménologique de Darcy généralisée à un fluide viscoélastique vérifiant l’approximation de Boussinesq. Ce problème admet une solution de conduction, et on trouve que deux types de structures sont susceptibles d’apparaître lorsque l’état de conduction perd sa stabilité : des structures stationnaires et des structures oscillatoires. Les seuils d’apparition de ces structures sont étudiés en fonction des paramètres adimensionnés du problème, à savoir le nombre de Rayleigh, les temps de relaxation et de retardation associés à l’élasticité du fluide. Une étude linéaire et non linéaire est donc menée. Il est intéressant de noter qu’une compétition entre les structures stationnaires et oscillatoires peut exister au voisinage d’un point appelé point de codimension 2. Une analyse non linéaire est donc menée au voisinage de ce point et est confrontée aux résultats issus des simulations numériques. Enfin, s’appuyant sur les propriétés de mélanges binaires des fluides viscoélastiques, une étude théorique est réalisée et nous montrons qu’il y a une compétition entre deux régimes : un régime où la viscoélasticité est dominante et un autre où l’aspect mélange binaire l’emporte. Ce résultat permet d’expliquer certaines observations expérimentales. / In this theoretical and numerical work, we study differents instabilities which can develop in a porous media saturated by viscoelastic fluid and heated from below. The mathematical formulation of the equations of this problem is based on phenomenological Darcy law generalized to a viscoelastic fluid verifying Boussinesq estimate. This problem admits a solution of conduction, and we find that two types of structures may appear when the conduction state loses his stability : stationary and oscillatory structures.The apparition thresholds of these thermo-convectives structures are studied and depend on the non-dimensionnalized parameters of the problem, Rayleigh number, relaxation and retardation time associated to the fluid elasticity. A linear and non linear stability is also realized. It is interesting to note that it can have a competition between stationary and oscillatory structures near a point named codimension 2 point. A linear analysis is also realized near this point and is compared to the numerical simulation results.Finally, by taking into account binary mixtures properties of the viscoelastic fluids, a theoretical study is realized and we show that there is a competition between two states : a state in which viscoelasticity is dominant and another state in which binaries properties are also dominants. This result can explain experimental observations.

Η περιγραφή των αλληλεπιδράσεων που αναπτύσσονται μεταξύ ενός συνόλου σωμάτων που είτε είναι ακίνητα είτε κινούνται και ενός ρευστού, ρέοντος ή ηρεμούντος, βρίσκεται στον πυρήνα ενός μεγάλου αριθμού επιστημονικών προσπαθειών. Η μοντελοποίηση των αλληλεπιδράσεων αυτών μπορεί να γίνει με πολλούς τρόπους, καθένας από τους οποίους εμφανίζει πλεονεκτήματα και μειονεκτήματα, ανάλογα με το είδος του συστήματος στο οποίο επιχειρείται να εφαρμοστεί. Σκοπός της εργασίας αυτής είναι η ανάδειξη ενός εξ αυτών που είναι η Μέθοδος των Προτύπων Ροών. / The description of interactions between a number of solid particles, which may be in motion or fixed in space and a fluid which in turn may either be flowing or resting, is placed at the core of a large number of scientific publications. Modelling such interactions may follow a variety of approaches. Each approach has its own pros and cons, regarding the nature and the specific features of the system which is aimed to describe. Current thesis goal is to describe one of these approaches, which is the approach of Singularities or the Method of Elementary Flows.

It is well-known that, increasing the flow rate in polymer extrusion, the flow becomes unstable and the smooth extrudate surface becomes wavy and disordered to an increasing degree. In order to investigate the mechanisms responsible for these instabilities we perform a linear stability analysis of the steady extrusion of a viscoelastic fluid flowing through a planar die under creeping flow conditions. We consider the Phan-Thien-Tanner (PTT) model to account for the viscoelasticity of the material. We employ the mixed finite element method combined with an elliptic grid generator to account for the deformable shape of the interface. The generalized eigenvalue problem is solved using Arnoldi’s algorithm. We perform a thorough parametric study in order to determine the effects of all material properties and rheological parameters. We investigate in detail the effect of interfacial tension and the presence of a deformable interface. It is found that the presence of a finite surface tension destabilizes the flow as compared to the case of the stick-slip flow. We recognize two modes which are found to become unstable beyond a critical value of the Weissenberg number and perform an energy analysis to examine the mechanisms responsible for the destabilization of the flow and compare against the mechanisms that have been suggested in the literature. / --

Στόχος της παρούσας διδακτορικής διατριβής είναι η μελέτη χρονικά μεταβαλλόμενων και ασταθών υπό συνθήκες ροών ιξωδοπλαστικών και ιξωδο-ελαστικών ρευστών λόγω της σημασίας που παρουσιάζουν τα υλικά αυτά σε διάφορες βιομηχανικές διεργασίες. Στα πλαίσια της εργασίας αυτής μελετήθηκε αρχικά το πρόβλημα της συμπίεσης ιξωδοπλαστικού ρευστού μεταξύ δύο δίσκων το οποίο είναι ένα βασικό πείραμα ρεολογικού χαρακτηρισμού μη Νευτωνικών υλικών. Η χρονικά μεταβαλλόμενη προσομοίωση επέτρεψε τον προσδιορισμό των σημαντικών διαφορών μεταξύ των δύο εκδοχών του συγκεκριμένου ρεολογικού πειράματος (δίσκοι κινούμενοι υπό σταθερή ταχύτητα ή υπό σταθερή δύναμη) που ήταν αδύνατος με τα προηγούμενα μοντέλα, ενώ παρουσιάζεται και πλήρης παραμετρική μελέτη της διεργασίας αυτής. Το υπόλοιπο μέρος της παρούσας εργασίας αφιερώθηκε στη μελέτη της διεργασίας εκβολής ιξωδοελαστικών υλικών η οποία συναντάται ευρέως στη βιομηχανία μορφοποίησης πολυμερών. Παρουσιάζεται μια πλήρης παραμετρική ανάλυση για την εκβολή ενός ιξωδοελαστικού ρευστού Phan-Tien & Tanner (PTT) από ένα δακτυλιοειδή αγωγό προκειμένου να εξετάσουμε την επίδραση της γεωμετρίας του αγωγού, και των ρεολογικών ιδιοτήτων του ρευστού. Επιπλέον ένα σημαντικό πρόβλημα που εμφανίζεται κατά την εκβολή από ένα κυλινδρικό, επίπεδο ή δακτυλιοειδή αγωγό είναι η εμφάνιση ασταθειών οι οποίες επηρεάζουν σημαντικά την ποιότητα του τελικού προϊόντος. Προκειμένου να διερευνηθεί ο μηχανισμός που προκαλεί την εμφάνιση αυτών των ασταθειών πραγματοποιήθηκε γραμμική ανάλυση ευστάθειας για τη αξονοσυμμετρική και καρτεσιανή δισδιάστατη ροή επικόλλησης-ολίσθησης (stick-slip flow) ενός ρευστού PTT. Γύρω από τη λύση μόνιμης κατάστασης πραγματοποιείται γραμμική ανάλυση ευστάθειας, και υπολογίζονται οι ιδιοτιμές του γενικευμένου προβλήματος με τη βοήθεια της μεθόδου Arnoldi. Έτσι προσδιορίζεται η επίδραση των ιδιοτήτων του υλικού στην ευστάθεια ή μη της ροής. / The purpose of this dissertation is to perform a study on time-dependent and conditionally unstable flows of viscoplastic and viscoelastic fluids because of their importance in various industrial processes. A typical experiment which is used widely for the rheological characterization of non-Newtonian fluids is squeeze flow which is the subject of the first chapter of this dissertation. The time-dependent simulation permitted the determination of the distinct differences between the two versions of this rheological experiment (disks moving with constant velocity or under constant force) which was impossible with the models used up to now, while a complete parametric analysis of this process is presented. The rest of this work focuses in the important problem of the extrusion process of a viscoelastic fluid, which is frequently encountered in the polymer industry. A complete parametric analysis is presented for the extrusion of a Phan-Thien Tanner (PTT) fluid from an annular die in order to examine the effect of the die geometry and the rheological properties of the fluid on this process. Moreover, it is widely known that during the extrusion from a cylindrical, planar or an annular die various flow instabilities may arise affecting significantly the quality of the final product. In order to investigate the mechanism, which causes the appearance of these instabilities, a linear stability analysis is performed for the cylindrical or planar stick-slip flow of a PTT fluid. The stability analysis is performed around the steady state solution and the eigenvalues of the generalized problem are calculated using the Arnoldi algorithm. With this method the effect of the various rheological properties of the fluid on the stability of the flow is determined.

La méthode de Boltzmann sur réseau est une formulation discrète particulière de l'équation de Boltzmann. Depuis ses débuts, il y a trente ans, cette méthode a gagné une certaine popularité, et elle est maintenant utilisée dans presque tous les problèmes habituellement rencontrés en mécanique des fluides notamment pour les écoulements multi-espèces. Dans le cadre de ce travail, une force de friction intermoléculaire est introduite pour modéliser les interactions entre les molécules de différent types causant principalement la diffusion entre les espèces. Les phénomènes de dissipation visqueuse (collision usuelle) et de diffusion moléculaire (force de friction intermoléculaire) sont séparés et peuvent être ajuster indépendamment. Le principal avantage de cette stratégie est sa compatibilité avec des optimisations de la collision usuelle et les opérateurs de collision avancés. Adapter un code mono-espèce pour aboutir à un code multi-espèces est aisé et demande beaucoup moins d’effort comparé aux précédentes tentatives. De plus, il n’ y a pas d’approximation du mélange, chaque espèce a ses propres coefficients de transport pouvant être calculés à l’aide de la théorie cinétique des gaz. En général, la diffusion et la convection sont vus comme deux mécanismes séparés : l’un agissant sur la masse d’une espèce, l’autre sur la quantité de mouvement du mélange. En utilisant une force de friction intermoléculaire, la diffusion et la convection sont couplés par l’intermédiaire la quantité de mouvement de chaque espèce. Les mécanismes de diffusion et de convection sont intimement liés dans de nombreux phénomènes physique tel que la digitation visqueuse.L’instabilité de digitation visqueuse est simulée en considérant dans un milieu poreux deux espèces dans des proportions différentes soit un mélange moins visqueux déplaçant un mélange plus visqueux. Les principaux moteurs de l’instabilité sont la diffusion et le contraste de viscosité entre les espèces. Deux stratégies sont envisagées pour simuler les effets d’un milieu poreux. Les méthodes de rebond partiel et de force de Brinkman bien que basées sur des approches fondamentalement différentes donnent dans notre cas des résultats identiques. Les taux de croissance de l’instabilité calculés à partir de la simulation coïncident avec ceux obtenus à partir d’analyses de stabilité linéaire. L’évolution de la longueur de mélange peut être divisée en deux étapes dominées d’abord par la diffusion puis par la convection. La physique de la digitation visqueuse est ainsi correctement simulée. Toutefois, les effets de diffusion multi-espèces ne sont généralement pas pris en compte lors de la digitation visqueuse de trois espèces et plus. Ces derniers ne sont pas négligeable puisque nous mettons en avant une configuration initialement stable qui se déstabilise. La diffusion inverse entraîne la digitation dont l’impact dépend de la diffusion entre les espèces. / The lattice Boltzmann method (LBM) is a specific discrete formulation of the Boltzmann equation. Since its first premises, thirty years ago, this method has gained some popularity and is now applied to almost all standard problems encountered in fluid mechanics including multi-component flows. In this work, we introduce the inter-molecular friction forces to take into account the interaction between molecules of different kinds resulting primarily in diffusion between components. Viscous dissipation (standard collision) and molecular diffusion (inter-molecular friction forces) phenomena are split, and both can be tuned distinctively. The main advantage of this strategy is optimizations of the collision and advanced collision operators are readily compatible. Adapting an existing code from single component to multiple miscible components is straightforward and required much less effort than the large modifications needed from previously available lattice Boltzmann models. Besides, there is no mixture approximation: each species has its own transport coefficients, which can be calculated from the kinetic theory of gases. In general, diffusion and convection are dealt with two separate mechanisms: one acting respectively on the species mass and the other acting on the mixture momentum. By employing an inter-molecular friction force, the diffusion and convection are coupled through the species momentum. Diffusion and convection mechanisms are closely related in several physical phenomena such as in the viscous fingering instability.A simulation of the viscous fingering instability is achieved by considering two species in different proportions in a porous medium: a less viscous mixture displacing a more viscous mixture. The core ingredients of the instability are the diffusion and the viscosity contrast between the components. Two strategies are investigated to mimic the effects of the porous medium. The gray lattice Boltzmann and Brinkman force models, although based on fundamentally different approaches, give in our case equivalent results. For early times, comparisons with linear stability analyses agree well with the growth rate calculated from the simulations. For intermediate times, the evolution of the mixing length can be divided into two stages dominated first by diffusion then by convection, as found in the literature. The whole physics of the viscous fingering is thus accurately simulated. Nevertheless, multi-component diffusion effects are usually not taken into account in the case of viscous fingering with three and more species. These effects are non-negligible as we showcase an initial stable configuration that becomes unstable. The reverse diffusion induces fingering whose impact depends on the diffusion between species.

Ce travail est consacré à l'étude de la viscosité de cisaillement par simulation numérique de dynamique moléculaire classique à l'équilibre avec une attention particulière à l'influence des hautes pressions sur cette propriété. La viscosité est obtenue à partir des trajectoires générées par ces simulations et en appliquant la formule de Green-Kubo. Un large panel de systèmes a été étudié, allant de fluides atomiques purs à un liquide moléculaire, en passant par des mélanges binaires. En premier lieu, nous nous sommes focalisés sur les métaux alcalins. La conclusion majeure de cette étude est que la viscosité des alcalins a un comportement universel sur une large plage du diagramme de phases. Par ailleurs, sur cet intervalle, la relation universelle que nous avons proposée permet de prédire la valeur de la viscosité de n'importe quel élément avec une incertitude inférieure à 10%. La validité de la relation de Stokes-Einstein, reliant le coefficient d'autodiffusion à la viscosité, a également été vérifiée. Une étude systématique a ensuite été menée sur des mélanges modèles de type Lennard-Jones an de tester l'influence des paramètres d'interaction sur le comportement de la viscosité. Une estimation théorique basée sur le modèle de fluide effectif pur a été proposée. D'autre part, la relation de Stokes-Einstein a été étendue aux mélanges avec succès. Ces observations ont été confrontées aux cas de deux alliages réels : K-Cs et Li-Bi. Pour finir, une étude préliminaire a été entreprise sur l'eau en modélisant les interactions par deux potentiels : SPC/E, non polarisable, et BK3, polarisable. L'effet de l'introduction de la polarisabilité sur le calcul de la viscosité a été étudié. La validité des relations de Stokes-Einstein et de Stokes-Einstein-Debye, faisant intervenir la rotation de la molécule, a été évaluée à très haute pression / This work is devoted to the study of the shear viscosity by numerical simulation of equilibrium classical molecular dynamics with a particular attention to the influence of high pressures on this property. From trajectories generated by these simulations and using the Grenn-Kubo formula, the viscosity is obtained. A broad range of systems has been studied, covering from pure atomic fluids to a molecular liquid, as well as binary mixtures. First, we focused on alkali metals. The main outcome of this study is that the viscosity of these metals has a universal behavior over a wide range of phase diagram. Furthermore, over this interval, the universal relation that we have proposed permits the prediction of the viscosity value of any elements with an uncertainty lower than 10%. The validity of the Stokes-Einstein relation, connecting the self-diffusion coefficient and the viscosity, has also been verified. Then, a systematic study has been carried out on model mixtures of Lennard-Jones fluids to test the influence of interaction parameters on the viscosity behavior. A theoretical estimation based on the effective one-component fluid model has been proposed. Moreover, the Stokes-Einstein relation has been successfully extended to mixtures. These observations have been compared with two real alloys: K-Cs and Li-Bi. Lastly, a preliminary study on water has been undertaken by modeling the interactions with two models: SPC/E, non-polarizable and BK3, polarizable. The effect of the introduction of the polarizability on the viscosity has been studied. The validity of Stokes-Einstein and Stokes-Einstein-Debye, involving the rotation of the molecule, has been evaluated under very high pressure