Cette thèse consiste à mettre en place un couplage externe aéro-thermo-mécanique, sur la base d'un schéma partitionné, entre les codes de recherche CEDRE (mécanique des fluides, volumes finis) et Z-set (modules indépendants pour la mécanique des structures et la thermique du solide, éléments finis). Les résultats numériques sont confrontés à ceux de l'expérience (une campagne de mesures a été menée dans le cadre de cette étude), dans le cas d'un problème complexe lié au domaine de l'aérospatial : l'interaction flamme-paroi. Ce phénomène est piloté par la thermique, à travers le flux de chaleur pariétal généré par la flamme. A cause de la disparité des temps caractéristiques thermiques entre les milieux fluide et solide, la partie aéro-thermique du couplage est traitée de façon simplifiée, en considérant le fluide comme une suite d'états stationnaires. L'échauffement de la plaque métallique provoque sa déformation (la loi de comportement mécanique du matériau est de type élasto-visco-plastique). Le déplacement de l'interface fluide-structure est propagé sur le maillage fluide. En se basant sur les similitudes entre jets non réactifs et réactifs (de type flamme) dans le cas de l'impact, des calculs couplés sont menés dans des configurations 2D et 3D de l'impact d'un jet chaud non réactif. / This thesis consists in setting up an external fluid-thermal-structural coupling, based on a partitionned scheme, between the research codes CEDRE (fluid mechanics, finite volumes) and Z-set (independent solvers for structural mechanics and heat transfer through the solid). The numerical results are compared with experimental data, to study a complex problem related to the aerospace certification process: the flame-wall interaction. This phenomenon is is driven by the heat flux generated by the flame close to the wall. Because of the disparity of thermal characteristic times between the fluid and the solid, the aero-thermal part of the coupling is simplified by considering the fluid as a sequence of steady states. The heating of the metallic plate causes its deformation (the material has a viscoplastic behavior). The displacement of the fluid-structure interface is propagated through the fluid mesh. Based on similitudes between impinging reacting jets (flames) and non-reacting jets, coupled computations are performed in 2D and 3D configurations with an equivalent non-reacting hot jet.

Le travail présenté intègre l'utilisation des vibrations mécaniques dans les procédés de mise en forme à froid des matériaux de comportement viscoplastique. Les avantages de procédé de mise en forme assisté par vibrations se mesurent par une réduction de l'intensité de la force de mise en forme, une diminution de la contrainte d'écoulement, la réduction du frottement entre l'outil et lopin, une distribution uniforme de la déformation viscoplastique et une amélioration de l'état de surface de la pièce finale. A cet effet, un modèle analytique avec une seule source de vibration a été développé pour le procédé de forgeage. Les vibrations sinusoïdale et triangulaire ont été intégrées dans le modèle et il se trouve que les vibrations triangulaire est plus avantageux comparer à les vibrations sinusoïdal a cause du contrôle sur le rapport cyclique. On observe également que les paramètres principaux qui contrôlent la réduction de la force de forgeage sont : le rapport de vitesse définie par l'amplitude de la vitesse de vibration et la vitesse moyenne de l'outillage et la sensibilité à la vitesse de déformation. Le modèle analytique a été validé au moyen des simulations éléments finis réalisés avec le code de calcul FORGE2011 ® et les essais expérimentaux en dispositif expérimentale de pot vibrant a été conçu et les essais ont été réalisé sur les machines Lloyd LR30K ou ZWICK / Roell 1200. L'étude a ensuite été enrichi avec les plusieurs sources de vibrations. L'étude cinématique montre que les vibrations multiples générés une onde progressive. Cette nouvelle cinématique a été intégrée dans un modèle numérique sur le code de simulation par éléments finis FORGE2011 ® pour étudier le forgeage assisté par plusieurs sources des vibrations. L'analyse de résultats montre une réduction considérable de la force de forgeage et de la contrainte de l'écoulement ainsi qu'une modification de la loi de frottement générer par le changement de la direction du vecteur de vitesse. / This work is focused on the application of mechanical vibrations in cold metal forming process for visco-plastic domain to improve it by obtaining mean forming load reduction, low stresses, less friction, uniform distribution of strains and good surface qualities. For this purpose, analytical model has been developed for vibration assisted forming process with single vibration source. Sinusoidal and triangular vibration waveforms have been integrated in these models and it is found that triangular vibration is more advantageous as compare to traditional sinusoidal due to the control on cycle duty time. It is also observed that the main parameters that control the amount of forging force reduction are: the speed ratio defined by the amplitude of the vibration speed over the mean speed of the die and the sensitivity to strain rate. Analytical model is verified with the help of Finite element simulations performed in FORGE2011 ® and experiments performed with the help of specifically design vibration pot and Lloyd LR30K Universal testing machine or ZWICK/Roell 1200. The scope of this work is further extended by the use multiple vibration sources in the metal forming process. Kinematic study is carried out and it is found that application of multiple vibrations generates progress wave. This kinematic is applied in FORGE2011 ® during the metal forging process with the help of generic press and the reduction in forging load has been observed along with reduction in stress in the material and the modification of the friction on the lower die due to the change in velocity vector direction.

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.

