Le transport des bruts paraffiniques, et tout particulièrement leur remise en écoulement après un arrêt, dans de longues conduites sous-marines soumises à de basses températures, peut être difficile du fait de l'augmentation de leur viscosité. Le comportement rhéologique d´une huile paraffinique modèle, possédant des propriétés macroscopiques d'écoulement analogues à celles des bruts paraffiniques, est d'abord analysé en utilisant la vélocimétrie par imagerie par résonance magnétique associée à des mesures de contrainte de cisaillement au sein d'une géométrie Couette. Nous montrons que lors d'un écoulement forcé à température constante le matériau subit une déstructuration irréversible qui dépend de l'intensité du cisaillement. Ainsi la contrainte apparente critique permettant l'écoulement du matériau dépend de l'histoire thermique et d'écoulement subie par le matériau. Nous étudions ensuite le comportement rhéologique complet de deux bruts réels à partir de différents types de tests rhéométriques (fluages, redémarrage, régime permanent, changement brusque de vitesse) pour différentes histoires d'écoulement, notamment pendant la période de refroidissement. Le comportement détaillé du matériau en régime transitoire ainsi observé peut alors être modélisé. De plus les variations du seuil de contrainte en fonction de l'histoire thermique et de l'écoulement sont aussi décrites, ce qui nous donne le champ de contrainte seuil dans la conduite à l'état initial. Le modèle dans son ensemble est finalement implémenté dans un code de calcul pour simuler le redémarrage de l´écoulement d'un brut paraffinique dans une conduite réelle / Transporting waxy crude oils through long pipelines at low temperatures may be challenging, particularly its flow restart after a pipeline shut-in, due to its viscosity increase. The rheological behavior of a model waxy oil with macroscopic flow properties analogous to waxy crude oils is first analyzed using Magnetic Resonance Imaging velocimetry associated to stress measurements in a Couette geometry. While flowing at constant temperature, major irreversible structure break depending on shear intensity are observed. Thus, the critical apparent shear stress beyond which the material flows depends on the thermal and flow histories of the oil. Next, the rheological behavior of two waxy crude oils is studied using rheometrical tests (creep tests, flow restarts, abrupt changes of shear rate and steady flow) after different flow histories, notably during the cooling process. Then, those experimentally observed trends are modeled. Additionally, a comprehensive study of the yield stress in function of flow and temperature histories is presented. It provides an approach for describing the yield stress field inside the pipeline at the flow restart moment. Finally, the entire rheological model is implemented in the computational code for simulating waxy crude oils flow restart of a real scale pipeline

