Cette thèse est consacrée à l'amélioration des capacités de prédiction des émissions d'oxydes d'azote (NO et NO2) des foyers de combustion aéronautiques. Les travaux, exclusivement numériques, consistent d'abord dans une étude de la cinétique chimique responsable des émissions polluantes. Cetteétude conduit à l'écriture d'un modèle, nommé NOMANI (pour Nitrogen Oxide emission model with one-dimensional MANIfold), basé sur l'approche PCM-FPI (pour Presumed Conditional Moments - Flame Prolongation of ILDM) avec une variable de progrès additionnelle afin calculer l'avancement de la chimie azotée une fois la chimie carbonée à l'équilibre. Différentes validations sur des configurations laminaires simples puis des flammes de laboratoire de Sandia sont présentées. Les résultats en terme de structure de flamme et d'émission de monoxyde d'azote sont confrontés aux mesures expérimentales. Le dernier volet de ces travaux, disponible uniquement dans la version confidentielle du manuscrit, consiste dans le développement d'un modèle de prédiction de polluants associé au modèle TF-LES (pour Thickening Flame for Large Eddy Simulation). Le modèle développé est ensuite appliqué à des calculs d'une chambre de combustion aéronautique.

Anticiper l’évolution de la microstructure d’un matériau en condition d’usage est d’une importance cruciale pour l’industrie. Cette maîtrise du vieillissement nécessite une compréhension claire des mécanismes sous-jacents, qui agissent sur une large gamme d’échelles spatiales et temporelles. Dans cette optique, ce travail de thèse a choisi d’appliquer la méthode de champ de phase qui, en raison du saut d’échelle qu’elle réalise naturellement, est un outil intensivement employé dans le domaine des matériaux, pour prédire l’évolution en temps long de la microstructure. L’enjeu de l’étude a été d’étendre cette méthode à un système porté loin de l’équilibre thermodynamique, en particulier en présence d’irradiation. Nous avons ainsi adopté le formalisme du mélange ionique, introduit par Gras-Marti pour décrire le mélange balistique au sein d’une cascade de déplacements. Par l’utilisation conjointe d’un schéma numérique et d’une approche analytique, il nous a été possible d’établir le diagramme de phase générique d’un matériau irradié. Nous avons ensuite étudié le vieillissement de l’alliage binaire test AgCu sous irradiation, par l’utilisation conjointe de la méthode du champ de phase et d’approches atomistiques, dans une démarche multi-échelles. En fixant les paramètres de contrôle que sont le flux d’irradiation et la température, il nous a ainsi été possible de prédire la taille,la concentration ainsi que la distribution spatiale des nodules de cuivre produits sous irradiation dans cet alliage. La connaissance de ces informations a permis de simuler un diagramme de diffraction en incidence rasante, directement comparable aux diagrammes expérimentaux. / It is of dramatic matter for industry to be able to predict the evolution of material microstructure under working conditions. This requires a clear understanding of the underlying mechanisms, which act on numerous space and time scales. Because it intrinsically performs a scale jump, we chose to use a phase field approach, which is widely used amidst the condensed matter community to study the aging of materials. The first challenge of this work was to extend this formalism beyond its thermodynamic scope and embrace the case of far from equilibrium systems when subjected to irradiation. For that purpose, we adopted the model of ion mixing, developed by Gras Marti to account for ballistic exchanges within a displacements cascade. Based on a numerical scheme and ananalytical method, we were able to describe the generic microstructure signature for materials under irradiation.We then applied this formalism to the particular case of the immiscible binary alloy AgCu.With the joined use of the phase field approach and atomistic methods, we managed to predict how the temperature and the irradiation flux tailor the main microstructure features such as the size, the concentration and the distribution of copper precipitates. This eventually allowed us to simulate a diffraction pattern in grazing incidence, which is directly comparable to experimental ones.

