Dans le cadre général de la simulation du procédé de frittage en phase solide, ce travail propose une approche numérique, à l'échelle des particules, de la consolidation d'un compact pulvérulent céramique. Le frittage est un procédé mettant en jeu plusieurs chemins de diffusion activés thermiquement. Parmi ces chemins de diffusion, les plus importants sont : la diffusion surfacique, la diffusion aux joints des grains et la diffusion volumique. La mise en place de cette physique dans un cadre de calcul intensif doit permettre de mieux comprendre ces mécanismes de diffusion ainsi que leur influence sur l'évolution de la microstructure. Le but de ce travail consiste à développer un modèle ainsi qu'une stratégie numérique capable d'intégrer les différents mécanismes de diffusion dans un cadre de calcul intensif. Le flux de matière est calculé en fonction du Laplacien de la courbure dans les cas de la diffusion surfacique, tandis que pour la diffusion volumique ce flux est proportionnel au gradient de la pression hydrostatique. Le modèle physique est tout d'abord présenté dans le cadre de la mécanique des milieux continus. Ensuite, la stratégie numérique développée pour la simulation du frittage d'un empilement granulaire est détaillée. Cette stratégie est basée sur une discrétisation du problème par des éléments finis stabilisés couplée avec une méthode Level-set pour décrire la surface libre des particules. Cette stratégie nous permet de faire des simulations avec un "grand" nombre de particules. Plusieurs simulations en 3D, menées dans un cadre de calcul parallèle, montrent l'évolution qui a lieu sur un empilement granulaire réaliste. / Within the general context of solid-state sintering process, this work presents a numerical modeling approach, at the particle scale, of ceramic particle packing consolidation. Typically, the sintering process triggers several mass transport paths that are thermally activated. Among those diffusion paths, the most important ones are: surface diffusion, grain boundary diffusion and volume diffusion. Including this physics into a high-performance computing framework would permit to gain precious insights about the driving mechanisms. The aim of the present work is to develop a model and a numerical strategy able to integrate the different diffusion mechanisms into continuum mechanics framework. In the cases of surface diffusion and volume diffusion, the mass flux is calculated as a function of the surface curvature Laplacian and the hydrostatic pressure gradient, respectively. The physical model describing these two transport mechanisms is first presented within the framework of continuum mechanics. Then the numerical strategy developed for the simulation of the sintering of many particles is detailed. This strategy is based on a discretization of the problem by using a finite element approach coupled with a Level-Set method used to describe the particles free surface. This versatile strategy allows us to perform simulations involving a relatively large number of particles. Furthermore, a mesh adaptation technique allows the particles surface description to be improved, while the number of mesh elements is kept reasonable. Several 3D simulations, performed in a parallel computing framework, show the changes occurring in the structure of 3D granular stacks.

