Ce mémoire présente un modèle paramétrique visant à illustrer l'effet du revenu sur la demande de biens ayant des caractéristiques environnementales supérieures (biens verts). Le modèle incorpore un paramètre reflétant la préférence environnementale et le revenu d'un consommateur, ainsi que deux biens substituts, l'un étant une version verte et l'autre, une version traditionnelle. En caractérisant la demande pour ces deux biens, on note un effet revenu dans le comportement de consommation à l'égard du bien vert. Même si un individu se préoccupe de l'environnement, ses choix de consommation sont limités par son niveau de revenu. Selon ce niveau, il peut consommer exclusivement du bien traditionnel, exclusivement du bien vert ou simultanément des deux. L'effet revenu observé avec le modèle va dans le même sens que les résultats obtenus empiriquement dans les travaux menés précédemment sur la demande de bien vert.

Ce travail a été réalisé dans le cadre du projet GRACY regroupant l’IEMN et d’autres laboratoires de recherche : CALISTO et IMS Bordeaux. Ce manuscrit fait état d’une synthèse exhaustive des études et avancées menées dans le cadre de ce travail de thèse au sein de l’Institut d’Électronique, de Microélectronique et de Nanotechnologie (IEMN) dans le groupe CARBON. Le principal axe de réflexion de ce travail repose sur la conception, la modélisation et la caractérisation des dispositifs actifs et passifs sur substrat souples et rigides en vue du développement de nouveaux composants et de circuits électroniques avec des critères de performances de plus en plus importants. Au cours de ce travail, l’accent a été principalement portée sur les étapes essentielles à la réalisation de circuit intégré en ondes millimétriques utilisant la technologie coplanaire en impression jet d’encre et les transistors à effet de champ à base de graphène (GFETs). Ce mémoire montre en particulier l’intérêt et les potentialités du graphène pour son intégration au sein des circuits électroniques. De plus, une attention particulière a été portée sur la modélisation et les techniques de caractérisations relatives aux dispositifs passifs sur substrat souple. Par conséquent, un banc de caractérisation de ces éléments sur substrat flexibles a été développé au cours de cette thèse afin de vérifier et consolider expérimentalement leurs comportements. / This work was carried out under the project involving GRACY IEMN and other research laboratories: CALISTO and IMS Bordeaux. This manuscript reports a comprehensive overview of studies and advanced conducted as part of this thesis in the Institute of Electronics, Microelectronics and Nanotechnology (IEMN) in CARBON group. The main reflection axis of this work is based on the design, modeling and characterization of active and passive devices on flexible and rigid substrates for the development of new components and electronic circuits with increasingly important performance criteria. During this work, the focus was mainly focused on the essential steps to achieving integrated circuit millimeter wave using coplanar technology by inkjet printing and field effect transistors based on graphene (GFETs). This memory in particular shows the importance and potential of graphene for integration into electronic circuits. In addition, special attention was paid on modeling and characterization techniques related to passive devices on flexible substrates. Therefore, a characterization bench of these elements on flexible substrate has been developed during this thesis to verify and consolidate their behavior experimentally.

Résumé : voir fichier joint

