Dans cette thèse, nous nous intéressons à l'étude d'un modèle de transport de sédiments en nous plaçant sous deux angles d'approche différents. L'un concerne la modélisation numérique du problème et propose une méthode de résolution numérique basée sur un solveur de Riemann approché pour le système de Saint-Venant-Exner qui reste un des systèmes les plus répandus pour traiter le transport sédimentaire par charriage. Ce dernier repose sur un couplage du modèle hydraulique de Saint-Venant et du modèle morphodynamique d'Exner. Le point essentiel de la méthode proposée se situe au niveau du traitement du couplage de ce système. En effet, il existe deux stratégies; la première consiste à découpler la résolution de la partie fluide de la partie solide et les faire interagir à des instants donnés alors que la seconde considère une résolution couplée du système en mettant à jour conjointement les grandeurs hydrauliques et solides aux mêmes instants. Se posera alors la question du choix de la stratégie de résolution pour laquelle nous apporterons des éléments de réponses en comparant les deux approches. L'autre se concentre sur la mise en place d'une méthodologie pour l'étude des incertitudes liées au même modèle. Pour ce faire, nous proposons une formulation stochastique du système de Saint-Venant-Exner et nous cherchons à caractériser la variabilité des sorties par rapport aux paramètres d'entrée naturellement aléatoires. Cette première étude révélera la nécessité de revenir à un système de Saint-Venant avec un fond bruité pour étudier la sensibilité des grandeurs hydrauliques par rapport aux perturbations topographiques. / In this thesis, we are interested on the study of a sediment transport model through two different approaches. One of them concerns the numerical modelling of the problem and proposes a numerical problem-solving method based on an approximate Riemann solver for the Saint-Venant-Exner system which is one of the most common model to deal with sedimentary bed-load transport. This last one is based on a coupling between the hydraulic model of Saint-Venant and the morphodynamic model of Exner. The key point of the proposed method is the treatment of the coupling issue. Indeed, there exists two strategies; the first one consists on decoupling the resolution of the fluid part from the solid part and making them interact at fixed times whereas the second one considers a coupled approach to solve the system by jointly updating the hydraulic and solid quantities at same times. We then raise the issue of the choice of the strategy for which we suggest answers by comparing both approaches. The other one focuses on the development of a methodology to study the uncertainties related to the model previously mentioned. To this end, we propose a stochastic formulation of the Saint-Venant-Exner system and we look for characterizing the variabilities of the outputs in relation to the naturally random input parameters. This first study reveals the need for a return to the Saint-Venant system with a perturbed bed to understand the sensitivity of the hydraulic quantities on the topographical perturbations.