Solides élastiques au repos, les fluides à seuil s’écoulent comme un liquide au-delà d’une certaine contrainte. Plusieurs applications industrielles concernent l’écoulement de ces fluides dans des milieux poreux. On peut citer par exemple les émulsions dans le processus de récupération du pétrole, les opérations de cimentation dans le sol, ou le nettoyage d’un sol contaminé par une boue. Pour ces applications, il est nécessaire de connaitre la pression nécessaire pour un débit voulu à la sortie du milieu poreux. Dans de tels cas, l’écoulement est perturbé par la complexité de la géométrie. Les modèles développés pour décrire la loi de Darcy supposent une loi rhéologique appliquée localement, mais ces modèles décrivent mal ce type d’écoulement. De plus, des effets complexes peuvent s’ajouter comme le glissement à la paroi ou la thixotropie. Dans cette thèse, nous étudions l’écoulement de carbopol (ETD 2050) à travers différentes géométries. Tout d’abord au rhéomètre, nous montrons que le fluide, sous certaines conditions, correspond bien à un fluide à seuil modèle. Nous démontrons que le protocole expérimental utilisé est très important et qu’un comportement thixotropique peut apparaitre s’il n’est pas respecté. Ce comportement apparait notamment lorsque le fluide reste sous le seuil, l’impact augmentant avec le temps d’attente. Ensuite, nous comparons la loi d’écoulement obtenue au rhéomètre à l’écoulement dans un canal droit obtenu par microfabrication. Nous montrons alors l’importance du glissement proche du seuil et ses conséquences sur la loi d’écoulement. Enfin nous étudions l’écoulement du carbopol dans un milieu poreux. Le milieu poreux de 5x5cm est obtenu par microfabrication. La largeur moyenne des canaux est égale à celle du canal droit. Nous avons développé une nouvelle méthode de mesure des champs de vitesse. Nous montrons l’apparition d’une chenalisation de l’écoulement à travers quelques canaux du milieu poreux. Nous comparons ensuite la loi d’écoulement du milieu poreux à celle obtenue dans le canal droit. On remarque que la vitesse d’écoulement est plus faible dans le milieu poreux que dans le canal droit. / Elastic solids at rest, yield stress fluids flow like a liquid beyond a certain stress. Many industrial applications required the flow of these fluids in porous media, for example: the emulsion flow in oil recovery processes, the cementing operations in the ground, or the cleaning of sludge in a contaminated soil. For many applications, it could be interesting to know the pressure required for a desired flow rate. In such cases, the flow behavior of the fluid is complicated by the complexity of the geometry. The models developed to describe Darcy's law assume a rheological law applied locally, but these models poorly describe this type of flow. Furthermore, complex effects can be added like the wall slip or the thixotropy. In this thesis, we study the flow of carbopol (ETD 2050) through different geometries. First we show that the fluid, for some conditions, corresponds to model yield stress fluids. The experimental protocol used is very important and a thixotropic behavior can appear if it is not respected. This behavior appears especially when the fluid remains below the yield stress, the impact increases with the waiting time. We then compare the flow law obtained by rheometer in a straight channel obtained by microfabrication. We show the importance of the wall slip near the yield stress and the impact on the flow law. Finally, using a new method to measure the velocity fields developed during this thesis, we study the flow of carbopol in a porous medium. This porous medium of 5x5cm is obtained by microfabrication. The mean width of the channels is equivalent to the one of the straight channel. We show the emergence of a channeling flow through some channels of the porous medium. We then compare the flow law of the porous medium to the one obtained in the straight channel. It can be observed that the flow rate is lower in the porous medium than in the straight channel.