Cette étude développe une plate-forme mathématique allant de la création d’une base de connaissances à la mise en place d’un modèle adaptatif. Une nouvelle approche en ce qui concerne la modélisation du four rotatif du clinker a été nécessaire pour mener à bien notre étude. Ainsi, un modèle d’état du système aux paramètres repartis, à la base des phénomènes physico-chimiques, a été conçu à l’aide d’équations aux dérivées partielles. La structure du modèle est basée sur trois variables d’état qui peuvent décrire correctement les principaux phénomènes (diffusions, déplacements, transferts et échanges thermiques, etc.) qui se déroulent dans le four rotatif de clinker. En outre, ces variables d’état sont accessibles à travers des paramètres mesurables tels que les températures de la surface extérieure du four et de l’air ambiant. Les paramètres du modèle proposé sont des fonctions de trois variables d’état, qui tout en remplaçant les paramètres mal connus et inaccessibles, décrivent les phénomènes physicochimiques qui caractérisent la complexité du système et sont définis comme des fonctions opératoires. Cette approche permet non seulement de minimiser les erreurs (hypothèse de simplification) mais aussi permet un accès facile et aisé des paramètres. Nous avons élaboré une procédure d’identification basée sur des analyses phénoménologique et dimensionnelle. Ainsi, à partir de données issues d’un complément d’instrumentation associé à la base de connaissances, l’identification des fonctions opératoires (paramètres) a été effectuée à partir d’un état stationnaire du système qui représente notre condition initiale. Cette opération est importante dans la mesure où elle donne aux paramètres du modèle les propriétés de la dynamique du système. Ces résultats sont corroborés par les estimations des variables d’état dont les valeurs et évolutions sont acceptables. Après l’évaluation de l’influence des fonctions opératoires sur les variables d’état, une modélisation des paramètres (fonctions opératoires) du modèle a été proposée à la base d’une décomposition physique à la lumière des connaissances fournies dans la littérature. Une procédure d’ajustement a été également développée afin de pouvoir fournir des variables d’état optimales. Ainsi, l’erreur peut être compensée par ajustement des fonctions opératoires. Par conséquent, à travers la rapidité de la procédure d’ajustement et l’adaptation des paramètres du modèle dans un horizon de temps bien déterminé, le modèle renferme les caractéristiques d’un modèle adaptatif qui pourra être un support essentiel dans la commande automatique du système. / This study develops a mathematical platform going from the establishment of a knowledge database to the setting up of an adaptive model. This has required a new approach of modeling of the clinker rotary kiln (CRK). Thus, a state model of distributed parameter systems, based on physico-chemical phenomena, was designed using partial differential equations. The model structure is based on three state variables which are: the gas, clinker temperatures and the clinker mass distributions, and are elaborated with the help of heat, pressure and mass balance equations. The model parameters are defined by the functions of three state variables. Moreover, the resulting state model, decomposed into five phenomenological zones of CRK, is used as a first step to define a set of Operating Functions (OFs). These OFs has also been decomposed into longitudinal distribution of CRK to replace the constant, unknown or unmeasured parameters. We develop an identification procedure based on phenomenological and dimensional analysis where the identification of operational functions (model parameters) was performed from a stationary state of the CRK. Once the restores state variables have been evaluated, the desired input (which is treated as the control of the CRK) can be more easily found by the proposed model than by simple trial and error. Moreover, the fact that the computation time, to estimate-calibrate the OFs above-mentioned, is very short, then this dynamic computation works faster than real-time. In summary, the cooperation and coordination in real-time between industrial computers and the CRK allows for an adaptable model, where each specific set of the OFs must be analyzed by its accuracy.

Cette thèse étudie les forces potentielles des mécanismes parallèles plans à deux degrés de liberté équipés d'embrayages de sécurité (limiteur de couple). Les forces potentielles sont étudiées sur la base des matrices jacobienne. La force maximale qui peut être appliquée à l'effecteur en fonction des limiteurs de couple ainsi que la force maximale isotrope sont déterminées. Le rapport entre ces deux forces est appelé l'efficacité de la force et peut être considéré ; comme un indice de performance. Enfin, les résultats numériques proposés donnent un aperçu sur la conception de robots coopératifs reposant sur des architectures parallèles. En isolant chaque lien, les modèles dynamiques approximatifs sont obtenus à partir de l'approche Newton-Euler et des équations de Lagrange pour du tripteron et du quadrupteron. La plage de l'accélération de l'effecteur et de la force externe autorisée peut être trouvée pour une plage donnée de forces d'actionnement. / This thesis investigates the force capabilities of two-degree-of-freedom planar parallel mechanisms that are equipped with safety clutches (torque limiters). The force capabilities are studied based on the Jacobian matrices. The maximum force that can be applied at the end-effector for given torque limits (safety index) is determined together with the maximum isotropic force that can be produced. The ratio between these two forces, referred to as the force effectiveness, can be considered as a performance index. Finally, some numerical results are proposed which can provide insight into the design of cooperation robots based on parallel architectures. Considering each link and slider system as a single body, approximate dynamic models are derived based on the Newton-Euler approach and Lagrange equations for the tripteron and the quadrupteron. The acceleration range or the external force range of the end-effector are determined and given as a safety consideration with the dynamic models.

