O presente trabalho apresenta um estudo numérico do escoamento laminar em cavidade quadrada aletada sob o efeito de convecção mista. O escoamento proposto é assumido bidimensional, laminar e permanente. Objetiva-se através do “Constructal Design” a obtenção de geometria ótima da aleta de forma a maximizar a transferência de calor entre o fluido que escoa no interior da cavidade e a aleta aquecida cuja base está localizada no centro da base da cavidade. Para isto é fixada a relação das dimensões externas da cavidade (H/L) = 1 e variada a relação entre altura e comprimento da aleta (H1/L1) para otimização da troca térmica. A área da aleta apresenta 5% da área total da cavidade e este valor é mantido fixo. O fluido que escoa no inteiror da cavidade possui as propriedades termofísicas do ar para Pr = 0,7. A variação das forças de empuxo no escoamento é realizada através do uso de diferentes números de Rayleigh no intervalo Ra= 10³ a 106. As diferentes magnitudes das forças inerciais serão aplicadas ao escoamento através da variação do número de Reynolds variando entre ReL = 10 e 1000. Para solução numérica das equações de conservação de massa, quantidade de movimento e energia é utilizado o método de volumes finitos (VFM), programa comercial Fluent®, sendo o acoplamento entre velocidade e pressão realizado através do algoritmo SIMPLEC e a discretização espacial pelo método upwind de primeira ordem. Os resultados apresentam um acréscimo significativo na transferência de calor entre a aleta e o fluido a medida que o número de Rayleigh aumenta. Considerando o caso de maior influência do mecanismo de transferência de calor por convecção mista houve um aumento de 779% em comparação com o mesmo caso considerando apenas convecção forçada, o que comprova a importância da convecção natural na maximização da transferência de calor entre cavidade e fluido para os casos analisados. / The present work shows a numerical study of laminar flow inside C-shaped lid-driven square cavity under mixed convection effect. The flow is assumed to be two-dimensional, laminar and permanent. The main objective of this work is by means of Constructal Design to maximize the heat transfer between the fluid and the heated central fin intruded in the bottom of the cavity. The aspect ratio of the cavity is fixed and the fin aspect ratio (H1/L1) varies from 0.1 to 10 ranges in order to maximize heat transfer. The ratio area between fin and cavity (H/L) = 1 is kept fixed at 5%. The thermophysical properties of fluid the air are set at Pr = 0,71. To vary the magnitude of buoyancy forces the Rayleigh number is ranged between Ra=10³ and 106.The inertial forces of flow are ranged by the use of different Reynolds numbers between ReL=10 and 1000. In order to solve the proposed problem, the commercial software Fluent® based on finite volume method was used to solve mass, momentum and energy equations, making the pressure-velocity couple using SIMPLEC method and the spatial discretization using first order upwind scheme. The results showed a significant increase of heat transfer between fin and fluid as consequence of Rayleigh number increase. Considering the mixed convection most influenced case, an increase of 779% was sense in comparison with the same case with forced convection mechanism only, which makes evident the importance of natural convection in the maximization of heat transfer inside cavity in the analized cases.

We study heat transfer in one dimension with and without convection, also called advection-diffusion. This is done using the Finite Element Method (FEM) to discretise the mathematical model, i.e. the heat equation. The results are compared to analytic Fourier series solutions. Our main result is that the FEM could be used to better model the heat transfer which allow for more accurate models than today's use of steady state models.

