Tableau d'honneur de la Faculté des études supérieures et postdoctorales, 2011-2012 / La radiochirurgie stéréotaxique est un mode de traitement délivrant la dose de prescription en une seule fraction. Pour ce faire des petits champs et des arcs non coplanaires sont utilisés de manière à diminuer la dose aux tissus sains. Afin de s'assurer que la dose planifiée est bien conforme à celle donnée lors du traitement, des mesures de contrôles de qualité s'imposent. Or, il n'existe pas encore de détecteur idéal permettant de connaître précisément la dose dans ces conditions de traitement. Afin de palier à cette lacune, nous avons développé un détecteur à fibre scintillante possédant des caractéristiques favorables pour la radiochirurgie stéréotaxique telle qu'une grande résolution spatiale, l'équivalence à l'eau et l'indépendance angulaire et énergétique. Pour ce faire, des mesures de facteurs de sortie brute (soit, plus rigoureusement, la multiplication du facteur collimateur avec le facteur fantôme) et de profils furent effectuées avec les collimateurs coniques utilisés en radiochirurgie stéréotaxique. Les films Gafchromic EBT2, la diode de radiochirurgie non blindée SFD, la diode blindée 60008, la micro-chambre à ionisation A16 (des détecteurs présents sur le marché) ainsi que le détecteur à fibre scintillante furent utilisés lors des mesures, de façon à comparer la réponse des différents dosimètres pour mettre à jour leurs forces et leurs faiblesses. Par la suite, un prototype composé de 49 détecteurs à fibre scintillante disposés en croix a été développé et construit afin de prendre des mesures de distributions de dose en deux dimensions. Finalement, une mesure d'un plan de traitement complet a été effectuée avec la fibre scintillante, la diode de radiochirurgie non blindée SFD et la micro-chambre à ionisation A16 afin d'évaluer la précision avec laquelle il était possible de faire une mesure de dose lors d'un contrôle de qualité avec ces différents détecteurs.

Le récent développement de la méthode de discrimination spectrale pour la calibration des dosimètres à fibres scintillantes plastiques (PSD) a conduit au développement du premier dosimètre commercial de ce type, soit l’Exradin W1 (STANDARD IMAGING INC., Middleton, WI, U.S.). Comme tout nouveau dosimètre, celui-ci se devait d’être caractérisé afin d’évaluer sa fiabilité en conditions standards d’irradiation. Le but premier du projet était d’effectuer cette validation. Les résultats obtenus ont démontré que le dosimètre à scintillation était un excellent candidat pour réaliser des mesures de dose hors-champ telles que rencontrées lors de l’évaluation de la réduction de dose aux appareils cardiaques implantables (ACI) par blindage. Les mesures effectuées avec l’Exradin W1 ont démontré une importante réduction de dose autant pour les simples champs antéro-postérieurs (de 40 % à 80 % de réduction de dose) que pour les traitements cliniques complexes (entre 5 % et 45 % de réduction de dose). / The recent development of the spectral method for correction of the Cˇ erenkov light in plastic scintillation detectors (PSD) have led to the first commercial dosimeter of this type, the Exradin W1 (STANDARD IMAGING INC., Middleton, WI, U.S.). As any new commercial dosimeter, it had to be characterized and validated in different situations of radiation. The first purpose of this project was to perform this validation. The results of characterization have shown that PSD is an efficient dosimeter to make out-of-field doses measurements as required in the evaluation of the feasibility of using a lead sheet to reduce dose to Cardiac Implantable Electronic Devices, CIED. Measurements performed with PSD shown an important reduction of doses by the shielding for square anterior-posterior field (40 % to 80 %) and also for complex clinical treatments (5 % to 45 %).