La considération adéquate des différentes sources d’incertitude est un aspect crucial de la prévision hydrométéorologique. La prévision d’ensemble, en fournissant des informations sur la probabilité d’occurrence des sorties du modèle, représente une alternative séduisante à la prévision déterministe traditionnelle. De plus, elle permet de faire face aux différentes sources d’incertitude qui se trouvent le long de la chaîne de modélisation hydrométéorologique en générant des ensembles là où ces incertitudes se situent. Le principal objectif de cette thèse est d’identifier un système qui soit capable d’appréhender les trois sources principales d’incertitude que sont la structure du modèle hydrologique, ses conditions initiales et le forçage météorologique, dans le but de fournir une prévision qui soit à la fois précise, fiable et économiquement attractive. L’accent est mis sur la cohérence avec laquelle les différentes incertitudes doivent être quantifiées et réduites. Notamment, celles-ci doivent être considérées explicitement avec une approche cohésive qui fasse en sorte que chacune d’entre elles soit traitée adéquatement, intégralement et sans redondance dans l’action des divers outils qui composent le système. Afin de répondre à cette attente, plusieurs sous-objectifs sont définis. Le premier se penche sur l’approche multimodèle pour évaluer ses bénéfices dans un contexte opérationnel. Dans un second temps, dans le but d’identifier une implémentation optimale du filtre d’ensemble de Kalman, différents aspects du filtre qui conditionnent ses performances sont étudiés en détail. L’étape suivante rassemble les connaissances acquises lors des deux premiers objectifs en réunissant leurs atouts et en y incluant la prévision météorologique d’ensemble pour construire un système qui puisse fournir des prévisions à la fois précises et fiables. Il est attendu que ce système soit en mesure de prendre en compte les différentes sources d’incertitude de façon cohérente tout en fournissant un cadre de travail pour étudier la contribution des différents outils hydrométéorologiques et leurs interactions. Enfin, le dernier volet porte sur l’identification des relations entre les différents systèmes de prévisions précédemment créés, leur valeur économique et leur qualité de la prévision. La combinaison du filtre d’ensemble de Kalman, de l’approche multimodèle et de la prévision météorologique d’ensemble se révèle être plus performante qu’aucun des outils utilisés séparément, à la fois en précision et fiabilité et ceci en raison d’unemeilleure prise en compte de l’incertitude que permet leur action complémentaire. L’ensemble multimodèle, composé par 20 modèles hydrologiques sélectionnés pour leurs différences structurelles, est capable de minimiser l’incertitude liée à la structure et à la conceptualisation, grâce au rôle spécifique que jouent les modèles au sein de l’ensemble. Cette approche, même si utilisée seule, peut conduire à des résultats supérieurs à ceux d’un modèle semi-distribué utilisé de façon opérationnelle. L’identification de la configuration optimale du filtre d’ensemble de Kalman afin de réduire l’incertitude sur les conditions initiales est complexe, notamment en raison de l’identification parfois contre-intuitive des hyper-paramètres et des variables d’état qui doivent être mises à jour, mais également des performances qui varient grandement en fonction du modèle hydrologique. Cependant, le filtre reste un outil de première importance car il participe efficacement à la réduction de l’incertitude sur les conditions initiales et contribue de façon importante à la dispersion de l’ensemble prévisionnel. Il doit être malgré tout assisté par l’approche multimodèle et la prévision météorologique d’ensemble pour pouvoir maintenir une dispersion adéquate pour des horizons dépassant le court terme. Il est également démontré que les systèmes qui sont plus précis et plus fiables fournissent en général une meilleure valeur économique, même si cette relation n’est pas définie précisément. Les différentes incertitudes inhérentes à la prévision hydrométéorologique ne sont pas totalement éliminées, mais en les traitant avec des outils spécifiques et adaptés, il est possible de fournir une prévision d’ensemble qui soit à la fois précise, fiable et économiquement attractive. / A proper consideration of the different sources of uncertainty is a key point in hydrometeorological forecasting. Ensembles are an attractive alternative to traditional deterministic forecasts that provide information about the likelihood of the outcomes. Moreover, ensembles can be generated wherever a source of uncertainty takes place in the hydrometeorological modeling chain. The global objective of this thesis is to identify a system that is able to decipher the three main sources of uncertainty in modeling, i.e. the model structure, the hydrological model initial conditions and the meteorological forcing uncertainty, to provide accurate, reliable, and valuable forecast. The different uncertainties should be quantified and reduced in a coherent way, that is to say that they should be addressed explicitly with a cohesive approach that ensures to handle them adequately without redundancy in the action of the different tools that compose the system. This motivated several sub-objectives, the first one of which focusing on the multimodel approach to identify its benefits in an operational framework. Secondly, the implementation and the features of the Ensemble Kalman Filter (EnKF) are put under scrutiny to identify an optimal implementation. The next step reunites the knowledge of the two first goals by merging their strengths and by adding the meteorological ensembles to build a framework that issues accurate and reliable forecasts. This system is expected to decipher the main sources of uncertainty in a coherent way and provides a framework to study the contribution of the different tools and their interactions. Finally, the focus is set on the forecast economic value and provides an attempt to relate the different systems that have been built to economic value and forecast quality. It is found that the combination of the EnKF, the multimodel, and ensemble forcing, allows to issue forecasts that are accurate and nearly reliable. The combination of the three tools outperforms any other used separately and the uncertainties that were considered are deciphered thanks to their complementary actions. The 20 dissimilar models that compose the multimodel ensemble are able to minimize the uncertainty related to the model structure, thanks to the particular role they play in the ensemble. Such approach has the capacity to outperform more a complex semi-distributed model used operationally. To deal optimally with the initial condition uncertainty, the EnKF implementation may be complex to reach because of the unintuitive specification of hyper-parameters and the selection of the state variable to update, and its varying compatibility with hydrological model. Nonetheless, the filter is a powerful tool to reduce initial condition uncertainty and contributes largely to the predictive ensemble spread. However, it needs to be supported by a multimodel approach and ensemble meteorological forcing to maintain adequate ensemble dispersion for longer lead times. Finally, it is shown that systems that exhibit better accuracy and reliability have generally higher economic value, even if this relation is loosely defined. The different uncertainties inherent to the forecasting process may not be eliminated, nonetheless by explicitly accounting for them with dedicated and suitable tools, an accurate, reliable, and valuable predictive ensemble can be issued.