O presente trabalho apresenta um estudo numérico do escoamento laminar em cavidade quadrada aletada sob o efeito de convecção mista. O escoamento proposto é assumido bidimensional, laminar e permanente. Objetiva-se através do “Constructal Design” a obtenção de geometria ótima da aleta de forma a maximizar a transferência de calor entre o fluido que escoa no interior da cavidade e a aleta aquecida cuja base está localizada no centro da base da cavidade. Para isto é fixada a relação das dimensões externas da cavidade (H/L) = 1 e variada a relação entre altura e comprimento da aleta (H1/L1) para otimização da troca térmica. A área da aleta apresenta 5% da área total da cavidade e este valor é mantido fixo. O fluido que escoa no inteiror da cavidade possui as propriedades termofísicas do ar para Pr = 0,7. A variação das forças de empuxo no escoamento é realizada através do uso de diferentes números de Rayleigh no intervalo Ra= 10³ a 106. As diferentes magnitudes das forças inerciais serão aplicadas ao escoamento através da variação do número de Reynolds variando entre ReL = 10 e 1000. Para solução numérica das equações de conservação de massa, quantidade de movimento e energia é utilizado o método de volumes finitos (VFM), programa comercial Fluent®, sendo o acoplamento entre velocidade e pressão realizado através do algoritmo SIMPLEC e a discretização espacial pelo método upwind de primeira ordem. Os resultados apresentam um acréscimo significativo na transferência de calor entre a aleta e o fluido a medida que o número de Rayleigh aumenta. Considerando o caso de maior influência do mecanismo de transferência de calor por convecção mista houve um aumento de 779% em comparação com o mesmo caso considerando apenas convecção forçada, o que comprova a importância da convecção natural na maximização da transferência de calor entre cavidade e fluido para os casos analisados. / The present work shows a numerical study of laminar flow inside C-shaped lid-driven square cavity under mixed convection effect. The flow is assumed to be two-dimensional, laminar and permanent. The main objective of this work is by means of Constructal Design to maximize the heat transfer between the fluid and the heated central fin intruded in the bottom of the cavity. The aspect ratio of the cavity is fixed and the fin aspect ratio (H1/L1) varies from 0.1 to 10 ranges in order to maximize heat transfer. The ratio area between fin and cavity (H/L) = 1 is kept fixed at 5%. The thermophysical properties of fluid the air are set at Pr = 0,71. To vary the magnitude of buoyancy forces the Rayleigh number is ranged between Ra=10³ and 106.The inertial forces of flow are ranged by the use of different Reynolds numbers between ReL=10 and 1000. In order to solve the proposed problem, the commercial software Fluent® based on finite volume method was used to solve mass, momentum and energy equations, making the pressure-velocity couple using SIMPLEC method and the spatial discretization using first order upwind scheme. The results showed a significant increase of heat transfer between fin and fluid as consequence of Rayleigh number increase. Considering the mixed convection most influenced case, an increase of 779% was sense in comparison with the same case with forced convection mechanism only, which makes evident the importance of natural convection in the maximization of heat transfer inside cavity in the analized cases.

Le secteur tertiaire représente une source potentielle d'économie incontournable pour parvenir à réduire la dépendance énergétique de la France. Le taux de renouvellement du parc immobilier Français étant relativement faible, un effort doit être porté sur l'existant. Dans ces travaux, nous étudions un procédé innovant de rénovation par l'intérieur, dont l'élément principal est une paroi ventilée multifonctionnelle, assurant l'isolation, l'émission de chaleur basse température ainsi que la finition des surfaces murales. Les premiers travaux sur cette paroi ventilée furent menés sur un prototype dimensionné _a l'aide d'un modèle numérique simplifié. Deux séries d'expériences menées dans une cellule climatique nous ont permis de quantifier les flux de chaleur à travers le système. Le bon fonctionnement de la paroi ventilée repose sur les mécanismes de convection naturelle dans un canal vertical. Les résultats issus du prototype ont montré la présence de phénomènes complexes intervenant au sein de l'écoulement. Nous avons donc choisi d'étudier plus en détails les phénomènes thermoconvectifs dans un système du type source chaude/cheminée avant de poursuivre l'étude sur le système global. Une étude théorique et une expérience ont été menées sur un cas académique du problème. A l'issue des résultats expérimentaux, nous avons observé plusieurs régimes d'écoulements, dépendants du rapport de forme du canal et du nombre de Richardson en sortie. Enfin, nous proposons un modèle analytique de la paroi ventilée comprenant l'ensemble des variables géométriques influentes. Ce modèle a été implémenté dans l'environnement de simulation Trnsys, dans la perspective d'effectuer des simulations annuelles à l'échelle du bâtiment. / In France, energy consumption due to buildings heating is an important part of the global primary energy consumption. The tertiary sector represents an unavoidable source of economy in order to reduce energy dependency of France. The turnover of French real estate being relatively low, an effort must be focused on the existing. In this work, we investigate on an innovative process of internal thermal renovation, whose main element is a multifunctional ventilated wall, providing insulation, low temperature heat emission and the wall surfaces finishing. The first works on this ventilated wall were conducted on a prototype designed using a simplified numerical model. Two series of experiments conducted in a climatic cell allowed us to quantify the heat flow through the system. The smooth functioning of the ventilated wall is based on the natural convection in a vertical channel and the results from the prototype showed the presence of complex phenomena within the flow. We therefore chose to study in detail the thermoconvective phenomena in a chimney-like system before continuing the study of the overall component. A theoretical study and PIV experiment were conducted on an academic case of the problem. At the end of the experimental results, we observed several flow regimes, depending on the channel aspect ratio and the outlet Richardson number. Finally, we propose an analytical model of the ventilated wall including all influential geometrical variables. This model has been implemented in the simulation environment Trnsys with the perspective to make annual simulations on a building scale.