Cette thèse se concentre sur le développement d'une méthode précise et efficace pour la simulation des grandes échelles (LES) des écoulements turbulents. Une approche de la LES basée sur la méthode variationnelle multi-échelles (VMS) est considérée. La VMS applique aux équations de la dynamique des fluides une séparation d'échelles a priori sans recours à des hypothèses sur les conditions aux limites ou sur l'uniformité du maillage. Afin d'assurer effectivement une séparation d'échelles dans l'espace des nombres d'onde associé, nous choisissons d'utiliser les ondelettes de deuxième génération (SGW), une base polynomiale qui présente des propriétés de localisation spatiale-fréquence optimales. A partir de la séparation d'échelles ainsi réalisée, l'action du modèle sous-maille est limitée à un intervalle de nombres d'onde proche de la coupure spectrale. Cette approche VMS-LES basée sur les ondelettes est désignée par WAVVMS-LES. Elle est incorporée dans un solveur d'ordre élevé pour la simulation des écoulements incompressibles sur la base d'une méthode de Galerkin discontinue (DG-FEM) stabilisée pour la pression. La méthode est évaluée par réalisation de LES sur des maillages fortement sous-résolus pour le cas test du tourbillon de Taylor-Green 3D à deux nombres de Reynolds différents. / This thesis focuses on the development of an accurate and efficient method for performing Large-Eddy Simulation (LES) of turbulent flows. An LES approach based upon the Variational Multiscale (VMS) method is considered. VMS produces an a priori scale-separation of the governing equations, in a manner which makes no assumptions on the boundary conditions and mesh uniformity. In order to ensure that scale-separation in wavenumber is achieved, we have chosen to make use of the Second Generation Wavelets (SGW), a polynomial basis which exhibits optimal space-frequency localisation properties. Once scale-separation has been achieved, the action of the subgrid model is restricted to the wavenumber band closest to the cutoff. We call this approach wavelet-based VMS-LES (WAV-VMS-LES). This approach has been incorporated within the framework of a high-order incompressible flow solver based upon pressure-stabilised discontinuous Galerkin FEM (DG-FEM). The method has been assessed by performing highly under-resolved LES upon the 3D Taylor-Green Vortex test case at two different Reynolds numbers.

En 1988, le Bassin du fleuve Sénégal au Mali a vu la mise en service d’un barrage par l’Organisation pour la Mise en Valeur du fleuve Sénégal (OMVS). L’objectif principale est la production d’énergie électrique partagée entre : le Mali (52%), le Sénégal (33%) et la Mauritanie (15%). Mais l’apparition de cet ouvrage constitue une perturbation profonde dans l’organisation et le fonctionnement des espaces riverains. Les territoires autour du Bafing, affluent où se localise le barrage, se trouvent ainsi cloisonnés en une partie amont et aval du lac de retenue. Ce dernier a insufflé une nouvelle dynamique spatiale avec l’installation récente de 25 campements de pêche.Cette nouvelle économie constitue un facteur d’attraction de pêcheurs professionnels venus des régions du Centre du Mali, particulièrement de Mopti et de Ségou. Le lac devient ainsi une immense réserve de poissons avec des tailles plus importantes par rapport à ceux capturés dans les affluents du Bafing, du Bakoye, ou même du fleuve Sénégal. En amont du barrage, la pêche devient la première activité économique poussant des jeunes agriculteurs et éleveurs «autochtones» à une reconversion professionnelle pas toujours aboutie. Les campements de pêcheurs occupent un ancien site des villages déplacés lors de la construction du barrage. Trente-trois villages sont actuellement réinstallés en aval du barrage dans le finage d’autres hameaux préexistants. Cette cohabitation bouleverse l’occupation de l’espace et entraine des tensions foncières. Certains sites comme Manantali à 5 km du barrage en sont les grands bénéficiaires. Ce village s’est transformé en un véritable « centre-rural », en accueillant les cadres et les ouvriers qualifiés et toute la main d’œuvre venus du Mali voire de l’étranger. Cet afflux d’habitants urbanisés a profondément changé la configuration du site et l’a surtout fortement ségrégué. Les bureaux de la société d’exploitation au pied du barrage, le vieux village de Manantali, les cités des ouvriers et celle des cadres sont autant de témoins d’un espace urbain en devenir, fonctionnel et très inégalitaire. Le principal objet de cette étude est l’étude des impacts du barrage de Manantali sur l’organisation socio-spatiale et physique des territoires riverains. L’entrée par une lecture des représentations sociales, consensuelles comme conflictuelles, par les habitants et les différents acteurs, est privilégiée. L’analyse des données quantitatives et des différents discours identifie les expressions tant des effets environnementaux physiques que socioéconomiques suscités par l’ouvrage. Le concept de représentations sociales est posé de la façon suivante: «les représentations forment des codes mémorisés par le cerveau, mobilisables de façon consciente et se prêtant à de multiples utilisations mentales. Ces codes servent en particulier à décrypter notre environnement géographique, mais aussi à communiquer avec autrui, à rêver, imaginer, planifier et orienter nos conduites ou nos pratiques les plus diverses» (DI MEO, 2008). Cette étude d’impact se positionne donc le domaine de la géographie sociale.Au cœur de cette étude se trouvent les acteurs, les responsables politiques à différentes échelles, et surtout l’habitant ordinaire qui vit à l’ombre du barrage de Manantali. Cette notion d’ombre doit être comprise dans la polysémie des impacts de l’ouvrage, tout autant néfastes que bénéfique, et par rapport à son aire d’influence. Les principales zones d’étude considérées se situent à « l’ombre proche » de l’ouvrage, ou à l’échelle locale, dans un rayon de 50 kilomètres autour du barrage. Il s’agit de 8 villages déplacés et anciens. Mais également des 25 campements de pêche autour du lac de retenue. / In 1988 the Senegal River basin in Mali has witnessed the building of a dam by the Organization for the Development of the Senegal River (OMVS) in French. The main objective is the production of electric power shared between: Mali (52%), Senegal (33%) and Mauritania (15%). But the birth of this dam has deeply disrupted the organization and functioning of all the waterside’s areas. The territories around the Bafing, the river where the dam is localized, are now strictlydivided by the reservoir between an upstream and downstream portions. This artificial lake has created new types of spatial organization with the recent settlements of 25 fishing camps.This new economy has created a pull factor for professional fishermen coming from the central regions of Mali, particularly Mopti and Segou. The lake hence becomes a huge fish reserve with larger sizes compared to those caught in the Bafing, the Bakoye or even in the Senegal River. In the upstream areas of the dam, fishing has become the first business activityencouraging local young farmers and herdsmen towards an unlikely professional retraining.Fishermen settlements occupy former site of displaced villages due to the dam construction. Thirty-three villages have been relocated downstream of the dam in the lands of existing hamlets. This cohabitation disruptstraditional land use and leads to social strains.

