Tableau d'honneur de la Faculté des études supérieures et postdoctorales, 2016-2017 / Introduction: En 2015, 65 pays avaient des programmes de vaccination contre les VPH. La modélisation mathématique a joué un rôle crucial dans leur implantation. Objectifs: Nous avons réalisé une revue systématique et analysé les prédictions de modèles mathématiques de l’efficacité populationnelle de la vaccination sur la prévalence des VPH-16/18/6/11 chez les femmes et les hommes, afin d’évaluer la robustesse/variabilité des prédictions concernant l’immunité de groupe, le bénéfice ajouté par la vaccination des garçons et l’élimination potentielle des VPH-16/18/6/11. Méthodes: Nous avons cherché dans Medline/Embase afin d’identifier les modèles dynamiques simulant l’impact populationnel de la vaccination sur les infections par les VPH-16/18/6/11 chez les femmes et les hommes. Les équipes participantes ont réalisé des prédictions pour 19 simulations standardisées. Nous avons calculé la réduction relative de la prévalence (RRprev) 70 ans après l’introduction de la vaccination. Les résultats présentés correspondent à la médiane(10ème; 90èmeperccentiles) des prédictions. Les cibles de la vaccination étaient les filles seulement ou les filles & garçons. Résultats: 16/19 équipes éligibles ont transmis leurs prédictions. Lorsque 40% des filles sont vaccinées, la RRprev du VPH-16 est 53%(46%; 68%) chez les femmes et 36%(28%; 61%) chez les hommes. Lorsque 80% des filles sont vaccinées, la RRprev est 93%(90%; 100%) chez les femmes et 83%(75%; 100%) chez les hommes. Vacciner aussi les garçons augmente la RRprev de 18%(13%; 32%) chez les femmes et 35%(27%; 39%) chez les hommes à 40% de couverture, et 7%(0%; 10%) et 16%(1%; 25%) à 80% de couverture. Les RRprev étaient plus élevées pour les VPH-18/6/11 (vs. VPH-16). Si 80% des filles & garçons sont vaccinés, les VPH-16/18/6/11 pourraient être éliminés. Interprétation: Même si les modèles diffèrent entre eux, les prédictions s’accordent sur: 1)immunité de groupe élevée même à basse couverture, 2)RRprev supérieures pour les VPH-18/6/11 (vs. VPH-16), 3)augmenter la couverture chez les filles a un meilleur impact qu’ajouter les garçons, 4)vacciner 80% des filles & garçons pourraient éliminer les VPH-16/18/6/11. / Background: As of 2015, 65 countries have introduced HPV vaccination programmes. Mathematical models have played a key role in the implementation of these programmes. Objectives: We conducted a systematic review and pooled-analysis of model predictions of population-level effectiveness of HPV vaccination against HPV-16/18/6/11 infection in women and men, to examine the robustness/variability of predicted populationnal effects, incremental benefit of vaccinating boys, and potential for HPV vaccine-type elimination. Methods: We searched Medline and Embase (2009-2015) for transmission-dynamic modeling studies predicting the population-level impact of vaccination on HPV-16/18/6/11 infections among women and men in high-income countries. Participating modeling teams produced predictions for 19 standardized scenarios. We derived pooled relative reduction in HPV prevalence (RRprev) 70 years after vaccination, using the median (10th; 90thpercentile) of model predictions. Strategies investigated were Girls-Only and Girls & Boys vaccination at 12 years of age. Findings: 16/19 eligible models, from ten high-income countries provided predictions. With 40% Girls-Only vaccination coverage, HPV-16 RRprev among women and men was 53%(46%; 68%) and 36%(28%; 61%), respectively. With 80% Girls-Only vaccination coverage, HPV-16 RRprev among women and men was 93%(90%; 100%) and 83%(75%; 100%), respectively. Vaccinating boys in addition to girls increased HPV-16 RRprev among women and men by 18%(13%; 32%) and 35%(27%; 39%) for 40% coverage, and 7%(0%; 10%) and 16%(1%; 25%) for 80% coverage, respectively. RRprev were greater for HPV-18/6/11 than HPV-16 for all scenarios. Finally at 80% coverage, most models predicted that Girls & Boys vaccination would eliminate HPV-16/18/6/11, with a median RRprev of 100% for women and men for all types. Interpretation: Although HPV models differ in structure, data used for calibration and setting, population-level predictions were generally concordant: 1) strong herd effects from vaccinating Girls-Only, even with low coverage, 2) greater post-vaccination reductions in HPV-18/6/11 infection (vs. HPV-16), 3) increasing coverage in girls provides greater impact than including boys, and 4) reaching 80% coverage in Girls would eliminate HPV-16/18/6/11, with a median RRprev of 100% for women and men for all types. Interpretation: Although HPV models differ in structure, data used for calibration and setting, population-level predictions were generally concordant: 1) strong herd effects from vaccinating Girls-Only, even with low coverage, 2) greater post-vaccination reductions in HPV-18/6/11 infection (vs. HPV-16), 3) increasing coverage in girls provides greater impact than including boys, and 4) reaching 80% coverage in Girls would eliminate HPV-16/18/6/11, with a median RRprev of 100% for women and men for all types. Interpretation: Although HPV models differ in structure, data used for calibration and setting, population-level predictions were generally concordant: 1) strong herd effects from vaccinating Girls-Only, even with low coverage, 2) greater post-vaccination reductions in HPV-18/6/11 infection (vs. HPV-16), 3) increasing coverage in girls provides greater impact than including boys, and 4) reaching 80% coverage in Girls & Boys could eliminate HPV-16/18/6/11.

