Depuis quelques années, l’industrie québécoise de production de bois d’œuvre opère dans un contexte économique difficile expliqué par plusieurs facteurs. Parmi ceux-ci, mentionnons notamment une plus forte concurrence sur les marchés internationaux, l’accès limité au marché américain résultant d’un conflit commercial, un dollar canadien fort, la présence de matériaux de substitution et la rareté croissante de la matière première ainsi que la diminution de sa qualité. Parmi les pistes de solutions les plus souvent évoquées, on trouve la mise en place de réseaux logistiques plus efficaces. L’industrie possède des outils permettant d’estimer l’augmentation de la production associée à l’achat d’un équipement. Cependant, ces outils ne considèrent généralement qu’une seule étape du procédé de production. Or, il est fréquent que la modification d’une étape du procédé de production affecte les autres également. De plus, la rentabilité globale de l’investissement dépendra de quelle façon la stratégie de production et de vente de l’entreprise sera ajustée suite à l’investissement. Par conséquent, nous proposons un cadre décisionnel permettant d’analyser l’impact combiné d’investissements au sein de plusieurs étapes de production, tout en définissant simultanément la meilleure stratégie de production et de vente. Alors que de manière classique, ces décisions sont prises de manière indépendantes, nous préconisons plutôt l’utilisation d’un modèle mathématique afin d’intégrer ces niveaux de décision. Le cadre proposé spécifie les étapes qui devraient être suivies par le décideur pour utiliser ce modèle mathématique; il comprend également modèle de données et des interfaces-utilisateur supportant l’analyse du décideur. Afin de s’assurer de l’applicabilité du cadre décisionnel proposé, deux études de cas ont été réalisées en industrie. Elles ont permis de montrer les avantages du cadre décisionnel et l’importance d’analyser les plans tactiques lors de la prise de décision d’investissement. De plus, ces études ont permis de mettre en évidence certaines embuches pouvant être rencontrées lors de l’application de la démarche.

Ce travail de mémoire présente la modélisation numérique et l'analyse de l'écoulement dans un seul passage de la roue d'une turbine hydraulique axiale de type hélice à l'aide du logiciel ANS YS CFX 11.0. Une approche Reynolds Averaged Navier-Stokes (RANS), basé sur le modèle Shear Stress Transport (SST), a été utilisée pour modéliser numériquement l'écoulement. Les simulations prennent en considération le jeu à l'extrémité des aubes. Souvent, le jeu est omis dans les simulations numériques dans les roues de turbines hydrauliques. Cependant, il donne lieu à l'écoulement de fuite interne et autres phénomènes turbulents, tel que le tourbillon de bout d'aube. La comparaison des simulations avec les mesures expérimentales réalisées dans la roue du modèle réduit de la turbine hélice démontrent que l'inclusion du jeu dans les simulations est essentielle pour mieux représenter les caractéristiques réelles de l'écoulement dans la roue de la turbine hélice et ainsi améliorer la prédiction de la performance de la roue d'une turbine hydraulique axiale.

Le comportement des individus est influencé par leur environnement, notamment le comportement de leurs pairs. Cette influence est particulièrement forte pour les enfants et les adolescents dans le milieu scolaire qui interagissent régulièrement avec le même petit ensemble de personnes. Dans ce mémoire, les effets de pairs sur la performance scolaire en français ont été estimés à l'aide d'une banque de données du ministère de l'Éducation, du Loisir et du Sport du Québec. Cette banque de données permet l'observation des résultats à l'examen commun de français des étudiants québécois de secondaire 5. Tout d'abord, un modèle naïf (linéaire simple) est estimé par MCO. Ce type de modèle ne permet pas la prise en compte de la simultanéité dans les décisions ainsi que les effets corrélés qui ne sont pas des effets sociaux. Ensuite, un modèle linéaire en moyenne avec effet fixe de groupe est estimé par la méthode de variables instrumentales généralisée. Ce modèle est d'abord estimé en ignorant l'un des deux effets sociaux, soit l'effet endogène (impact de la performance scolaire des pairs) ou les effets exogènes (impact des caractéristiques des pairs). La méthode d'estimation par variables instrumentales généralisée permet l'obtention d'estimateur asymptotiquement efficace. Grâce à la non-linéarité introduite par la variation dans les tailles de groupe, l'estimation simultanée de l'effet endogène et des effets exogènes est possible, même en prenant en compte les effets corrélés. Les résultats de l'effet endogène se sont avérés quelque peu surprenants. En effet, les estimations ont révélé un effet endogène négatif dans les spécifications pour lesquelles il est significatif. Un tel résultat pourrait s'expliquer par un effet de découragement ou par le fait que le professeur porte plus d'attention aux meilleurs élèves qu'aux élèves en difficultés.