The rheological properties of cement-based grouts play a key role in determining the final spread in grouted rock formations. Rheologically, cement grouts are known to be complex thixotropic fluids, but their steady flow behavior is often described by fitting the simple Bingham constitutive law to flow curve data. The resultant Bingham parameters are then used in grouting design of e.g. tunnels, to estimate the penetration length. Since cement grouts are thixotropic suspensions, the interpretation of their flow curves as obtained from flow sweeps in concentric cylinder rotational rheometers is often complicated by: the presence of wall slip, sedimentation and unstable flow at low shear rates. A systematic approach to study these effects within the constraints of the concentric cylinder geometry (Couette) and for different cement grout concentrations was carried out as part of the Licentiate research work. Of particular interest was the influence of geometry and flow sweep measurement interval on flow curves, including the characteristic unstable flow branch that appears at applied shear rates that are below the critical shear rate. The unstable flow branch observed below the critical shear rate has been described as a characteristic feature in the flow curves of thixotropic suspensions, e.g. cement grouts, laponite. From a practical standpoint, this information can then be readily used to improve rheological measurements of cement grouts. The existence of the critical shear rate below which no stable flow occurs, plus the complex wall slip phenomenon are then discussed by considering how they affect actual spread in rough and smooth rock fractures. Another major part of the research presented in this thesis relates to the measurement of model yield stress fluid (YSF), i.e. Carbopol, velocity profiles within the radial flow geometry. Radial flow between parallel plates, is an idealized fundamental flow configuration that is often used as a basis for grout spread estimation in planar rock fractures. Compared to other flow configurations with YSFs, e.g. channels, only a limited amount of work has presented analytical solutions, numerical models and especially experimental work for radial flow. Thus, as a first step towards more systematic studies of the plug flow region of YSFs in radial flow the current work presents the design, manufacture and for the first time velocity profile measurements that were conducted by using the pulsed Ultrasound Velocity Profiling (UVP) technique. The current observations for tests carried out with different disk spacings and flow rates show a distinct plug region, coupled with wall slip effects for the Carbopol model YSF fluid that was used. The theoretically predicted velocity profiles and the measured ones agree reasonably well, and the main discrepancies are discussed. Future studies, would then be targeted at improving the current experimental setup, for detailed measurements of the plug flow region along the radial length, which remains a challenging issue for studies on YSFs and engineering applications such as rock grouting design. / Cementbaserade injekteringsmedels reologiska egenskaper har en stor påverkan på strömning och inträngningslängd i sprickigt berg. Medlens reologi är komplex, inklusive tixotropi, men strömningen beskrivs ändå oftast med den enkla linjära Bingham modellen i injekteringssammanhang. De två parametrarna från denna modell, flytgräns och viskositet, används sedan inom injekteringsdesign, för t.ex. tunnlar och dammar, för att bedöma inträngningen. Eftersom cementbaserade medel är tixoptropa suspensioner försvåras utvärderingen vid mätning med konventionella rotationsviskometrar på grund av glidning vid fasta begränsningsytor, sedimentation/separation av partiklarna och instabila flöden vid låga deformationshastigheter. En systematisk mätprocedur för att studera ovanstående problem med rotationsviskometer och koncentriska cylindrar samt olika vanliga vattencementtal, har utförts inom ramen för detta licentiatarbete. Av särskilt intresse har varit att studera effekten av olika geometrier och tidsintervallet mellan mätningarna, inklusive den instabila delen av flödeskurvan då deformationshastigheten är lägre än ett kritiskt värde. Denna del av kurvan har i litteraturen beskrivits som karakteristisk för tixotropa suspensioner, som t.ex. cementbaserade injekteringsmedel. Praktiskt kan ovanstående kunskap användas för att förbättra mätningen av de reologiska egenskaperna. Existensen av en kritisk deformationshastighet under vilken det inte finns något stabilt flöde, i kombination med glidning vid fasta begränsningsytor, diskuteras särskilt med hänsyn till dess påverkan på faktisk inträngning i släta och råa bergsprickor. Ett annat fokus i licentiatarbetet har varit att studera icke-Newtonska modellvätskors (Carbopol) radiella strömning mellan parallella plattor. Denna typ av strömningsgeometri används ofta som en idealiserad konfiguration för strömning i bergsprickor. I jämförelse med andra enklare geometrier, finns endast en begränsad forskning utförd för denna geometri både då det gäller analytiska och numeriska beräkningar men framförallt då det gäller experiment. Som ett första steg inför en mer systematisk undersökning av icke-Newtonsk radiella strömning presenteras i detta arbete framtagandet av en fysisk laboratoriemodell där hastighetsprofilerna mellan plattorna för första gången visualiserats med hjälp av ultraljud. De utförda mätningarna med tre olika öppningar mellan plattorna sam tre olika värden på det konstanta flödet, visar på en distinkt plugg som är ett resultat av vätskans flytgräns samt glidning i gränsskiktet mellan vätskan och plattornas fasta begränsningsytor. En jämförelse mellan uppmätta hastighetsprofiler och analytiskt beräknade diskuteras där resultaten överensstämmer relativt väl, med beaktande av de långtgående förenklade antaganden som krävs för beräkningarna. Fortsatta studier kommer att fokuseras på att förbättra laboratoriemodellen för en mer detaljerad studie av icke-Newtonska vätskors strömning och hur pluggen utvecklas under den radiella inträngningen, vilket fortsättningsvis är av betydelse för design av injektering i bergsprickor. / <p>QC 20190521</p>