Thèse en cotutelle, doctorat en génie civil : Université Laval, Québec, Canada, Philosophiæ doctor (Ph. D.) et Université Toulouse III - Paul Sabatier, Toulouse, France / Un modèle de comportement mécanique des bétons fibrés est développé et implanté dans le logiciel éléments finis Cast3m. Il vient compléter le modèle de fissuration orthotrope du béton Fluendo3D, en ajoutant la capacité de traiter des matériaux à fibres courtes, cylindriques et rectilignes. Ce modèle, qui s'appuie sur des données d'essais et certains mécanismes modélisés issus de la littérature, permet d'apporter des éléments de compréhension du comportement de ce type de matériau, et notamment concernant le phénomène de multi-fissuration. La caractérisation des phénomènes mis enjeu durant l'extraction des fibres est le point de départ de cette étude. Les effets de l'inclinaison des fibres par rapport à la direction d'extraction sont pris en compte et interviennent dans le comportement du modèle qui présente la capacité d'utiliser des orientations préférentielles de fibres. Le phénomène de multi-fissuration est représenté grâce à une loi d'effet d'échelle de Weibull qui permet de tenir compte de la dispersion des résistances à la traction du béton et d'expliquer le développement de la multi-fissuration. Cette représentation permet d'obtenir des distributions d'ouvertures de fissures dans des macro-éléments finis, et apporte donc une nouvelle précision dans le calcul des ouvertures de fissures dans des structures de grandes dimensions. / A model for the mechanical behavior of fiber-reinforced concrete is developed and implemented in Cast3m finite element software. It completes the orthotropic concrete cracking model Fluendo3D, by adding the ability to treat materials with short, cylindrical and straight fibers. This model, which is based on tests results and assumptions of mechanisms modelled from litterature, provides elements of understanding of the behavior of this type of material, especially concerning the phenomenon of multi-cracking. The characterization of the phenomena involved during the extraction of the fibers is the starting point of this study. The effects of the inclination of the fibers with respect to the direction of extraction are taken into account and intervene in the behavior of the model which presents the capacity to use preferential orientations of fibers. The multi-cracking phenomenon is represented by a Weibull scaling law that allows to take into account the dispersion of the concrete tensile strengths and to explain the development of multi-cracking process. This representation makes it possible to obtain crack opening distributions in macro-finite elements, and thus brings a new precision in the calculation of cracks openings in large structures.