Depuis de nombreuses années, le givrage a été identifié comme un danger dans le domaine de l’aéronautique.L’accrétion de givre se produit lorsque des gouttelettes d’eau surfondue se déposent sur une surface, enparticulier le bord d’attaque d’une aile ou la lèvre d’entrée d’air moteur, et gèlent après l’impact. Ceci peutensuite engendrer une dégradation des performances aérodynamiques, un dysfonctionnement des sondes ouencore un endommagement du moteur. C’est pourquoi cette problématique est étudiée avec attention. Lesessais en vol et en soufflerie étant longs et coûteux, la simulation numérique de l’accrétion de givre est devenueun outil nécessaire dans le processus de conception et de certification des avions. Cette thèse s’inscrit dans le contexte de la modélisation 3D de l’accrétion de givre, et plus particulièrement des couches limites dynamique et thermique qui se développent autour du corps givré. Les outils numériques devant être rapides et robustes, l’approche proposée dans cette thèse pour le calcul aérodynamique est une méthode couplée Euler/couche limite intégrale. Ainsi, un modèle intégral est développé pour représenterle développement de la couche limite dynamique. La partie thermique est modélisée soit par une méthodesimplifiée basée sur des approches algébriques, soit par une méthode intégrale. Cette modélisation des coucheslimites dynamique et thermique est valable sur paroi lisse ou rugueuse et permet de fournir notamment lecoefficient de frottement et le coefficient d’échange thermique qui sont nécessaires pour un calcul d’accrétion degivre. Les équations intégrales de couche limite, associées à leurs relations de fermeture, sont ensuite résoluespar une méthode Volumes-Finis sur maillage surfacique non structuré, qui est bien adaptée pour les géométriescomplexes. De plus, des traitements numériques spécifiques sont mis en œuvre pour améliorer la précision dela méthode au voisinage du point d’arrêt et pour rendre le code robuste au passage du décollement.Après la validation de la méthode de couche limite, le code est utilisé dans les chaînes de givrage 2D et 3Dde l’ONERA pour des applications d’accrétion de givre. Ceci permet de montrer l’intérêt de la méthode entermes de robustesse et de précision par rapport aux codes de couche limite habituellement utilisés dans lescodes de givrage actuels. / Icing has since long been identified as a serious issue in the aeronautical world. Ice accretion occurs whensupercooled water droplets impinge on a surface, particularly the leading edge of a wing or an engine inlet, andfreeze after the impingement. This can lead to degradation of aerodynamic performances, sensor malfunctionor engine damage. This is why this issue is being carefully studied. The lengthy and costly flight and windtunnel tests have made numerical simulation of ice accretion a necessary tool in the aircraft design andcertification process. The present work deals with the 3D numerical modeling of ice accretion, and more particularly the modeling of the dynamic and thermal boundary layers that develop around an iced body. Since numerical tools must befast and robust, the approach proposed in this thesis for aerodynamic computation is a coupled Euler/integralboundary layer method. Thus, an integral model is developed to represent the development of the dynamicboundary layer. The thermal part is modeled either by a simplified method based on algebraic approaches,or by an integral method. This modeling of the dynamic and thermal boundary layers is valid on smoothor rough walls and provides the friction coefficient and heat exchange coefficient that are necessary for thecalculation of ice accretion. The integral boundary layer equations, associated with their closure relations,are then solved by a Finite-Volume method on unstructured surface mesh, that is well suited for complexgeometries. In addition, specific numerical treatments are implemented to improve the accuracy of the methodin the vicinity of the stagnation point and to make the code robust to separated boundary layers.After validation of the boundary layer method, the code is used in ONERA’s 2D and 3D icing tools foricing applications. This demonstrates the value of the method in terms of robustness and accuracy comparedto the boundary layer codes more commonly used in current icing tools.

