Le modèle de gestion de la qualité de l'eau SWAT pourra assister les acteurs de l'eau dans la mise en application de politiques de gestion de l'eau à l'échelle du bassin versant au Canada. Cependant, il y a peu d'applications du modèle SWAT sous des conditions climatiques caractéristiques du sud-est du Canada. Un protocole de mise en place du modèle SWAT est proposé. La performance des fonctions hydrologiques du modèle a fait l'objet d'une évaluation en distinguant les saisons hydrologiques hiver et été. Les sources d'incertitude dans les données de calibration pour la saison hivernale et dans la structure du modèle de fonte de SWAT ont fait l'objet d'une analyse détaillée. La capacité des fonctions hydrologiques de SWAT à converger vers un jeu unique de paramètres et à assurer que les voies d'entraînement des indicateurs de la qualité des eaux sont adéquatement modélisées est discutée. / The river basin model SWAT could assist stakeholders in the setting up of integrated watershed management policies in Canada. However, there is a lack of SWAT's applications under climatic conditions characteristic of Southeastern Canada. An implementation protocol of the SWAT model is proposed. The performance of the hydrologic component has been evaluated in distinguishing the hydrologic seasons winter and summer. The sources of uncertainty in the calibration data of the winter season, and in the structure of the SWAT snowmelt algorithm have been the object of a detailed analysis. The capabilities of the SWAT hydrologic component to converge towards a unique parameter set and to guarantee that the flow paths of water quality indicators are correctly modeled is addressed.

Le barrage anti-marée Joseph-Samson, situé à l'embouchure de la rivière Saint-Charles, crée des conditions idéales pour l'accumulation de sédiments qui menacent les écosystèmes et la sécurité des gens qui habitent les abords de la partie aval de la rivière Saint-Charles. Des observations sur le terrain ont permis de constater l'étendue des effets de la sédimentation. Des simulations effectuées à l'aide des modèles River2D et HEC-RAS dans différentes conditions de débits et de marées ont permis de constater les avantages d'ouvrir la vanne de fond du barrage lors des cycles de marnage normaux tant au niveau de la qualité des écosystèmes qu'au niveau de la sécurité et de l'évacuation des sédiments. Dans la grande majorité des scénarios, l'ouverture de la vanne offre des vitesses d'écoulement bonnes pour l'évacuation du limon accumulé en amont du barrage et un échange important avec les eaux mieux oxygénées du Fleuve Saint-Laurent.

On étudie la propagation d'une flamme dans une chambre prolongée par un tube. C'est une configuration classique rencontrée lorsqu'on veut limiter le niveau de pression généré par une explosion interne. Mais au changement de section une explosion secondaire peut se produire dans la canalisation, ce qui conduit à une surpression supplémentaire dans la chambre.<br />Les résultats expérimentaux obtenus pour un prémélange propane-air ont été confrontés à des simulations numériques. Pour des petits diamètres de tube, on observe l'apparition de l'effet de flamme tulipe dans la chambre, qui contribue à limiter le débit d'évacuation des gaz. Pour les grands tubes, l'explosion secondaire à l'entrée du tube conduit à l'intensification de la combustion retardée des poches de gaz frais situés dans les coins de la chambre, causant ainsi une remontée brutale de la pression.<br />Pour éviter la formation de l'explosion secondaire une solution consiste à placer un insert (réseau de fils métalliques) à l'entrée du tube afin d'obtenir une extinction totale ou partielle de la flamme. Les expériences montrent qu'une optimisation du dispositif est nécessaire pour obtenir l'atténuation maximale de la supression. Un modèle thermique d'absorption de la chaleur par l'insert a été conçu, s'appuyant sur une méthode nodale. Le critère de transmission de la flamme au delà de l'insert est basé sur la comparaison entre la température d'auto-inflammation du mélange avec une température moyenne des gaz calculée en sortie de l'insert. Les résultats numériques, en bon accord avec les résultats expérimentaux, confirment le rôle prépondérant des pertes thermiques provoquées par l'insert sur le coincement de la flamme.