Cette thèse porte sur le développement et la validation d’outils afin d’améliorer l’efficacité des traitements de curiethérapie à haut débit de dose (HDR) du sein. Le but du projet vise donc à concevoir et valider un nouveau système de planification et de guidage en temps réel des traitements de curiethérapie HDR du sein, basé principalement sur l’échographie 3D (3DUS). Dans un premier temps, une étude clinique a été réalisée à l’aide d’un premier prototype 3DUS linéaire développé par le groupe. L’étude a permis de constater que le volume mesuré par échographie 3D est environ trois fois plus petit que celui mesuré par tomodensitomètre et que le système possède plusieurs limitations (ex : volume d’acquisition trop petit, guidage des cathéters impossible). Dans un deuxième temps, un algorithme d’optimisation de la position et du nombre de cathéters a été développé. Il s’est avéré que la méthode développée est robuste aux erreurs d’implantation et qu’il est possible de réduire significativement le nombre de cathéters sans avoir un impact important sur la dosimétrie. Dans un troisième temps, une étude a comparé cet algorithme au seul algorithme commercialement disponible, HIPO. Les résultats ont démontré que l’algorithme HIPO produit des plans significativement inférieur, en terme de dosimétrie, à l’algorithme présenté dans cette thèse. Dans un quatrième temps, deux méthodes ont été développées afin de réaliser une planification personnalisée, en temps réel des traitements de curiethérapie HDR du sein. Les deux méthodes se sont avérées efficaces pour réduire le nombre de cathéter. Une preuve de concept a été validée et elle a démontré le potentiel d’une planification personnalisée, en temps réel des traitements de curiethérapie HDR du sein. À l’aide des connaissances acquises dans l’étude clinique, un nouveau prototype 3DUS a été développé. Il inclut un nouveau mode d’acquisition hybride et un module de guidage de la position des cathéters. Les résultats présentés dans cette étude ont mis en évidence le fait que le système 3DUS hybride a la capacité de mesurer avec justesse les dimensions linéaires ainsi que les volumes. Il reconstruit également avec justesse la trajectoire des cathéters. Pour terminer, dans le but d’effectuer une reconstruction dynamique des cathéters, un système de guidage électromagnétique a été validé. Les résultats ont démontré que la reconstruction des cathéters, en curiethérapie HDR, est significativement plus juste et plus précise à l’aide d’un système de guidage électromagnétique qu’en utilisant les méthodes conventionnelles. / This thesis focuses on the development and validation of different tools to increase the efficacy of high dose rate (HDR) breast brachytherapy treatments. The project goal aim at designing and validating a new system for real-time guidance and planning of HDR breast brachytherapy treatments, based mainly on 3D ultrasound (3DUS). As a first step, a clinical study was performed using the first linear 3DUS prototype developed by our group. This study has shown the limitations of the current system (ex : small acquisition volume, no catheter tracking possibility) and that 3DUS volumes are three times smaller than computed tomography volumes. As a second step, a catheter optimization algorithm was developed. The algorithm was shown to be robust to catheter implantation errors and it was possible to reduce significantly the number of catheters without having a significant negative impact on the dosimetry. As a third step, a study was designed to compare the present catheter optimization algorithm to the only commercially available algorithm, HIPO. The results demonstrated that the HIPO algorithm produce significantly worse plan, in term of dosimetry, than the algorithm that was developed in the present thesis. As a fourth step, two methods were developed for personalized, real-time planning of breast HDR brachytherapy treatments. The two methods were efficient and they were able to reduce the number of catheters. A proof-of-concept was validated and it demonstrated the potential of a personalized, real-time planning approach for breast HDR brachytherapy. Using the experience acquired during the clinical study, a new 3DUS system was developed. The system includes a new hybrid acquisition approach and a module for catheter tracking. The results presented in this study have shown the ability of the hybrid 3DUS system to accurately measure linear dimensions and volumes. Furthermore, it allows the reconstruction of the catheters trajectory with accuracy as well as track them in real-time. Finally, in order to dynamically reconstruct catheters, an electromagnetic tracking system was validated. This study has shown that the reconstruction of catheters, in HDR brachytherapy, is significantly more accurate and precise with an electromagnetic tracking system than with the conventional methods.

