[pt] A dissertação apresenta um estudo da resolução de problemas de acoplamento termo-hidráulico através do uso técnicas numéricas. O estudo tem por objetivo a validação de implementações de efeitos de gradientes de temperatura no fluxo de umidade. Para isto utiliza-se o programa UNSATCHEM- 2D que simula o fluxo bidimensional de umidade em meios de saturação variável, transporte de calor, transporte de CO2 e transporte de solutos. O programa desenvolvido por Simunek e Suarez (1993), possui código aberto (escrito em Fortran 77), onde as equações governantes de fluxo e transporte de calor são resolvidas por elementos finitos. A partir do programa UNSATCHEM-2D, foram realizadas implementações no código a fim de permitir o acoplamento parcial (staggered) do fluxo de umidade e o fluxo de calor que reproduz o fenômeno descrito por de Vries (1958). Paralelamente, desenvolveram-se mudanças no programa MTool(Tecgraf/PUC- Rio), utilizado como pré- e pós-processador. Posteriormente, foram validadas as implementações realizadas no programa comparando-se os resultados da solução numérica com resultados experimentais a partir de um ensaio de coluna (Bach, 1989), com resultados considerados satisfatórios. Finalizando, são apresentadas simulações com fontes de calor interna e externa paralelamente com o emprego de condições de contorno atmosféricas que simulam coberturas de resíduos de minério. / [en] The dissertation presents a study of the resolution of thermo-hydraulic coupling problems, through the use of numerical techniques. The study has the objective the implementations validation of the temperature gradients effects in the moisture flux. For this, the program UNSATCHEM-2D is used to simulate the flow of water in mediums of variable saturation to two dimensions, transport of heat, transport of CO2 and solute transport. The program developed by Simunek and Suarez (1993), has open code (written in Fortran 77), where the governing equations of flow and transport of heat are solved for the method of the finite elements. Starting from the original program UNSATCHEM-2D, alterations were carried out in the code in order to allow the partial coupling (staggered) of moisture flow and heat flow as described by de Vries (1958). In parallel, changes were developed in the software MTool (Tecgraf/PUC-Rio), used as pre- and posprocessor. Later, the modifications carried out in the program were validated by the comparison of the numeric solution results with experimental results starting from a column test (Bach, 1989), with acceptable results. Concluding, simulations are presented with internal and external sources of heat together with the use of atmospheric boundary conditions, one that simulates covers for mining wastes.

This thesis reports on an experimental study of slow viscoelastic flow through models of fibrous porous media. The models were square arrays of parallel cylinders, with solid volume fractions or ‘solidities’ of 2.5%, 5.0%, and 10%. An initial study using a Newtonian fluid provided a baseline for comparison with results for two Boger fluids, so that the effects of fluid elasticity could be determined. Boger fluids are elastic fluids that have near constant viscosities and can be used in experiments without having to account for shear-thinning effects. The experimental approach involved measurements of pressure loss through the three arrays and interior velocity measurements using particle image velocimetry (PIV). For the Newtonian flows, pressure loss measurements were in good agreement with the analytical predictions of Sangani and Acrivos (1982). PIV measurements showed velocity profiles which were symmetrical and independent of flow rate. Pressure loss measurements for the Boger fluid flows revealed that the onset of elastic effects occurred at a Deborah number of approximately 0.5, for both fluids and the three arrays. Flow resistance data collapsed for the two Boger fluids, and increased with solidity. For all three models, the flow resistance increased monotonically with Deborah number, reaching values up to four times the Newtonian resistance for the 10% model. PIV measurements showed that the transverse velocity profiles for the Newtonian and Boger fluids were the same at Deborah numbers below the elastic onset. Above onset, the profiles became skewed. The skewness, like the flow resistance, was observed to increase with both Deborah number and solidity. In the wake regions between cylinders in a column, periodic flow structures formed in the spanwise direction. The structures were staggered from column to column, consistent with the skewing. As either Deborah number or solidity increased, the flow structures became increasingly three-dimensional, and the stagger became more symmetric. An analysis of fluid stresses reveals that the elastic flow resistance is attributed to additional normal stresses caused by shearing, and not by extension.