Le but de ce mémoire de recherche est d'étudier les convections naturelle et mixte d'écoulements fluides, et de développer et valider des méthodes numériques pour le suivi d'interfaces afin de traiter plus tard des écoulements incompressibles de fluides immiscibles. Dans une première étape, une méthode numérique originale, basée sur des discrétisations Volumes Finis, est développée pour modéliser les écoulements à faible nombre de Mach et grands écarts de température. Trois applications physiques, portant sur l'écoulement d'air à travers des plaques verticales parallèles chauffées, sont étudiées. Nous avons montré que l'espacement optimal, correspondant au pic de flux de chaleur transféré d'un réseau de plaques parallèles isothermes refroidies par convection mixte, est plus faible que ceux obtenus en convections naturelle ou forcée lorsque la chute de pression à la sortie est constante. Nous avons également prouvé que les écoulements de convection mixte à débit imposé peuvent présenter des solutions physiques inattendues ; un modèle alternatif basé sur une pression totale imposée à l'entrée et une pression fixée à la sortie donne de meilleurs résultats. Pour des canaux soumis un flux de chaleur sur une paroi seule, le rayonnement de surface tend à supprimer l'apparition des recirculations à la sortie et à uniformiser les températures des parois. Dans une seconde étape, le modèle mathématique couplant les équations de Navier-Stokes incompressibles et la méthode Level-Set pour le suivi d'interfaces est développé. Des améliorations de la conservation du volume fluide par l'utilisation de schémas de discrétisation d'ordres élevés (ENO-WENO) pour l'équation de transport et des variantes de l'équation de la distance signée sont discutées / The aim of this research dissertation is at studying natural and mixed convections of fluid flows, and to develop and validate numerical schemes for interface tracking in order to treat incompressible and immiscible fluid flows, later. In a first step, an original numerical method, based on Finite Volume discretizations, is developed for modeling low Mach number flows with large temperature gaps. Three physical applications on air flowing through vertical heated parallel plates were investigated. We showed that the optimum spacing corresponding to the peak heat flux transferred from an array of isothermal parallel plates cooled by mixed convection is smaller than those for natural or forced convections when the pressure drop at the outlet keeps constant. We also proved that mixed convection flows resulting from an imposed flow rate may exhibit unexpected physical solutions; alternative model based on prescribed total pressure at inlet and fixed pressure at outlet sections gives more realistic results. For channels heated by heat flux on one wall only, surface radiation tends to suppress the onset of recirculations at the outlet and to unify the walls temperature. In a second step, the mathematical model coupling the incompressible Navier-Stokes equations and the Level-Set method for interface tracking is derived. Improvements in fluid volume conservation by using high order discretization (ENO-WENO) schemes for the transport equation and variants of the signed distance equation are discussed

Les alliages de titane α/β ont beaucoup d’applications dans des domaines industriels divers comme l’aéronautique. Le phénomène de globularisation qui se produit lors de traitements thermomécaniques est un phénomène important dans la mesure où une microstructure globulaire présente une tenue mécanique et une ductilité accrues.Les microstructures lamellaires sont constituées de colonies de lamelles de phase α parallèles qui se développent dans les grains β au cours de leur refroidissement. La globularisation se fait en deux étapes : les lamelles se subdivisent d’abord en segments plus courts, puis ces segments globularisent, au cours de la déformation à chaud et des traitements thermiques. La subdivision des lamelles se fait au niveau des sous-joints formés au cours de la déformation à chaud. Lors de traitements thermiques prolongés, les plus gros globules de phase α grossissent au détriment des plus petits. La formation des sous-joints et le grossissement des globules n’ont pas été étudiés en détail dans ce travail. L’accent a été mis sur les mécanismes de migration des interfaces α/β et α/α qui conduisent à la subdivision des lamelles et à leur globularisation.Des essais de compression à chaud et des traitements thermiques appliqués à des échantillons de Ti-6Al-4V ont permis de confirmer que l’épaisseur des lamelles et de leur orientation sont des facteurs importants pour la globularisation. Les lamelles plus fines et celles orientées parallèlement à l’axe de compression globularisent plus facilement. Ce travail expérimental a permis de mieux appréhender la complexité du phénomène de globularisation et d’introduire un cadre numérique adapté pour sa simulation. Une méthode à champ complet, basée sur la méthode Level-Set dans un cadre Eléments Finis, a ainsi été testée pour simuler les mécanismes physiques de migration interfaciale menant à la subdivision des lamelles et à l'évolution vers une forme globulaire. Les premiers résultats sont très prometteurs et illustrent le potentiel du cadre numérique proposé. / Α/β titanium alloys have many industrial applications in various fields such as aeronautics. Spheroidization is a phenomenon that occurs in initially lamellar α/β titanium alloys during thermomechanical processing and receives considerable attention as spheroidized microstructures exhibit enhanced strength and ductility.Lamellar microstructures are made of colonies of parallel α lamellae developed inside β grains while they are cooled down. Spheroidization actually proceeds in two successive steps: the lamellae first split into smaller α laths, which subsequently undergo spheroidization. This occurs during hot-deformation and subsequent annealing. Lamella splitting occurs where subboundaries were formed inside lamellae during hot-deformation. Over long term annealing the spheroidized α phase particles undergo coarsening. The formation of subboundaries and coarsening were not addressed in this work. The focus has been placed on the interfacial kinetics mechanisms leading α lamellae splitting during the first stages of spheroidization.Hot compression tests and subsequent annealings carried out on Ti-6Al-4V samples confirmed that the thickness and the orientation of the lamellae are important factors with regards to spheroidization. Thinner lamellae and lamellae oriented parallel to the compression axis spheroidize faster. Those experiments contributed to a better understanding of the phenomenon and allowed to introduce a suitable numerical framework to simulate the early stages of spheroidization. A full field method in a Finite Element/ Level Set framework has thus been tested for simulating the involved physical mechanisms of interface migration that lead to lamellae splitting and the subsequent shape evolution of the α laths towards a spheroidal shape. First results are promising and illustrate the potential of this numerical framework