La réglementation thermique sur les bâtiments en France impose le recours aux énergies renouvelables. L'énergie solaire, à la fois thermique et photovoltaïque, présente un grand potentiel en termes de disponibilité. Néanmoins, du fait de son caractère intermittent, des systèmes de stockage associés sont nécessaires. Pour ce qui est des systèmes solaires thermiques, les capteurs solaires intégrant le stockage (CSIS) ont pour avantages une grande simplicité ainsi que des coûts réduits comparés à des systèmes plus classiques. Ils permettent également d'éviter la problématique de la mise en place du stockage, les rendant ainsi adaptés à la rénovation énergétique.Le présent travail de thèse étudie les phénomènes de convection naturelle dans une cavité fermée à haut rapport de forme (H/L = 13) représentant le réservoir de stockage d'un CSIS. Dans un premier temps, un état de l'art des systèmes existants ainsi qu'une étude bibliographique de l'étude de ces systèmes et des phénomènes s'y déroulant sont réalisés. Il y est mis en évidence l'importance de la stratification des températures au sein du réservoir pour l'amélioration des performances globales du système. Dans un second temps, une étude expérimentale de la cavité est effectuée pour la compréhension des écoulements liés à la mise en place d'une stratification thermique satisfaisante, ainsi que des paramètres pouvant l'influencer. Dans cette étude, une plaque de stratification est ajoutée à l'intérieur de la cavité pour canaliser le fluide chaud ascendant, limitant ainsi le brassage. La zone de chauffe représentant l'échange entre la boucle solaire et le réservoir est placée au milieu de la face supérieure du système, résultant en une bonne stratification. Un écoulement inverse est observé au niveau de la sortie de ce canal. Ce phénomène est principalement responsable de l'écart de température obtenu entre le haut et le bas du système. Une étude numérique CFD du système a donc été conduite pour étudier l'influence de différents paramètres sur le comportement thermique et de l'écoulement. La diminution de l'écoulement inverse permet une meilleure circulation du fluide dans la globalité du système. Afin de déterminer si la stratification thermique observée dans le système de stockage conduit à une performance énergétique annuelle satisfaisante du capteur auto-stockeur, une modélisation nodale a été effectuée. Différents paramètres tels que la position de la zone d'échange entre le collecteur et le réservoir de stockage, ou encore l'épaisseur de l'isolant ont été étudiés pour quantifier leur influence sur les performances du système complet. / Thermal regulation for dwellings in France makes the resort to renewable energy mandatory. Solar energy has an important potential regarding its availability, for both thermal and photovoltaic applications. However, due to its intermittency, there is a need for a storing device. Concerning solar thermal systems, Integrated Collector Storage Solar Water Heaters (ICSSWH) have the advantage of a simple design and use combined with reduced costs compared to more classical devices. These systems also enable to avoid the problem of finding space for the installation of a storage system in an existing building, making them more suitable for energy renovation of houses.This present work studies the phenomena of natural convection in an enclosed cavity with high aspect ratio (H/L = 13) representing the storage of an ICSSWH. Firstly, a state of the art of existing systems as well as a literature review of these systems and of the phenomena occurring in them are conducted. The importance of thermal stratification inside the storage tank for the increase of global efficiency for the whole system is highlighted. Secondly, an experimental study of the cavity is performed to understand the fluid flow linked to the satisfying stratification obtained in the system, and find the parameters influencing it. In the experimental setup, a stratification plate is added inside the cavity to harness the ascending hot flow, hence limiting the mixing of fluid which is responsible for the destratification. The heated zone, which represents the exchange between the solar collector and the storage, is placed in the middle of the upper wall of the system, resulting in a good thermal stratification. A reverse flow can be observed at the channel outlet, which is the main parameter for the high temperature difference obtained between the top and the bottom of our system. A CFD study of the experimented system is then conducted to determine the influence of several parameters on the thermal and fluid flow behavior of the system. The decrease of reverse flow allows a better loop of the fluid in the whole cavity. In order to discover whether the thermal stratification observed in the storage leads to a better annual efficiency for the whole ICSSWH, a nodal model is used. Several parameters such as the position of the heat exchange zone between the collector and the storage or the insulation thickness are studied to quantify their influence on thermal stratification, and hence on the overall efficiency of the system.