Le travail de thèse présenté dans ce manuscrit porte sur l'étude de la fragmentation du 12C à 95 MeV par nucléon par collision nucléaire avec des cibles épaisses de PMMA. Cette étude est motivée par le développement d'une nouvelle technique d'irradiation des tumeurs cancéreuses : la thérapie par ions carbone ou hadronthérapie. L'objectif est de comparer des mesures expérimentales de fragmentation avec des simulations GEANT4 afin de tester différents modèles nucléaires. Le but recherché est de déterminer si les modèles sont suffisamment prédictifs au regard des critères exigés en hadronthérapie : 3% d'incertitude sur le dépôt de dose en fin de parcours. Pour atteindre cet objectif, une première expérience a été réalisée au GANIL avec un faisceau de 12C à 95 MeV par nucléon et des cibles épaisses de PMMA. Cette expérience a permis d'obtenir les taux de production, les distributions angulaires et énergétiques des différents fragments produits lors des collisions nucléaires. Des comparaisons entre ces données expérimentales et les résultats simulés obtenus par le modèle de cascade binaire intra-nucléaire (BIC) et le modèle de dynamique moléculaire quantique (QMD) disponibles dans GEANT4 ont été effectuées. Ces comparaisons font ressortir l'inaptitude des modèles testés à reproduire la fragmentation du carbone à 95 MeV par nucléon avec la précision requise en hadronthérapie.

Le travail présenté dans ce mémoire de thèse concerne l'étude, la réalisation et la caractérisation de scintillateurs plastiques dopés au Bore. Leur mise en œuvre au sein de dispositifs de détection de neutrons rapides et thermiques y est examinée.<br />Dans une première partie, le processus de fabrication et la méthode de dopage en éléments borés sont présentés de façon détaillée. Différentes analyses physico-chimiques (I.R., UV-vis et RMN) ont été utilisées pour caractériser le matériau obtenu et pour confirmer sa composition et sa structure. De plus, une série d'expériences sous rayonnements alpha et sous faisceaux de protons a permis de mesurer les caractéristiques de scintillation. L'émission lumineuse peut ainsi être déterminée grâce à la relation semi-empirique de Birks.<br />Une seconde partie rend compte de la modélisation complète de l'ensemble du détecteur : interaction rayonnement - matière, émission et parcours de la lumière, formation du signal. Le code de calcul MCNP (Monte Carlo N-Particles) a été couplé à un code de génération et de propagation de lumière développé au cours de ce travail. Ces moyens de simulation ont permis d'optimiser les dimensions du détecteur pour la détection neutronique et de cerner l'influence de la géométrie sur la collecte des photons de scintillation. L'efficacité de détection et la durée de vie des neutrons dans ce type de détecteurs ont été estimées. Les mesures confirment la validité des simulations.<br />Une troisième partie fait état des méthodes de discrimination neutron - gamma analogique (méthode du passage par zéro) et numérique (méthode de la comparaison de charges) mises en œuvre pour ce scintillateur. Leurs performances respectives sont comparées et discutées.<br />La dernière partie concerne les applications qui peuvent bénéficier de l'utilisation de scintillateurs plastiques dopés au Bore, comme les fibres scintillantes dopées, la spectrométrie neutronique ou encore les systèmes de mesures des multiplicités neutroniques.

Ce travail de thèse s'inscrit dans le cadre du projet SPIRAL (Système de Production d'Ions Radioactifs Accélérés en Ligne) dont la mise en service débutera en septembre 2001 au GANIL (Grand Accélérateur National d'Ions Lourds) de Caen. Cette thèse a essentiellement porté sur le développement de systèmes de production d'ions radioactifs (ensemble cible source) par l'étude approfondie de chaque étape de la production par la méthode ISOL (Isotopic Separation On Line) : la production par fragmentation de la cible et/ou du projectile, la diffusion hors du matériau cible, l'effusion jusqu'à la source d'ionisation et enfin l'ionisation des atomes radioactifs. Une étude bibliographique et des simulations thermiques ont permis d'optimiser au mieux les matériaux et la forme des cibles de production et de diffusion. Une première cible a été optimisée et fiabilisée pour la production de gaz rares radioactifs (Argon, Néon,...). Une deuxième cible dédiée à la production d'Héliums radioactifs a entièrement été conçue (du cahier des charges jusqu'aux tests “hors ligne” et “en ligne”). Enfin, un troisième ensemble cible source a été défini pour la production d'alcalins radioactifs monochargés. Les intensités de faisceaux secondaires prévues pour SPIRAL y sont présentées. Une étude détaillée sur les efficacités de diffusion effusion pour ces différentes cibles a montré que l'utilisation d'un carbone dont la microstructure est fine (grain de 1µm) améliore la diffusion ainsi que l'importance de l'épaisseur des lamelles dans l'effusion des isotopes de courts temps de vie.