Ce travail porte sur la simulation numérique et la modélisation du comportement thermo-mécanique des matériaux composites renforcés par des fibres. Spécifiquement les matériaux SMC (Sheet Moulding Compound) sont utilisés dans le processus de moulage par compression pour construire des pièces automobiles de haute performance. Ce travail est divisé en quatre chapitres, décrivant tout d’abord un modèle thermo-mécanique entièrement couplé pour les matériaux SMC standards et innovants à haute concentration en fibres (> 25% en volume). Le SMC est traité comme un mélange incompressible de fibre et de résine complété éventuellement par une phase de porosité compressible. Son anisotropie est modélisée au moyen de tenseurs structurels. La cinétique de réaction et de consolidation de la pièce est également modélisée et étudiée. Les données expérimentales mécaniques et thermiques enregistrées sur des échantillons de matériaux SMC sont comparées au modèle et à la solution numérique fournie par ce travail. D’un point de vue numérique, nous utilisons la méthode des domaines immergées o`u chaque phase est distinguée par une fonction distance signée. Nous décrivons le procédé de moulage par compression en proposant une résolution compressible anisotrope unifiée capable de décrire la transition compressible / incompressible du matériau SMC sous déformation. Cela permet de décrire la réponse mécanique du SMC et de prédire localement la consolidation (durcissement) de la pièce le long du cycle thermique. / This work deals with the numerical simulation and modeling of thermomechanical analysis of fiber reinforcedcomposites materials. Specifically for SMC (Sheet Molding Compound) materials that are used in compression molding processes to build automotive high performance parts. The work is divided into fourchapters, firstly describing a fully coupled thermo-mechanical model for standard SMC materials and for innovative SMC with high fiber concentration (> 25% in volume). The SMC is treated as an incompressible mixtureof fibers and paste complemented by a compressible porosity phase. Its anisotropy is modeled by means of structural tensors. Kinetic of reaction and consolidation of the part is also modeled and studied. Mechanicaland thermal experimental data recorded on samples of SMC materials are compared to the model and numerical solution provided in this work. A numerical framework, we use the immersed boundary method and the level set method. We describe the compression molding process by proposing an unified anisotropic compressible resolution able to describe the transition between compressible/ incompressible of SMC materials under deformation. We are able to describe the mechanical response of the SMC and to predict locally the consolidation (curing) of thepart throughout the thermal cycle.

L’objectif de cette étude est de simuler l’érosion d’un sol cohésif sous l’effet d’un écoulement incompressible. Le modèle élaboré décrit une vitesse d’érosion interfaciale qui dépend de la contrainte de cisaillement de l’écoulement. La modélisation numérique proposée est une approche eulérienne, où une méthode de pénalisation de domaines est utilisée pour résoudre les équations de Navier-Stokes autour d’un obstacle. L’interface eau/sol est décrite par une fonction Level Set couplée à une loi d’érosion à seuil.L’approximation numérique est basée sur un schéma DDFV (Discrete Duality Finite Volume) autorisant des raffinements locaux sur maillages non-conformes et non-structurés. L’approche par pénalisation a mis en évidence une couche limite d'inconsistance à l'interface fluide/solide lors du calcul de la contrainte de cisaillement. Deux approches sont proposées pour estimer précisément la contrainte de ce problème à frontière libre. La pertinence du modèle à prédire l’érosion interfaciale du sol est confirmée par la présentation de plusieurs résultats de simulation, qui offrent une meilleure évaluation et compréhension des phénomènes d'érosion / This study focuses on the numerical modelling of the interfacial erosion occurring at a cohesive soil undergoing an incompressible flow process. The model assumes that the erosion velocity is driven by a fluid shear stress at the water/soil interface. The numerical modelling is based on the eulerian approach: a penalization procedure is used to compute Navier-Stokes equations around soil obstacle, with a fictitious domain method, in order to avoid body- fitted unstructured meshes. The water/soil interface’s evolution is described by a Level Set function coupled to a threshold erosion law.Because we use adaptive mesh refinement, we develop a Discrete Duality Finite Volume scheme (DDFV), which allows non-conforming and non-structured meshes. The penalization method, used to take into account a free velocity in the soil with non-body-fitted mesh, introduces an inaccurate shear stress at the interface. We propose two approaches to compute accurately the erosion velocity of this free boundary problem. The ability of the model to predict the interfacial erosion of soils is confirmed by presenting several simulations that provide better evaluation and comprehension of erosion phenomena.

Ce mémoire présente les travaux réalisés en aérodynamique afin de pouvoir disposer d’une chaîne de conception complète nécessaire au traitement des différentes problématiques sur les pales éoliennes, qui ne peuvent être pris en compte à l’aide d’une méthode unique. Afin de pouvoir faire de l’optimisation de pales en fonction de différents critères, un outil numérique simplifié (VALDAG) a été développé. Le module de simulation utilise la méthode du Disque Actif Généralisé, qui s’appuie sur la résolution des équations Navier-Stokes, complété par des corrections empiriques. Il respecte un compromis entre la précision et le coût de calcul. Cet outil pour lequel une interface web a été développée pour le rendre adapté à des travaux d’ingénierie est susceptible de se calibrer automatiquement sur une nouvelle géométrie de pale grâce à la méthode Nelder-Mead Simplex. On peut ensuite réaliser une optimisation de la performance de la pale en modifiant cette géométrie, et générer les fichiers 3D qui serviront d’entrée à la simulation 3D et de visualisation des optimisations de forme. Ces designs optimisés sont ensuite validés par des simulations DNS à l’aide de l’outil NaSCar 3D développé à l’INRIA. Cet outil résout des équations Navier-Stokes sur un maillage cartésien 3D et prend en compte des obstacles immergés via la fonction Level-Set et la pénalisation. Après des adaptations sur le traitement de la géométrie de pales, un compromis de configuration CFD est trouvé pour simuler un rotor éolien. En conclusion, ce paquet d’outils VALDAG est peu coûteux, facile à utiliser et efficace. En associant VALDAG aux simulations 3D, une chaîne de conception est complétée. / This Ph.D. thesis presents some research work on aerodynamics of wind turbine blades, in order to dispose a conception chain necessary for different problems, which cannot be treated by one unique method. A simplified numerical toolkit (VALDAG) has been developed to optimize the performance of blades in different creteria. The simulation module use the Generalized Actuator Disc model, which relies on the solution of Navier-Stokes equations and completed with empiric corrections. This tool respects a reasonable compromise between model complexity and computational reliability. An automatic calibration mechanism was implemented using the Nelder-Mead Simplex algorithm. A web users interface (WUI) is also available to adapt VALDAG for industrial engineers. Optimization is then carried on by modifying the blades’ geometry parameters and the designs optimized is stocked in files which can be used for 3D simulation and/or visualization. The blade designed with VALDAG are then simulated by a 3D numerical tool to validate previous predictions. This 3D tool called NaSCar 3D is developed in INRIA and resolve Navier-Stokes equations on to a cartesian mesh, in which the immersed obstacles are considered with the Level-Set function and the penalization method. After some necessary adaptation for the treatment of blades’ geometry, a compromise on CFD configuration is found to simulate a wind turbine rotor. To conclude, this design toolkit VALDAG is low time-costly, user-friendly and efficient. Associated with 3D simulations, a conception chain has been established.

Background Camera-based photoplethysmography (cbPPG) is a measurement technique which enables remote vital sign monitoring by using cameras. To obtain valid plethysmograms, proper regions of interest (ROIs) have to be selected in the video data. Most automated selection methods rely on specific spatial or temporal features limiting a broader application. In this work, we present a new method which overcomes those drawbacks and, therefore, allows cbPPG to be applied in an intraoperative environment. Methods We recorded 41 patients during surgery using an RGB and a near-infrared (NIR) camera. A Bayesian skin classifier was employed to detect suitable regions, and a level set segmentation approach to define and track ROIs based on spatial homogeneity. Results The results show stable and homogeneously illuminated ROIs. We further evaluated their quality with regards to extracted cbPPG signals. The green channel provided the best results where heart rates could be correctly estimated in 95.6% of cases. The NIR channel yielded the highest contribution in compensating false estimations. Conclusions The proposed method proved that cbPPG is applicable in intraoperative environments. It can be easily transferred to other settings regardless of which body site is considered.

Cette thèse est consacrée au développement et à la comparaison des méthodes de suivi d'interface pour les écoulements diphasiques incompressibles. Elle s'intéresse à la sélection de méthodes robustes de suivi d'interface, puis à leur couplage avec le solveur des équations de Navier-Stokes. La méthode level-set est en premier lieu étudiée, en particulier l'influence du schéma d'advection et de l'étape de réinitialisation sur la qualité des résultats du suivi d'interface. Il a été montré que la méthode de réinitialisation avec contrainte de volume est robuste et précise en combinaison avec des schémas conservatifs WENO d'ordre 5 pour l'advection. Il a été constaté que les erreurs du suivi d'interface augmentent de manière abrupte lorsque la condition CFL est trop petite. Comme remède, la réinitialisation du champ level-set effectuée moins souvent réduit la diffusion numérique et le déplacement non-physique de l'interface. La conservation de la masse n'est pas assurée avec les méthodes level-set. Les méthodes VOF (volume-of-fluid) qui conservent naturellement la masse du fluide de référence sont alors étudiées. Une résolution géométrique avec un schéma consistent et conservatif est alors adoptée, ainsi qu'une autre technique alternative plus aisément extensible en 3D. Il a été trouvé que ces deux dernières méthodes donnent des résultats très proches. La méthode MOF (moment-of-fluid), qui reconstruit l'interface en utilisant le centre de masse du fluide de référence, est plus précise que les méthodes VOF. Différentes méthodes couplées entre level-set et VOF sont alors étudiées, notamment: CLSVOF, MCLS, VOSET et CLSMOF. Il a été observé que la méthode level-set tend à épaissir les filaments minces, tandis que VOF et les méthodes couplées les fragmentent en petites particules. Finalement, on a couplé les méthodes level-set et VOF avec le solveur incompressible des équations de Navier-Stokes. On a comparé différentes manières de prise en compte des conditions de saut à l'interface (lisse et raide). Il a été montré que les méthodes VOF sont plus robustes, et donnent d'excellents résultats pour quasiment toutes les simulations. Deux méthodes level-set donnant de très bons résultats, comparables à ceux de VOF, sont aussi identifiées. / This thesis is devoted to the development and comparison of interface methods for incompressible two-phase flows. It focuses on the selection of robust interface capturing methods, then on the manner of their coupling with the Navier-stokes solver. The level-set method is first investigated, in particular the influence of the advection scheme and the reinitialization step on the accuracy of the interface capturing. It is shown that the volume constraint method for reinitialization is robust and accurate in combination with the conservative fifth-order WENO schemes for the advection. It is found that interface errors increase drastically when the CFL number is very small. As a remedy, reinitializing the level-set field less often reduces the amount of numerical diffusion and non-physical interface displacement. Mass conservation is, however, not guaranteed with the level-set methods. The volume-of-fluid (VOF) method is then investigated, which naturally conserves the mass of the reference fluid. A geometrical consistent and conservative scheme is adopted, then an alternative technique more easily extended to 3D. It is found that both methods give very similar results. The moment-of-fluid (MOF) method, which reconstructs the interface using the reference fluid centroid, is found to be more accurate than the VOF methods. Different coupled level-set and VOF methods are then investigated, namely: CLSVOF, MCLS, VOSET and CLSMOF. It is observed that the level-set method tends to thicken thin filaments, whereas the VOF and coupled methods break up thin structures in small fluid particles. Finally, we coupled the level-set and volume-of-fluid methods with the incompressible Navier-Stokes solver. We compared different manners (sharp and smoothed) of treating the interface jump conditions. It is shown that the VOF methods are more robust, and provide excellent results for almost all the performed simulations. Two level-set methods are also identified that give very good results, comparable to those obtained with the VOF methods.

Cette thèse porte sur l’étude des écoulements de Taylor (ou poche/bouchon) gazliquide en microcanal. Ces écoulements où les effets de tension de surface sont prépondérants ont été étudiées expérimentalement et numériquement pour des géométries rectangulaires avec divers rapports d’aspects. Une première partie expérimentale a consisté à caractériser la formation de bulles (taille, fréquence) en fonction des conditions opératoires, des propriétés des fluides (notamment à travers le nombre capillaire) et du mode de mise en contact des fluides. La dynamique de l’écoulement établi a par la suite été étudiée à l’aide du code JADIM. La simulation de ces écoulements dominés par la tension de surface a nécessité de lever les limitations liées à la prise en compte de la force capillaire. En effet des courants parasites numériques sont créés à proximité de l’interface lors de la simulation d’écoulements capillaires. Une méthode Level Set a été implémentée et comparée à la méthode Volume of Fluid d’origine en termes de courants parasites. Des simulations numériques 3D ont permis l’étude des effets du nombre capillaire et de la géométrie sur la dynamique des bulles de Taylor (vitesse, pression et formes de bulles). Les effets inertiels souvent négligés ont été considérés et leur influence, notamment sur les sauts de pression à l’interface, a été mise en évidence. Le mélange dans le bouchon liquide a également été étudié. / This thesis focuses on the hydrodynamics of gas-liquid Taylor flow (or slug flow) in microchannels. These flows, which are generally dominated by surface tension forces, have been investigated in rectangular channels of various cross-sectional aspect ratios by means of both experimental visualizations and numerical simulations. The first experimental part aims at characterizing the bubble generation process (bubble length and frequency of break-up) depending on the operating conditions, the fluid properties, as well as the junction where both fluids merge. Numerical simulations of fully developed Taylor flow have been carried out with the JADIM code. The computation of such surface tension dominated flows requires an accurate calculation of the surface tension force. Some limitations of the Volume of Fluid method have been highlighted and a Level Set method has been developed in order to improve the calculation of capillary effects. Both methods have been compared in detail in terms of spurious currents. 3D numerical simulations have been performed and the influence of the capillary number, as well as the effects of geometry have been highlighted. Inertial effects have been taken into account and their influence on the pressure drop has been shown to be non-negligible. Mixing in the liquid slug has also been studied.

Ce travail porte sur la mise en forme des matériaux composites à matrice polymère renforcée par des fibres courtes dont les performances physiques et mécaniques sont directement reliées à la distribution spatiale et à l’orientation des renforts employés. Il se focalise sur l’étude des mécanismes de déformation se produisant au cours de l’écoulement de ces systèmes qui se comportent comme des suspensions de fibres non-newtoniennes. Le problème est abordé par une procédure originale combinant images 3D acquises en temps réel et simulations numériques avancées, réalisées à l’échelle des fibres. Dans le premier cas, des suspensions modèles avec fluide suspensif non-newtonien ont été déformées en compression dans des conditions confinées dans un microtomographe à rayons X synchrotron. Cette technique a permis l’acquisition en temps réel de clichés 3D à forte résolution spatiale de l’écoulement des suspensions. Dans le deuxième cas, un code de calculs éléments finis 3D a été utilisé, celui-ci étant capable de décrire finement des objets immergés dans des fluides non-newtoniens, par des level-sets et des techniques de remaillage anisotrope. La pertinence des simulations numériques dans les régimes de concentration dilués à semi-dilués a été jaugée par une comparaison expériences-simulations avancée.De là, dans le régime de concentration dilué, nous montrons que le confinement de l’écoulement et le comportement rhéofluidifiant du fluide suspensif ont une influence mineure sur la cinématique des fibres, si ces dernières sont suffisamment éloignées des plateaux de compression. Si ce prérequis n’est pas respecté, l’effet du confinement devient important. Des modifications au modèle heuristique d’haltère de la littérature ont été proposées pour corriger la cinématique de fibres. Dans le régime semi-dilué, des déviations de la cinématique de fibres sont également observées au cœur des suspensions. Ces déviations sont principalement liées aux interactions hydrodynamiques entre fibres suffisamment voisines. La cinématique des fibres prédite par le modèle de Jeffery et les approximations de champ affine sont mises en défaut. Dans le régime concentré, si l’évolution de l’orientation globale de la suspension est étonnamment bien décrite par l’équation de Jeffery, de très importantes fluctuations des champs de translation et de rotation des fibres sont observées à l’échelle des fibres. Celles-ci sont induites par les nombreux contacts entre fibres qui peuvent par ailleurs être correctement prédits par le modèle de tube. / This study focuses on the processing of short fibre-reinforced polymer composites. The physical and mechanical properties of these materials are mainly affected by the position and orientation distribution of fibres induced during their forming. Thus, we analysed the flow-induced micro-mechanisms that arose at the fibre scale during the forming stage of these complex systems which behave as non-Newtonian fibre suspensions. For that purpose, an original approach was developed by combining 3D imaging technique and direct numerical simulation, both performed at the fibre scale. Hence, several model fibre suspensions with a non-Newtonian suspending fluid and with a concentration regime that ranged from dilute to concentrated were prepared . They were subjected to confined lubricated compression loadings using a rheometer mounted on a synchrotron X-ray microtomograph. Thanks to very short scanning times, 3D images of the evolving fibrous microstructures at high spatial resolution were recorded in real-time. These experiments were also simulated using a dedicated Finite Element library enabling an accurate description of fibre kinematics in complex suspending fluids thanks to high performance computation, level sets and adaptive anisotropic meshing. The efficiency of the numerical simulation from the dilute to semi-dilute concentration regimes was assessed through experimental and numerical comparisons.Then, we showed that the confinement effect and the non-Newtonian rheology of the suspending fluid had a weak effect on the fibre kinematics, if the fibres were sufficiently far from the compression platens, typically the fibre-platen distance should be larger than twice the fibre diameter. Otherwise, confinement effects occurred. Some extensions of the dumbbell model were proposed to correct the fibre kinematics in this flow conditions. In semi-dilute concentration, deviations of the fibre kinematics compared to the Jeffery’s predictions were also observed and related to hydrodynamic interactions between fibres. In this case, the predictions of Jeffery’s model and the related assumption of affine fibre motions are less relevant. In the concentrated regime, even if the overall orientation of fibre suspension could be astonishingly well described by using the Jeffery’s model, strong fluctuations on each fibre motion and rotation were observed. These deviations were induced by the numerous fibre-fibre contacts, which could be correctly predicted by the tube model.