Le Disjoncteur Basse-Tension (DBT) est un appareil classique de la distribution électrique depuis plus de cinquante ans. Mais aujourd'hui, avec l'arrivée de produits bas-coût fabriqués dans les pays émergents, les industriels sont soumis à une forte pression pour développer de nouveaux systèmes moins encombrants, utilisant d'autres matériaux, ou incorporant davantage de fonctionnalités. Cette recherche est très compliquée dans la mesure où le DBT est un système hautement multi-physique (mécanique, thermique, physique des matériaux, physique des plasmas, ...). De fait, le développement de nouveaux produits passe par un processus empirique long et coûteux. Cet effort pourrait être réduit par l'utilisation de modèles prédictifs permettant d'arriver plus vite à un système fonctionnel. De nos jours, avec l'augmentation des moyens de résolution numérique, de plus en plus de travaux portent sur la description multi-physique en 3D du DBT et notamment sur la chambre de coupure ou l'arc électrique est amorcé, se déplace et doit être éteint, l'objet de nôtre étude. Le travail de cette thèse se divise en deux axes complémentaires : le développement d'un modèle fluide 3D en méthode des volumes finis simulant l'arc électrique et son déplacement dans la chambre de coupure; ainsi que la mise en place d'un dispositif expérimental permettant d'analyser le phénomène physique en œuvre. Pour ces deux points la problématique est abordée dans une configuration simplifiée de DBT où l'arc se déplace entre deux rails parallèles dans une chambre parallélépipédique. Basé sur le savoir-faire du groupe AEPPT, un modèle numérique est établi pour simuler le plasma thermique. Les particularités de ce modèle, du fait de l'application, sont la nécessité d'une résolution précise du champ magnétique en utilisant le calcul de Biot & Savart pour les conditions limites ainsi que l'utilisation de méthodes permettant le déplacement et la commutation de l'arc. La validation de ce modèle se fera à géométrie similaire par confrontation avec l'expérience. En s'inspirant de précédents travaux nous avons réalisé une maquette expérimentale composée d'un réacteur faisant office de chambre de coupure et d'un mécanisme permettant l'amorçage de l'arc dans le réacteur par ouverture rotative du contact à vitesse contrôlée. D'autres paramètres modifiables sont la taille du réacteur ainsi que les matériaux qui le constituent. Les diagnostiques disponibles en plus de la mesure de courant et de tension sont l'imagerie rapide et la mesure de pression en différents points de la chambre de coupure. Notre expérience est utile pour la réalisation d'études paramétriques en découplant facilement les paramètres. En outre, par la mise en évidence des phénomènes prépondérants, notre maquette aide à la mise en place du modèle en plus de permettre sa validation expérimentale. Cette thèse est donc une étape cruciale vers la mise en place d'un modèle prédictif. / Low-Voltage Circuit Breakers (LVCB) are classical apparatuses of electrical distribution since more than fifty years. But nowadays, with the outbreak of low-cost products from the developing countries, industry is under a strong stress in order to improve their devices by making them more compact, using different materials or to implement new functionalities. This research is harsh since LVCB are highly multiphysics systems (mechanics, thermal properties, materials, plasma physics, ...). Therefore, developing new products goes through a long and expensive empirical process. Those efforts could be reduced by using predictive models allowing to get faster to a functional device. With the improvements of the numerical solution capacity, there are more and more works toward the 3D multiphysical description of the LVCB, especially on the extinction chamber where the electrical arc is ignited, moved and must be quenched. This is the subject of our work. The study described here is divided in two complementary parts: development of a 3D fluid model with finite volume method simulating the electrical arc and its movement inside the arc extinction chamber; and the set-up of experimental means to analyse this physical phenomenon. For those two points, we use a simplified LVCB configuration with an arc moving between two parallel rails inside a rectangular box chamber. Based on AEPPT's know-how, a numerical model is established to simulate thermal plasma. Particularities on this model, due to the application, are the resolution of Biot&Savart law to calculate precisely the magnetic field for the boundary condition and the development of methods to model the arc roots movement and commutation of the arc from the moving contact to the rail. Validation of this model will be done with the same geometry by confrontation with the experiment. Inspired by precedent works we designed a test apparatus with a reactor representing the extinction chamber of the LVCB and an opening mechanism allowing arc ignition by contact opening at a specified speed. Other parameters such as size of the chamber and materials can be modified. Measurements will include high speed imaging and pressure acquisition in several points of the reactor in addition to the classical current and voltage measurements. This experiment is useful for parametric studies with its easy uncoupling of the parameters. Moreover, by highlighting the dominating phenomena for arc movement, this set-up helps in the build-up of the model over and above the experimental validation.