This thesis reports on an experimental study of slow viscoelastic flow through models of fibrous porous media. The models were square arrays of parallel cylinders, with solid volume fractions or ‘solidities’ of 2.5%, 5.0%, and 10%. An initial study using a Newtonian fluid provided a baseline for comparison with results for two Boger fluids, so that the effects of fluid elasticity could be determined. Boger fluids are elastic fluids that have near constant viscosities and can be used in experiments without having to account for shear-thinning effects. The experimental approach involved measurements of pressure loss through the three arrays and interior velocity measurements using particle image velocimetry (PIV). For the Newtonian flows, pressure loss measurements were in good agreement with the analytical predictions of Sangani and Acrivos (1982). PIV measurements showed velocity profiles which were symmetrical and independent of flow rate. Pressure loss measurements for the Boger fluid flows revealed that the onset of elastic effects occurred at a Deborah number of approximately 0.5, for both fluids and the three arrays. Flow resistance data collapsed for the two Boger fluids, and increased with solidity. For all three models, the flow resistance increased monotonically with Deborah number, reaching values up to four times the Newtonian resistance for the 10% model. PIV measurements showed that the transverse velocity profiles for the Newtonian and Boger fluids were the same at Deborah numbers below the elastic onset. Above onset, the profiles became skewed. The skewness, like the flow resistance, was observed to increase with both Deborah number and solidity. In the wake regions between cylinders in a column, periodic flow structures formed in the spanwise direction. The structures were staggered from column to column, consistent with the skewing. As either Deborah number or solidity increased, the flow structures became increasingly three-dimensional, and the stagger became more symmetric. An analysis of fluid stresses reveals that the elastic flow resistance is attributed to additional normal stresses caused by shearing, and not by extension.

Utilizing the analytical solution of the steady state sharp interface saltwater intrusion model in coastal aquifers, a multi-objective optimization formulation of pumping rates and well locations in a coastal aquifer is formulated to solve problems in water management practice. The proposed optimization problem uses progressive genetic algorithm technique and the method developed is applied to the previous work of Cheng et al. [2000]. Through this analysis, several other applications are provided to demonstrate the use of the model in practical applications. This work is the first to optimize pumping rates as well as well locations simultaneously in coastal aquifer management. Known the limitation of the analytical solution, the work is expanded to cover the physics of saltwater intrusion in a more realistic way. This is variable density flow in a variably saturated porous medium. In this method, mixing between two fluids such as saltwater and freshwater can be described and the porous medium is also expanded to cover saturated and unsaturated zones together. One of the objectives is to develop a three dimensional physical model, verify the model, and apply to various applications in coastal aquifers. The developed model, TechFlow, is used to investigate instability issues associated with the numerical solution of the Elder problem in the perspective that includes physical instability issues associated with density differences used in numerical solutions, sensitivity of the solution to idealization irregularity, and the importance of accurate estimation of the velocity field and its association to the grid density levels that is necessary to solve the problem accurately. Saltwater intrusion hydrodynamics in a beach under the influence of tidal effects is also investigated using TechFlow. Based on the results of TechFlow with the use of various boundary conditions for the transport equation, the saltwater intrusion hydrodynamics in a beach under the influence of tidal effects shows unique dynamics. These solutions are primarily affected by density differences, tidal effects on a mild slope, variably saturated porous medium and finite domain solution condition. TechFlow is also used to investigate saltwater upconing beneath pumping wells both two- and three-dimensional applications.