Dans cette thèse, un modèle numérique a été proposé pour simuler la combustion liquide des carburants conventionnels et non-conventionnels, en particulier le mélange de biodiesel B20. La matrice de test numérique constitue de quatre cas d’écoulement réactifs c.à.d. avec combustion et d’un cinquième avec injection liquide sans combustion (écoulement non-réactif). Les modèles sont calculés à l’aide du logiciel FLUENT™ v.14 en 3D et a l’état stationnaire. Les flammes de diffusion turbulentes sont modélisées en utilisant l’approche de flammelette laminaire stable, avec une fonction de densité de probabilité jointe (PDF). La Validation est effectuée en comparant les mesures expérimentales disponibles avec les résultats obtenus de la CFD. L’aérodynamique de la chambre de combustion, ainsi que les températures de parois extérieures sont captures avec un degré de précision satisfaisant. La validation des principaux produits de combustion, tels que : CO2, H2O et O2, montre des résultats satisfaisants pour tous les cas d'écoulement réactifs, mais certaines incohérences ont été relevées pour les émissions de CO. On pense que le banc d'essai (la géométrie de la chambre de combustion et son état de fonctionnement) n'est pas suffisamment adéquat pour la combustion de combustibles liquides. D’autre part, et d’un point de vue numérique, l’approche de flammelette laminaire stable a été trouvé raisonnablement hors mesure de saisir les effets profonds du non-équilibre chimique qui sont souvent associés au processus de lente formation d’un polluant, comme le CO. / In this thesis, a CFD model was proposed to simulate the liquid combustion of conventional and non-conventional biodiesel fuels, in particularly the B20 biodiesel blend. The numerical test matrix consists of four reacting flow cases, and one non-reacting liquid fuel injection case. The models are computed using FLUENT™ v.14 in a 3D steady-state fashion. The turbulent non-premixed diffusion flames are modeled using the steady laminar flamelet approach; with a joint presumed Probability density function (PDF) distribution. Validation is achieved by comparing available experimental measurements with the obtained CFD results. Combustor aerodynamics and the outer wall temperatures are captured with a satisfactory degree of accuracy. Validation of the main combustion products, such as: CO2, H2O, and O2, shows satisfactory results for all the reacting flow cases; however, some inconsistencies were found for the CO emissions. It is believed that the test rig (combustor geometry and operating condition) is not sufficiently adequate for burning liquid fuels. On the other hand, from a numerical combustion point of view, the steady laminar flamelet approach was found not reasonably able to capture the deep non-equilibrium effects associated with the slow formation process of a pollutant, such as CO.