La simulation numérique fiable des performances de turbines hydrauliques suppose : i) de pouvoir inclure dans les calculs RANS (Reynolds-Averaged Navier-Stokes) traditionnellement mis en œuvre l'effet des incertitudes qui existent en pratique sur les conditions d'entrée de l'écoulement; ii) de pouvoir faire appel à une stratégie de type SGE (Simulation des Grandes Echelles) pour améliorer la description des effets de la turbulence lorsque des écarts subsistent entre calculs RANS et résultats d'essai de référence même après prise en compte des incertitudes. Les présents travaux mettent en oeuvre une démarche non intrusive de quantification d'incertitude (NISP pour Non-Intrusive Spectral Projection) pour deux configurations d'intérêt pratique : un distributeur de turbine Francis avec débit et angle d'entrée incertains et un aspirateur de turbine bulbe avec conditions d'entrée (profils de vitesse,en particulier en proche paroi, et grandeurs turbulentes) incertaines. L'approche NISP est utilisée non seulement pour estimer la valeur moyenne et la variance de quantités d'intérêt mais également pour disposer d'une analyse de la variance qui permet d'identifier les incertitudes les plus influentes. Les simulations RANS, vérifiées par une démarche de convergence en maillage, ne permettent pas pour la plupart des configurations analysées d'expliquer les écarts calcul / expérience grâce à la prise en compte des incertitudes d'entrée.Nous mettons donc également en ouvre des simulations SGE en faisant appel à une stratégie originale d'évaluation de la qualité des maillages utilisés dans le cadre d'une démarche de vérification des calculs SGE. Pour une majorité des configurations analysées, la combinaison d'une stratégie SGE et d'une démarche de quantification des incertitudes permet de produire des résultats numériques fiables. La prise en compte des incertitudes d'entrée permet également de proposer une démarche d'optimisation robuste du distributeur de turbine Francis étudié. / The reliable numerical simulation of hydraulic turbines performance requires : i) to includeinto the conventional RANS computations the effect of the uncertainties existing in practiceon the inflow conditions; ii) to rely on a LES (Large Eddy Simulation) strategy to improve thedescription of turbulence effects when discrepancies between RANS computations and experimentskeep arising even though uncertainties are taken into account. The present workapplies a non-intrusive Uncertainty Quantification strategy (NISP for Non-Intrusive SpectralProjection) to two configurations of practical interest : a Francis turbine distributor, with uncertaininlet flow rate and angle, and a draft-tube of a bulb-type turbine with uncertain inflowconditions (velocity distributions, in particular close to the wall boundaries, and turbulentquantities). The NISP method is not only used to compute the mean value and variance ofquantities of interest, it is also applied to perform an analysis of the variance and identify inthis way the most influential uncertainties. The RANS simulations, verified through a gridconvergence approach, are such the discrepancies between computation and experimentcannot be explained by taking into account the inflow uncertainties for most of the configurationsunder study. Therefore, LES simulations are also performed and these simulations areverified using an original methodology for assessing the quality of the computational grids(since the grid-convergence concept is not relevant for LES). For most of the flows understudy, combining a SGE strategy with a UQ approach yields reliable numerical results. Takinginto account inflow uncertainties also allows to propose a robust optimization strategy forthe Francis turbine distributor under study.