L'amélioration de la sécurité incendie au sein de l'habitat est un des principaux objectifs de la recherche actuelle. En effet, chaque année, un grand nombre de feux sont déclarés, générant la perte de nombreuses vies humaines, de fortes pertes financières, l'endommagement des structures et la pollution de l'environnement. Face à cette problématique, on remarque qu'un grand nombre de pays d'Europe possèdent une législation très pauvre vis-à-vis de la protection incendie dans l'habitat. Historiquement, les bâtiments ont été dessinés suivant des obligations prescriptives. La tendance de l'ingénierie de la sécurité incendie (Fire Safety Engineering, FSE selon le sigle Anglais) a changé amplement pendant la dernière décennie : des groupes de recherche dans le domaine de l'incendie ont mis au point les principes du design fondé sur la performance (Performance Building Design, PBD en Anglais). Le PBD a permis une approche de la sécurité incendie fondée sur la prédiction du comportement d'un incendie dans des scénarios donnés, en utilisant des outils numériques d'ingénierie. L'approche PBD de FSE est une méthodologie qui a été initialement développée pour les établissements recevant du public, toutefois peu à peu cette approche commence à être utilisée dans tout type d'habitat. La prédiction du comportement d'un incendie nécessite le calcul du débit calorifique (Heat Release Rate, HRR en Anglais) qui est la grandeur physique utilisée pour la mesure de la puissance d'un feu. En ingénierie, le HRR est indispensable à l'estimation de la sévérité du sinistre et des possibles endommagements causés dans un scénario donné. Sa détermination dépend des combustibles présents lors de l'incendie ainsi que de l'environnement du sinistre. La prédiction du HRR est réalisée à l'aide des codes de simulation numérique de l'incendie. Ceux-ci sont un assemblage de plusieurs sous modèles dont chacun calcule un ensemble des phénomènes présents dans la combustion p. ex. la pyrolyse, le rayonnement, la turbulence, etc. 8 La capacité à prédire correctement le HRR est limitée par les calculs très simplifiés du processus de décomposition thermique des solides. La décomposition est notamment dépendante des processus diffusifs et chimiques mis en jeu dans la zone comprise entre le solide et la flamme, lesquels ne sont pas modélisés de façon rigoureuse. Par le passé, plusieurs études expérimentales ont permis de mesurer le HRR d'un certain nombre de produits, cependant, ils ne contribuent pas à la compréhension de la physique du processus de décomposition de la matrice solide, donnée pourtant essentielle car source des espèces volatiles et du débit massique du combustible. En effet, un grand nombre de simulations trouvées dans la littérature font une approche empirique de la production de fuel ou considèrent une seule étape de décomposition. C'est dans ce contexte que prend place la présente étude qui vise à caractériser la cinétique de décomposition de combustibles solides et de formation des espèces volatiles : les changements survenus dans la phase solide sont pris en compte ensemble avec ceux de la phase gazeuse (dégagement d'espèces). La détermination du mécanisme de décomposition est une tâche fondamentale de l'analyse thermique. Le mécanisme doit considérer la succession des transformations de la matière pendant la gazéification des solides. Cette succession inclus les échantillons vierges ainsi que ceux qui ont déjà souffert des attaques thermiques (sous produits des étapes de décomposition). Le mécanisme de décomposition constitue une des principales données d'entrée de la grande majorité de modèles de décomposition thermique. Cette recherche tient compte de la décomposition thermique d'une mousse polyéther polyuréthane (PPUF) à trois échelles différentes. Chaque échelle caractérise le comportement au feu d'une masse différente de mousse et est concentrée sur l'étude de phénomènes particuliers : · L'échelle matière permet l'analyse du comportement d'échantillons avec des masses proches d'un milligramme. À l'échelle matière, les effets de transfert de chaleur et des espèces sont minimisés et l'effet de l'augmentation de la température du solide peut être étudié précisément. L'échantillon est considéré comme une particule de masse et de dimension négligeables, de sorte que sa température soit homogène. · La petite échelle permet l'analyse des échantillons avec des masses proches de dix grammes. À l'échelle matière des gradients de transfert de chaleur et d'espèces existent. L'échantillon est irradié seulement par une des surfaces, produisant ainsi 9 le déplacement du front de décomposition. La combustion de matériaux polymériques est complexe et concerne souvent des processus simultanés tels que la pyrolyse, la décomposition oxydative et le processus de combustion avec présence de flamme. · L'échelle produit concerne des échantillons avec des masses proches d'un kilogramme. À cette échelle, la géométrie et le positionnement d'un produit ont un rôle fondamental dans la croissance du feu. La ventilation (la disponibilité d'oxygène et la turbulence) affecte également le processus de combustion. L'échelle produit montre le comportement au feu d'une mousse dans des conditions d'utilisation proches de celles de la réalité. Les résultats obtenus dans cette recherche vérifient que le mécanisme de décomposition reste inchangé indépendamment de l'échelle. Dans la littérature, ces trois échelles n'ont jamais été considérées ensemble. Généralement, chaque échelle est considérée indépendamment et les chercheurs restent concentrés sur les phénomènes observés à l'échelle étudiée. De plus, les résultats de l'échelle matière sont souvent extrapolés à l'échelle produit. Toutefois, les phénomènes supplémentaires qui apparaissent entre une échelle et l'autre ne sont pas pris en compte, engendrant une grande incertitude dans la prédiction des résultats. Cette recherche propose une contribution vis-à-vis de l'intégration verticale des résultats obtenus dans les trois échelles. L'intégration verticale signifie explorer la possibilité d'identifier quelles propriétés de la matière doivent être mesurées et fournies en tant que données d'entrée des codes de simulation incendie afin de pouvoir prédire la décomposition thermique des solides. Ces travaux constituent un pas dans une vision globale de la science des matériaux qui permettrait une prédiction très juste du comportement au feu des solides à diverses échelles tout en utilisant principalement des mesures menées à l'échelle matière et la petite échelle. La cinétique de la décomposition a été étudié à la petite échelle grâce à des analyses thermogravimétriques (TGA). Cette technique a permis de mettre en évidence le nombre d'étapes, les espèces qui entrent en réaction et de détailler le mécanisme de réaction. En outre, des algorithmes génétiques ont été utilisés pour calculer les paramètres cinétiques optimum qui permettent de prédire le changement de la masse d'un échantillon en fonction de la température. Selon la démarche à échelle croissante décrite ci-dessus, les propriétés thermiques ainsi que les paramètres cinétiques de la 10 décomposition du PPUF ont été utilisés comme données d'entrée dans un code de simulation incendie. Les simulations ont été réalisées avec le code de calcul le plus amplement utilisé dans le monde, Fire Dynamics Simulator (FDS V 5.3). Les simulations tentent de prédire le comportement du PPUF en cône calorimètre (petite échelle). Un faible calage entre les courbes de changement de la masse expérimentales et numériques a été observé. Une grande incertitude vis-à-vis de la façon d'introduire les données d'entrée a été identifiée ainsi que de leur interprétation. Des possibles voies d'amélioration des modèles de pyrolyse ont été proposées.