Le contrôle thermique des composants électroniques embarqués dans les engins spatiaux est souvent assuré par des boucles fluides diphasiques à pompage capillaire (Loop Heat Pipe (LHP) ou Capillary Pumped Loop (CPL)). La présente étude est centrée sur les évaporateurs des LHP. Ils sont composés principalement d’un bâti métallique, d’une mèche poreuse et de cannelures. Le milieu poreux est initialement saturé en liquide. La charge thermique à évacuer est appliquée sur le bâti entraînant la vaporisation du liquide au sein de la mèche. La vapeur est ensuite récoltée au sein des cannelures pour être évacuée. L’étude est effectuée sur une cellule unitaire de l’évaporateur. Dans le but d’étudier les transferts de masse et de chaleur, un modèle de réseau de pores 3D dit mixte a été développé. Les champs de pression et de température sont calculés à partir des équations macroscopiques tandis que la capillarité est gérée à l’aide d’une approche réseau de pore classique. L’un des avantages d’une telle formulation est de pouvoir accéder à la répartition des phases liquide et vapeur au sein de l’espace poral du milieu poreux. Il a ainsi été mis en évidence qu’une zone diphasique (zone où le liquide et la vapeur coexistent) se met en place pour une large gamme de flux lorsque la vapeur apparait dans la structure capillaire. Cette zone diphasique est localisée sous le bâti métallique et est corrélée avec les meilleures performances thermiques de l’évaporateur. Cette observation diffère fortement de l’hypothèse souvent considérée de la présence d’une zone sèche dans cette région. Trois positions différentes de cannelures ont été étudiées. Il a ainsi pu être mis en évidence que la plus large gamme de flux, pour laquelle les performances de l’évaporateur sont les meilleures, est obtenue lorsque les cannelures sont usinées à la surface extérieure de la mèche. Toujours dans le but d’améliorer les performances thermiques de l’évaporateur, une étude paramétrique a été menée pour mettre en évidence les paramètres qui influencent positivement la conductance de l’évaporateur. Finalement, l’étude de l’influence d’une mèche biporeuse/bidispersée, c’est-à-dire d’un milieu poreux caractérisé par deux tailles de pores/liens différentes, a été menée. La distribution des phases liquide et vapeur au sein de la structure capillaire bidispersée est différente de celle d’un milieu mono-poreux du fait des chemins préférentiels créés par les larges pores. Par ailleurs, l’analyse thermique a montré qu’un tel milieu poreux permet de réduire considérablement la température du bâti ainsi que d’augmenter les performances thermiques de l’évaporateur. Un deuxième modèle basé sur une approche continue a été développé. Cette méthode utilise l’algorithme IMPES (IMplicit Pressure Explicit Saturation) et est couplé à la résolution du champ de température avec changement de phase. Ce type de résolution permet d’accéder à un champ de saturation. Les résultats ainsi obtenus sont en bon accord avec ceux prédits par le modèle réseau de pores mixte. Le modèle continu, moins gourmand en temps de calcul, permet d’envisager des simulations sur une plus grande partie de l’évaporateur. / The thermal control of electronic devices embedded in spacecraft is often carried out by capillary twophase loop systems (Loop Heat Pipe (LHP) or Capillary Pumped Loop (CPL)). This thesis focuses on the LHP evaporators. They mostly consist of a metallic casing, a porous wick and vapour grooves. The porous medium is initially saturated with liquid. The heat load is applied at the external surface of the casing inducing the vaporisation of the liquid within the wick. The vapour is then evacuated thanks to the vapour grooves. A unit cell of the evaporator is studied and corresponds to our computational domain. A so-called 3D mixed pore network model has been developed in order to study the heat and mass transfers. Pressure and temperature fields are computed from macroscopic equations, while the capillarity is managed using the classical pore network approach. The main advantage of such formulation is to obtain the liquid-vapour phase distribution within the porous medium pore space. The work highlights that a two-phase zone (characterized by the coexistence of the liquid and the vapour) exists for a large range of fluxes when vaporisation takes place within the capillary structure. This twophase zone is located right under the casing and is positively correlated with the best evaporator thermal performances. This result differs from the often made assumption of a dry region under the casing. Three different groove locations are tested. This investigation highlights that evaporator thermal performances are the best over a large range of fluxes for grooves manufactured at the external surface of the wick. In complement, a parametric study is performed to highlight parameters which impact positively the evaporator thermal performances. Finally, a biporous/bidispersed wick, i.e. a wick with a bimodal pore/throat size distribution, is studied. The liquidvapour phase distribution within the capillary structure is different from the one for a monoporous structure due to preferential vapour paths created by the large throats and pores-network. Moreover, the thermal analysis shows that such a porous medium permits to reduce considerably the evaporator wall temperature and to increase the evaporator thermal performances. A second model is developed based on a continuum approach. This method uses the IMPES (IMplicit Pressure Explicit Saturation) algorithm coupled with the heat transfer with phase change. Results are in good agreement with those predicted by the mixed pore network model. The continuum model, requiring less computing time, should allow considering larger sub domains of the evaporator.