Cette thèse, composée de deux parties, porte sur l’étude d’instabilités au sein d’un milieugranulaire.La première partie de cette thèse est consacrée à l’instabilité de tôle ondulée, c’est-àdirel’apparition d’un motif de rides sur une piste soumise au passage répété de véhicules.Nos travaux expérimentaux ainsi que des simulations numériques de dynamique moléculairenous ont permis de progresser dans la compréhension de ce phénomène. Pour uneroue rigide tirée à vitesse constante, nous avons étudié quantitativement l’impact de lacohésion du milieu granulaire sur les caractéristiques de l’instabilité : vitesse critique, longueurd’onde et taux de croissance. Nous avons enfin entrepris des études préliminairessur un fluide à seuil simple : un microgel de carbopol.Dans la seconde partie, nous étudions la ségrégation granulaire au sein des astéroïdeslâches. Plusieurs indices montrent qu’un grand nombre d’entre eux sont des empilementsde grains, sans cohésion interne et liés par la gravité. Leur répartition n’est pas homogène,avec des zones lisses recouvertes de fines poussières et d’autres où s’accumulent les grosblocs. Nous avons simulé numériquement un empilement granulaire et nous l’avons soumisà des secousses répétées : nous avons réussi à reproduire la ségrégation. Nous avonsétudié l’influence des paramètres physiques et numériques sur le niveau de ségrégation etla dynamique du phénomène. Nous nous sommes aussi intéressés aux phénomènes physiquespouvant causer la ségrégation : convection granulaire, tamisage cinétique, pressionde déplétion, etc. Enfin, nous avons mis en évidence une instabilité azimutale qui peutexpliquer les observations faites pendant les missions spatiales. / This thesis, consisting of two parts, is focusing on instabilities into a granular assembly.The first part deals with washboard road instability, which is the growth of a ripplepattern on a track subjected to repeated passages of vehicles. Our experimental work aswell as soft spheres numerical simulations provide us a better understanding of this phenomenon.In the case of a rigid wheel dragged at constant velocity, we quantitatively studiedthe impact of the cohesion into the granular media on the main features of the instability :critical velocity, wavelength and growth rates. Finally, we have begun preliminary studieson a yield stress fluid, namely a carbopol microgel.The second part is devoted to the granular segregation in asteroids. Several keys showthat a significant part of them are rubble-piles, without internal cohesion and held togetherby their weak self-gravity. The distribution of grains is heterogeneous : while some regionsconsist in fine sand or powder, large boulders seem to accumulate in other parts. Wenumerically simulated a granular pile and subjected it to repeated quakes : we managedto reproduce segregation. We studied the influence of physical and numerical parameterson the segregation level and the dynamics of the phenomenon. We also investigate variouscauses of the segregation : granular convection, kinetic sieving, depletion pressure, etc.Finally, we highlighted an azimuthal segregation which can explain observations duringspatial missions.

Cette thèse présente une étude expérimentale de l'impact de gouttes de fluides à seuil. Au-delà des applications (impression à jet d’encre solide, modélisation d’impact solide à grandes vitesses), cette étude permet de sonder le rôle de l'élasticité sur le comportement à temps court de ces fluides complexes. D'abord, nous nous sommes intéressés aux impacts sur une surface rigide. L'utilisation de fluides à seuil modèles (solutions concentrées d’argiles, micro-gel de Carbopol) et de surfaces d'impact variées (partiellement mouillante ou super-hydrophobe), révèle une grande variété de comportements, allant de l'étalement viscoplastique irréversible jusqu'à des déformations élastiques géantes. Un modèle minimal d'étalement inertiel, incluant une rhéologie élasto/viscoplastique, permet de décrire dans un cadre unique les principaux régimes observés. Au cours de cette étude, nous avons mis en évidence un phénomène spécifique avec le Carbopol : pour des grandes vitesses d'impact, on observe un étalement beaucoup plus grand sur des surfaces rugueuses hydrophobes que sur des surfaces lisses. Cette réduction apparente du frottement basal est discutée en termes de longueur de glissement et d'instabilité de « splash ». Enfin, nous avons étudié l’impact d'une goutte de fluide sur un sol constitué du même fluide, en utilisant un fluide à seuil transparent (Carbopol). La combinaison de lois d'échelle, d’expériences en « micro-gravité » et de mesures locales du champ de déformation montre que la dynamique du cratère transitoire est dominée par l’élasticité, même au-delà du seuil d’écoulement. Ces résultats pourraient avoir des implications dans le contexte des impacts de météorites en astrophysique. / This thesis presents an experimental study on the drop impact of yield-stress fluids. Beyond applications (solid ink-jet printing, lab modelling of high-speed collision of solids), this study offers a mean to probe the role of the elasticity on the short-time behaviour of these complex fluids. We have first studied drop impacts on solid rigid surfaces. Using different model yield-stress fluids (clay suspensions, Carbopol micro-gel) and impacted surfaces (partially wettable, super-hydrophobic), we have observed a rich variety of behaviours ranging from irreversible viscoplastic coating to giant elastic spreading and recoil. A minimal model of inertial spreading, including an elasto/viscoplastic rheology, allows explaining in a single framework the different regimes and scaling laws. In this study, we identified a specific phenomenon with Carbopol: for large impact velocities, the drop spreads much more on rough hydrophobic surfaces than on smooth surfaces. This apparent reduction of the basal friction is discussed in terms of slip length and splash instability. Endly, we investigated the impact of a drop onto a pool of the same fluid, using a transparent yield-stress fluid (Carbopol). The combination of scaling laws, micro-gravity experiments and local deformation measurements shows that the transient crater is dominated by elasticity, even beyond the flow threshold. These results could have implications for impact cratering in Planetary Sciences.