Les effets des transferts radiatifs sur les écoulements de convection naturelle sont étudiés en régimes transitionnel et turbulent. On considère des mélanges air/H2O/CO2 confinés dans des cavités cubiques différentiellement chauffées. Des simulations numériques de référence sont entreprises jusqu'à Ra=3x108 en couplant une méthode spectrale de collocation pour l'écoulement et une méthode de lancer de rayons, associée à un modèle ADF, pour le rayonnement. Pour l'étude du régime turbulent, une modélisation des transferts radiatifs basée sur un filtrage spatial est proposée : les contributions filtrées sont résolues par la méthode de lancer de rayons sur un maillage lâche et les contributions de sous-maille sont résolues de manière analytique dans l'espace de Fourier. Ce modèle est combiné à la simulation numérique directe de l'écoulement à Ra=3x109. Les transferts radiatifs ont pour effet de diminuer la stratification thermique verticale et d’augmenter la circulation générale. Lorsque les six parois de la cavité sont noires et le gaz transparent, deux zones de stratification thermique instable apparaissent en amont des couches limites verticales. Dès Ra=5x106, une instabilité de type Rayleigh-Bénard se développe dans ces zones, induisant des écoulements instationnaires. Lorsque les parois adiabatiques sont parfaitement réfléchissantes, les parois isothermes noires et le gaz rayonnant, des écoulements instationnaires chaotiques sont obtenus à partir de Ra=3x107. Des rouleaux contra-rotatifs à la sortie des couches limites verticales sont observés, ce qui suggère qu'une instabilité de force centrifuge soit responsable de la transition. / Radiative transfer effects on natural convection flows are investigated in transitional and turbulent regimes. Air/H2O/CO2 mixtures contained in cubical differentially heated cavities are considered. Benchmark numerical simulations are carried out up to Ra=3x108 by coupling a spectral collocation method for the flow and a ray tracing method, associated with an ADF model, for radiation. In order to study the turbulent regime, a radiative transfer model based on spatial filtering is proposed: filtered contributions are solved with the ray tracing method on a coarse grid and sub-grid contributions are obtained analytically in Fourier space. This model is combined with the direct numerical simulation of the flow at Ra=3x109. The effects of radiative transfer are a decrease of the vertical thermal stratification and an increase of the flow driven in the cavity. When the six cavity walls are black and the gas is transparent, two unstably stratified zones appear upstream the vertical boundary layers. From Ra=5x106, a Rayleigh-Bénard type instability in these zones triggers the unsteadiness. When the adiabatic walls are perfectly reflecting, the isothermal walls are black and the gas is participating, unsteady chaotic flows are obtained in this case from Ra=3x107. Counter rotating rolls at the exit of the vertical boundary layers are observed, which suggests that transition to unsteadiness is due to centrifugal forces.