Pour la conception des cœurs de réacteurs de 4ème génération, une précision accrue est requise pour les calculs des différents paramètres neutroniques. Les ressources mémoire et le temps de calcul étant limités, une solution consiste à utiliser des méthodes de raffinement de maillage afin de résoudre l'équation de transport des neutrons. Le flux neutronique, solution de cette équation, dépend de l'énergie, l'angle et l'espace. Les différentes variables sont discrétisées de manière successive. L'énergie avec une approche multigroupe, considérant les différentes grandeurs constantes sur chaque groupe, l'angle par une méthode de collocation, dite approximation Sn. Après discrétisation énergétique et angulaire, un système d'équations hyperboliques couplées ne dépendant plus que de la variable d'espace doit être résolu. Des éléments finis discontinus sont alors utilisés afin de permettre la mise en place de méthodes de raffinement dite hp. La précision de la solution peut alors être améliorée via un raffinement en espace (h-raffinement), consistant à subdiviser une cellule en sous-cellules, ou en ordre (p-raffinement) en augmentant l'ordre de la base de polynômes utilisée.Dans cette thèse, les propriétés de ces méthodes sont analysées et montrent l'importance de la régularité de la solution dans le choix du type de raffinement. Ainsi deux estimateurs d'erreurs permettant de mener le raffinement ont été utilisés. Le premier, suppose des hypothèses de régularité très fortes (solution analytique) alors que le second utilise seulement le fait que la solution est à variations bornées. La comparaison de ces deux estimateurs est faite sur des benchmarks dont on connaît la solution exacte grâce à des méthodes de solutions manufacturées. On peut ainsi analyser le comportement des estimateurs au regard de la régularité de la solution. Grâce à cette étude, une stratégie de raffinement hp utilisant ces deux estimateurs est proposée et comparée à d'autres méthodes rencontrées dans la littérature. L'ensemble des comparaisons est réalisé tant sur des cas simplifiés où l'on connaît la solution exacte que sur des cas réalistes issus de la physique des réacteurs.Ces méthodes adaptatives permettent de réduire considérablement l'empreinte mémoire et le temps de calcul. Afin d'essayer d'améliorer encore ces deux aspects, on propose d'utiliser des maillages différents par groupe d'énergie. En effet, l'allure spatiale du flux étant très dépendante du domaine énergétique, il n'y a a priori aucune raison d'utiliser la même décomposition spatiale. Une telle approche nous oblige à modifier les estimateurs initiaux afin de prendre en compte le couplage entre les différentes énergies. L'étude de ce couplage est réalisé de manière théorique et des solutions numériques sont proposées puis testées.