Le schéma traditionnel d'utilisation de l'énergie est désormais considéré comme un problème sérieux en raison de sa relation directe avec le changement climatique. Actuellement, notre dépendance vis-à-vis des combustibles fossiles augmente de façon spectaculaire, ce qui peut être attribué à la croissance de la population mondiale et à la forte demande d'énergie pour le développement économique. Ce modèle semble être préférable uniquement pour une économie florissante, mais ses perspectives pour les générations futures seront sans aucun doute décevantes. Dans ce scénario, un volume gigantesque de CO₂ produit par la combustion des combustibles fossiles dans les industries chimiques, les cimenteries et les centrales électriques, est rejeté de manière irresponsable dans l'atmosphère. Il ne fait aucun doute qu'une telle exploitation des combustibles fossiles nous conduit à des catastrophes environnementales sans précédent en ce qui concerne l'émission de CO₂, qui est le principal contributeur aux gaz à effet de serre (GES). L'une des solutions disponibles pour faire face à cette situation critique est de moderniser les centrales existantes qui émettent du CO₂ avec des technologies de capture et de stockage du carbone (CSC) afin de lutter systématiquement contre le changement climatique. Toutefois, les technologies actuelles de CSC présentent encore des problèmes techniques et des limites opérationnelles qui entraînent un surcoût pour les dépenses d'une usine et une augmentation de sa consommation d'énergie. La technologie membranaire est actuellement considérée comme une méthode de séparation prometteuse pour la séparation des gaz en raison de la simplicité de son procédé et de son mécanisme écologique. Par rapport aux autres méthodes de séparation, cette technologie est encore en cours de développement. Actuellement, la recherche se concentre sur l'amélioration des caractéristiques des membranes afin de faire face à un compromis bien connu entre la perméabilité et la sélectivité décrit par les graphiques de Robeson. Cette approche pourrait viser à commercialiser cette technologie plus efficacement dans le domaine de la séparation des gaz, tandis qu'une technologie d'absorption à base d'amines sera encore utilisée de manière dominante à cette fin pendant plusieurs années. Malgré cela, il est également nécessaire d'évaluer la performance des membranes fabriquées pour la séparation de différents mélanges de gaz avant de les utiliser pour des projets industriels réalistes. Pour ce faire, un outil de simulation est nécessaire pour prédire la composition des composants gazeux dans les flux de produits du rétentat et du perméat dans différentes conditions de fonctionnement. Ainsi, au chapitre 1, un modèle fiable est développé pour la simulation de la séparation des gaz à l'aide de modules de membranes à fibres creuses. Ensuite, ce modèle permet d'identifier les propriétés requises de la membrane, ce qui permet d'obtenir des performances intéressantes pour le module. Un procédé membranaire à plusieurs étapes est nécessaire pour atteindre les spécifications du produit qui sont une pureté et une récupération élevées du CO₂ dans le cas de projets de capture du CO₂. Dans ce cas, au chapitre 2, un modèle d'optimisation est proposé pour déterminer les valeurs optimales des paramètres de fonctionnement pour chaque étape et surtout pour déterminer une disposition optimisée à différents taux de récupération tout en minimisant le coût de la capture du CO₂. Dans le chapitre 3, nous comparons les performances de séparation de la technologie membranaire et du procédé d'absorption enzymatique en effectuant plusieurs analyses technico-économiques. Cette approche vise à démontrer la viabilité technique et l'efficacité économique de ces méthodes pour la modernisation d'une centrale électrique de 600 MWe par rapport aux procédés traditionnels d'absorption à base d'amines. Enfin, au chapitre 4, un système hybride est présenté en combinant les procédés d'absorption membranaire et enzymatique pour capturer le CO₂ des gaz de combustion d'une centrale électrique de 600 MWe. Ce système hybride est ensuite évalué pour révéler la faisabilité du procédé et pour étudier les performances de séparation en partageant la capture partielle du CO₂ entre ces deux unités de séparation. Dans l'ensemble, cette thèse de doctorat contribue à tirer parti de la fusion de la technologie membranaire avec d'autres méthodes de séparation conventionnelles telles que le procédé d'absorption enzymatique pour faciliter plus rapidement son intégration industrielle et sa commercialisation sur le marché de la séparation des gaz. / The traditional pattern of energy use is now regarded as a serious problem due to its direct relationship to the climate change. Currently, our dependency on fossil fuels is dramatically increasing which can be attributed to the world population growth and heavy energy demand for economic development. This model appears to be preferable only for flourishing economy but undoubtedly its outlook for the future generations will be disappointing. Under this scenario, a gigantic volume of CO₂ produced by burning the fossil fuels in chemical industries, cement manufactures, and power plants, is recklessly released in the atmosphere. Undoubtedly, such exploitation of the fossil fuels is bringing us further to unprecedented environmental disasters pertaining to the emission of CO₂ which is the major contributor to the greenhouse gases (GHGs). One of the available solutions to deal with this critical situation is to retrofit existing CO₂ emitter plants with carbon capture and storage (CCS) technologies in order to systematically combat with the climate change. However, the current CCS technologies still have technical issues and operational limitations resulting in incurring extra cost to a plant’s expenditures and increasing its energy consumption. Membrane technology is currently regarded as a promising separation method for gas separation due to its process simplicity and eco-friendly mechanism. In comparison to other separation methods, this technology is still under progress. Currently, the research focus is on the enhancement of membrane characteristics in order to deal with a well-known trade-off between permeability and selectivity described by Robeson plots. This approach might aim at commercializing this technology more efficiently in the gas separation area while an amine-based absorption technology will still be dominantly utilized for this purpose for several years. Despite this, it is also needed to evaluate the performance of fabricated membranes for the separation of different gas mixtures prior to utilizing for realistic industrial projects. To do so, a simulation tool is required to predict the composition of gas components in retentate and permeate product streams under different operating conditions. Thus, in Chapter 1, a reliable model is developed for the simulation of gas separation using hollow fiber membrane modules. Later, this model allows identifying the required membrane properties hence, resulting in module performances of interest. A multi-stage membrane process is required to hit product specifications which are high CO₂ purity and recovery in the case of CO₂ capture projects. In this case, an optimization model is proposed in Chapter 2 to determine the optimal values of operating parameters for each stage and more importantly to determine an optimized layout at different recovery rates while CO₂ capture cost is minimized. In Chapter 3, we compare the separation performance of membrane technology and the enzymatic-absorption process through performing several techno-economic analyses. This approach aims at demonstrating the technical viability and economic efficiency of these methods for retrofitting to a 600 MWe power plant compared to traditional amine-based absorption processes. Finally, a hybrid system is introduced in Chapter 4 by combining membrane and enzymatic-absorption processes to capture CO₂ from flue-gas of a 600 MWe power plant. This hybrid system is then assessed to reveal the process feasibility and to investigate separation performance through sharing partial CO₂ capture between these two separation units. Overall, this PhD thesis contributes to leverage the merge of membrane technology with other conventional separation methods such as the enzymatic-absorption process to more rapidly facilitate its industrial integration and commercialization in the gas separation market.