Les observations suggèrent que la barre est une structure très présente parmi les galaxies spirales, se retrouvant dans près de la moitié des galaxies spirales de l’Univers proche. La présence d’une barre a un impact important sur l’évolution séculaire de la galaxie hôte, puisque cette dernière modifie la distribution spatiale de la matière et crée un flot de gaz vers le centre de la galaxie, provoquant ainsi une rétroaction positive sur la formation stellaire centrale de la galaxie. Une autre structure très présente dans les galaxies est un noyau actif (NAG), soit un trou noir supermassif dans le centre qui interagit avec la matière autour de celui-ci, et la rétroaction de cette interaction forme une région très brillante dans le centre de la galaxie. La rétroaction du NAG réchauffera le gaz environnant par radiation, et pourra aussi appliquer une poussée cinétique sur celui-ci, le poussant vers l’extérieur du centre de la galaxie. Cependant, l’effet de la rétroaction du NAG sur l’évolution de la galaxie et de sa formation stellaire centrale est encore un sujet mal compris, et on ne sait toujours pas si l’effet de cette rétroaction sur la formation stellaire est positif ou négatif. Dans le cadre de ce travail, j’ai effectué un total de 12 simulations de galaxies spirales de masse égale, dont 10 sont barrées. Le but est de vérifier quel est l’effet du NAG en faisant varier la rétroaction du trou noir, la quantité de gaz disponible et l’impact de la barre sur la formation stellaire centrale et globale de la galaxie. Les résultats de mes simulations montrent que la présence d’un NAG aurait un effet plutôt positif sur la formation stellaire lorsqu’il s’agit d’une galaxie spirale barrée, et principalement lorsque la galaxie est encore à un jeune stade de son évolution. Dans le cas d’une galaxie spirale non barrée, l’effet serait au contraire négatif en son centre et négligeable globalement. / Observations suggest that the bar is a very common structure among spiral galaxies, finding itself in near half of the spiral galaxies of the nearby Universe. The presence of a bar has an important impact on the secular evolution of the host galaxy because the latter modifies the spatial distribution of the material and creates a stream of gas toward the center of the galaxy, causing a positive feedback on the central star formation in the galaxy. Another very common structure in galaxies is an active galactic nucleus (AGN), a supermassive black hole in the center of the galaxy, which interacts with matter around it and the feedback from this interaction forms a very bright region in the center of the galaxy. Feedback from the AGN will warm up the surrounding gas by radiation and it can also apply a kinetic kick on it, pushing it outward from the center of the galaxy. However, the effect from the AGN on the evolution of the galaxy and its central star formation is still a misunderstood subject and we still do not know if the effect of the feedback on stellar formation is positive or negative. As part of this work, I made a total of 12 simulations of spiral galaxies of equal mass, which 10 of them were barred. The aim is to verify what is the effect of AGN by varying strength, available gas fraction, and the impact of the bar on the central and global star formation rate in the galaxy. The results of my simulations show that the presence of AGN would have a rather positive effect on star formation in barred spiral galaxy, and especially when the galaxy is still at an early stage of its evolution. However, in the case of a non-barred spiral galaxy, the effect would be opposite, negative in the center and generally negligible.