Les capacités installées d’énergie éolienne continuent à croître rapidement et prennent une place de plus en plus significative dans le monde. Au fur et à mesure, les études menées sur la conception, la sureté de fonctionnement et la supervision de la chaîne éolienne ont pris progressivement de l’importance. Deux axes de recherche ont été privilégiés dans cette thèse. Le premier concerne la continuité de service d'une éolienne connectée au réseau en présence de défaut de court-circuit entre spires dans une phase du stator de la génératrice asynchrone à cage d'écureuil. L'analyse du défaut ainsi que son impact sur le système éolien et notamment sur la qualité de la puissance produite souligne l'intérêt de développement d'un algorithme de détection et d'isolation rapide, dédié par la suite à la reconfiguration de la commande. Ainsi, une commande tolérante au défaut (CTD) a été conçue de manière à éviter l'arrêt de la production, compenser l'impact de défaut et garder des performances acceptables de la qualité d'énergie produite. Le travail effectué s'est articulé sur les observateurs à mode glissant (OMG), communément connus comme outil puissant pour la supervision et la commande à la fois. Le deuxième axe porte sur la sécurité structurale et la stabilité du système éolien sous contraintes vibratoires. Les travaux se répartissent en deux parties complémentaires : L'établissement d'un modèle numérique tridimensionnel (3-D) sous un logiciel d’analyse par éléments finis (ANSYS) et la réalisation des essais vibratoires sous différentes excitations au sein d'une plateforme vibratoire (TREVISE). Dans ce cadre, un modèle numérique (3-D) d'une éolienne à axe horizontal couplée à un mât et une fondation adéquats a été développé en utilisant la méthode de volumes finis (FVM) afin d'appréhender son comportement vibratoire. Les essais vibratoires expérimentaux valident le modèle numérique et permettent l’identification de la réponse dynamique de la structure d'une manière fine. De plus, nous avons élaboré un modèle expérimental de la tenue de l’éolienne aux contraintes vibratoires de formes aléatoire, sinusoïdale et impulsionnelle. / The wind energy capacity carries on growing quickly and taking an increasingly significant place in the world. Progressively, research studies dealing with designing and supervising wind turbines have become more important. Two areas of research were developed in this thesis. The first one concerns the continuity of service of a wind turbine connected to the grid while an inter-turn short-circuit fault is present in the stator phase of the induction squirrel cage generator. The analysis of the fault as well as its impact on the wind turbine system and mainly on the quality of the produced power highlights the interest of development of a fast detection and isolation algorithm, dedicated to the reconfiguration of the control law. Hence, a fault tolerant control scheme has been established in order to avoid stopping production, compensate the fault impact and maintain acceptable performances of the quality of the produced energy. The carried out work was based on sliding mode observers, commonly known as robust tools for monitoring and controlling at the same time. The second axis concerns the structural modeling and stability checking of the wind system under vibratory stresses. The work is divided into two complementary parts: The establishment of a three-dimensional (3-D) numerical model using a finite element analysis software (ANSYS) and the realization of vibratory tests under different excitations within the platform (TREVISE). In this framework, a numerical (3-D) model of a horizontal axis wind turbine coupled to a suitable tower and foundation was developed basing on the finite volume method (FVM) in order to analyze its vibratory behavior. The experimental vibratory tests validate the numerical model and allow the identification of the dynamic response of the structure in a precise way. In addition, we have developed an experimental model of the behavior of the wind turbine under vibratory stresses of random, sinusoidal and impulse shapes.

Dans le but d’améliorer le comportement thermo-mécanique des aubes de turbine présentes dans les moteurs d’avion développés par SAFRAN, il est nécessaire de mettre au point un nouveau matériau composite. Le procédé de fabrication de ce matériau est complexe et requiert une densification par voie liquide divisée en deux étapes. Cette thèse s’intéresse à la modélisation numérique de la première étape appelée Slurry Cast/APS.Celle-ci consiste en l’injection et le confinement, dans la préforme fibreuse, de particules préalablement mises en suspension. Pour cela, nous avons développé à l’échelle des fibres,un modèle qui utilise les équations de Navier-Stokes incompressibles et monophasiques ,l’équation de Phillips [Phillips et al., 1992] et une loi rhéologique [Krieger, 1972]. Après validation des résultats numériques par comparaison avec des résultats expérimentaux [Hampton et al., 1997] et théoriques [Belfort et al., 1994], le modèle est utilisé pour simuler l’écoulement autour de géométries de tissage proches du matériau étudié. / In order to improve thermo-mechanical behavior of tubine blades in SAFRAN engines plane, a new composite material is necessary. The manufacturing process to obtain this composite is intricate and requires a two steps fluid densification process. This thesis focuses on numerical simulation of the first one called Slurry Cats/APS. In this step, suspended particles are introducted and captured in the reinforcement. For that purpose,we carry out a model at fiber scale, using Navier-Stokes equations in incompressible and monophasic formulation, Phillips equations [Phillips et al., 1992] and a rheological law [Krieger, 1972]. After validation step consisting in a comparison of computational results with experiments [Hampton et al., 1997] and theorical law [Belfort et al., 1994], this model has been used to simulate flow around geometries similar to those encountered in our composite material.