Asphalt materials exhibit self-healing characteristics. In order to evaluate the healing capabilities of these materials, numerous investigations are carried out by researchers. These investigations include various definitions of healing and mechanical test methods to quantify the extent of healing. Irrespective of these extensive studies, there is no unique method to evaluate the healing in asphalt materials. The biasedness to the test conditions and complexity due to the involvement of other phenomena such as stress relaxation and strain recovery make the characterization of healing a challenging process. The work presented in this thesis includes the development of experimental and analytical approaches to contribute to the understanding of healing in asphalt materials. The thesis focuses on both fracture mechanics and damage mechanics-based approaches to evaluate the healing in asphalt materials. In the first case, three-point bending tests are carried out to characterize the healing following the fracture. The experimental protocol is carefully designed to avoid the effect due to low temperature physical hardening of the material during the rest period. Different healing indices appealing to linear elastic fracture mechanics and viscoelastic fracture mechanics are compared. While healing is generally defined based on the recovery of certain parameters following the rest period, the comparison across healing indices shows that the quantitative interpretations of healing are dependent on the post-processing methods. The damage mechanics-based approach discussed in this study includes creep and recovery tests in shear. Here, unlike the first case, the samples are not fractured during the test. The damage is considered to be a part of the viscoplastic deformation and the recovery of viscoplastic deformation is defined as the healing. The method proposed in this study is useful to isolate the viscoelastic effects to quantify the healing of damage. / Asfaltmaterial uppvisar självläkande egenskaper. För att utvärdera dess läkande förmågor utförs många undersökningar av forskare. Dessa undersökningar inkluderar olika definitioner av läkning och mekaniska testmetoder för att kvantifiera graden av läkningen. Trots dessa omfattande studier finns det ingen unik metod för att utvärdera läkning i asfaltsmaterial. Påverkan av testförhållandena och komplexiteten från inverkan av andra fenomen, såsom spänningsrelaxation och töjningsåterhämtning, gör karakteriseringen av läkning en utmanande process. Arbetet som presenteras i denna uppsats inkluderar utvecklingen av experimentella och analytiska metoder för att bidra till förståelsen av läkning i asfaltsmaterial. Uppsatsen fokuserar på både brottmekaniks- och skademekaniksbaserade metoder för att utvärdera läkning i asfaltsmaterial. I det första fallet utförs trepunktsböjningstest för att karakterisera läkningen efter brott. Det experimentella protokollet är noggrant designat för att undvika effekten av hårdnande vid låg temperatur under viloperioden. Dessutom jämförs olika läkningsindex baserade på linjärelastisk brottmekanik och viskoelastisk brottmekanik. Trots att läkning generellt definieras baserat på återhämtningen av vissa parametrar efter viloperioden, så visar jämförelsen av läkningsindexen att den kvantitativa tolkningen av läkning beror på metoderna för efterbehandling.  Den skademekaniksbaserade metoden som diskuteras i denna studie inkluderar kryp- och återhämtningstest i skjuvning. Till skillnad mot det första fallet så spricker här inte proven. Skadan anses istället vara del av den viskoplastiska deformationen och återhämtningen av viskoplastisk deformation är definierad som läkning. Metoden som föreslås i denna studie är användbar för att isolera de viskoelastiska effekterna för att kvantifiera läkningen av skada.