L’objectif principal de cette thèse est le développement d’un outil numérique, en utilisant l’approche XFEM, permettant la simulation des problèmes transitoires non linéaires thermoélectriques dans un milieu fissuré en deux dimensions, avec prise en compte des échanges thermiques et électriques entre les lèvres de la fissure. La simulation numérique de la propagation de fissures présente un grand intérêt pour de nombreux secteurs industriels (production d’aluminium, aéronautique, nucléaire, etc.). De plus, c’est un problème complexe sur le plan numérique. La méthode d’éléments finis classiques présentent des contraintes importantes de raffinement de maillage en fond de fissure, de remaillage pendant la propagation de la fissure avec la projection des champs, ce qui a pour effet d’augmenter le temps de calcul et de dégrader la précision des résultats. D’autre part, la méthode des éléments finis étendue XFEM, a reçu un succès grandissant pour le traitement de problèmes avec fissures durant la dernière quinzaine d’années. Elle permet d’utiliser un maillage qui ne se conforme pas à la géométrie des fissures, ceci grâce à un enrichissement de l’approximation éléments finis. Dans cette thèse, on s’intéresse à étendre le champ d’application de la méthode XFEM pour les problèmes non linéaires thermoélectriques avec fissures. En effet, le problème thermique transitoire est couplé avec le problème électrique par la génération de la chaleur dans le solide, et la génération de chaleur à la fissure à cause de la résistance de l’interface. Les échanges thermiques et électriques entre les lèvres de la fissure sont aussi considérés, et dépendent, respectivement, du saut de la température et du potentiel électrique à la fissure. En raison de la génération de la chaleur dans le solide et aux lèvres L’objectif principal de cette thèse est le développement d’un outil numérique, en utilisant l’approche XFEM, permettant la simulation des problèmes transitoires non linéaires thermoélectriques dans un milieu fissuré en deux dimensions, avec prise en compte des échanges thermiques et électriques entre les lèvres de la fissure. La simulation numérique de la propagation de fissures présente un grand intérêt pour de nombreux secteurs industriels (production d’aluminium, aéronautique, nucléaire, etc.). De plus, c’est un problème complexe sur le plan numérique. La méthode d’éléments finis classiques présentent des contraintes importantes de raffinement de maillage en fond de fissure, de remaillage pendant la propagation de la fissure avec la projection des champs, ce qui a pour effet d’augmenter le temps de calcul et de dégrader la précision des résultats. D’autre part, la méthode des éléments finis étendue XFEM, a reçu un succès grandissant pour le traitement de problèmes avec fissures durant la dernière quinzaine d’années. Elle permet d’utiliser un maillage qui ne se conforme pas à la géométrie des fissures, ceci grâce à un enrichissement de l’approximation éléments finis. Dans cette thèse, on s’intéresse à étendre le champ d’application de la méthode XFEM pour les problèmes non linéaires thermoélectriques avec fissures. En effet, le problème thermique transitoire est couplé avec le problème électrique par la génération de la chaleur dans le solide, et la génération de chaleur à la fissure à cause de la résistance de l’interface. Les échanges thermiques et électriques entre les lèvres de la fissure sont aussi considérés, et dépendent, respectivement, du saut de la température et du potentiel électrique à la fissure. En raison de la génération de la chaleur dans le solide et aux lèvres / The main objective of this thesis is the development of a numerical tool, using the XFEM approach, for the simulation of transient nonlinear thermoelectrical problems in fractured media in two dimensions, taking into account thermal and electrical exchanges between the crack’s lips. Numerical simulations of crack propagation are of great interest for many industrial sectors (aluminum production, aerospace, nuclear, etc.). In addition, this is a numerically complex problem. The classical finite element method has important constraints of mesh refinement at the crack tip, remeshing during crack propagation and field projections, which has the effect of increasing the computation time and degrading the accuracy. On the other hand, the eXtended Finite Element Method (XFEM), has received a growing success for the treatment of the problems containing cracks in the last fifteen years. It allows using a mesh that does not conform to the geometry of the crack; this is possible by the enrichment of the finite element approximation. In this thesis, we are interested in extending application field of the XFEM method to the nonlinear thermoelectrical problems with cracks. Indeed, the transient thermal problem is coupled to the electrical problem by the heat generation in the solid, and the heat generation at the crack’s lips due to the interface resistance. The heat and electrical exchanges between the crack’s lips are also considered, and depend, respectively, on the temperature and the voltage jump at the crack. Due to the heat generation in the solid and in crack’s lips (Joule effect), and the temperature dependence of the physical parameters of the material, the problem is nonlinear and fully coupled. The discretized nonlinear system by the XFEM method is solved using the Newton-Raphson algorithm. The robustness of the proposed technique is demonstrated through the simulation of different examples, and the results shows an excellent agreement with the analytical solution, or with the finite element solution using a refined mesh.