Tableau d'honneur de la Faculté des études supérieures et postdoctorales, 2011-2012 / Les cellules thalamocorticales (TC) du noyau ventro-postéro-latéral (VPL) du thalamus relayent l'information du système somatosensoriel (synapses excitatrices lemniscales aux dendrites proximaux) à la région correspondante du cortex, mais reçoivent également en rétro-propagation des projections du cortex (synapses excitatrices corticothalamiques aux dendrites distaux). Afin d'étudier l'intégration synaptique aux différentes parties de la cellule TC, nous avons bâti un modèle multi-compartimental à partir de reconstructions tridimensionnelles de cellules du noyau VPL, ce qui consiste en une discrétisation spatiale des dendrites en une multitude de segments associés à des circuits RC interconnectés. Nous avons pu dégager quantitativement l'impact de la géométrie cellulaire (taille d'arborisation et diamètre dendritique) sur l'amplitude et sur la durée des réponses au soma. Nous avons par la suite comparé l'intégration synaptique pour différentes distributions des entrées aux dendrites proximaux et distaux et sous différentes conditions de courants intrinsèques et de potentiel membranaire. Dans tous les cas, la sommation des entrées proximales induisait une réponse indépendante de la distribution, alors que la réponse aux entrées distales saturait lorsqu'elles étaient localisées aux mêmes branches. Nos résultats ont permis d'apporter une explication physiologique au patron d'organisation synaptique chez les cellules TC. / Thalamocortical (TC) cells from the ventroposterolateral (VPL) nucleus of the thalamus relay the somatosensory inputs (excitatory lemniscal synapses at proximal dendrites) to the corresponding cortical area, but also receive feedback excitatory inputs from the cortex (corticothalamic synapses at distal dendrites). The goal of this study was to compare the synaptic integration of inputs coming to proximal vs. distal dendrites. A multicompartmental model was drawn from fully reconstructed cells of the VPL nucleus. Dendrites were spatially discretized in multiple segments associated to interconnected RC circuits. We were able to characterize the impact of neuronal size and dendritic diameter on the amplitude and on the time course of the somatic response. We also compared the synaptic integration for different distributions of proximal or distal inputs under different conditions of membrane potential and active properties. In all cases, the summation of proximal inputs was independent of their distribution, while the response induced by distal inputs saturated when those inputs were located at the same branches. The results obtained in this study suggest a physiological explanation of the synaptic pattern at TC cells.