O presente trabalho apresenta um estudo numérico do escoamento laminar em cavidade quadrada aletada sob o efeito de convecção mista. O escoamento proposto é assumido bidimensional, laminar e permanente. Objetiva-se através do “Constructal Design” a obtenção de geometria ótima da aleta de forma a maximizar a transferência de calor entre o fluido que escoa no interior da cavidade e a aleta aquecida cuja base está localizada no centro da base da cavidade. Para isto é fixada a relação das dimensões externas da cavidade (H/L) = 1 e variada a relação entre altura e comprimento da aleta (H1/L1) para otimização da troca térmica. A área da aleta apresenta 5% da área total da cavidade e este valor é mantido fixo. O fluido que escoa no inteiror da cavidade possui as propriedades termofísicas do ar para Pr = 0,7. A variação das forças de empuxo no escoamento é realizada através do uso de diferentes números de Rayleigh no intervalo Ra= 10³ a 106. As diferentes magnitudes das forças inerciais serão aplicadas ao escoamento através da variação do número de Reynolds variando entre ReL = 10 e 1000. Para solução numérica das equações de conservação de massa, quantidade de movimento e energia é utilizado o método de volumes finitos (VFM), programa comercial Fluent®, sendo o acoplamento entre velocidade e pressão realizado através do algoritmo SIMPLEC e a discretização espacial pelo método upwind de primeira ordem. Os resultados apresentam um acréscimo significativo na transferência de calor entre a aleta e o fluido a medida que o número de Rayleigh aumenta. Considerando o caso de maior influência do mecanismo de transferência de calor por convecção mista houve um aumento de 779% em comparação com o mesmo caso considerando apenas convecção forçada, o que comprova a importância da convecção natural na maximização da transferência de calor entre cavidade e fluido para os casos analisados. / The present work shows a numerical study of laminar flow inside C-shaped lid-driven square cavity under mixed convection effect. The flow is assumed to be two-dimensional, laminar and permanent. The main objective of this work is by means of Constructal Design to maximize the heat transfer between the fluid and the heated central fin intruded in the bottom of the cavity. The aspect ratio of the cavity is fixed and the fin aspect ratio (H1/L1) varies from 0.1 to 10 ranges in order to maximize heat transfer. The ratio area between fin and cavity (H/L) = 1 is kept fixed at 5%. The thermophysical properties of fluid the air are set at Pr = 0,71. To vary the magnitude of buoyancy forces the Rayleigh number is ranged between Ra=10³ and 106.The inertial forces of flow are ranged by the use of different Reynolds numbers between ReL=10 and 1000. In order to solve the proposed problem, the commercial software Fluent® based on finite volume method was used to solve mass, momentum and energy equations, making the pressure-velocity couple using SIMPLEC method and the spatial discretization using first order upwind scheme. The results showed a significant increase of heat transfer between fin and fluid as consequence of Rayleigh number increase. Considering the mixed convection most influenced case, an increase of 779% was sense in comparison with the same case with forced convection mechanism only, which makes evident the importance of natural convection in the maximization of heat transfer inside cavity in the analized cases.

Dans le contexte du stockage des déchets radioactifs en milieu poreux, nous considérons l’équation de convection-diffusion instationnaire et sa discrétisation par des méthodes numériques. La discontinuité des paramètres physiques et la variabilité des échelles d’espace et de temps conduisent à utiliser des discrétisations différentes en temps et en espace dans différentes régions du domaine. Nous choisissons dans cette thèse le schéma volumes finis en dualité discrète (DDFV) et le schéma de Galerkin Discontinu en temps couplés à une méthode de décomposition de domaine de Schwarz de type relaxation d’ondes optimisées (OSWR), ce qui permet de traiter des maillages espace-temps non conformes. La principale difficulté réside dans l’obtention d’une discrétisation amont du flux convectif qui reste locale à un sous-domaine et telle que le schéma monodomaine soit équivalent au schéma multidomaine. Ces difficultés sont appréhendées d’abord en une dimension d’espace où différentes discrétisations sont étudiées. Le schéma retenu introduit une inconnue hybride sur les interfaces entre cellules. L’idée du décentrage amont par rapport à cette inconnue hybride est reprise en dimension deux d’espace, et adaptée au schéma DDFV. Le caractère bien posé de ce schéma et d’un schéma multidomaine équivalent est montré. Ce dernier est résolu par un algorithme OSWR dont la convergence est prouvée. Les paramètres optimisés des conditions de Robin sont obtenus par l'étude du taux de convergence continu ou discret. Différents cas-tests, dont l’un est inspiré du stockage des déchets nucléaires, illustrent ces résultats. / In the context of nuclear waste repositories, we consider the numerical discretization of the non stationary convection diffusion equation. Discontinuous physical parameters and heterogeneous space and time scales lead us to use different space and time discretizations in different parts of the domain. In this work, we choose the discrete duality finite volume (DDFV) scheme and the discontinuous Galerkin scheme in time, coupled by an optimized Scwharz waveform relaxation (OSWR) domain decomposition method, because this allows the use of non-conforming space-time meshes. The main difficulty lies in finding an upwind discretization of the convective flux which remains local to a sub-domain and such that the multidomain scheme is equivalent to the monodomain one. These difficulties are first dealt with in the one-dimensional context, where different discretizations are studied. The chosen scheme introduces a hybrid unknown on the cell interfaces. The idea of upwinding with respect to this hybrid unknown is extended to the DDFV scheme in the two-dimensional setting. The well-posedness of the scheme and of an equivalent multidomain scheme is shown. The latter is solved by an OSWR algorithm, the convergence of which is proved. The optimized parameters in the Robin transmission conditions are obtained by studying the continuous or discrete convergence rates. Several test-cases, one of which inspired by nuclear waste repositories, illustrate these results.