[pt] Este trabalho investiga o impacto de suspensões de grafeno funcionalizado sobre reologia de uma dispersão aquosa de Carbopol. Os derivados de grafeno usados foram óxido de grafeno (GO) e óxido de grafeno amino funcionalizado (AFGO). A variação da concentração e do grau de oxidação do grafeno funcionalizado foram avaliadas e relacionadas à reologia das suspensões. O GO foi produzido a partir da síntese de óxido de grafite pelo método de Hummers modificado, e caracterizado pelas técnicas de DRX, Raman, TGA, FTIR, XPS, TEM e AFM. O GO foi funcionalizado com trietilenotetramina por reação assistida por microondas para produzir AFGO, caracterizado por TGA e XPS. As suspensões foram caracterizadas por experimentos de potencial zeta e testes reológicos de escoamento em estado estacionário e oscilatório. A caracterização do GO mostra que grupos funcionais oxigenados foram incorporados em sua superfície grafítica. O GO oxidado por 96 horas (GO 96 h) apresentou maior distância interplanar e também apresentou menos camadas quando comparado com GO oxidado por 2 horas (GO 2 h). A caracterização do AFGO aponta que grupos aminados foram covalentemente ligados às nanofolhas GO e o GO com maior grau de oxidação produziu AFGO com maior nível de aminação. Todas as nanofolhas sintetizadas são coloidalmente estáveis quando suspensas em meio neutro. Todas as suspensões foram bem modeladas pela equação de Hershel-Bulkley. O aumento da concentração de nanofolhas nas suspensões prejudica a microestrutura do fluido e leva a uma diminuição da viscosidade, tensão de escoamento e elasticidade. O GO 96 h promoveu menor decréscimo na viscosidade, tensão de escoamento e elasticidade do que a suspensão GO 2 h. No caso do AFGO, o maior grau de aminação pode levar a uma queda mais pronunciada nas propriedades reológicas da suspensão. Para a suspensão com maior concentração de GO 96 h, observou-se o aparecimento de histerese em baixas taxas de cisalhamento. Esses resultados mostram que pequenas mudanças na superfície das nanofolhas de grafeno funcionalizado podem influenciar as respostas reológicas de um fluido não newtoniano. / [en] This work investigates the impact of functionalized graphene suspensions on the rheology of a Carbopol aqueous dispersion. The graphene derivatives used were graphene oxide (GO) and amino-functionalized graphene oxide (AFGO). The variation of the functionalized graphene concentration and oxidation was evaluated and related with the suspensions rheology. GO nanosheets were produced from synthesis of graphite oxide by modified Hummers method, and they were characterized by XRD, Raman,TGA, FTIR, XPS, TEM and AFM techniques. The GO was functionalized with triethylenetetramine by microwave assisted reaction to produce the AFGO, which was characterized by TGA and XPS techiniques. The suspensions were characterized by zeta potential experiments and rheological tests through steady-state and oscillatory flow. The GO characterization shows that oxygenated functional groups were incorporated in its graphitical surface. GO oxidized for 96 hours (GO 96 h) showed higher interplanar distance and also presented fewer layers when compared with GO oxidized for 2 hours (GO 2 h). The AFGO characterization points out that aminated groups were covalently attached to the GO nanosheets and the GO with a higher oxidation degree produced an AFGO with higher amination level. All synthesized nanosheets are colloidally stable when suspended at neutral media. All suspensions were well modeled by the Hershel-Bulkley equation. The increase of the nanosheets concentration in the suspensions impairs the level of fluid structure and leads to a decrease in viscosity, yield stress, and elasticity. The GO 96 h promoted a lower decrease in viscosity, yield stress and elasticity than the GO 2 h suspension. In the case of AFGO, the greater amination degree can lead to a more pronounced drop in the suspension rheological properties. For the suspension with a higher concentration of GO 96 h, it was observed the appearance of hysteresis at low shear rates. These results show that small changes in the graphene functionalized nanosheets surface can influence the rheological responses of a non-Newtonian fluid.

[pt] É uma prática usual em reologia medir o tensão limite de escoamento. Nessas medidas, a tensão limite de escoamento é definida como o máximo valor absoluto de tensão ao qual abaixo não ocorrem escoamentos irreversíveis. Sendo assim, tensão limite de escoamento aparente estimada é usada em conjunto com o critério de von Mises em qualquer escoamento complexo. Este critério compara esta medida a intensidade do segundo invariante do tensor deviatórico das tensões. Acontece que, para escoamento simples de cisalhamento, o mesmo é composto por tensões cisalhantes e diferenças de tensão normais, mas a contribuição do último nunca foi considerada na determinação experimental da tensão limite de escoamento. Em vista de avaliar a importância da contribuição das diferenças de tensões normais na tensão limite de escoamento aparente, foram realizadas uma sequência de testes de creep para cada material, estimando a tensão crítica que representa o valor médio obtido entre os valores das curvas de tensão nas quais o material escoa e não escoa com uma tolerância considerável. Depois disso, foram propostos testes para avaliar os valores de N1 − N2 e apenas N1 no nível de tensão crítica. E em seguida avaliando-se adequadamente a tensão limite de escoamento. Observou-se que, para alguns materiais, a contribuição das diferenças de tensões normais é muito maior do que a contribuiçõ da tensão cisalhante. Por fim, a validade do critério de von Mises para materiais elasto-viscoplásticos foi avaliada. Para este fim, com o intuito de generalizar o estudo, ensaios de compressão a volume constante e de tração foram realizados avaliando-se as correspondentes tensões limites de escoamento. Como conclusão mais importante, o critério de von Mises não foi considerado adequado como critério de falha para os materiais elasto-viscoplásticos analisados. / [en] It is usual practice in rheology to measure the yield stress in a simple shear flow. In these measurements, the yield stress is identified as the maximum value of the shear stress below which no irreversible flow occurs. Then, the thus determined yield stress is used in conjunction with the von Mises criterion in any complex flow. The latter compares it with the intensity of the deviatoric stress tensor. It happens that for simple shear flow the intensity of the deviatoric stress is composed of both the shear stress and the normal stress differences, but the contribution of the latter is never considered in the experimental determination of the yield stress. In view of assess the importance of the contribution of the normal stresses to the yield stress, a sequence o standard constant shear stress tests were performed for each material, estimating the critical stress which represents the mean value obtained between the stress values of the curves in which the material flows and does not flow with an accurate tolerance. After that, proposed tests were performed in order to obtain the values of N1 − N2 and solely N1 at the critical stresses. Following the appropriate yield stress evaluation. It was observed that for some materials the normal stress contribution is much larger than the shear stress contribution. Furthermore, the validity of the von Mises yielding criterion for elasto-viscoplastic materials was evaluated. For this purpose, in order to generalize the study for different flow conditions, constant volume squeeze flow and traction tests were performed evaluating the corresponding yield stresses. As the most important conclusion, the von Mises yielding criterion was considered not to be accurate representing yielding for the elastoviscoplastic materials analyzed.