Les vents galactiques sont des éléments importants à considérer dans les simulations numériques à grande échelle car ils ont des impacts sur la formation des galaxies environnantes. Puisque les galaxies sont mal résolues dans de telles simulations, les vents galactiques sont habituellement générés par des méthodes semi-analytiques. Dans le cadre de ce projet, un modèle galactique a été développé afin d'améliorer le modèle semi-analytique de Pieri et al. (2007). Ce nouveau modèle permet de suivre de manière consistante l'évolution de l'enrichissement des galaxies en tenant compte des vents stellaires, des supernovae et de différents scénarios de formation stellaire. Les vents galactiques sont générés par l'énergie thermique provenant des supernovae et des vents stellaires à l'intérieur des galaxies. Avec ce formalisme, seules les galaxies ayant une masse inférieure ou égale à 10¹⁰ MQ risquent de contribuer à l'enrichissement du milieu intergalactique. La distribution des vents galactiques dans ce milieu est calculée en respectant l'ordre chronologique des éjectas. De plus, la composition de ce vent peut désormais être décomposée en 31 éléments chimiques. Pour la même quantité d'étoiles formées durant l'évolution galactique, un taux de formation stellaire de longue durée produit un plus long vent galactique qu'un taux de formation stellaire de courte durée. Cependant, ce vent est alors moins dense et moins concentré en métaux. En augmentant l'efficacité de formation stellaire, la portée et la métallicité du vent galactique augmentent également. Par contre, dans certains cas, une trop grande quantité d'étoiles peut complètement balayer le milieu interstellaire de son gaz, ce qui altère l'évolution du vent galactique. Pour respecter la quantité de métaux observée dans le milieu intergalactique, les vents galactiques doivent provenir des galaxies ayant possédé une métallicité initiale différente de zéro au moment de leur formation. Dans ce cas et lors d'une collision galactique, les vents stellaires peuvent contribuer de manière significative au bilan énergétique et à la quantité de carbone et d'azote éjectée dans le milieu intergalactique.