Cette thèse constitue une contribution à la modélisation numérique diphasique du transport sédimentaire en milieux côtiers et estuariens. L'approche diphasique diffère de l'approche dite classique en résolvant des équations de conservation de la masse, de la quantité de mouvement et de transport pour l'énergie cinétique turbulente séparément pour chaque phase. Les interactions entre le fluide et les particules sont alors représentées par des termes de transfert entre les deux phases. Le travail présenté dans cette thèse porte essentiellement sur la modélisation de la turbulence dans les écoulements diphasiques fluide particule.<br /><br />En terme de résultats, nous montrons que l'approche développée est capable de reproduire quantitativement les principaux processus mis en jeu dans le transport sédimentaire de particules non cohésives : la sédimentation et la dispersion turbulente des particules en milieu dilué. Le modèle développé confirme l'existence d'une différence de vitesse horizontale entre les particules et l'eau. Il simule les effets de dispersion des particules par le mouvement turbulent du fluide et l'atténuation de l'énergie cinétique turbulente du fluide due à la présence des particules. Une autre originalité de ce travail est de proposer un modèle diphasique à surface libre, bidimensionnel vertical, pour la simulation du transport sédimentaire. Nous avons identifié des lacunes entre les théories et les expériences notamment pour la simulation de la turbulence en écoulement dense. Nous proposons des solutions pour améliorer la simulation du comportement de matériaux cohésifs. Une tentative de simulation hydrosédimentaire sur l'estuaire de la Seine est présentée. Le phénomène de bouchon vaseux est qualitativement reproduit par le modèle sans qu'aucune loi d'érosion ou de dépôt ne soit imposée.

Reproduction de : Thèse de doctorat : Systèmes industriels : Toulouse, INPT : 2007. / Titre provenant de l'écran-titre. Bibliogr. 111 réf.