Cette étude porte sur le développement d’un outil de contrôle qualité basé sur l’expérience, dérivé du concept de frontière stochastique en économie et s’appuyant sur des connaissances géométriques spécifiques au patient pour améliorer la qualité des traitements de curiethérapie à haut débit de dose pour le cancer de la prostate. Cent plans cliniques de curiethérapie à haut débit de dose de la prostate ont été utilisés dans cette étude, dans laquelle l’échographie transrectale était la seule modalité d’imagerie. Une fraction unique de 15 Gy a était prescrite à tous ces patients. Un algorithme de recuit simulé de planification inverse a été appliqué pour réaliser tous les plans et Oncentra Prostate a été employé comme système d’imagerie et de planification du traitement en temps réel. Les recommandations relatives aux paramètres de dose de la société américaine de curiethérapie pour la cible et les organes à risque ont été suivies. Les relations entre les paramètres géométriques et les paramètres dosimétriques d’intérêt sont examinées. Les paramètres géométriques sont liés aux dimensions anatomiques des patients et ceux associés aux cathéters. Pour déterminer les paramètres géométriques dominants dans un modèle de frontière stochastique donné, les relations monotones entre les paramètres géométriqueset les paramètres dosimétriques d’intérêt sont mesurées avec une approche non paramétrique, à savoir le coefficient de corrélation de Spearman. Ensuite, une recherche de force brute est effectuée pour un modèle donné dans lequel différents modèles, incluant toutes les combinaisons possibles des paramètres géométriques dominantes, sont optimisés. L’optimisation est accomplie en utilisant une méthode de vraisemblance maximale implémentée dans le progiciel de calcul statistique R, avec son algorithme de recuit simulée généralisée. Le test du rapport de vraisemblance et sa valeur-p correspondante sont utilisés pour comparer la signification statistique de l’ajout de nouveaux paramètres géométriques aux modèles. Un modèle de production pour la cible et un modèle de coût pour chacun des organes à risque sont développés pour le traitement par curiethérapie à haut débit de dose guidé par l’échographie transrectale. De plus, pour valider si chacun des modèles développés est universel, nous l’appliquons à une autre catégorie de traitement de la curiethérapie à haut débit de dose, dans laquelle la tomodensitométrie était utilisée comme modalité d’imagerie plutôt que de l’échographie transrectale. Ainsi, une nouvelle cohorte de cent plans cliniques curiethérapie à haut débit de dose guidés par la tomodensitométrie est prise en compte. Un modèle de frontière stochastique de production pour la cible et trois modèles de coût pour les organes à risque basés sur la tomodensitométrie sont développés. Enfin, les modèles intégrés de la tomodensitométrie et de l’échographie transrectale sont comparés. / This thesis focuses on developing an experience-based quality control (QC) tool, derived from the concept of stochastic frontier (SF) analysis in economics and based on patient-specific geometric knowledge to improve the quality of the high-dose-rate brachytherapy (HDR-BT) treatment for prostate cancer. One hundred clinical HDR prostate BT plans, using the transrectal ultrasound (TRUS) asthe only imaging modality, all treated with a single fraction of 15 Gy, and made using Inverse PlanningSimulated Annealing (IPSA) algorithm, are studied. Also, Oncentra Prostate system is employed as the real-time 3D prostate imaging and treatment planning system. American Brachytherapy Society dose parameter recommendations for target and organs at risk (OARs) were followed. Relationships between all the different geometric parameters (GPs) and the four dosimetric parameters (DPs) V100 of the prostate, V75 of the bladder and rectum, and D10 of the urethra were examined. Geometricinformation of the patients and catheters are considered as different GPs. To find the dominant GPs in a given SF model, monotonic relationships between the GPs and DPs of interest are measured using a nonparametric approach: the Spearman correlation coefficient. Then, to determine the optimal SF model for each of the target production SF, and the OARs cost SF models, brute-force searches are performed. Different SF models including all the possible combinations of the dominant GPs in the SF model under study are optimized. Optimization is done using a maximum likelihood method implemented in the statistical computing package R, along with its Generalized Simulated Annealing algorithm. The likelihood ratio test and its corresponding p-value are used to compare the statistical significance of adding new GPs to SF models. A production SF (PSF) model for the target, and a costSF (CSF) model for each of the bladder, rectum, and urethra are developed for TRUS-guided HDR-BTtreatment. The difference between the dose value of a plan obtained by IPSA and the one predicted by an SF model is explored. Additionally, to verify if each of the models developed for the TRUS-guided category of the HDR-BT treatment for prostate is universal, we apply it on another category of HDR-BT treatment, in which computed tomography (CT) was used as the imaging modality. So, a different cohort of one hundred clinical CT-guided HDR-BT plans is taken into consideration. A target production SF and three OARs cost SF models are developed for the CT-based plans. Subsequently, the built-in SF models for the TRUS-based and CT-based plans are compared.