Le but de cette thèse est de faire l'étude de méthodes de volumes finis pour des équations aux dérivées partielles déterministes et stochastiques; nous effectuons des simulations numériques et démontrons des résultats de convergence d'algorithmes.Au Chapitre 1, nous appliquons un schéma semi-implicite en temps combiné avec la méthode de volumes finis généralisés SUSHI pour la simulation d'écoulements à densité variable en milieu poreux; il vient à résoudre une équation de convection-diffusion parabolique pour la concentration couplée à une équation elliptique en pression. Nous présentons ensuite une méthode de simulation numérique pour un problème d'écoulements à densité variable couplé à un transfert de chaleur.Au Chapitre 2, nous effectuons une étude numérique de l'équation de Burgers non visqueuse en dimension un d'espace, avec des conditions aux limites périodiques, un terme source stochastique de moyenne spatiale nulle et une condition initiale déterministe. Nous utilisons un schéma de volumes finis combinant une intégration en temps de type Euler-Maruyama avec le flux numérique de Godunov. Nous effectuons des simulations par la méthode de Monte-Carlo et analysons les résultats pour différentes régularités du terme source. Il apparaît que la moyenne empirique des réalisations converge vers la moyenne en espace de la condition initiale déterministe quand t → ∞. Par ailleurs, la variance empirique converge elle aussi en temps long, vers une valeur qui dépend de la régularité et de l'amplitude du terme stochastique.Au Chapitre 3, nous démontrons la convergence d'une méthode de volumes finis pour une loi de conservation du premier ordre avec une fonction de flux monotone et un terme source multiplicatif faisant intervenir un processus Q-Wiener. Le terme de convection est discrétisé à l'aide d'un schéma amont. Nous présentons des estimations a priori pour la solution discrète dont en particulier une estimation de type BV faible. A l'aide d'une interpolation en temps, nous démontrons deux inégalité entropiques vérifiées par la solution discrète, ce qui nous permet de prouver que la solution discrète converge selon une sous-suite vers une solution stochastique faible entropique à valeurs mesures de la loi de conservation.Au Chapitre 4, nous obtenons des résultats similaires à ceux du Chapitre 3 dans le cas où la fonction flux n'est pas monotone; le terme de convection est discrétisé à l'aide d'un schéma monotone. / This thesis bears on numerical methods for deterministic and stochastic partial differential equations; we perform numerical simulations by means of finite volume methods and prove convergence results.In Chapter 1, we apply a semi-implicit time scheme together with the generalized finite volume method SUSHI for the numerical simulation of density driven flows in porous media; it amounts to solve a nonlinear convection-diffusion parabolic equation for the concentration coupled with an elliptic equation for the pressure. We then propose a numerical scheme to simulate density driven flows in porous media coupled to heat transfer. We use adaptive meshes, based upon square or cubic volume elements.In Chapter 2, We perform Monte-Carlo simulations in the one-dimensional torus for the first order Burgers equation forced by a stochastic source term with zero spatial integral. We suppose that this source term is a white noise in time, and consider various regularities in space. We apply a finite volume scheme combining the Godunov numerical flux with the Euler-Maruyama integrator in time. It turns out that the empirical mean converges to the space-average of the deterministic initial condition as t → ∞. The empirical variance also stabilizes for large time, towards a limit which depends on the space regularity and on the intensity of the noise.In Chapter 3, we study a time explicit finite volume method with an upwind scheme for a first order conservation law with a monotone flux function and a multiplicative source term involving a Q-Wiener process. We present some a priori estimates including a weak BV estimate. After performing a time interpolation, we prove two entropy inequalities for the discrete solution and show that it converges up to a subsequence to a stochastic measure-valued entropy solution of the conservation law in the sense of Young measures.In Chapter 4, we obtain similar results as in Chapter 3, in the case that the flux function is non-monotone, and that the convection term is discretized by means of a monotone scheme.