No description available.

Ette thèse traite d’écoulements de fluides à seuil (ou viscoplastiques) en milieu confiné. Les difficultés analytiques et numériques sont dues à la multivaluation du tenseur des contraintes dans les zones plastiques ainsi qu’à la non-différentiabilité du problème de minimisation associé. Cette thèse s’articule en deux parties.Dans un premier temps, des simulations numériques parallèles très précises à l’aide d’algorithmes de dualité ont été effectuées. Elles ont permis de retrouver des résultats observés expérimentalement dont l’existence d’une ligne de glissement pour l’écoulement au dessus d’un obstacle et le caractère quasi-Poiseuille de la vitesse au-delà de cette ligne. Par ailleurs, la théorie de couche limite viscoplastique définie par Oldroyd (1947, à nombre de Bingham asymptotiquement grand) a été revisitée à nombre de Bingham modéré en milieu confiné. L’étude a mis en œuvre des allers-retours entre ces simulations et les expériences physiques de Luu et al. d’IRSTEA ainsi qu’une dérivation théorique. L’approximation de couche limite est vérifiée dans une certaine mesure à l’intérieur de la cavité. Une adaptation de la notion de couche limite viscoplastique est alors exhibée et permet d’étendre les scalings dérivés par Oldroyd (1947) et Balmforth et al. (J. of Fluid Mech, 2017). Ces scalings sont aussi généralisés au cas de la loi d’Herschel-Bulkley. Dans un second temps, on présente une analyse asymptotique des champs de vitesses et de contraintes pour des écoulements en faible épaisseur (ε). Un développement à l’ordre ε2 de la vitesse permet de trouver une équation de Reynolds à la même précision. Cette équation de Reynolds prolonge les résultats déjà existants dans le cadre newtonien, d’une part et dans le cadre fluide à seuil avec une surface libre, d’autre part. / This thesis is devoted to the flow of yield stress (or viscoplastic) fluids in pipes.Analytical and numerical difficulties lie in the multivaluation of the stress tensor in the plastic regions and in the non-differentiability of the associated minimization problem. This manuscript is organized following two main axes.First, very accurate numerical simulations were carried out using duality methods and parallel multifrontal solvers. Thus, experimental observations were recovered, namely the existence of a slip line for the flow over an obstacle and the Poiseuille-like behaviour of the velocity above this line. Moreover, the viscoplastic boundary layer theory defined by Oldroyd (1947 at high Bingham numbers) was revisited at moderate Bingham numbers in confined areas. This study provided an opportunity to go back and forth between these simulations and the physical measures of Luu et al. from IRSTEA and to perform a theoretical derivation. The boundary layer approximation is valid up to a certain extent in the cavity. An adaptation of the viscoplastic boundary layer definition is then given and allows to generalize the scalings shown by Oldroyd (1947) and Balmforth et al. (JFM 2017). These scalings are also generalized to the Herschel-Bulkley case. Then, an asymptotic analysis of the velocity and stress fields for thin layer (ε) flows is presented. A velocity development up to ε2 lets us find a Reynolds equation of same accuracy. This Reynolds equation extends the already existing results, on the one hand in the newtonian case and on the second hand for free surface flows.

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.