Nous avons utilisé un algorithme P³M afin de simuler la formation de galaxies dans une boîte cosmologique de (15 Mpc)³ dans un univers ΛCDM d'un décalage spectral de 24 à 2. Nous avons ensuite simulé des vents galactiques anisotropes et suivi l'effet de rétroaction de ces vents sur l'évolution des galaxies naines simulées. Les vents sont modélisés comme étant bipolaires et les angles d'ouverture étudiés sont de α=60°, 90°, 120°, 150° et 180°. Ces vents modélisés ont tendance à se propager dans une direction perpendiculaire à la structure à grande échelle (filament, crêpe) hébergeant la galaxie. Nous incluons l'effet de la suppression de galaxies par la photoionisation à partir d'un décalage spectral z=6 dans cinq de nos simulations et l'ignorons dans les cinq autres afin de comparer cette dernière série de résultats avec des travaux antérieurs. Nous incluons aussi les interactions entre les halos (accrétion, fusion et fission) pour modéliser l'apparition de sursauts d'étoiles produisant des supernovae. Les vents anisotropes sont les plus susceptibles d'enrichir les régions de faibles densités. De plus, ils ont moins de chances de se superposer, ce qui augmente la fraction de milieu intergalactique enrichie en métaux. L'anisotropie grandissante des vents diminue aussi la probabilité qu'un halo soit frappé et dépossédé de son gaz, empêchant ainsi la formation d'une galaxie. Lorsque l'on diminue l'angle d'ouverture des vents de α=180° (un vent isotrope) à un angle α=60°, le nombre de galaxies créées double. Cela a pour effet de produire le double de vents. Ces vents étant plus anisotropes, ils se superposent moins, ce qui a pour résultat que la fraction de volume enrichie passe de 8%, dans le cas de vents isotropes, à 28% pour des vents ayant un angle d'ouverture de 60°. L'anisotropie des vents augmente l'enrichissement des régions de toutes densités, ce qui est dû en partie au fait que les vents de hautes anisotropies sont plus nombreux. Lorsque nous ne tenons pas compte de cet effet, nous remarquons que l'anisotropie des vents crée une augmentation de l'enrichissement pour des densités allant jusqu'à 10p̄, où p̄ est la densité moyenne. Nous attribuons cet effet à l'évolution dynamique de nos simulations. Le gaz situé dans les régions de faible densité est fortement enrichi par les vents anisotropes, mais une partie de ce gaz est ensuite accrétée par les structures à grandes échelles. Les régions de hautes densités sont plus efficacement enrichies que les régions de faibles densités (~80% comparé à ~20%), mais celles-ci sont privilégiées (un maximum d'enrichissement par volume à 0.3p̄). L'effet de photoionisation diminue grandement le nombre de galaxies de faibles masses formées à z<3, ce qui produit une baisse de la fraction de volume enrichie après z=3 puisque l'effet d'accrétion l'emporte sur le nombre de vents créés, particulièrement dans les régions de faibles densités.

Tableau d'honneur de la Faculté des études supérieures et postdoctorales, 2011-2012 / Les réseaux sont des structures mathématiques abstraites qui sont utilisées pour modéliser une grande variété de systèmes complexes. Ils trouvent une application importante en épidémiologie, où on les utilise pour modéliser les contacts entre les individus d'une population afin d'étudier comment un agent viral pourrait s'y propager. Jusqu'à maintenant, la majorité des modèles mathématiques en épidémiologie sur réseaux ont été développés dans le cadre du paradigme de V épidémie isolée, dans lequel le seul processus étudié est la propagation d'un agent viral. En réalité, plusieurs processus peuvent toutefois avoir lieu simultanément et interagir au sein d'une même population. On parle alors de processus coévolutifs. Ce mémoire est consacré au développement des outils nécessaires à la modélisation de processus coévolutifs en épidémiologie sur réseaux. On s'intéresse plus particulièrement à deux types de coévolution, soit la coévolution d'un agent viral avec la topologie du réseau et la coévolution de deux agents viraux.

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.

SITELLE est le nouveau spectromètre imageur à transformée de Fourier qui doit être installé au Télescope Canada-France-Hawaii. Le développement de son modèle de performance a permis d’évaluer les paramètres instrumentaux qui sont critiques par rapport aux exigences scientifiques. D’une part, une configuration hors-axe à miroirs plans a été sélectionnée pour répondre plus facilement à l’exigence de haute efficacité dans le proche ultraviolet. D’autre part, les critères d’asservissement ont été définis afin de concevoir un instrument limité par le bruit de photons pour une scène pertinente donnée. Ces performances doivent être maintenues pendant une acquisition totale de 4h sous des vibrations opérationnelles et toute vibration extérieure comme des bourrasques de vent sur le télescope. Enfin, le modèle de performance de SITELLE est au cœur du simulateur de rapport signal sur bruit qui permettra aux astrophysiciens d’évaluer les bénéfices potentiels de l’utilisation de ce spectromètre imageur pour leurs sujets d’études. / SITELLE is a new imaging Fourier transform spectrometer to be installed at the Canada-France- Hawaii Telescope. The development of its dedicated performance model drives the design of the instrument and the flow down of the science cases requirements into system requirements. First, the selected configuration with off-axis flat mirrors makes the achievement of a high efficiency in the near ultra violet easier. Secondly, servomechanism’s desirable performances were defined in order to design a photon noise limited instrument, based on a relevant scene. These performances should be maintained during a 4 hours data-cube acquisition, under operational vibrations and external effects such as wind gust hitting the telescope. Ultimately, this instrument performance model is the core of the signal to noise ratio simulator that will help observers to evaluate the potential benefits of SITELLE for their target.