Dans ce travail, des méthodes de diagnostics et des techniques de développement de modèles sous-maille sont proposées pour la simulation des grandes échelles (SGE) du mélange turbulent. Plusieurs modèles sous-maille issus de ces stratégies sont ainsi présentés pour illustrer ces méthodes.Le principe de la SGE est de résoudre les grandes échelles de l'écoulement responsables des transferts principaux et de modéliser l'action des petites échelles de l'écoulement sur les échelles résolues. Au cours de ce travail, nous nous sommes appuyés sur le classement des modèles sous-maille en deux catégories. Les modèles "fonctionnels" qui s'attachent à reproduire les transferts énergétiques entre les échelles résolues et les échelles modélisées et les modèles "structurels" qui cherchent à bien reproduire le terme sous-maille. Le premier enjeu important a été d'évaluer la performance des modèles sous-maille en prenant en compte leur comportement à la fois fonctionnel (capacité à reproduire les transferts d'énergie) et structurel (capacité à reproduire le terme sous-maille exact). Des diagnosctics des modèles sous-maille ont pu être conduits avec l'utilisation de la notion d'estimateur optimal ce qui permet de connaitre le potentiel d'amélioration structurelle des modèles. Ces principes ont dans un premier temps servi au développement d'une première famille de modèles sous-maille algébrique appelée DRGM pour "Dynamic Regularized Gradient Model". Cette famille de modèles s'appuie sur le diagnostic structurel des termes issus de la régularisation des modèles de la famille du gradient. D'après les tests menés, cette nouvelle famille de modèle structurel a de meilleures performances fonctionnelles et structurelles que les modèles de la famille du gradient. L'amélioration des performances fonctionnelles consiste à supprimer la prédiction excessive de transferts inverses d'énergie (backscatter) observés dans les modèles de la famille du gradient. Cela permet ainsi de supprimer le comportement instable classiquement observé pour cette famille de modèles. La suite de ce travail propose ensuite d'utiliser l'estimateur optimal directement comme modèle sous-maille. Comme l'estimateur optimal fournit le modèle ayant la meilleure performance structurelle pour un jeu de variables donné, nous avons recherché le jeu de variable optimisant cette performance. Puisque ce jeu comporte un nombre élevé de variables, nous avons utilisé les fonctions d'approximation de type réseaux de neurones pour estimer cet estimateur optimal. Ce travail a mené au nouveau modèle substitut ANNM pour "Artificial Neural Network Model". Ces fonctions de substitution se construisent à partir de bases de données servant à émuler les termes exacts nécessaire à la détermination de l'estimateur optimal. Les tests de ce modèle ont montré qu'il avait de très bonnes perfomances pour des configurations de simulation peu éloignées de la base de données servant à son apprentissage, mais qu'il pouvait manquer d'universalité. Pour lever ce dernier verrou, nous avons proposé une utilisation hybride des modèles algébriques et des modèles de substitution à base de réseaux de neurones. La base de cette nouvelle famille de modèles ACM pour "Adaptative Coefficient Model" s'appuie sur les décompositions vectorielles et tensorielles des termes sous-maille exacts. Ces décompositions nécessitent le calcul de coefficients dynamiques qui sont modélisés par les réseaux de neurones. Ces réseaux bénéficient d'une méthode d'apprentissage permettant d'optimiser directement les performances structurelles et fonctionnelles des modèles ACM. Ces modèles hybrides allient l'universalité des modèles algébriques avec la performance élevée mais spécialisée des fonctions de substitution. Le résultat conduit à des modèles plus universels que l'ANNM. / This work develops subgrid model techniques and proposes methods of diagnosis for Large Eddy Simulation (LES) of turbulent mixing.Several models from these strategies are thus presented to illustrate these methods.The principle of LES is to solve the largest scales of the turbulent flow responsible for major transfers and to model the action of small scales of flowon the resolved scales. Formally, this operation leads to filter equations describing turbulent mixing. Subgrid terms then appear and must bemodeled to close the equations. In this work, we rely on the classification of subgrid models into two categories. "Functional" models whichreproduces the energy transfers between the resolved scales and modeled scales and "Structural" models that seek to reproduce the exact subgrid termitself. The first major challenge is to evaluate the performance of subgrid models taking into account their functional behavior (ability to reproduce theenergy transfers) and structural behaviour (ability to reproduce the term subgrid exactly). Diagnostics of subgrid models have been enabled with theuse of the optimal estimator theory which allows the potential of structural improvement of the model to be evaluated.These methods were initially involved for the development of a first family of models called algebraic subgrid $DRGM$ for "Dynamic Regularized GradientModel". This family of models is based on the structural diagnostic of terms given by the regularization of the gradient model family.According to the tests performed, this new structural model's family has better functional and structural performance than original model's family of thegradient. The improved functional performance is due to the vanishing of inverse energy transfer (backscatter) observed in models of thegradient family. This allows the removal of the unstable behavior typically observed for this family of models.In this work, we then propose the use of the optimal estimator directly as a subgrid scale model. Since the optimal estimator provides the modelwith the best structural performance for a given set of variables, we looked for the set of variables which optimize that performance. Since this set of variablesis large, we use surrogate functions of artificial neural networks type to estimate the optimal estimator. This leads to the "Artificial Neural Network Model"(ANNM). These alternative functions are built from databases in order to emulate the exact terms needed to determine the optimal estimator. The tests of this modelshow that he it has very good performance for simulation configurations not very far from its database used for learning, so these findings may fail thetest of universality.To overcome this difficulty, we propose a hybrid method using an algebraic model and a surrogate model based on artificial neural networks. Thebasis of this new model family $ACM$ for "Adaptive Coefficient Model" is based on vector and tensor decomposition of the exact subgrid terms. Thesedecompositions require the calculation of dynamic coefficients which are modeled by artificial neural networks. These networks have a learning method designedto directlyoptimize the structural and functional performances of $ACM$. These hybrids models combine the universality of algebraic model with high performance butvery specialized performance of surrogate models. The result give models which are more universal than ANNM.