Les approches 1D sont valides pour la modélisation des écoulements en rivière tant qu'ils peuvent être négligés dans la plaine d'inondation. Quand ces approximations sont invalides, il est fait appel aux modèles 1D à casiers ou aux modèles 2D. Dans le premier cas, les transferts de quantité de mouvement sont toutefois négligés. Pour la modélisation 2D des écoulements en rivière, un maillage précis du lit mineur doit être réalisé pour une prise en compte correcte de sa topographie. Le maillage est donc composé de mailles de petite taille et le pas de temps de calcul doit donc être réduit pour garantir la stabilité du schéma numérique. Une méthode alternative consiste à coupler les modèles 1D et 2D. Ces modèles 1D-2D ne sont généralement pas complètement satisfaisant dans la mesure où seuls les transferts de masse entre les modèles sont pris en compte. Dans le logiciel proposé SW12D, le couplage est réalisé en incluant le modèle 1D dans le 2D et en prenant en compte les transferts de masse et de quantité de mouvement. Ce formalisme permet une réduction considérable du nombre de mailles, en comparaison à un maillage 2D classique, et par conséquent de la durée de simulation. Des tests numériques ainsi que des exemples sur topographie réelle ont été réalisés avec un modèle 2D classique et SW12D. L'approche proposée représente correctement les pertes de charge dues à la sinuosité de l'écoulement ainsi que les court-circuits de méandre. Par ailleurs, les résultats sont suffisamment précis pour des applications d'ingénierie courante et obtenus avec un temps de calcul réduit. Cela constitue une amélioration significative par rapport aux modèles 1D et 2D.

L'objectif de cette thèse est l'analyse d'erreur a posteriori et la proposition de stratégies d'adaptivité basées sur des critères d'arrêt et de raffinement local de maillage. Nous traitons une classe d'équations paraboliques dégénér ées multidimensionnelles modélisant des problèmes importants pour l'industrie. Au chapitre 1 nous considérons le problème de Stefan instationaire a deux phases qui modélise un processus de changement de phase régi par la loi de Fourier. Nous régularisons la relation entre l'enthalpie et la température et nous discrétisons le problème par la méthode d'Euler implicite en temps et un schéma numérique conforme en espace tel que les élément finis conformes, ou les volumes finis centrés aux sommets du maillage. Nous démontrons une borne supérieure de la norme duale du résidu, de l'erreur sur l'enthalpie dans L2(0; T;H-1) et de l'erreur sur la température dans L2(0; T;L2), par des estimateurs d'erreur entièrement calculables. Ces estimateurs comprennent : un estimateur associé à l'erreur de régularisation, un estimateur associé à l'erreur d'une méthode de linéarisation (par exemple, la méthode de Newton), un estimateur associé à l'erreur en temps et un estimateur associé à l'erreur du schéma en espace. Par conséquent, ces estimateurs permettent de formuler un algorithme adaptatif de résolution où les erreurs associées peuvent être équilibrées. Nous proposons également une stratégie de raffinement local de maillages. En fin, nous prouvons l'efficacité de nos estimations d'erreur a posteriori. Un test numérique illustre l'efficacité de nos estimateurs et la performance de l'algorithme adaptatif. En particulier, des indices d'efficacité proches de la valeur optimale de 1 sont obtenus. Au chapitre 2 nous développons des estimations d'erreur a posteriori pour l'écoulement de Darcy polyphasique et isothermique, décrit par un système couplé d'équations aux dérivées partielles non linéaires et d'équations algébriques non linéaires. Ce système est discrétisé en espace par une méthode de volume finis centrés par maille et la méthode d'Euler implicite en temps. Nous etablissons une borne supérieure d'une norme duale du résidu augmentée d'un terme qui tiens compte de la non-conformité des volumes finis par des estimateurs d'erreur a posteriori entièrement calculables. Dans ce chapitre, nous nous concentrons sur la formulation d'un critère d'arrêt de l'algorithme de linéarisation du problème discrète (tel que la méthode de Newton) avec un critère d'arrêt du solveur algébrique de résolution du système linéarité (par exemple la méthode GMRes), de sort que les contributions des estimateurs d'erreur correspondant n'affectent plus la somme globale des estimateurs d'erreur de manière significative. Nous appliquons notre analyse sur des exemples réalistes d'ingénierie de réservoir pour confirmer qu'en général notre ajustement des critères d'arrêt apporte une économie significative (jusqu'au un ordre de magnitude en termes du nombre total des itérations du solveur algébrique), déjà sur des maillages fixes, et ceci sans perte notable de précision. Au chapitre 3 nous complétons le modèle décrit au chapitre 2 en considérant une condition non-isothermique pour l'écoulement a fin de traiter le modèle général d'écoulement polyphasique thermique dans les milieux poreux. Pour ce problème, nous développons des estimateurs d'erreur analogues a ceux du chapitre 2 pour lesquels nous établissons une borne supérieure d'erreur entièrement calculable, pour une norme duale du résidu complétée par un terme d'évaluation de la non-conformité. Nous montrons ensuite comment estimer séparément chaque composante d'erreur, ce qui nous permet d'ajuster les critères d'arrêt et d'équilibrer les contributions des différents estimateurs d'erreur : erreur d'approximation en temps, erreur d'approximation en espace, erreur de linéarisation et erreur du solveur algébrique. Ce chapitre se termine par une application des estimateurs au modèle d'huile morte. La preuve de l'efficacité de notre estimation a postiriori est egalement fournie. Finalement, au chapitre 4 nous considérons les procédés de récupération assistée d'huile. Plus précisément, nous étudions une technique de récupération thermique d'huile de type huile morte par injection de vapeur destinée a augmenter la mobilité des hydrocarbures. Dans ce chapitre, nous appliquons l'analyse a posteriori des chapitres 2 et 3, nous proposons une formule de quadrature pour simplifier l'évaluation des estimateurs, nous proposons un algorithme adaptatif de raffinement de maillages en espace et en temps basé sur les estimateurs et nous illustrons pas des essais numériques sur des exemples réalistes la performance de cette stratégie de raffinement. Notamment, des gains significatifs sont réalisés en terme du nombre de mailles nécessaires pour la simulation sur des exemples en dimension trois.

