Thermo-mechanical failure of components in planar-type solid oxide fuel cells (SOFCs) is a major obstacle on the path to bringing this technology to commercial viability. The probability of material degradation and failure in SOFCs depends strongly on the local temperature gradients at the interfaces of different materials. Therefore, it is of paramount importance to accurately predict and manage the temperature fields within the stack, especially near the interfaces. In this work we consider three effects in detail. First, we analyze radiative heat transfer effects within the semi-transparent solid electrolyte and compared them to thermal conduction. We also, present the modeling approach for calculation of surface-to-surface exchange within the flow channels and from the stack to the environment. The simplifying assumptions are identified and their carefully justified range of applicability to the problem at hand is established. This allows thermal radiation effects to be properly included in overall thermal modeling efforts with the minimum computational expense requirement. Second, we developed a series of reduced-order models for the transient heating and cooling of a cell, leading to a framework for optimization of these processes. The optimal design is one that minimizes heating time while maintaining thermal gradients below an allowable threshold. To this end, we formulated reduced order models (validated by rigorous CFD simulations) that yield simple algebraic design rules for predicting maximum thermal gradients and heating time requirements. Several governing dimensionless parameters and time scales were identified that shed light on the essential physics of the process. Finally, an analysis was performed to assess the degree of local thermal non-equilibrium (LTNE) within porous SOFC electrodes, and through a simple scaling analysis we discovered the parameter that gives an estimate of the magnitude of LTNE effects. We conclude that because of efficient heat transfer between the solid and gas in the microscale pores of the electrodes, the temperature difference between gas and solid is often negligible. However, if local variations in current density are significant, the LTNE effects may become significant and should be considered.

No description available.

Un modèle constitutif complètement couplé est présenté pour l'analyse rigoureuse de la déformation, de l'écoulement de fluides et de transfert de chaleur dans les milieux poreux saturés à double porosité soumis à des chargements thermo-hydro-mécaniques, y compris ceux induisant un non-équilibre thermique local. La phase solide contient deux cavités distinctes: le bloc poreux et le réseau des fissures. Les équations de champs sont obtenues à partir des équations de conservation de la masse, du mouvement et de l'énergie et sont résolues par une approche par élément finis. Le modèle est utilisé pour deux types d'applications: la stabilité d'un puits de forage stimulée thermiquement pour la récupération de pétrole et l'extraction de chaleur dans un réservoir géothermique fracturé. Les différences substantielles, particulièrement de la contrainte effective, soulignent l'influence majeure de la double porosité et du non-équilibre thermique pour prédire le comportement des milieux fracturés. / A fully coupled constitutive model is presented for a rigorous analysis of deformation, hydraulic and heat flows in saturated dual porosity media subject to thermo-hydro-mechanical loadings including those able to cause local thermal non-equilibrium. The solid phase is assumed to contain two distinct cavities: the porous blocks and the fissure network. The governing equations are derived based on the equations of conservation of mass, momentum and energy. Solution to the governing equations is obtained numerically using the finite element approach. The capabilities of the model address two energy applications: the stability of a borehole in a thermally enhanced oil recovery context and the heat extraction of enhanced geothermal systems. Substantial differences, particularly in the effective stress response, highlight the major influence of the dual porosity model and the importance of the local thermal non-equilibrium assumption to predict the behaviour of fractured media.

Ce travail concerne l'utilisation de la technique de Fracturation Hydraulique (FH) pour exploiter l'énergie géothermique des réservoirs profonds de roches sèches chaudes (HDR). La fracturation hydraulique est réalisée par injection de fluides géothermiques dans des réservoirs partiellement fracturés de faible perméabilité. Les fluides à haute pression sont destinés à faire évoluer les fissures et leur connectivité. Les valeurs de débit/pression auxquelles les fluides géothermiques doivent être pompés, ainsi que le calendrier de pompage pour initier la fracturation hydraulique, dépendent principalement des conditions géostatiques existantes (contraintes géostatiques, pression fluide et température initiales de l'HDR) ainsi que des propriétés des fissures de l'HDR (longueur, épaisseur, densité et distribution directionnelle initiales moyennes de fissures). Tous ces éléments, en sus de leurs effets sur la stabilité des forages, sont analysés dans cette recherche. Des modèles de fracturation, qui sont capables de suivre l'évolution des fissures dans toutes les orientations spatiales possibles, sont utilisés pour obtenir le tenseur anisotrope de perméabilité. Ces modèles sont intégrés dans un code domestique d'éléments finis qui est développé pour résoudre des problèmes aux limites thermo-poroélastiques. Pour supprimer/diminuer les oscillations qui accompagnent les solutions paraboliques et/ou hyperboliques lors de la convection forcée, plusieurs techniques de stabilisation ont dû être implémentées. / The application of the Hydraulic Fracturing (HF) technology to exploit geothermal energy from Hot Dry Rocks (HDR) reservoirs is addressed. HF is achieved by extensively pumping geothermal fluids to already existing fractured HDR reservoirs of low permeability. High fluid pressures are expected to drive cracks to evolve and connect. The newly created burgeoning hydraulic conduits should supposedly enhance the permeability of the existing HDR reservoirs. The flow rate/pressure values at which geothermal fluids should be pumped, as well as the pumping schedule to initiate HF, depend primarily on the existing geostatic conditions (geostatic stresses, initial HDR pressure and temperature) as well as on HDR fracture properties (initial mean fracture length, mean fracture aperture, density and orientational distribution of fractures). While these components, in addition to their effects on borehole stability, are scrutinized in this research, focus is on the evolution during circulation processes of the fracture properties. A fracturing model that is capable of tracking fracture evolution in all possible spatial orientations is used to obtain the time course of the anisotropic permeability tensor. This evolving property is integrated into a domestic finite element code which is developed to solve thermo–poroelastic BVPs: emphasis is laid on the efficiency of the doublet flow technique where a fluid gains thermal energy by circulating through the HDR reservoir from the injection well to the production well. The spurious oscillations in the hyperbolic solutions of the approximated finite element approach that are commensal with the phenomenon of forced heat convection are healed/mitigated through several stabilization approaches.