Ce travail de thèse vise à développer un procédé de synthèse d'esters dans un réacteur enzymatique membranaire fonctionnant en milieu CO2 supercritique. Un tel procédé représente une alternative intéressante à la synthèse chimique classiquement utilisée en industrie car il permet, d'une part, d'utiliser le CO2 supercritique comme solvant et de substituer ainsi les solvants organiques généralement utilisés et, d'autre part, d'avoir un produit final doté d'un label naturel grâce à l'utilisation d'un catalyseur biologique. Dans cette étude, une membrane enzymatique de taille industrielle a été développée. Un pilote spécifique permettant la conduite de réactions enzymatiques en milieu supercritique a été conçu. Le procédé a été testé à l'aide d'une réaction modèle : la synthèse d'acétate d'anisyle à partir d'alcool et d'acétate de vinyle. L'influence de divers paramètres opératoires tels que la température, la pression, ou le débit sur les performances du procédé a été évaluée. Les performances du procédé ont été également comparées à celles d'un réacteur à lit fixe. / This study deals with the development of an enzymatic membrane reactor working in supercritical carbon dioxide for ester synthesis. This process is an alternative to the chemical synthesis classically used in industry because it allows, on the one hand, the use of supercritical carbon dioxide as a solvent instead of organic solvents and on the other hand it leads to natural label ester thanks to the use of a biological catalyst. In this work, an industrial size enzymatic membrane was developed. A special pilot plant was designed to achieve enzymatic reactions in supercritical carbon dioxide medium. The process was studied with a model reaction : the anisyl acetate synthesis from anisic alcohol and vinyl acetate. The impact of several operating conditions like temperature, pressure and flow rate on the process performances was studied. The enzymatic membrane developed in this study was active and showed an interesting conversion rate. The performances were compared to those obtained with a packed bed reactor.

Ce travail de thèse de doctorat a eu pour objectif d’étudier la cristallisation d’un principe actif pharmaceutique, le sulfathiazole, par procédé Supercritique Anti-Solvant. Nous nous sommes plus particulièrement intéressés à l’influence des conditions opératoires sur les caractéristiques des cristaux générés en termes de taille, nature polymorphique et faciès. Lors d’une campagne expérimentale de cristallisation, nous avons fait varier la température (313 K et 328 K), la pression (10 et 20 MPa), le rapport molaire solvant/CO2 (2,5 à 15%) et la concentration massique de la solution organique (0,5 à 1,8%). La durée de l’étape de cristallisation a été aussi variée de 4, 5 et 8 h. Nous avons également étudié au préalable les équilibres de phases et mesuré les solubilités du principe actif dans la phase fluide afin de connaître l’état de sursaturation dans le milieu de cristallisation pour toutes les conditions étudiées. Les résultats ont montré que dans certaines conditions la forme polymorphique IV (la plus stable à pression ambiante) est obtenue pure alors que dans d’autres conditions elle est obtenue en mélange avec la forme instable I. Nous avons par ailleurs observé une variation des caractéristiques des cristaux en fonction de la durée de l’étape de cristallisation. Dans les conditions correspondant à de faibles sursaturations et à des phénomènes de transfert limités, la forme IV est initialement obtenue alors que la formation de la forme I est ensuite favorisée au cours du temps. Ces résultats démontrent clairement l’importance du contrôle de la durée de l’étape de cristallisation par le procédé SAS. / This Ph.D. work aimed to study the crystallization of an active pharmaceutical ingredient, sulfathiazole, by Supercritical Anti-Solvent process. We have particularly focused on the influence of the operating conditions on the characteristics of the formed crystals in terms of size, polymorphism concerns more than 80% of active pharmaceutical ingredients and the formation of an undesired polymorphic form can have dramatic consequences on the drug bioavailability or on its stability over time. During an experimental crystallization campaign, we varied the temperature (313 K and 328 K), the pressure (10 and 20 MPa), the solvent / CO2 molar ratio (2.5 to 15%) and the mass concentration of the organic solution (0.5 to 1.8%). The duration of the crystallization step was also varied to 4, 5 and 8 hours. We have also studied the phase equilibria before hand and measured the solubilities of the active principle in the fluid phase in order to evaluate the supersaturation level in the crystallization medium for all the studied conditions.The results showed that under certain conditions, the polymorphic form IV (the most stable at ambient pressure) is obtained pure whereas under other conditions, it is obtained in mixture with the unstable form I. We also observed a variation of the crystal characteristics depending on the duration of the crystallization step. For conditions corresponding to low supersaturations and limited transfer phenomena, form IV is initially obtained whereas the formation of form I is favored over time. These results clearly demonstrate the importance of controlling the duration of the SAS crystallization step.