La détérioration des ressources en eau et la grande quantité d'eau polluée générée dans les sociétés industrialisées donnent une importance fondamentale aux procédés de traitement des eaux usées pour préserver les ressources, conformément à l'objectif 6 des 17 objectifs de développement durable des Nations Unies. Le rejet de nutriments tels que l'ammoniac par les eaux usées est un problème important, l'élimination de l'azote (N) est donc l'un des processus critiques de toute station de récupération des ressources en eau (StaRRE). L'objectif de ce projet de recherche doctoral est d'améliorer la compréhension des mécanismes d'élimination de l'azote dans le traitement biologique des eaux usées grâce à la modélisation, et d'optimiser les StaRRE existantes pour réduire la consommation d'énergie et de ressources. Dans ce cadre, 3 études différentes ont été réalisées. Tout d'abord, un modèle de décanteur réactif unidimensionnel a été développé. Celui-ci prédit le comportement de décantation de boues à des concentrations élevées de boues ainsi que les conversions biocinétiques dans le processus de décantation secondaire (DS). Il a été constaté qu'une description précise des réactions biocinétiques dans la DS impose des défis de calibration élevés pour le modèle de décantation, car ce dernier doit capturer les profils de concentration complets de la biomasse active dans la couverture de boues. Le modèle calibré a pu prédire avec précision les profils de concentration des effluents et du lit de boues dans la DS. Le modèle développé peut être utilisé pour le contrôle et la simulation des StaRRE afin d'obtenir de meilleures prédictions des concentrations d'effluents et des boues de retour, et aussi de calculer correctement le bilan massique d'azote d'une StaRRE. Deuxièmement, un modèle à l'échelle de l'usine a été mis en place pour un système de pré-dénitrification conventionnel pour la StaRRE pilEAUte à l'échelle pilote. Une méthodologie de calibration du modèle par étapes a été adoptée en fusionnant les principaux protocoles de calibration de modèle, tout en mettant l'accent sur le modèle biocinétique. Le modèle de la StaRRE pilEAUte, y compris le décanteur réactif développé, a été calibré et validé pour simuler les variables de modèle sélectionnées, puis utilisé pour une analyse de scénarios plus approfondie de l'optimisation de la consommation d'énergie et des ressources. Les résultats de l'analyse des scénarios ont montré le potentiel d'optimisation du système conventionnel d'élimination d'azote grâce à la réduction de l'aération et du retour interne des nitrates. Ils ont également démontré que la dénitrification dans le décanteur secondaire peut avoir une contribution significative à la capacité globale d'élimination d'azote d'une StaRRE lorsque la liqueur mixte peut traverser le lit de boues. Troisièmement, l'étude visait à évaluer l'applicabilité des stratégies de commande continu et intermittent du rapport de l'ammoniac par rapport aux NOX-N (commande AvN) sur la StaRRE pilEAUte. Les stratégies de commande de l'aération par AvN sont appliquées en amont d'un réacteur de désammonification, qui est un processus d'élimination efficace d'azote avec un besoin de ressources réduit (en termes d'aération et carbone) par rapport aux systèmes conventionnels. Les deux stratégies de commande testées pourraient être réalisées grâce à une commande automatique. Cependant, le maintien du rapport AvN dans l'effluent à la valeur souhaitée (1) dépend fortement des conditions opérationnelles telles que les variations de l'affluent, le temps de rétention des boues et la fiabilité des capteurs. Même si la recherche est guidée par les études de StaRRE à l'échelle pilote, les méthodologies développées pour démontrer et modéliser les processus et les conditions opérationnelles économes en énergie et en ressources sont applicables et transférables à d'autres études de cas à plein échelle. / Deterioration of water resources and the large amount of polluted water generated in industrialized societies gives fundamental importance to waste water treatment processes to preserve resources in accordance with goal 6 of the 17 sustainable development goals of the United Nations. Discharge of nutrients such as ammonia with waste water is a significant issue, thus nitrogen (N) removal is one of the critical processes of any water resource recovery facilities (WRRF). The objective of this PhD research project was to improve the understanding of N-removal mechanisms in biological treatment of wastewater through modelling and to optimize existing WRRFs to reduce energy and resource consumption. Within this context, 3 different studies were carried out. First, a one dimensional reactive settler model was developed that predicts the settling behaviour at high sludge concentrations together with biokinetic conversions in the secondary settling process. It was found that an accurate description of biokinetic reactions in the SST puts high calibration requirements on the settling model as it must properly capture the full concentration profiles of active biomass in the sludge blanket. The calibrated model was able to accurately predict the effluent and sludge blanket concentration profiles in the SST. The developed model can be used for control and simulation of WRRFs for better predictions of SST effluent and underflow concentrations and also properly calculate the nitrogen mass balance of a WRRF. Second, a plant-wide model was set up for a conventional pre-denitrification system for the pilot-scale pilEAUte WRRF. A step-wise model calibration methodology was adopted by merging main existing model calibration protocols while placing emphasis on the biokinetic model. The pilEAUte model, including the developed reactive settler, was calibrated and validated to simulate the selected model variables and used for further scenario analysis for energy and resource optimization. The scenario analysis results showed the optimization potential of conventional N removal systems through application of reduced aeration and internal nitrate recycling. It also demonstrated that denitrification in the secondary settler can contribute significantly to the overall N removal capacity of the WRRF when mixed liquor can pass through the sludge blanket. Third, it was aimed to evaluate the applicability of continuous and intermittent Ammonia vs NOₓ-N (AvN) control strategies on the pilEAUte WRRF. The AvN aeration control strategies are applied prior to a deammonification stage which is a short-cut N removal process with reduced resource (aeration and carbon) requirements in comparison to conventional systems. Both strategies could be achieved through automatic control. However, keeping the AvN ratio in the effluent on the desired value highly depends on operational conditions such asinfluent variations, sludge retention time and the sensor's measurement reliability.