Tableau d'honneur de la Faculté des études supérieures et postdoctorales, 2011-2012 / Ce mémoire présente la méthodologie et les résultats d'une étude numérique portant sur la modélisation et l'optimisation d'un système de chauffage par panneaux radiants hydroniques à faible masse thermique. L'objectif principal du projet est de fournir des données et de proposer des solutions au problème d'optimisation des performances énergétiques et du confort de tels systèmes. L'évolution du processus d'optimisation entrepris passe d'abord par une analyse détaillée de la modélisation du système, au terme de laquelle sont créés des modèles représentant la géométrie réelle des panneaux. Entre autres, un modèle semi-analytique est développé pour les panneaux avec configuration de tubes en serpentin et un modèle de volumes finis est développé pour les configurations de tubes en spirale. Les modèles choisis offrent la possibilité d'imposer des propriétés et des taux de transfert de chaleur non uniformes afin de représenter adéquatement la réalité lors d'un couplage avec l'environnement. Ces modèles servent ensuite à identifier l'effet de plusieurs paramètres de conception sur les performances énergétiques des panneaux radiants, d'abord dans un contexte découplé de l'environnement. Finalement, le problème complet du confort et de la demande énergétique de ces systèmes en climat hivernal québécois est traité dans la dernière partie du projet. Cette étape présente le couplage d'un modèle de dynamique des fluides numériques (DFN) avec le modèle semi-analytique de panneaux radiants pour tubes en serpentin développé préalablement. Le tout est intégré dans une pièce typique d'un bâtiment résidentiel afin d'analyser l'effet de la température d'entrée du fluide caloporteur ainsi que du positionnement et du dimensionnement des panneaux. Le modèle de DFN calcule les champs de vitesse de l'air, de température et de rayonnement dans la pièce pour fournir des valeurs locales du coefficient de transfert de chaleur au modèle de panneau radiant. Cette méthodologie permet de considérer le système de panneau radiant intégré à son environnement pour en optimiser les caractéristiques en regard du confort et de la performance énergétique.

Les travaux de recherche effectués au cours de cette thèse ont mené au développement d'un modèle numérique permettant de décrire le déplacement simultané de plusieurs ions dans le réseau poreux de matériaux cimentaires. L'algorithme sépare le transport des ions des réactions chimiques selon une approche SNIA (Sequential Non Iterative Approach). Dans la partie consacrée au transport, le déplacement des ions dans un milieu non saturé est basé sur la résolution de l'équation de Nernst-Planck étendue tenant compte de la diffusion, du couplage électrique entre les ions et de l'activité chimique. Un terme d'advection est ajouté à cette équation afin de considérer le transport des ions sous l'effet du mouvement du fluide dans le réseau poreux en présence de gradients de teneur en eau. Ces équations de transport sont discrétisés dans l'espace par la méthode des éléments finis. Elles sont résolues pou: sept espèces ioniques de manière couplée en utilisant la méthode de Newton-Raphson. Un module d'équilibre chimique corrige ensuite les profils de concentration obtenus à l'étape de transport. Dans ce module chimique, les phases solides considérées sont la portlandite, les C-S-H, les monosulfoaluminates, l'ettringite, le gypse et les chloroaluminates. Certaines de ces phases solides sont soient dissoutes, soient précipitées afin de retrouver l'équilibre entre la pâte de ciment hydraté et la solution dans les pores du matériau. L'effet de la dissolution ou de la précipitation de ces phases solides sur les propriétés de transport du matériau est également considéré, en fonction des variations de porosité. Le modèle donne comme résultats les profils de concentrations en solution et la teneur en tout point du matériau de chacune des phases solides. Les applications présentées à la fin du mémoire montrent l'utilisation du modèle pour des cas d'attaque aux sulfates. Une méthode permettant de déterminer le coefficient de diffusion des ions dans les matériaux cimentaires suite à une analyse des résultats de l'essai de migration a également été développée dans le cadre de ces travaux de recherche.