L'importance du phénomène entrepreneurial, et plus largement de la place des TPE, dans le développement durable des économies, conduit à s'interroger sur les éléments qui peuvent amener un individu à créer une entreprise. Les recherches indiquent clairement que certaines personnes réussissent mieux que d’autres à identifier et à exploiter des opportunités d’affaires qui se présentent. De nombreuses explications sont avancées : la perception que les individus peuvent avoir des ressources disponibles, l’accessibilité des ressources, la valeur attendue de l’exploitation des ressources, les coûts d’opportunité associés à différents choix, l’optimisme, l’orientation vers l’action et une grande variété de dispositions personnelles comme la tolérance à l’ambiguïté ou le besoin d’accomplissement. De ce foisonnement d’études, deux principaux ensembles d’explications permettent de comprendre l’acte d’entreprendre : d’une part, la création d’entreprise répond à des logiques individuelles qui s’organisent autour du concept central de motivation ; d’autre part, les dynamiques personnelles s’inscrivent dans des environnements particuliers, qui peuvent faciliter ou parfois limiter les actions des porteurs de projet. La recherche que nous nous proposons de mener a pour but de montrer que la motivation entrepreneuriale a pour origine un ensemble de besoins dont l'intensité varie d'un entrepreneur à un autre. L'intensité du besoin s'explique par des antécédents personnels et professionnels. Le modèle d'affaires à la base du projet (innovation, reproduction, etc.) est en partie déterminé par les caractéristiques de la motivation entrepreneuriale. La thèse s'appuie sur une grande variété de références bibliographiques afin de repérer les besoins à l'origine de la motivation entrepreneuriale, les antécédents et l'impact de ces variables sur les caractéristiques du modèle d'affaires du porteur de projet. Notre intention est de tester les relations à partir d'un échantillon d'entrepreneurs présents au Sénégal et au Mali, et de confronter nos résultats aux études antérieures sur le sujet. / The importance of entrepreneurial phenomena and, more generally, of the place occupied by VSBs in the sustainable development of economies, leads us to consider the factors that lead an individual to create a business. Research clearly illustrates that some people are more successful than others when it comes to identifying and exploiting the business opportunities that present themselves. Many explanations have been put forward for this: the way individuals may perceive the resources available, the accessibility of these resources, the value expected from exploiting these resources, the opportunity costs related to the various choices made, optimism, the move towards action and the level of personal propensity towards factors such as tolerating ambiguity or the need for accomplishment. This profusion of studies has brought two main sets of explanations to light that help us better understand the act of creating business: first of all, creating a business is based on individual considerations, organized around the core concept of motivation; secondly, personal dynamics fall within specific environments, which may make the actions of project owners easier or, may occasionally limit them.The aim of the research that we propose to undertake is to illustrate that the origin of entrepreneurial motivation is a range of needs with an intensity that varies from one entrepreneur to another. This intensity is explained through personal and professional antecedents. The business model established as the base of the project (innovation, reproduction, etc.) is determined in part by entrepreneurial motivation characteristics. The argument is supported by a wide range of bibliographic references, used to identify the needs that underlie this entrepreneurial motivation, the antecedents and the impact of these variables on the specificities of the project owner’s business model. We wish to test the links (relationships) from a sample group of entrepreneurs located in Senegal and Mali, and to compare our results with previous studies on the subject.