L’objectif de cette thèse est d’appréhender la complexité des mécanismes biologiques de l’acupuncture afin de construire un modèle mathématique multi-échelle. Ce modèle est étudié théoriquement et numériquement. L’acupuncture est une des plus vielles pratiques de l’histoire de la médecine et une partie intégrante de la médecine traditionnelle chinoise. Dans sa pratique la plus classique, une ou plusieurs aiguilles sont placées à des endroits spécifiques, nommés points d’acupuncture. L’aiguille est ensuite manipulée en utilisant des mouvements de rotation et de translation de façon à stimuler le point d’acupuncture. Les effets cliniques de l’acupuncture pourraient être le résultat d’effet de cascades de réactions produites par les interactions entre l’hypoderme et les systèmes nerveux, endocrinien et immunitaire. Le travail présenté s’articule sur la modélisation de l’insertion d’une aiguille dans le tissu conjonctif de l’hypoderme. Un modèle d’écoulement en milieu poreux du liquide interstitiel de l’hypoderme a permis d’étudier numériquement les composantes de contrainte qui agissent sur les récepteurs à la surface des cellules du tissu et notamment des mastocytes.Un modèle mathématique de la réponse chimiotactique des mastocytes à une contrainte physique créée par le traitement d’acupuncture est développé. Ce modèle prend en compte les mécanismes de signalisation cellulaire. La contrainte physique induit la libération rapide et continue, grâce au recrutement chimotactique de mastocytes, d’attractants et de médiateurs chimiques. Le modèle est basé sur le modèle de chimiotaxie de type Keller-Segel. / The objective of this thesis is to comprehend the complexity of the underlying basis of acupuncture. Acupuncture needling is investigated in order to establish a multiscale model that takes into account the complexity of biology but is mathematically simple enough to run simulations.Acupuncture is one of the oldest practices in the history of medicine and is the core of Traditional Chinese Medicine. Once needles are inserted in the right locations, called acupoints, they are manipulated via manual needling to stimulate the acupoint. The physiological reactions of acupuncture needling lead to therapeutic effects which can be explained by a series of interactions between the skin and the nervous, the endocrine, and the immune systems.In the present work, the thrusting and lifting of an acupuncture needle inserted in subcutaneous connective tissue is modeled. A porous media model is used to run simulations and compute the pressure and shear stress affecting the organization of fibers and of isolated cells in their matrix. A mathematical model was conceived to take into account cell signaling. There is ample evidence that needle manipulation in acupuncture can cause degranulation of mastocytes directly through a physical stress to occur. Activated mastocytes rapidly release granules containing chemical mediators. These chemical mediators play a key role recruiting mastocytes in their environment and are known to affect the excitability of nerve endings as well as local microcirculation permeability and size for the appropriate transfer of long-term acting endocrine signals. The process is sustained by the recruitment of mastocytes through chemotaxis.