Ce mémoire rend compte de l'expertise développée par l'auteur au Centre de recherchede l'Institut universitaire en santé mentale de Québec (CRIUSMQ) en endoscopie fibrée. Il décrit la construction d'un nouveau type de microscope optique, le MicroscopeInterstitiel Panoramique (PIM). Par la juxtaposition d'un court morceau de fibre à gradientd'indice et d'un prisme à l'extrémité d'une fibre monomode, la lumière laser estfocalisée sur le côté de la sonde. Pour former une image, cette dernière est rapidementtournée autour de son axe pendant qu'elle est tirée verticalement par un actuateurpiézo-électrique. Ce design de système rotatif d'imagerie interstitielle peu invasif est uneffort pour limiter les dégâts causés par la sonde tout en imageant la plus grande régionpossible en imagerie optique cérébrale profonde. / This thesis documents the expertise developed by the author at the Centre de recherchede l'Institut universitaire en santé mentale de Québec (CRIUSMQ) in fibered endoscopy, particularly the design and construction of a new kind of optical microscope: ThePanoramic Interstitial Microscope (PIM). Through the juxtaposition of a short piece ofGraded-Index fibre and a prism at the end of a single-mode fibre, laser light is focussedon the side of the probe. To form an image, the latter is quickly spun around its axiswhile it is being pulled vertically by a piezoelectric actuator. This minimally invasivefluorescence rotary interstitial imaging system is an endeavor to limit the damage causedby the probe while imaging enough tissue to provide good context to the user in deep brain optical imaging.

Le multidétecteur HERACLES a été utilisé pour étudier les réactions ²⁵Na; ²⁵Mg + ¹²C à 9.23 MeV par nucléon à l’accélérateur ISAC-2 de TRIUMF. Les données brutes recueillies par HERACLES doivent être traitées pour en tirer des résultats physiques. Nous avons conçu un programme d’analyse codé en C++ pour traiter ces données. Le développement d’algorithmes d’identification a permis d’accélérer l’identification des particules dans les spectres des nombreux détecteurs d’HERACLES. Les résultats de l’identification montrent une identification en charge des particules pour les détecteurs des 4 premiers anneaux couvrant des angles de 4.8° à 24° par rapport au faisceau. Les particules légères sont identifiées en masse jusqu’à l’hélium dans les 2 autres anneaux couvrant les angles de 24° à 46°. L’étalonnage en énergie des détecteurs est déterminé à l’aide de la formule de Parlog. Nous avons réalisé une simulation hybride AMD et GEMINI et nous avons comparé les résultats aux données expérimentales pour les systèmes ²⁵Na; ²⁵Mg + ¹²C C à 9.23 MeV par nucléon.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2007-2008. / L’étude de l’interférence entre un plasmon polariton de surface (SPP) avec le faisceau plan lumineux incident sur une nano-fente métallique de dimension sous-longueur d’onde est devenue récemment un sujet fondamental de recherche dans le domaine de la plasmonique afin de mieux comprendre le mécanisme du rehaussement de la transmission optique de la nano-fente. Par la méthode de différences-finies dans le domaine de temps (FDTD), nous avons étudié le mécanisme de couplage du SPP dans la fente nanométrique. L’objectif de ce projet consiste dans un premier temps à la conception et l’implémentation d’un algorithme FDTD flexible et robuste, capable de simuler l’interaction de sources électromagnétiques avec des nanostructures de métal. L’algorithme développé sur la plate-forme Matlab permet de modéliser dans un espace bi-dimensionel des structures diélectriques dispersives. Dans un deuxième temps, nous avons employé ce simulateur FDTD pour étudier les mécanismes de couplage entre un SPP et une fente nanométrique percée dans une mince couche d’argent. Notre analyse démontre que le SPP incident à la fente est diffusé par le rebord de la fente et que les charges électriques induites sur les rebords ré-irradient l’énergie électromagnétique à l’intérieur de la fente. De plus, ces charges électriques génèrent des nouveaux SPPs sur les parois de la fente, qui contribuent à la formation de modes Fabry-Pérot le long de l’axe central de la fente. Nous démontrons aussi la formation de modes Fabry-Pérot créés par la multi-réflexion des ondes diffusées entre les parois de la fente. La combinaison de ces deux modes Fabry-Pérot produit une distribution de champ asymétrique dans la fente. Nous démontrons que la phase du SPP, relativement au faisceau incident normal, détermine les conditions d’interférence constructive et destructive correspondant respectivement au rehaussement et à la diminution de la transmission à travers la nano-fente. Finalement, nous avons confirmé la théorie d’interférence entre les champs induits par le SPP incident et le faisceau incident normal par l’addition de leurs amplitudes instantanées. / The study of the interference of the surface plasmon polariton (SPP) with the incident plane wave on a subwavelength metallic slit has become recently a fundamental subject of research in the domain of plasmonics. One of the objectives is to better understand the phenomenon of enhanced optical transmission through the nano-slit. Using the numerical method of finite-difference time-domain (FDTD), we have investigated the coupling mechanisms of the SPP inside the nano-slit. The objective of this project first consists in the conception and implementation of a flexible and robust FDTD algorithm capable of simulating the interaction between electromagnetic sources and metallic nanostructures. The algorithm that we have developed on the Matlab platform is able to model two-dimensional dispersive dielectric structures. The second and main objective is to use this FDTD simulator to investigate the coupling mechanisms between the SPP and a nano-slit pierced into a thin silver film. Our analysis demonstrates that the incident SPP is scattered by the inner edge of the slit and that the scattered waves induce oscillating electric charges on the slit edges which reradiate electromagnetic energy inside the slit. Moreover, these electric charges generate new SPPs on the slit walls which contribute to the formation of Fabry-Pérot modes along the central slit axis. We also show the formation of Fabry-Pérot modes set up by the multi-reflections of scattered waves between the slit walls. The combination of these two Fabry-Pérot modes produces an asymmetric field distribution inside the slit. We demonstrate that the phase of the SPP relatively to the normal incident beam determine the conditions of constructive and destructive interference which correspond, respectively, to the enhancement and suppression of the optical transmission through the nano-slit. Finally, we have confirmed the theory of interference between the induced fields by the SPP and incident normal beam via the superposition of their instantaneous fields.