Le contrôle thermique des réactions dans les réacteurs à lits fixes nécessite la maîtrise du transfert thermique en proche paroi. Afin d’optimiser leur conception et en particulier de maintenir un chemin réactionnel le plus proche possible de l’isothermicité, un modèle de transfert thermique pertinent dans un tel milieu est indispensable. Deux premières thèses au Laboratoire d'Energétique et de Mécanique Théorique et Appliquée ont permis de mettre au point un modèle dispersif utilisable à cœur de réacteur. La présente thèse met en évidence un effet de paroi non négligeable causé par la variation de porosité du lit à l’approche de la paroi. Ce travail propose un raccordement du cœur, où le modèle thermique précédent est légitime, à la paroi, par une couche homogénéisée. Cette approche est validée par une caractérisation expérimentale des paramètres du modèle sur un dispositif de laboratoire en utilisant une méthode d’inversion originale reposant sur une approche Bayesienne. Une validation de cette caractérisation dans une géométrie plus proche d’un réacteur industriel est également effectuée / Thermal control for chemical reactions conducted in fixed-bed reactors requires a fine knowledge about heat transfer in the near-wall region of the bed. In order to optimize the process design, to minimize the mass of catalyst that is needed, one must attempt to maintain the operational path of the reactor the closest of the isotherm functioning. A pertinent and reliable heat transfer model is then required. Two thesis were realized in the Laboratoire d'Energétique et de Mécanique Théorique et Appliquée in order to construct a corresponding dispersive model. This model can be used at the core of the reactor, where the porous medium can be homogenized. This thesis shows a wall effect that cannot be neglected. This wall effect is directly caused by the important variation of the porosity distribution near the wall of the bed. This work proposes a junction between a core layer, where the previous model is still valid, and the wall through an homogenized near-wall layer. This approach is validated by the mean of an experimental characterization of the model parameters, using an original inversion technique based on a Bayesian approach. A validation of this characterization in another geometry is also done. This geometry is closer to an industrial reactor

Parmi la multitude des données pétrophysiques utilisées pour décrire une formationgéologique, certaines permettent spécifiquement d'en prédire la capacité de production,à savoir : la porosité, la perméabilité intrinsèque et le coefficient de Klinkenberg.Dans le cas particulier des gas shales, ces trois propriétés essentielles sont extrêmementdifficiles à mesurer précisément, compte tenu de la complexité de cesroches. Cette thèse s'inscrit dans la continuité de travaux menés au laboratoire I2M départementTREFLE, sur l'analyse et l'amélioration de la méthode Pulse Decay, quiconstitue la méthode de mesure transitoire classiquement utilisée dans l'industriepétrolière pour identifier une ou plusieurs des propriétés d'intérêt. Les multiplespoints faibles de la méthode Pulse Decay sont ici présentés, de même que les pointsforts de la nouvelle méthode issue du perfectionnement de la méthode Pulse Decay,à savoir la méthode Step Decay, développée au laboratoire I2M et brevetée en partenariatavec TOTAL. Plus précisément, les performances de la méthode Step Decaysont ici étudiée aussi bien numériquement qu'expérimentalement, en condition homogène,comme en condition hétérogène. Ce manuscrit fournit également les résultatsd'une analyse portant sur la méthode Pulse Decay sur broyat, qui forme une alternativepossible à la méthode Pulse Decay sur carotte mais dont la fiabilité est fortementremise en question. / Among the multitude of petrophysical data used to describe a geological formation,some of them allow specifically to predict the production capacity, namely: the porosity,the intrinsic permeability and the Klinkenberg coefficient. In the particular case ofgas shales, these three essential properties are extremely difficult to measure precisely,because of the complexicity of these rocks. This thesis is the continuity ofworks led in I2M laboratory-TREFLE department, on the analysis and the improvementof the Pulse Decay method, which is the classical transient method of measurementused in the oil&gas industry to identify one or several of the properties ofinterest. The numerous weaknesses of the Pulse Decay method are here presented,as the strengths of the new method derived from the improvement of the Pulse Decaymethod, namely the Step Decay method, developed in I2M laboratory and patentedin partnership with TOTAL. More exactly, the performances of the Step Decaymethod are here studied numerically as well as experimentally, in homogeneouscondition, as in heterogeneous condition. This manuscript provides also the results ofan analysis dealing with the Pulse Decay method on cuttings, which forms a possiblealternative to the Pulse Decay method on plug but which reliability is highly questioned.