Les gestionnaires des entreprises sont souvent des personnes différentes des actionnaires de ces entreprises. Une littérature récente (Danthine et Donaldson, 2004) explore les conséquences potentielles de cette séparation entre gestion et propriété, habituellement ignorée, pour la modélisation macroéconomique. Toutefois, les travaux de Danthine et Donaldson se basent sur un modèle simple et sans frictions et il est donc difficile de mesurer leur impact sur la modélisation macroéconomique récente, dans laquelle interagissent plusieurs frictions et chocs. Le présent mémoire adapte le modèle de Danthine et Donaldson à un environnement avec chocs monétaires et frictions nominales (contrats de salaires); à ce chapitre, le mémoire représente donc un premier pas dans la généralisation de ce modèle.

Dans ce mémoire, nous adressons l'efficacité du marché financier de l'électricité en Amérique du Nord. Nous regardons spécifiquement la complexité des modèles à terme et la qualité de leurs prévisions. Nous concentrons notre étude sur deux modèles à terme, le modèle simple de Black et Scholes et un plus complexe, défini par Pilipovic. Nous regardons six marchés dans les États-Unis comme la MidColombia, le NP15, le ComEd, le Cinergy, le PJM et le NePool. La plage des données quotidiennes utilisées est de 1997 à 2002, soit pré et post Enron. Nous avons incorporé dans les modèles les propriétés stochastiques fondamentales associées aux prix à terme comme la diffusion de saut et moyenne inversée. Nous employons le facteur de corrélation pour étudier les modèles contre des prix passés et les données projetés. Nos résultats montrent que Black-Scholes est légèrement supérieur pour prévoir le prix court à terme. Les deux modèles ne sont pas appropriés pour l'évaluation à long terme. / In this thesis, we address the efficiency of the financial market of electricity in North America. We specifically look at the complexity of future models and the quality of the forecast. We focus our study on two future models, from a basic model, the Black-Scholes to a more complex, one defined by Pilipovic. We look at six markets in United-States as Mid-Colombia, NP15, ComEd, Cinergy, PJM and NePool. The daily data range from 1997 to 2 002, as pre and post Enron. We incorporated fundamental stochastic properties associated to spot and futures prices as jump diffusion and mean reversion. We use a correlation factor to define the fitness of the models against past and real prices. Our findings show that Black-Scholes is slightly better to foresee short term price. Both models are not appropriate for long term pricing.

Dans cette thèse, on étudie deux modèles décrivant l'évolution d'un bassin sédimentaire sous une contrainte sur le taux d'érosion. Ces modèles sont obtenus par l'application de la loi de conservation de masse sur le flux de matières, ce qui conduit à l'équation de Darcy ou à l'équation de Darcy-Barenblatt selon l'expression du flux choisie parmi deux expressions possibles d'après les géologues. L'équation de Darcy-Barenblatt est obtenue de celle de Darcy en ajoutant un terme de diffusion. En outre, la contrainte d'érosion maximale est implicitement contenue dans la formulation du modèle de Darcy-Barenblatt mais pas dans celle de Darcy en dimension 2. Après la présentation de ces modèles dans l'introduction de la thèse, la première partie est consacrée au modèle de Darcy-Barenblatt. On a obtenu un résultat d'existence d'une solution par une méthode de point fixe de Schauder-Tikhonov. Ensuite, on a montré un résultat de régularité en utilisant des résultats de Meyers et de Necas sur les équations elliptiques à coefficients höldériens, ce résultat de régularité est propre à une dimension inférieur ou égale à 2. La première partie se termine par la démonstration d'un résultat d'unicité de la solution. Le modèle de Darcy est étudié dans la deuxième partie de la thèse, on a obtenu une solution du problème discrétisé en temps, mais en dimension 2 d'espace le passage à la formulation continue fait apparaitre des produits de deux convergences faibles et soulève des difficultés théoriques non surmontées. En dimension 1, on a obtenu une solution continue pour le cas de la sédimentation marine en résolvant un problème à frontière libre de type Bernoulli d'évolution.