Le dioxyde de vanadium (VO₂) est un matériau thermochrome de la famille des oxydes de métaux de transition (OMT). Sa température de transition (Tt), de 68 degrés Celsius, est parmi tous les OMT la plus proche de la température ambiante. Le VO₂ est métallique au-dessus de Tt et semi-conducteur en dessous de celle-ci. Le vanadium est un métal de transition avec plusieurs degrés d'oxydation. Lors de la fabrication des couches minces de VO₂, les différents degrés d'oxydation du vanadium entrent en compétition, ce qui rend très difficile la fabrication de couches minces stœchiométriques. On étudie les couches minces de VO₂ déposées par pulvérisation cathodique réactive ac double magnétrons assistée de faisceaux d'ions oxygène/azote pour des applications en optique et photonique. On dépose des couches minces de vanadium qu'on oxyde par la suite en dioxyde de vanadium. Les films ainsi fabriqués sont étudiés par des techniques de caractérisation optique, électrique, structurale et surfacique. Les propriétés de surface et de transport dépendent fortement de la proportion d'ions réactifs d'oxygène dans le faisceau ionique et du taux de conversion du V en VO₂. Le bombardement ionique diminue la température de transition et la résistance électrique des couches minces de VO2. Il mène à la formation de la phase monoclinique M2 du VO₂. Celle-ci est stabilisée par l'oxydation du V en VO₂ par la méthode "Rapid Thermal Annealing and Cooling" (RTAC). Ce processus permet un contrôle fin et aisé de la microstructure et de la texture des couches minces menant à de nouvelles propriétés optiques et électriques des couches minces de VO₂. Les applications du VO₂ proposées dans cette thèse vont des revêtements intelligents actifs et passifs aux lentilles minces nanométriques planes en passant par des dispositifs en couches minces thermo-électrochromes.