No description available.

Nous présentons une contribution autour de la modélisation des écoulements viscoplastiques. En vue d'applications réalistes telle que la simulation numérique des coulées de lave volcanique, le travail se concentre particulièrement sur les fluides complexes dont la rhéologie dépend fortement de grandeurs physiques telle que la température ou la concentration en particule. Nous développons un nouvel algorithme de résolution numérique des équations de Herschel-Bulkley combinant une méthode de Lagrangien augmenté à paramètre d'augmentation variable, une méthode des caractéristiques d'ordre 2 et une adaptation de maillage automatique. Sur des problèmes stationnaires ou en évolution tel que le problème test de la cavité entraînée, il apporte une solution efficace pour garantir à la fois une précision numérique élevée et un temps de calcul raisonnable. Cet algorithme est ensuite étendue et adapté au cas des rhéologies non-isothermes et aux suspensions. Concernant la simulation numérique des coulées de lave volcanique, nous détaillons une méthode de réduction par analyse asymptotique des équations de Herschel-Bulkley pour des écoulements de faible épaisseur sur une topographie arbitraire. Elle permet alors de décrire ces écoulements tridimensionnels de fluides viscoplastiques à surface libre par des équations bidimensionnelles surfaciques. Cette approche est ensuite étendue au cas non-isotherme en y ajoutant l'équation de la chaleur et des dépendances thermiques sur la rhéologie. Par intégration verticale de l'équation de la chaleur, on retrouve un modèle bidimensionnel. Le modèle non-isotherme est validé sur une expérience de dôme réalisée en laboratoire et une simulation numérique est réalisée autour d'une coulée qui a eu lieu sur le volcan du Piton de la Fournaise à la Réunion, en décembre 2010. La comparaison donne des résultats qui sont de notre point de vue satisfaisants et encourageants. / We present a contribution about modeling of viscoplastic flows. For realistic applications such as numerical simulation of volcanic lava flows, the work focuses particularly on complex fluids whose rheology strongly depends on physical quantities such as temperature or the particle concentration. We develop a new numerical resolution algorithm of Herschel-Bulkley's equations combining an augmented Lagrangian method with variable augmentation parameter, a second order characteristic method and an auto-adaptive mesh procedure. On stationary or evolving problems as the lid-driven cavity flow benchmark, it provides an effective solution to ensure both a high numerical accuracy within a reasonable computing time. This algorithm is then extended and adapted to the case of non-isothermal rheological and suspensions. On the numerical simulation of volcanic lava flows, we describe a method of reducing by asymptotic analysis of the Herschel-Bulkley's equations for thin flows on arbitrary topography. It allows to describe the three-dimensional flows of viscoplastic fluid with free surface by bidimensional surface equations. This approach is then extended to the non-isothermal case by adding the heat equation and thermal dependencies on rheology. By vertical integration of the heat equation, a two-dimensional model is maintained . The non-isothermal model is validated on a laboratory experiment of dome and a numerical simulation is performed on a December 2010 Piton de la Fournaise lava flow from La Réunion island. In our view, the comparison gives satisfactory and encouraging results.