Actuellement, les mousses de polymères thermoplastiques (ou plastiques expansés) sont utilisées dans de nombreuses applications de la vie courante. On peut les retrouver dans l’isolation, l’absorption de choc, l’emballage ou encore dans des applications spécifiques du secteur biomédical. Un des inconvénients majeurs est l’utilisation de matières synthétiques. Pour pallier cette problématique, l’utilisation de polymères biosourcés et biodégradables représente une alternative intéressante. L’objectif de ce travail est d’élaborer des mousses de PLA (poly-(acide lactique)), un polymère biosourcé, biodégradable et biocompatible. Pour réaliser ces matériaux poreux, la technique d’extrusion couplée à l’injection de CO2 supercritique est une solution performante. En effet, selon les conditions opératoires et la solubilité du CO2 dans le biopolymère utilisé, son injection dans le fourreau d’une extrudeuse va modifier les propriétés rhéologiques du polymère. Il va, de plus, jouer le rôle d’agent moussant lors de la dépressurisation que subit le polymère au cours de son passage dans la filière. L’avantage de cette technique par rapport aux techniques classiques est, qu’ici, l’agent physique moussant (le CO2) n’est pas toxique, est acceptable d’un point de vue environnemental et ne laisse aucun résidu dans le produit. Ce procédé a montré un grand potentiel pour la création de matériaux poreux avec différentes morphologies. Des mousses avec une porosité supérieure à 95 % ont pu être obtenues. À ces très hautes porosités, deux types de structure ont été observés en fonction de la température : une mousse fortement expansée radialement avec des cellules ouvertes de grande taille, ou une mousse fortement expansée longitudinalement présentant une faible teneur en cellules ouvertes de petite taille. Une étude paramétrique (T, P, % CO2), a été menée avec le PLA seul et en mélange à de l’amidon sous différentes formes. Ce dernier permet d’améliorer la biodégradabilité du PLA et de modifier son comportement durant le moussage. Afin de parvenir à une meilleure compréhension et donc une meilleure maîtrise du procédé, un modèle a été développé. Il s’appuie sur le couplage entre l’écoulement du mélange dans la filière avec les phénomènes de nucléation et croissance des bulles dans ce mélange. Il a été validé grâce à la comparaison avec les résultats expérimentaux préalablement obtenus. / Nowadays, foams of thermoplastic polymers (or expanded plastics) are used in many applications of day life. They can be found in insulation, shock absorption, packaging or biomedical specific applications. One of the major drawbacks lies in the use of synthetic materials. To overcome this problem, the use of biobased and biodegradable polymers represents an interesting alternative.The aim of this work is to elaborate foams with the biopolymer poly(lactic acid), PLA. To manufacture such a porous material, the process of extrusion assisted by supercritical CO2 is an efficient solution. Indeed, depending on operating conditions and CO2 solubility in the biopolymer, its injection in the barrel of an extruder modifies the rheological behaviour of the polymer. It also plays the role of an expansion agent during the depressurization undergone by the polymer while flowing through the die. The main advantages in comparison with traditional methods, are that, here, the physical blowing agent (CO2) is not toxic, environmentally friendly while leaving no residue in the final product. This process shows great potential for the creation of porous materials with different morphologies. Foams with porosity as high as 95 % have been produced. For such a high porosity, two structures have been observed: foams with a large radial expansion with only open cells of large size, or foams with a large longitudinal expansion with low content of open cell having small size. A parametric study (T, P, CO2 %) was conducted with PLA both pure and mixed with starch in various forms. The latter makes it possible to improve the biodegradability of the PLA and to modify its behaviour during the foaming. Modelling of the process could allow a better understanding. A model based on the coupling between the polymer flow in the die with the nucleation and growth phenomena, has been developed. It has been validated by comparison with the experimental results previously obtained.

Cette thèse porte sur l’étude de la propagation d’ondes mécaniques dans des milieux granulaires secs ou mouillés, avec pour objectif de relier les phénomènes de l’échelle microscopique (dynamique des grains, potentiels d’interactions entre particules, rhéologie du fluide interstitiel) aux propriétés de l’échelle macroscopique (relation de dispersion, vitesse et atténuation des ondes dans l’approximation des grandes longueurs d’ondes). Les systèmes étudiés sont soit des milieux granulaires unidimensionnels de grande taille, analogues des chemins de plus forts contacts entre particules (les chaînes de force) dans les empilements de grains réels, soit les milieux granulaires réels eux-mêmes. Dans un premier temps, nous étudions expérimentalement la transmission d’ondes au travers d’un alignement de sphères centimétriques sèches, que nous modélisons via le potentiel de Hertz. Nous montrons que le couplage élastofrictionnel entre les grains et un substrat (le support des sphères) engendre un potentiel élastique local, qui induit à son tour une bande interdite a fréquence nulle dans la fonction de transfert. Dans un deuxième temps, nous montrons que la présence d'une quantité infime de fluide visqueux au contact entre chaque particule induit une interaction élasto-hydrodynamique (EHD). Ce dernier induit une modification de l’atténuation des ondes et une augmentation très significative de la vitesse de propagation, qui dans ce cas dépendent de manière non-triviale de l’élasticité des particules, de la viscosité du fluide et de la fréquence. Dans un troisième temps, nous vérifions la fiabilité de notre analyse pour décrire la propagation d'ondes ultrasonores dans des milieux granulaires réels, tel que le sable mouillé ou non ; les particules sont ici des sphères millimétriques. Dans le cas sec, nos résultats sont en accord avec un modèle connu de milieux effectifs (EMT) qui relève de l’interaction de Hertz-Mindlin dans l'approximation des grandes longueurs d'ondes. Dans le cas mouille, le modèle EMT combiné à un mécanisme EHD reproduit de manière acceptable nos observations préliminaires. / This thesis deals with the study of mechanical wave propagation in dry or wet granular media, with the aim of relating the phenomena at the microscopic scale (particles dynamics, interaction potentials between grains, rheology of the interstitial fluid) to the features at the macroscopic scale (dispersion relation, wave speed and attenuation in the long wavelength approximation). The systems under study are either large one-dimensional granular media, as the analogs of the paths of the most compressed grains (the force chains) in real granular packings, or the real granular media themselves. In a first place, we study experimentally the wave transmission through alignments of dry centimetric spheres, which we model via the Hertz potential. We show that the elasto-frictional coupling between the grains and a substrate (the spheres’ support) induces an on-site elastic potential, which in turn induces a band gap at zero frequency in the transfer function. In a second place, we show that the presence of an infinitesimal amount of viscous fluid at the contact between every particle induces an elasto-hydrodynamic (EHD) interaction. The later affects the attenuation of waves in addition to a significant increase of the wave speed, which in this case both non-trivially depend on the elasticity of the particles, on the viscosity of the fluid and on the frequency. In a third place, we check the reliability of our analysis to describe ultrasonic wave propagation in real granular materials such as dry or wet sand; our particles are here millimetric spheres. In the dry configuration, our results are consistent with an effective medium theory (EMT) which relies on the Hertz-Mindlin interaction in the long wavelength approximation. In the wet configuration, the EMT model combined with an EHD mechanism fairly reproduces our preliminary observations.