Les lasers à fibre optique dopée à l’ytterbium en régime continu de haute puissance ontune part de marché grandissante pour des applications d’usinage de métaux. Malgré qu’il s’agisse d’une technologie relativement répandue, un important problème subsiste dans laplupart de ces lasers. À haute puissance, la diffusion Raman stimulée transfère une partie de la puissance de la longueur d’onde principale d’émission des lasers autour de 1070 nm à la première bande de Stokes Raman autour de 1120 nm, ce qu’on appelle l’instabilité Raman. Cette puissance transférée est inutilisable et peut même être dangereuse pour le système laser et ses utilisateurs. Malgré les travaux théoriques et expérimentaux effectués sur ces lasers, très peu d’explications ont été fournies sur les liens entre les paramètres de la fibre optique et des réseaux de Bragg formant la cavité laser et l’instabilité Raman. Le but du projet de maîtrise présenté dans ce mémoire est donc de développer un modèle de simulation numérique de ces lasers, afin de comprendre et d’identifier les mécanismes dominants qui favorisent l’instabilité Raman et de trouver des configurations de montage la minimisant. Ce mémoire présente les deux modèles de simulation développés dans le cadre de ce projet. Le premier traite la propagation du signal laser comme étant unidirectionnelle, alors que le second la traite comme étant bidirectionnelle, ce qui se rapproche plus de la situation expérimentale. Le montage typique simulé est constitué d’une fibre optique à double gaine dopée à l’ytterbium ayant une grande aire modale effective, d’un réseau de Bragg à haute réflectivité et d’un réseau de Bragg à faible réflectivité servant de coupleur de sortie. Les simulations ont permis d’identifier cinq paramètres de la cavité laser ayant un impact important sur l’instabilité Raman. Une faible puissance moyenne du signal, une courte fibre optique de gain, une configuration de pompage en contrapropagation, c’est-à-dire par le côté du coupleur de sortie, ainsi qu’une plus faible réflectivité et une large bande réfléchissante du réseau de Bragg à faible réflectivité permettent de limiter la génération de l’instabilité Raman. L’optimisation de ces paramètres permet d’obtenir une cavité laser ayant extrêmement peu d’instabilité Ra-man. Ce faible niveau d’instabilité Raman semble être causé par une plus faible puissance intracavité, une courte distance de propagation et des modulations rapides de la puissance du signal. Des montages simulés incluant un filtre dans la cavité à la longueur d’onde de Stokes Raman, un réflecteur non linéaire ou une cavité de basse puissance amplifiée ont également montré une réduction significative de l’intensité de l’onde de Stokes Raman. / Continuous high-power ytterbium-doped fiber lasers have an increasing market share formetal processing applications. Despite their widespread use, these lasers still suffer a ma-jor problem. At high power, stimulated Raman scattering shifts the power from the main emission wavelength around 1070 nm to the first Raman Stokes sideband around 1120 nm. This process is called Raman instability. The shifted power becomes useless and can even be dangerous for both the laser system and its users. Previous experimental and theoretical analyses have failed to provide clear explanations on the link between the Raman instability and the parameters of the ytterbium-doped optical fiber and the fiber Bragg gratings forming the laser cavity. The goal of this master’s degree project was to develop a simulation model for continuous high-power ytterbium-doped fiber lasers in order to identify and understand how the parameters of the laser cavity affect the Raman instability and to find cavity configurations that reduce it. This master’s thesis presents the two simulation models developed during this project. The first model considers unidirectionnal propagation of the laser signal while the second one considers bidirectionnal propagation. The latter is thus a more realistic model of such lasers.The typical simulated setup is made of a double-clad ytterbium-doped fiber with a large mode area, a high reflectivity Bragg grating and a low reflectivity Bragg grating that isused as output coupler. The simulations allowed to identify five cavity parameters having an impact on the Raman instability. A low average power, a short gain fiber, a counter-propagation pumping setup as well as a low reflectivity and a large reflective bandwidth for the fiber Bragg grating used as the output coupler help minimizing the Raman instability.The optimisation of these parameters creates a laser cavity with an extremely low power shift to the Raman Stokes sideband. The low Raman instability seems to be caused by a lower intra-cavity power, a shorter propagation distance and fast power modulations in thesignal. Incorporating a filter in the cavity, using a nonlinear reflector as output coupler or using a setup that includes a low-power master oscillator in combination with a high-power amplifier have also been simulated and show a reduction of the Raman instability.