Cette thèse est centrée autour de l'étude théorique et de l'analyse numérique des équations paraboliques non linéaires avec divers conditions aux limites. La première partie est consacrée aux équations paraboliques dégénérées mêlant des phénomènes non-linéaires de diffusion et de transport. Nous définissons des notions de solutions entropiques adaptées pour chacune des conditions aux limites (flux nul, Robin, Dirichlet). La difficulté principale dans l'étude de ces problèmes est due au manque de régularité du flux pariétal pour traiter les termes de bords. Ceci pose un problème pour la preuve d'unicité. Pour y remédier, nous tirons profit du fait que ces résultats de régularités sur le bord sont plus faciles à obtenir pour le problème stationnaire et particulièrement en dimension un d'espace. Ainsi par la méthode de comparaison "fort-faible" nous arrivons à déduire l'unicité avec le choix d'une fonction test non symétrique et en utilisant la théorie des semi-groupes non linéaires. L'existence de solution se démontre en deux étapes, combinant la méthode de régularisation parabolique et les approximations de Galerkin. Nous développons ensuite une approche directe en construisant des solutions approchées par un schéma de volumes finis implicite en temps. Dans les deux cas, on combine les estimations dans les espaces fonctionnels bien choisis avec des arguments de compacité faible ou forte et diverses astuces permettant de passer à la limite dans des termes non linéaires. Notamment, nous introduisons une nouvelle notion de solution appelée solution processus intégrale dont l'objectif, dans le cadre de notre étude, est de pallier à la difficulté de prouver la convergence vers une solution entropique d'un schéma volumes finis pour le problème de flux nul au bord. La deuxième partie de cette thèse traite d'un problème à frontière libre décrivant la propagation d'un front de combustion et l'évolution de la température dans un milieu hétérogène. Il s'agit d'un système d'équations couplées constitué de l'équation de la chaleur bidimensionnelle et d'une équation de type Hamilton-Jacobi. L'objectif de cette partie est de construire un schéma numérique pour ce problème en combinant des discrétisations du type éléments finis avec les différences finies. Ceci nous permet notamment de vérifier la convergence de la solution numérique vers une solution onde pour un temps long. Dans un premier temps, nous nous intéressons à l'étude d'un problème unidimensionnel. Très vite, nous nous heurtons à un problème de stabilité du schéma. Cela est dû au problème de prise en compte de la condition de Neumann au bord. Par une technique de changement d'inconnue et d'approximation nous remédions à ce problème. Ensuite, nous adaptons cette technique pour la résolution du problème bidimensionnel. A l'aide d'un changement de variables, nous obtenons un domaine fixe facile pour la discrétisation. La monotonie du schéma obtenu est prouvée sous une hypothèse supplémentaire de propagation monotone qui exige que la frontière libre se déplace dans les directions d'un cône prescrit à l'avance.