Aujourd'hui, dans de nombreux champs d'applications, de plus en plus de problèmes d'ingénierie demandent d'avoir une évaluation précise et efficace de quantités d'intérêt.Très souvent, ces quantités dépendent de la solution d'une équation aux dérivées partielles (EDP) paramétrée où les paramètres -- physiques ou géométriques -- sont les entrées du modèle et les quantités d'intérêt -- valeurs moyennes -- en sont les sorties.Les techniques de réduction d'ordre, notamment la méthode des bases réduites qui est la méthode utilisée tout au long de ces travaux,permettent de répondre à ces demandes.Dans cette thèse nous nous intéressons à la mise en place d'un framework en C++, supportant le calcul parallèle, permettant d'appliquer la méthode des bases réduites à des problèmes multi-physiques non-linéaires tels queles problèmes de convection naturelle (couplage fluide-thermique), ou encore la modélisation d'aimants de type résistifs à hauts champs (nous nous limitons au couplage thermo-electrique) aboutissant à une étude sur la quantification d'incertitude.La méthode des bases réduites s'appuie naturellement sur une approximation obtenue via la discrétisation élément fini du problème à traiter. Pour cela nous utilisons la librairie de calcul Feel++, spécialisée dans la résolution d'EDPs.Nous nous intéressons également aux problèmes de type multi-échelles.La particularité de ces problèmes est de manipuler un ensemble de phénomènes mettant en jeu des échelles différentes, comme c'est le cas par exemple lorsque nous considérons un écoulement en milieu poreux.La méthode des éléments finis multi-échelles permet d'avoir le comportement "global", associé aux grandes échelles, de la solution du problème sans devoir le résoudre sur les petites échelles.Nous proposons une nouvelle construction des fonctions de base élément fini multi-échelles basée sur la méthode des bases réduites. / Today, in many fields of applications, more and more engineering problems require to have an accurate and efficient evaluation of quantities of interest.Often, these quantities depend on a partial differential equation (PDE) parameterized solution -- physical or geometrical -- are the model inputs and the quantities of interest -- average values ​​-- are the outputs.The order reduction techniques, including reduced basis method which is the method used throughout this work, can meet these demands.In this thesis, we focus on the establishment of a framework in C ++ supporting parallel computing, which applies the reduced basis method to nonlinear multiphysics problems such as problems with natural convection (fluid-thermal coupling) or the high field resistive magnet modeling (we limit ourselves to thermo-electric coupling) leading to a study on the uncertainty quantification.The reduced basis method naturally relies on an approximation obtained using the finite element discretization of the problem being treated. For this, we use the Feel ++ computation library specialized in PDE resolution.We are also interested by multiscale problems.The particularity of these problems is to manipulate a set of phenomena involving different scales, as this is the case for example when we consider a flow in porous media.The multiscale finite element method allows having a "global" behavior, linked with large scales, of the problem solution without solving it on small scales.We propose a new construction of multiscale finite element basis functions based on the reduced basis method.