Dans l’industrie nucléaire, la simulation de transitoires accidentels à l’échelle d’un réacteur devient une composante d’importance croissante de la démonstration de sûreté à destination des agences de surveillance nationales. Elle permet ainsi de limiter le recours à des expériences complexes et coûteuses tout en facilitant l’évaluation des stratégies de mitigation. Cependant, les modèles mis en jeu sont inévitablement volumineux et construits avec une finesse de modélisation rendant difficile la prise en compte de détails géométriques locaux pourtant susceptibles d’influencer significativement la solution globale. Dans ce travail de thèse, on propose ainsi des approches multi-modèles pour l’intégration de tels détails dans un modèle global sans modification du maillage initial (on parle aussi de zoom numérique). Des techniques sont proposées aussi bien pour les structures que pour les fluides, avec un souci de démontrer la précision et la stabilité de la solution multi-modèles couplée comparée à une solution de référence à une seule échelle. Ce travail intègre deux spécificités propres, à savoir son adéquation avec les contraintes de la dynamique rapide avec intégration temporelle explicite et l’objectif de traiter simultanément la superposition de modèles et l’interaction fluide-structure. / In nuclear industry, simulating accidental transient sequences at full reactor scale is becoming an increasingly important feature of the safety demonstration towards national agencies. It thus allows limiting the number of complex and costly experiments, while simplifying and accelerating the evaluation of mitigation strategies. However, the implemented numerical models are inevitably heavy to build and maintain, with a global modelling scale making it difficult to account for local geometric details yet able to significantly influence the physical solution. To provide an answer to these problems, this PhD work is dedicated to multi-model approaches designed to integrate such details into bigger models with no modification at the global level (techniques often designated as numerical zoom techniques). Some methods are proposed for both structures and fluids, with special care given to the accuracy and stability of the coupled multi-scale solution compared to a single-scale reference solution. This work handles two very specific topics, namely its compatibility with numerical features imposed by fast transient dynamics with explicit time integration, and the general objective of simultaneously dealing with superimposed models and fluid-structure interaction.