Ce travail de recherche porte sur le comportement des milieux poreux (triphasiques), plus particulièrement les sols non saturés sous sollicitations hydro-mécaniques. Un modèle constitutif élastoplastique couplé est développé. Ce modèle original est formulé selon les principes suivants: une loi constitutive est développée pour décrire le comportement de chaque phase (squelette solide, liquide, et gaz). Ensuite, des relations de couplage sont ajoutées entre chacune des phases. Pour le comportement du squelette solide, une loi élastoplastique non associée est adoptée, avec deux surfaces de charges, en cisaillement et en compression. La partie hydrique est décrite par une formulation qui permet de prendre en compte l’effet d’hystérésis. Ce modèle a été enrichi par une relation de couplage hydromécanique qui permet d’exprimer la pression d’entrée d’air en fonction de la porosité. Ensuite, le couplage complet se fait avec la contrainte effective de Bishop en utilisant une nouvelle définition du paramètre de succion χ grâce à laquelle, les différents phénomènes présents dans la réponse des milieux poreux sous différentes sollicitations peuvent être reproduits. Ce modèle est validé par une confrontation à des données expérimentales issues de la littérature sur différents types de sol (sable, limon,…). Le modèle est implanté dans le code aux éléments finis Cast3M. L’analyse de problèmes particuliers, tels que la mise en œuvre d’un cas test d’un sol d’assise soumis à un cycle pluvial, ainsi que l’étude de la stabilité d’une pente, permette de montrer la capacité du modèle à reproduire le comportement des milieux poreux non saturés. / This work focuses on the behaviour of porous triphasic media, particularly on unsaturated soils subjected to hydromechanical loading. A coupled elastoplastic constitutive model has been developed. This original model is formulated according to the following principles: (1) a constitutive law describing the behaviour of different phases (solid skeleton, liquid and gas). (2) coupling relationships between each phase. For the behaviour of the solid skeleton, a non associated elastoplastic constitutive law is adopted, with two loading surfaces: shear surface and compression cap surface. The hydric part is discribed using a formulation which allows to take into account the hysteresis effect. This model has been extended using a hydromechanical coupling relation between the air entry value and the porosity. Then the coupling is completed with the Bishop effective stress, using a new definition for the suction parameter χ. Using this formulation, the various phenomena present in the porous media behaviour under different loading can be reproduced. The developed model has been validated through a comparison with experimental data on different types of soil (sand, silt,…). This model is implemented in the french finite element code Cast3M. The analysis of specific problems, such as (1) the study of shallow foundation subjected to cyclic rain event, as well as (2) the study of slope stability, show the model capacity to reproduce the behaviour of unsaturated porous media.

L'étude du transfert de chaleur et de masse dans les milieux poreux saturés et insaturés fortement chauffés à leur surface possèdent de nombreuses applications, notamment en archéologie, en agriculture et en géothermie. La première partie de ce travail concerne l'amélioration de la méthode AHC (Accumulation de chaleur latente) qui permet de traiter le changement de phase, dans un milieu homogène : l'intervalle de changement de température au moment du changement de phase apparaît comme un paramètre important, et il doit être choisi proportionnel à la taille des mailles. Des résultats à la fois précis et lisses sont obtenus grâce à un raffinement du maillage localisé près de l'interface de changement de phase. La deuxième partie se rapporte à l'estimation des propriétés thermophysiques du sol par problème inverse à l'aide de données à la fois synthétiques et expérimentales. La méthode de Gauss-Newton avec relaxation et l'algorithme de Levenberg-Marquardt sont utilisés pour résoudre le problème inverse. Le choix de l'intervalle de température de la méthode AHC apparaît crucial : la convergence n'est obtenue parfois qu'au prix d'un enchaînement de plusieurs problèmes inverses. La troisième partie présente un modèle simple pour calculer la conductivité thermique effective d'un milieu granulaire contenant une faible quantité d'eau liquide. La forme exacte de ces ménisques est calculée à l'équilibre. Les résultats montrent un phénomène très net d'hystérésis quand on étudie la variation de la conductivité thermique effective en fonction de la quantité d'eau liquide ; un futur travail concernant un nouveau modèle insaturé, limité au cas du régime pendulaire et présenté à la fin de cette thèse, devrait pouvoir utiliser ces résultats. / In this work, we are interested in studying heat and mass transfer in water saturated and unsaturated porous medium with a strong heating at the surface. Applications concerned are archaeology, agriculture and geothermal engineering. The first part of this work concerns the improvement of the AHC (Apparent Heat Capacity) method used in the numerical resolution of phase change problem in a homogeneous medium: the phase change temperature interval, over which the heat capacity varies, appears as a key parameter which must be chosen proportional to the mesh size. Accurate and smooth results are obtained thanks to a local refinement of the mesh near the phase change interface. The second part is about the estimation of the thermophysical properties of the soil by inverse problem using both synthetic and experimental data. The Damped Gauss-Newton and the Levenberg-Marquardt algorithms are used to solve the problem. In relation with the AHC method, the choice of the phase change temperature interval caused convergence problems which have been fixed by chaining many inverse problems. The obtained results show good convergence to the desired solution. The third part presents a simple model to calculate the effective thermal conductivity of a granular medium which contains a small quantity of liquid water. The exact shape of the liquid menisci between the grains is calculated at equilibrium. The effective thermal conductivity experiences a hysteresis behavior with respect to the liquid volume. A future work that concerns a new unsaturated model, restricted to the pendular regime and detailed at the end of this thesis, should be able to use this result.