Cette thèse porte sur l’analyse thermique et énergétique d’une fondation géothermique ventilée. A l’instar des échangeurs air-sol classiques (EAHE), celle-ci permet de rafraichir ou préchauffer selon la saison l’air destiné au renouvellement sanitaire des bâtiments. Face aux contraintes de rationalisation des consommations et aux exigences de confort thermique croissantes, ces systèmes passifs apparaissent comme étant prometteurs. Le principe de cette fondation est simple et similaire à celui des EAHE : faire circuler de l’air dans une conduite enterrée dans le sol (un à trois mètres) pour qu’il bénéficie - via convection - de l’inertie thermique du sol. La différence réside dans le fait que le canal dans lequel circule l’air n’est pas un tube en PVC ou aluminium mais fait partie intégrante de la structure du bâtiment, à savoir la fondation en béton armé. Ceci présente comme avantage majeur le gain de place lié à l’espace requis pour l’enfouissement des tuyaux. D’un point de vue thermique, la fondation échange non seulement de la chaleur avec le sol exposé aux sollicitations météorologiques mais aussi, et simultanément, aux sollicitations venant du bâtiment. De plus, la profondeur de la fondation – imposée par des raisons structurelles et économiques – est moindre que pour un EAHE traditionnel. Additionné au fait que le béton est poreux, la présence d’humidité peut fortement influencer la performance thermique de la fondation. Le présent travail propose donc d’étudier le comportement thermique complexe de cette fondation par deux approches. La première est expérimentale : un EHPAD équipé de deux fondations a été lourdement instrumenté et des données ont été accumulées sur plus d’un an. L’autre est numérique : deux modèles validés par comparaison avec les données expérimentales ont été développés. Le premier a vocation d’outil de dimensionnement, l’autre de compréhension fine des phénomènes physiques et prends en compte les transferts couplés de chaleur et de masse. / This thesis deals with the thermal and energy analysis of a geothermal ventilated fonudation. Similarly to earth-to-air heat exchangers (EAHE) this foundation enables, according to the season, to preheat or to cool down the air for the hygienic air change. Considering the energy consumption constraints and the buildings users thermal comfort desire, these systems appears to be relevant. The principle of this foundation is simple: to force the air to circulate in a hollowed beam buried into the ground (1 to 3m depth) so that it takes advantage - via convection - to the thermal inertia of the ground. The difference lays on the fact that the channel is not a plastic or aluminium pipe but it a part of the building structure, namely the reinforced concrete foundation. This induces a significant space gain, usually devoted to the pipe burying. From a thermal point of view, the foundation exchanges heat with both the soil beneath the building, and with the soil exposed to the weather thermal loads. Furthermore, the depth - imposed by structural and economical purposees - is lower than that of traditional EAHE. In addition to the fact that concrete is a porous material, the humidity content may strongly influence the thermal performance of the foundation. The current work thus proposes to study the complex thermal behaviour of this foundation in two ways. The first is experimental: an retirement home equipped with two foundation has been intensively instrumented and data recorded over more than one year. The other is numerical: two models validated against the experimental data have been developed. The first is intended to be a designing tool, the second a tool to allow a fine comprehension of the physical phenomenon and take into account coupled heat and moisture transfers.