L'identification des systèmes dynamiques non linéaires à partir d'un ensemble de données entrée/sortie est d'une importance fondamentale pour les applications pratiques puisque beaucoup de systèmes physiques possèdent des caractéristiques non linéaires. La structure du modèle de Volterra peut être utilisée pour représenter une classe générale de systèmes non linéaires. Cependant, l'usage pratique d'une telle représentation est souvent limité à cause du grand nombre de paramètres associé à une telle structure. Pour pallier à cet inconvénient, plusieurs solutions sont proposées dans cette thèse. La première utilise des développements en série des différents noyaux sur des bases de fonctions orthogonales. La deuxième est basée sur l'utilisation de techniques faisant appel à des décompositions d'ordre réduit des tenseurs relatifs aux noyaux d'ordre supérieur ou égal à trois. Diverses bases de fonctions (Laguerre, Kautz et Bases Orthogonales Généralisées (BOG)) sont tout d'abord étudiées en vue de leur utilisation pour la modélisation des systèmes linéaires puis pour la représentation des noyaux de modèle de Volterra. Le problème d'identification comporte plusieurs volets : détermination des pôles caractéristiques des bases de fonctions orthogonales, de l'ordre des développements des différents noyaux, des coefficients de Fourier du développement et de l'incertitude relative à ces coefficients. Une représentation d'état associée à un développement sur une base de fonctions orthogonales généralisées est développée puis utilisée pour la construction de prédicteurs de la sortie du système ainsi modélisé. Ensuite, plusieurs décompositions tensorielles sont étudiées. La décomposition PARAFAC est plus particulièrement considérée. Des modèles de Volterra à complexité réduite inspirés de cette technique sont proposés. En considérant le noyau quadratique de Volterra comme une matrice et les autres noyaux comme des tenseurs d'ordres supérieurs à deux, nous utilisons une décomposition à l'aide des valeurs singulières (SVD) pour le noyau quadratique et la décomposition PARAFAC pour les noyaux d'ordres supérieurs à deux afin de construire le modèle réduit de Volterra appelé SVD-PARAFAC-Volterra. Un nouvel algorithme appelé ARLS (Alternating Recursive Least Squares) est présenté. Cet algorithme essentiellement basé sur la technique RLS appliquée d'une manière alternée estime les paramètres de tels modèles de Volterra. Enfin, de nouvelles méthodes d'identification robuste dites à erreur bornée sont présentées. Elles sont utilisées pour l'identification de modèles linéaires issus des BOG, travail qui vise à étendre au cas des systèmes non linéaires incertains des résultats obtenus récemment pour des systèmes linéaires incertains. Parmi les techniques d'identification à erreur bornée présentées, l'approche polytopique est plus particulièrement considérée. Cette approche nous permet d'estimer les intervalles d'incertitude des coefficients de Fourier du développement sur les différentes bases orthogonales étudiées. Ces mêmes méthodes d'identification sont utilisées aussi afin d'identifier les intervalles d'incertitude des paramètres du modèle SVD-PARAFAC-Volterra. Les méthodes proposées permettent de réaliser une importante réduction de complexité numérique et un gain en temps de calcul considérables.

Les problèmes aérothermiques typiquement rencontrés en ingénierie imposent un recours à des simulations numériques (CFD), qui, pour être réalistes, doivent nécessairement prendre en compte les phénomènes transitoires. L'emploi de modèles réduits ou de surfaces de réponse offre une alternative rapide et peu coûteuse à ces calculs exigeants en temps de calcul et en espace mémoire. L'objet de cette thèse est d'étoffer l'applicabilité de la méthode d'identification modale (MIM) dans le contexte de la mécanique des fluides. Il s'agit d'une méthode de réduction de modèle qui s'appuie sur la minimisation d'un écart entre les sorties d'un modèle de référence et celles d'un modèle réduit. Les paramètres du modèle réduit sont vus comme les variables d'un problème d'optimisation. Nous développons la méthode dans le cadre d'écoulements turbulents. L'écoulement autour d'un cylindre circulaire constitue le cas d'application principal de ces travaux. Il nous permet de tester nos méthodologies dans diverses configurations incompressibles : stationnaires, instationnaires, écoulements isothermes ou de convection forcée thermique. Ces modèles produisent des champs de vitesse, de température, des densités de flux de chaleur . . . Une autre application concerne un profil d'aile dans un écoulement sub à transsonique. Pour celle-ci, les sorties des modèles sont des coefficients de traînée et de portance, ainsi que des répartitions de coefficients de pression le long de la surface portante. Cet exemple nous permet de mettre en évidence l'impact de l'échantillonnage sur la construction du modèle.