En forêt boréale, le flux total de CO₂ du sol (F[indice CO2]) est le deuxième flux de carbone le plus important après la photosynthèse et peut représenter de 48 à 71% de la respiration totale de l'écosystème. Le F[indice CO2] dépend de l'état du sol à petite échelle, c'est-à-dire de sa température, de sa teneur en eau, de la quantité et de la qualité de la matière organique et de l'activité racinaire et microbienne. D'autres facteurs, appelés « facteurs externes », tels quel a topographie, les événements météorologiques et la hauteur de végétation influencent l'état du sol et donc indirectement le F[indice CO2]. Prédire le F[indice CO2] est complexe en raison du nombre important de facteurs à prendre en compte et des différents patrons de variation temporelle à modéliser. L'objectif de l'étude était de quantifier l'effet des facteurs externes susceptibles d'influencer la variabilité saisonnière, intrasaisonnière et spatiale des F[indice CO2] dans la forêt boréale. Des mesures de F[indice CO2] furent prises à la Forêt Montmorency (Québec, Canada) entre juin et octobre 2020. Un modèle de type Random forest regression fut utilisé pour prédire le F[indice CO2] en fonction de facteurs externes liés à la topographie, à la végétation et à la météo récente. Le modèle obtenu a permis d'expliquer 79 % de la variation des mesures de F[indice CO2]. Nos résultats ont démontré l'importance de la variation spatiale et intrasaisonnière (44 %) dans les mesures de F[indice CO2] par rapport à la variation saisonnière (35 %), ce qui a des implications pour la mesure et la modélisation de F[indice CO2]. La hauteur des arbres et la température de l'air furent les deux facteurs les plus importants pour prédire la valeur du F[indice CO2]. La méthodologie proposée permettrait de simuler le F[indice CO2] à partir de facteurs externes qui peuvent être obtenus à grande échelle à partir de relevés LiDAR et des données climatiques. / In the boreal forest, total soil CO₂ efflux (F[indice CO2]) is the second most important carbon flux after photosynthesis and can account for 48-71% of total ecosystem respiration. F[indice CO2]depends on the small-scale condition of the soil, i.e., its temperature, water content, quantity and quality of organic matter and root and microbial activity. Other factors, called "external factors", such as topography, weather events and vegetation height influence the soil condition and thus indirectly F[indice CO2]. Predicting F[indice CO2] is complex due to the large number of factors to be considered and the different patterns of temporal variation to be modeled. The objective of the study was to quantify the effect of external factors that may influence the seasonal, intra-seasonal and spatial variability of F[indice CO2] in the boreal forest. F[indice CO2] measurements were taken in the Montmorency Forest (Quebec, Canada) between June and October 2020. A random forest regression model was used to predict F[indice CO2] as a function of external factors such as topography, vegetation and recent weather. The developed model was able to explain 79% of the variation in F[indice CO2] measurements. Our results demonstrated the importance of spatial and intra-seasonal variation (44%) in F[indice CO2] measurements compared to seasonal variation (35%), which has implications for F[indice CO2] measurement and modeling. Tree height and air temperature were the two most important factors in predicting F[indice CO2] value. The proposed methodology would allow F[indice CO2] to be predicted from external factors that can be obtained on a large scale from LiDAR surveys and climate data.