Nous présentons une contribution autour de la modélisation des écoulements viscoplastiques. En vue d'applications réalistes telle que la simulation numérique des coulées de lave volcanique, le travail se concentre particulièrement sur les fluides complexes dont la rhéologie dépend fortement de grandeurs physiques telle que la température ou la concentration en particule. Nous développons un nouvel algorithme de résolution numérique des équations de Herschel-Bulkley combinant une méthode de Lagrangien augmenté à paramètre d'augmentation variable, une méthode des caractéristiques d'ordre 2 et une adaptation de maillage automatique. Sur des problèmes stationnaires ou en évolution tel que le problème test de la cavité entraînée, il apporte une solution efficace pour garantir à la fois une précision numérique élevée et un temps de calcul raisonnable. Cet algorithme est ensuite étendue et adapté au cas des rhéologies non-isothermes et aux suspensions. Concernant la simulation numérique des coulées de lave volcanique, nous détaillons une méthode de réduction par analyse asymptotique des équations de Herschel-Bulkley pour des écoulements de faible épaisseur sur une topographie arbitraire. Elle permet alors de décrire ces écoulements tridimensionnels de fluides viscoplastiques à surface libre par des équations bidimensionnelles surfaciques. Cette approche est ensuite étendue au cas non-isotherme en y ajoutant l'équation de la chaleur et des dépendances thermiques sur la rhéologie. Par intégration verticale de l'équation de la chaleur, on retrouve un modèle bidimensionnel. Le modèle non-isotherme est validé sur une expérience de dôme réalisée en laboratoire et une simulation numérique est réalisée autour d'une coulée qui a eu lieu sur le volcan du Piton de la Fournaise à la Réunion, en décembre 2010. La comparaison donne des résultats qui sont de notre point de vue satisfaisants et encourageants. / We present a contribution about modeling of viscoplastic flows. For realistic applications such as numerical simulation of volcanic lava flows, the work focuses particularly on complex fluids whose rheology strongly depends on physical quantities such as temperature or the particle concentration. We develop a new numerical resolution algorithm of Herschel-Bulkley's equations combining an augmented Lagrangian method with variable augmentation parameter, a second order characteristic method and an auto-adaptive mesh procedure. On stationary or evolving problems as the lid-driven cavity flow benchmark, it provides an effective solution to ensure both a high numerical accuracy within a reasonable computing time. This algorithm is then extended and adapted to the case of non-isothermal rheological and suspensions. On the numerical simulation of volcanic lava flows, we describe a method of reducing by asymptotic analysis of the Herschel-Bulkley's equations for thin flows on arbitrary topography. It allows to describe the three-dimensional flows of viscoplastic fluid with free surface by bidimensional surface equations. This approach is then extended to the non-isothermal case by adding the heat equation and thermal dependencies on rheology. By vertical integration of the heat equation, a two-dimensional model is maintained . The non-isothermal model is validated on a laboratory experiment of dome and a numerical simulation is performed on a December 2010 Piton de la Fournaise lava flow from La Réunion island. In our view, the comparison gives satisfactory and encouraging results.

Ce travail de thèse concerne les écoulements de fluides à seuil de contrainte autour d‘un obstacle cylindrique en milieu confiné avec une configuration d‘écoulement de Poiseuille.Expérimentalement, un banc d‘essai permettant d‘obtenir un écoulement en continu dans un circuit fermé a été mis en place. Les régimes d‘écoulement rampant, recirculant et instationnaire périodique ont été étudiés. De nouveaux résultats ont été obtenus avec un fluide newtonien et des solutions de Carbopol, polymère permettant de réaliser des fluides à seuil modèles utilisés en recherche et développement et dans l‘industrie. Une caméra rapide et un éclairage plan laser a servi pour l‘établissement d‘images qui ont ensuite été traitées par PIV. Les champs de vitesses cinématiques, les morphologies d‘écoulement et les paramètres critiques de transitions de régimes ont été déterminés.Numériquement, un modèle viscoplastique basé sur la loi de Herschel-Bulkley régularisée a été utilisé. Des résultats comme les morphologies d‘écoulement, la localisation des zones rigides, les champs de vitesses ont été obtenus. Ceci a permis de comparer les différences entre les effets liés à la nature des gels de Carbopol et la modélisation viscoplastiques. Une étude spécifique sur le glissement à l‘interface fluide-structure a également été réalisée avec l‘utilisation d‘un modèle de lubrification élasto-hydrodynamique. / The flow of yield stress fluids around a circular cylinder in a confined geometry has been investigated with a Poiseuille flow configuration.Experimentally, a test set-up was built which provides a continuous flow in a closed loop. We studied creeping, recirculating and vortex shedding flow regimes. New results has been realised with a Newtonian fluid and Carbopol solutions, models for yield stress behaviour in laboratory experiments and in industry. A high speed camera and a laser sheet have been used to perform images which are treated by PIV. Kinematic fields, flow morphologies and critical transition parameters have been determined.Numerically, a viscoplastic model based on the regularised Herschel-Bulkley law has been used. Results as flow morphologies, rigid areas and local flow parameters fields have been performed. That allowed us to compare the intrinsic effects of Carbopol solutions and the viscoplastic numerical model. A specific study on the wall slip has also been considered with an elasto-hydrodynamic lubrication model.