Cette thèse concerne la modélisation des interactions fluide-structure et les méthodes numériques qui s’y rattachent. De ce fait, la thèse est divisée en deux parties. La première partie concerne l’étude des interactions fluide-structure par la méthode des domaines fictifs. Dans cette contribution, le fluide est incompressible et laminaire et la structure est considérée rigide, qu’elle soit immobile ou en mouvement. Les outils que nous avons développés comportent la mise en oeuvre d’un algorithme fiable de résolution qui intégrera les deux domaines (fluide et solide) dans une formulation mixte. L’algorithme est basé sur des techniques de raffinement local adaptatif des maillages utilisés permettant de mieux séparer les éléments du milieu fluide de ceux du solide que ce soit en 2D ou en 3D. La seconde partie est l’étude des interactions mécaniques entre une structure flexible et un fluide incompressible. Dans cette contribution, nous proposons et analysons des méthodes numériques partitionnées pour la simulation de phénomènes d’interaction fluide-structure (IFS). Nous avons adopté à cet effet, la méthode dite «arbitrary Lagrangian-Eulerian» (ALE). La résolution fluide est effectuée itérativement à l’aide d’un schéma de type projection et la structure est modélisée par des modèles hyper élastiques en grandes déformations. Nous avons développé de nouvelles méthodes de mouvement de maillages pour aboutir à de grandes déformations de la structure. Enfin, une stratégie de complexification du problème d’IFS a été définie. La modélisation de la turbulence et des écoulements à surfaces libres ont été introduites et couplées à la résolution des équations de Navier-Stokes. Différentes simulations numériques sont présentées pour illustrer l’efficacité et la robustesse de l’algorithme. Les résultats numériques présentés attestent de la validité et l’efficacité des méthodes numériques développées. / This thesis is concerned with the modeling of fluid-structure interactions (FSI) and the corresponding specific numerical methods. The thesis is divided into two principal parts. The first part concerns the study of fluid-structure interactions using the fictitious domain method. In this contribution, the fluid is incompressible and laminar and the structure is considered rigid, whether stationary or moving. The tools we have developed include the implementation of a reliable resolution algorithm that incorporates both domains (fluid and solid) in a common mixed formulation. The algorithm is based on adaptive local mesh refinement techniques used to distinguish the elements in the fluid from those of the solid either in 2D or 3D. The second part is the study of the mechanical interactions between a flexible structure and an incompressible fluid. In this context, we propose and analyze partitioned numerical methods for simulating fluid-structure interaction phenomena (FSI). We adopt an "arbitrary Lagrangian-Eulerian" (ALE) formulation for this purpose. The fluid resolution is performed iteratively by means of a projection scheme and the structure is modeled by hyperelastic models in large deformations. We have introduced new mesh update methods to achieve large deformation of the structure. Finally, a more complex strategy for FSI problem is proposed. The turbulence and two-phase flows modelling are introduced and coupled to the resolution of the Navier-Stokes equations for studying FSI problems. The numerical results presented attest the validity and effiency of the proposed numerical methods developed.

Cette thèse concerne l’étude des ailes oscillantes flexibles et des méthodes numériques qui s’y rattachent. De ce fait, la thèse est divisée en deux parties. La première contribution concerne le développement d’un algorithme de couplage fluide-structure qui prend en charge les interactions entre un solide élastique en grands déplacements et un fluide incompressible. L’algorithme est basé sur une approche partitionnée et permet d’utiliser des codes numériques de mécanique des fluides et de mécanique des solides existants. L’utilisation d’un terme de compressibilité artificiel dans l’équation de continuité du fluide combinée à des choix algorithmiques judicieux permet d’utiliser cette méthode de couplage efficacement avec un code de mécanique des fluides utilisant une méthode de projection de type SIMPLE ou PISO. La seconde contribution est l’étude de l’effet de flexibilité des ailes sur le vol à ailes battantes. Deux principaux régimes de vol sont mis en évidence concernant la déformation de l’aile : déformation causée par la pression et déformation causée par l’inertie. Les effets de ces régimes sur la topologie de l’écoulement et sur les performance de l’aile en propulsion sont discutés. Il est montré que les cas avec des déformations causées par la pression présentent généralement des efficacités plus élevées avec une flexibilité modérée. Il en est de même pour la force de poussée lorsque l’amplitude de tangage est faible. D’autre part, lorsque les déformations sont causées par l’inertie, les performances de l’aile sont généralement réduites. Certains cas montrent une augmentation marginale des performances lorsque le synchronisme des déformations est optimal, mais ces cas représentent davantage une exception que la norme. Il est également démontré que la flexibilité peut être utilisée comme mécanisme de tangage passif tout en conservant des performances intéressantes. Enfin, un modèle d’aile oscillante flexible non contraint est présenté. Il est démontré que le mouvement de déviation observé dans la nature est une conséquence d’un phénomène aérodynamique de mise en drapeau. / This thesis concerns the study of flexible flapping wings and the related numerical methods. It thus contains two distinct themes. The first contribution is the implementation of an efficient fluid-structure interaction algorithm that handles the interaction of an elastic solid undergoing large displacement with an incompressible fluid. The algorithm is based on the partitioned approach and allows state-of-the-art fluid and structural solvers to be used. Stabilization with artificial compressibility in the fluid continuity equation along with judicious algorithmic choices make the method suitable to be used with SIMPLE or PISO projection fluid solvers. The second contribution is the study of the effects of wing flexibility in flapping flight. The different regimes, namely inertia-driven and pressure-driven wing deformations are presented along with their effects on the topology of the flow and, eventually, on the performance of the flapping wing in propulsion regime. It is found that pressure-driven deformations can increase the thrust efficiency if a suitable amount of flexibility is used. Thrust increases are also observed when small pitching amplitude cases are considered. On the other hand, inertia-driven deformations generally deteriorate aerodynamic performances of flapping wings unless meticulous timing is respected, making them less practical. It is also shown that wing flexibility can act as a passive pitching mechanism while keeping decent thrust and efficiency. Lastly, a freely-moving flexible flapping wing model is presented. It is shown that the deviation motion found in natural flyers is a consequence of a feathering mechanism.