Cette thèse traite du problème d’optimisation de formes dans le contexte de la mécanique des solides en contact. Le modèle physique considéré est celui de solides linéaires élastiques en petites déformations, en contact (glissant ou avec frottement de Tresca) avec un corps rigide. Les formulations mathématiques étudiées sont deux versions régularisées de l’inéquation variationnelle décrivant le système d’origine : la formulation pénalisée et la formulation Lagrangien augmenté. Comme ces deux formulations présentent des non différentiabilités, nous proposons une approche par dérivées directionnelles pour obtenir les dérivées de forme associées. En particulier, pour chacune des formulations, nous exprimons des conditions suffisantes pour que la solution soit dérivable par rapport à la forme. Ceci nous permet de construire un algorithme d’optimisation topologique de type gradient, s’appuyant sur les dérivées obtenues et une représentation des formes par des ensembles de niveau (level-sets). L’algorithme bénéficie en outre d’une technique de découpage de maillage qui permet d’obtenir une représentation explicite de la forme à chaque itération, et ainsi d’appliquer fortement les conditions aux limites sur la zone de contact. Après avoir détaillé les différentes étapes de la méthode, nous présentons des résultats numériques en deux et trois dimensions pour en tester la validité. / This thesis deals with shape optimization for contact mechanics. More specifically, the linear elasticity model is considered under the small deformations hypothesis, and the elastic bodyis assumed to be in contact (sliding or with Tresca friction) with a rigid foundation. The mathematical formulations studied are two regularized versions of the original variational inequality: the penalty formulation and the augmented Lagrangian formulation. In order to get the shape derivatives associated to those two non-differentiable formulations, we suggest an approach based on directional derivatives. Especially, we derive sufficient conditions for the solution to be shape differentiable. This allows to develop a gradient-based topology optimization algorithm, built on these derivatives and a level-set representation of shapes. The algorithm also benefits from a mesh-cutting technique, which gives an explicit representation of the shape at each iteration, and enables to apply the boundary conditions strongly on the contact zone. The different steps of the method are detailed. Then, to validate the approach,some numerical results on two-dimensional and three-dimensional benchmarks are presented.