La problématique des échanges turbulents en convection thermique est abordée au travers de deux écoulements complémentaires en eau, Prandtl compris entre 2.5 et 6.5.Le premier est une cellule de Rayleigh-Bénard asymétrique. La plaque chaude située au bas de la cellule présente des rugosités périodiques de taille contrôlée. La plaque froide en haut de la cellule est, au contraire, lisse. Le flux de chaleur est imposé à la plaque chaude. Dans la gamme de forçage explorée dans ce travail, les deux plaques sont indépendantes ce qui permet une comparaison in-situ entre une plaque lisse et une plaque rugueuse. Par comparaison avec des résultats de référence, nous mettons en évidence une augmentation importante du nombre de Nusselt (flux de chaleur adimensionné) associé à la plaque rugueuse alors que celui de la plaque lisse demeure inchangé. Cette augmentation semble pilotée par la hauteur des rugosités choisies.Le second écoulement étudié ici est un écoulement en canal thermique. La géométrie du dispositif permet de s'affranchir des conditions aux limites à l'injection du flux de chaleur. L'expérience étant inclinable, il est possible d'étudier les effets de la stratification sur l'écoulement, puisqu'elle est de plus en plus importante avec l'augmentation de l'inclinaison. En fonction du couple (flux de chaleur, inclinaison), quatre régimes d'écoulements distincts peuvent apparaître. Nous proposons une modélisation pour trois d'entre eux, dont un modèle d'écoulement turbulent. Ce modèle est ensuite validé à plus haut Reynolds dans une nouvelle expérience à grande échelle, qui permet d'explorer le régime de turbulence inertielle induit par l'écoulement. / The question of turbulent exchanges in thermal convection is studied with two kind of complementary flows in water, Prandtl number is kept between 2.5 and 6.5.The first one is an asymmetrical Rayleigh-Bénard cell. The hot bottom plate of the cell is covered with periodical obstacles of controlled size, whereas the cold top plate is kept smooth. The thermal flux is imposed at the hot plate. Considering the range of thermal forcing imposed to the cell, the two plates can be considered independent, which allows in-situ comparison between the smooth plate and the rough one. Comparison with reference results underlines a huge increase of the Nusslet number (non dimensional thermal transfer) of the rough plate whereas the transfer of the smooth plate remains unchanged. This increase seems to be controlled by the height of the obstacles.The second flow studied here is a thermal channel. The geometry of the apparatus allows the channel to be free of the influence of boundary conditions. By tilting the experiment, it is possible to study the influence of stratification on the flow. The larger the inclination is, the larger is the stratification. Considering the (heat flux, inclination) couple, four regimes of flow can be identified. We detailed a model for three of them, including a turbulent flow model. This model is then tested at higher Reynolds number in a new experiment at large scale. In this case the inertial range of turbulence is investigated.