Afin de répondre aux besoins du marché de l’acier, le procédé de cokéfaction du charbon doit s’adapter. Cependant, changer les paramètres de cokéfaction du charbon, tels que la température du four, le temps de cuisson ou la composition de la pâte à coke enfournée, peut générer un endommagement prématuré de la maçonnerie des fours. En effet, la transformation du charbon en coke s’accompagne d’une poussée sur les parois du four fortement dépendantes d’un grand nombre de paramètres. Afin d’anticiper ce problème, un projet européen nommé « Swelling PRessure In a coke oven, Transmission on oven walls and COnsequences on wall » a été mis en place. Cette thèse s’inscrit dans ce programme et vise à déterminer la poussée maximale pouvant être admise par un piédroit de cokerie lors de la pyrolyse du charbon. Pour modéliser ces structures composées de plusieurs centaines de milliers de briques, le point de vue macroscopique est le plus approprié. La maçonnerie est remplacée par un matériau homogène équivalent dont le comportement varie en fonction de l’état d’endommagement de la maçonnerie, ramené localement à un état d’ouverture des joints de mortier. Afin de détecter ces ouvertures, un critère de type Mohr-Coulomb en contraintes est utilisé. Il repose sur la comparaison des limites à rupture d’un sandwich brique-Mortier déterminé expérimentalement à haute température avec les contraintes mésoscopiques issues de la simulation. Un protocole expérimental novateur a été développé pour caractériser la tenue en traction du sandwich brique / mortier / brique jusqu’à 1000°C. Les limites à rupture issues de cet essai de traction directe ont été comparées à celles obtenues par des essais de fendage réalisés à l’université de Leoben. Les résultats sont concordants et confirment l’importance de l’état de surface avant maçonnage. Selon l’état d’endommagement considéré, les contraintes mésoscopiques sont obtenues grâce à un tenseur de localisation ou grâce à une sous-Structuration. Cette étape de sous-Structuration consiste à simuler localement une cellule à l’échelle mésoscopique en lui appliquant le champ de déplacement macroscopique obtenu grâce à la simulation. L’outil numérique a été validé par confrontation avec un cas test de référence. Pour finir, l’outil numérique développé a été utilisé pour caractériser l’influence de différents paramètres tels que la prise en compte de la thermique, la mise en compression de la structure…. Enfin, la simulation de cuissons sur des piédroits complets (sains ou initialement endommagés) a été réalisée. L’importance de l’endommagement initial est clairement soulignée par les résultats. Enfin, un nouveau modèle, appelé « deux carneaux avec poutres», est proposé pour réduire le coût de calcul. Plus complet que le modèle « deux carneaux » utilisé au CPM, il donne accès à de très bons résultats pour un coût nettement moindre que celui du piédroit complet avec homogénéisation et sous-Structuration. / To face coke and steel market requirements, the coking process has to be more flexible. Changing process parameters such as coking temperature, blend composition and cooking time can damaged coke oven battery heating wall. Indeed, the coking process generates a swelling pressure on wall which depends on a lot of parameters. To study this point, a European project named « Swelling PRessure In a coke oven, Transmission on oven walls and COnsequences on wall » has been set up. This work is a part of it and aims to determine the admissible pushing pressure for the coke oven heating walls to prevent crack formation. To model large masonries composed of numerous bricks, a mesoscopic point of view is more appropriate. Bricks and mortar are replaced by a Homogeneous Equivalent Material (HEM) whose behaviour depends on the joint state. In order to represent joint opening mechanism, a Mohr-Coulomb criterion in stress is used. This criterion compares the level of stress to the ultimate tensile or shear stress at mesoscopic scale. Ultimate stresses are obtained thanks to an experimental campaign using a new protocol developed at PRISME Laboratory. The brick-Mortar behaviour is experimentally characterised at high temperature (20°C to 1000°C). To validate the tensile test developed, a second experimental campaign using “wedge splitting tests” has been done at Leoben University. Results are similar and confirm the importance of the brick surface state. Depending on the initial damage of the structures, mesoscopic stresses are obtained by localization tensor or by sub-Modelling. The sub-Modelling step aims to simulate a local part of the masonry at the mesoscopic scale. This step aims to simulate with a mesoscopic point of view a local part of the global model. This numerical tool has been validated thanks to a literature test. Finally, the numerical tool has been used to characterise the influence of some parameters (thermal, force due to the cross tie rod,..). Next, the simulation of the whole coke oven heating wall has been performed (undamaged or initially damaged masonry). These FE simulations show the influence of initial damage on the final masonry damage. Finally, a two flues model with beams is proposed to take into account compression due to cross tie rod and to limit computational cost. It permits to obtain better results than the existing two flues model used at CPM with a lower cost compared to the whole coke oven heating wall model.