L’utilisation des matériaux composites devient de plus en plus importante dans les structures aéronautiques de nouvelle génération. Le gain de masse engendré, et donc de carburant, pousse les constructeurs aéronautiques à les employer de façon optimale. Néanmoins, ces matériaux se dégradent rapidement lorsqu’ils sont soumis à des flux de chaleur importants, entraînant une perte de leur résistance mécanique. Ce problème peut être dramatique pour la sécurité des passagers car la tenue de ces nouvelles structures peut ne plus être assurée dans le cas d’un incendie. Les méthodes actuelles de certification de la tenue au feu des matériaux composites aéronautiques reposent principalement sur l’utilisation de moyens expérimentaux, dont les résultats ne sont représentatifs que des conditions particulières dans lesquelles les essais ont été réalisés. La compréhension des différents phénomènes thermiques, chimiques et mécaniques intervenant lors de la dégradation de ces matériaux, avec l’appui de simulations numériques et d’expériences, peut permettre d’améliorer les méthodes existantes et donc d’optimiser les futures structures aéronautiques dès la phase de conception.Cette étude s’est attachée à développer et valider un modèle thermo-chimique de dégradation des matériaux composites multi-dimensionnel et multi-constituants. Ce modèle permet de traiter des cinétiques de dégradation complexes suivant plusieurs réactions de décompositions et de prendre en compte le transport des gaz produits depuis leur formation jusqu’à leur évacuation hors du matériau. L’utilisation de lois d’homogénéisation avancées est proposée afin de rendre compte des effets des transformations sur les transferts de chaleur et de masse se produisant au sein du matériau. L’application du modèle thermo-chimique à un cas de dégradation sous flux thermique connu mais non-uniforme dans un environnement contrôlé permet de confronter les résultats de simulation aux mesures expérimentales et ainsi de valider l’approche multi-constituants adoptée. Enfin, l’étude numérique de la dégradation d’un composite soumis à une flamme met en avant l’effet des gaz de décomposition éjectés à l’interface sur le flux thermique pariétal échangé. / Composite materials are increasingly used in new generation aircraft structures. Mass and as a consequence fuel savingsencourage aircraft manufacturers to use them optimally. However, these materials can degrade quickly when exposed tosignificant heat fluxes, resulting in a loss of mechanical strength. This problem can be dramatic for passenger safety asmechanical resistance of such innovative structures can not be ensured in case of fire events. Current certification methodsof fire resistance of aeronautical composite materials are mainly based on experiments, that are only representative of thespecific conditions under which they were carried out. The understanding of thermal, chemical and mechanical phenomenaoccurring during the decomposition of these materials, with the support of numerical simulations and experiments, can helpimproving existing methods and optimizing the future aeronautical structures from the design chain. This study deals withthe development and validation of a multi-components and multi-dimensional thermo-chemical model of decomposing compositematerials. It can deal with complex degradations following several decomposition reactions as well as transport ofpyrolysis gases from their formation up to their ejection out of the material. The use of advanced homogenization laws isproposed to account for the chemical transformations on heat and mass transfers occurring in the material. The applicationof the thermo-chemical model to a laser degradation study under known but non-uniform heat flux in a controlled environmentallows to confront the simulation results with experimental measurements and thus validate the multi-componentsapproach. Finally, the numerical analysis of a decomposing composite material submitted to a flame highlights the effectof emitted decomposition gases on the exchanged parietal heat flux.