Depuis sa formation, la Terre subit un refroidissement lent. La chaleur provenant du noyau et de la désintégration des éléments radioactifs présents dans le manteau est évacuée vers la surface par convection. L’évolution des structures thermiques ainsi créées contrôle de nombreux phénomènes de surface tels que le mouvement des continents et le niveau marin. L’étude présentée ici s’attache à déterminer quelles structures convectives terrestres peuvent être reconstruites, sur quelle période de temps et avec quelle précision. La chaoticité de la convection implique que les incertitudes initialement présentes sur le champ de température croissent exponentiellement au cours du temps et peuvent créer des structures convectives artificielles dans les modèles. A l’aide de la méthode des expériences jumelles initialement développée par Lorenz [1965] en météorologie, le temps de Lyapunov et l’horizon de prédiction sont calculés pour la première fois en géodynamique mantellique. Différentes rhéologies sont étudiées. La valeur du temps de Lyapunov pour notre modèle le plus proche de la Terre suggère qu’une erreur de 5% sur les conditions initiales limite l’horizon de prédiction à 95 millions d’années. D’autre part, la qualité de la prédiction des structures thermiques dépend de notre capacité à décrire de façon réaliste les propriétés rhéologiques du manteau. L’utilisation d’une rhéologie pseudo-plastique dans les modélisations de convection en 3D sphérique, permet aujourd’hui de générer une tectonique de plaques compatible au premier ordre avec les caractéristiques cinématiques de la surface terrestre. Une stratégie cohérente de reconstruction peut alors être élaborée. L’état thermique actuel du manteau est reconstruit en imposant les vitesses de surface de ces 200 derniers millions d’années [Seton et al., 2012; Shephard et al., 2013] sur un modèle de convection généré par le code StagYY [Tackley, 2008]. La morphologie et la position des slabs reconstruits varient considérablement avec le contraste de viscosité et la pseudo-plasticité. L’erreur introduite par l’utilisation de rhéologies différentes lors des reconstructions est ainsi plus importante que les erreurs liées aux incertitudes sur les conditions initiales et les vitesses de surface. Ces résultats montrent l’importance du choix la rhéologie sur la qualité des prédictions réalisées. Ils mettent également en évidence rôle clé du contraste de viscosité et de la pseudo-plasticité pour reconstruire des slabs cohérents et des subductions plates, structures propres à la convection terrestre. / Since its formation, the Earth is slowly cooling. The heat produced by the core and the radioactive decay in the mantle is evacuated toward the surface by convection. The evolving convective structures thereby created control a diversity of surface phenomena such as vertical motion of continents or sea level variation. The study presented here attempts to determine which convective structures can be predicted, to what extent and over what timescale. Because of the chaotic nature of convection in the Earth’s mantle, uncertainties in initial conditions grow exponentially with time and limit forecasting and hindcasting abilities. Following the twin experiments method initially developed by Lorenz [1965] in weather forecast, we estimate for the first time the Lyapunov time and the limit of predictability of Earth’s mantle convection. Our numerical solutions for 3D spherical convection in the fully chaotic regime, with diverse rheologies, suggest that a 5% error on initial conditions limits the prediction of Earth’s mantle convection to 95 million years. The quality of the forecast of convective structures also depends on our ability to describe the mantle properties in a realistic way. In 3D numerical convection experiments, pseudo plastic rheology can generate self-consistent plate tectonics compatible at first order with Earth surface behavior [Tackley, 2008]. We assessed the role of the temperature dependence of viscosity and the pseudo plasticity on reconstructing slab evolution, studying a variety of mantle thermal states obtained by imposing 200 million years of surface velocities extracted form tectonic reconstructions [Seton et al., 2012; Shephard et al., 2013]. The morphology and position of the reconstructed slabs largely vary when the viscosity contrast increases and when pseudo plasticity is introduced. The errors introduced by the choices in the rheological description of the mantle are even larger than the errors created by the uncertainties in initial conditions and surface velocities. This work shows the significant role of initial conditions and rheology on the quality of predicted convective structures, and identifies pseudo plasticity and large viscosity contrast as key ingredients to produce coherent and flat slabs, notable features of Earth’s mantle convection.