La théorie de couplage de mode (MCT) est l'une des méthodes les plus utilisées pour étudier les transitions vitreuses dans les fluides classiques. Ses prédictions sont en général en accord semi quantitatif avec les simulations. Sa mise en oeuvre nécessite la détermination de la structure statique, généralement par résolution des équations d'Orsntein-Zernike avec une fermeture adéquate. Partant de fermetures utilisant des fonctions « bridges » déduites de la fonctionnelle de référence du mélange de sphères dures, notre travail a consisté d'abord à étudier l'influence de la qualité de cette structure statique sur les prédictions relatives aux états non ergodiques dans des mélanges binaires dissymétriques. Nous avons ensuite considéré les résultats de la théorie du couplage de modes dans sa version naïf (NMCT ) et complète, afin d'analyser les mécanismes d'arrêt, les comparer au fluide effectif et aux approches stochastiques (équations de Langevin généralisées). Enfin, nous proposons une version pragmatique de cette méthode qui fournit des prédictions en meilleur accord quantitatif avec les résultats des simulations pour une variété de potentiels d'interaction / The mode coupling theory (MCT) is one of the most widely used methods for studying the glass transition in classical fluids. Its predictions are usually in semi-quantitative agreement with simulation. Its implementation requires the determination of the static structure usually from the Ornstein-Zernike equations with a suitable closure. Starting from closures that use bridge functions deduced from the hard-sphere reference functional, our work consisted first in studying the influence of the quality of this static structure on the predictions concerning the non-ergodic states in asymmetric binary mixtures. We next considered the results of the mode coupling theory in its naive and full versions, in order to analyze the arrest mechanisms and compare them to the effective fluid and the stochastic approaches (generalized Langevin equations). Finally, we propose a pragmatic version of this method that provides predictions in better quantitative agreement with simulations for a variety of interaction potentials

L’objectif de cette thèse est de développer un modèle numérique 3D afin d’étudier le phénomène de laves torrentielles ; écoulements visqueux fortement chargés en matière solide, surgissant en montagne lors d’orages violents. Aujourd’hui, la prévision des zones de vulnérabilité s’appuie sur des outils de calcul 0D, 1D ou 2D. Or ces outils ne peuvent représenter pleinement le comportement à surface libre des écoulements du fait de nombreuses approximations et hypothèses. Ainsi cette thèse met en œuvre un code numérique 3D pour étudier ces écoulements. Ce travail aboutit à la construction d’un modèle 3D à partir de données réelles de terrain. Plusieurs scénarios ont été étudiés et comparés à des résultats issus d’un modèle 2D. Les résultats mettent en évidence les apports non négligeables de la modélisation 3D : zones d’étalement et de dépôt, phénomènes d’encombrement, modélisation fine des écoulements dans les zones chenalisées. Préalablement, le modèle 3D a été validé en comparant les résultats numériques à des données expérimentales issues de la littérature, pour des typologies d’écoulement représentatives de celles observées sur des sites grandeur nature. / The objective of this thesis is to develop a 3D numerical model to assess debris flow. These viscous flows, heavily loaded with solid matter, form when heavy rain occurs in mountains. Today, forecasts of potentially impacted areas are based on 0D, 1D and 2D numerical tools. However, these tools cannot fully represent the free surface behaviour of debris flows due to the approximations and assumptions on which they are based. Thus, this work utilises a 3D numerical code to study this phenomenon. A specific model is built with real field data. Several flow scenarios are studied and compared with a 2D numerical model. The results highlight the significant benefits of a 3D approach by providing information on the fine representation of flow dynamics over the catchment area. The model also predicts the impact of debris flow (overflowing on a road bridge) and the zones of deposition and spreading. It highlights possible congestion phenomena and reproduces flows in the channels by fully accounting for parietal friction, capabilities not provided by 2D models. Prior to this application, the 3D model was evaluated with five sets of experimental data to validate its ability to represent viscoplastic flows. Different types of flows are studied and are representative of those observed on real sites when debris flow occur.