Le contexte de cette étude est l’optimisation de la combustion turbulente prémélangée de syngaz pour la production propre d’énergie. Un brûleur CH4/air de type bec Bunsen avec forçage turbulent multi-échelles produit par un système de trois grilles, est simulé numériquement par différentes techniques de simulation des grandes échelles (SGE), et les résultats sont comparés à l’expérience. On a développé et appliqué une formulation bas-Mach du solveur Navier-Stokes basé sur différents schémas numériques, allant des différences finies centrées d’ordre 4 à des versions avancées des schémas WENO d’ordre 5. La méthodologie est évaluée sur une série de cas-tests classiques (flamme laminaire 1D prémélangée, turbulence homogène et isotrope en auto-amortissement), et sur des simulations 2D de la flamme turbulente prémélangée expérimentale. Les SGE implicites (ILES), i.e. sans aucune modélisation sous-maille, et explicites avec le modèle de flamme épaissie et un modèle de plissement sous-maille nouvellement élaboré (TFLES), sont appliquées à la simulation 3D du brûleur expérimental. Les résultats montrent que l’approche TFLES avec un schéma d’ordre élevé à faible dissipation numérique prédit correctement la longueur de la flamme et la densité de surface de flamme. La SGE implicite avec un schéma WENO avancé produit une flamme trop courte mais réaliste à condition que la taille de la maille soit de l’ordre de l’épaisseur de flamme laminaire. La représentation des interactions flamme/turbulence est néanmoins très différente entre TFLES et ILES. / The context of this study is the optimization of premixed turbulent combustion of syngas for clean energy production. A Bunsen-type CH4/air turbulent premixed burner with a multi-scale grid generator is simulated with different Large Eddy Simulation (LES) strategies and compared to experimental results. A low-Mach formulation of a compressible Navier-Stokes solver based on different numerical methods, ranging from 4th order central finite difference to 5th order advanced WENO schemes, is developed and applied. Classical test cases (1D laminar premixed flame, decaying HIT), and 2D simulations of the turbulent premixed flame are performed to assess the numerical methodology. Implicit LES (ILES), i.e. LES without any explicit subgrid modeling, and explicit LES with the Thickened Flame model and subgrid scale flame wrinkling modelling (TFLES) are applied to simulate numerically the 3D experimental burner. Results show that TFLES with a high-order low dissipation scheme predicts quite well the experimental flame length and flame surface density. ILES with advanced WENO schemes produces a slightly shorter although realistic flame provided the grid spacing is of order of the laminar flame thickness. The representation of flame/turbulence interactions in TFLES and ILES are however quite different.