Le déminage mécanique consiste à scarifier le sol avec une " charrue ", outil formé d'une lame en V munie de plusieurs dents, poussée par un engin : les dents déstructurent le sol et le font remonter devant la lame qui l'évacue sur le côté avec les mines qu'il contient. L'objectif de nos travaux est de mettre en œuvre la simulation numérique par éléments finis du problème fortement non linéaire issu de la modélisation de la scarification du sol. Le code d'éléments finis implicite Forge3®, dédié à la mise en forme des métaux, a été choisi comme support numérique. Il permet de prendre en compte les grandes déformations, en particulier à l'aide de son remaillage automatique. Nous avons dans un premier temps implanté dans ce code deux modèles élastoplastiques de comportement, l'un incompressible réservé aux sols fins saturés, purement cohésifs, l'autre compressible, fondé sur la notion d'état critique, pour les matériaux purement frottants ou cohésifs et frottants. Ces modèles adoucissants sont intégrés par un schéma de retour radial généralisé, au sein d'une formulation implicite du problème aux limites. Nous montrons que la matrice de raideur est non symétrique dans le cas compressible et que la symétrisation du système ne conduit pas à une approche robuste. Nous avons donc importé et testé un solveur itératif non-symétrique : Bi-CGSTAB. Nous avons validé la programmation de ces modèles sur des essais triaxiaux. Pour les comportements adoucissants, on constate des oscillations dans la relation contrainte/déformation passé le pic de contrainte. Ces difficultés numériques sont traitées par linéarisation et régularisation. Dans un deuxième temps, nous avons mis en œuvre des simulations de scarification pour différents niveaux de complexité : une dent seule, une dent + une tranche de lame, plusieurs dents, en faisant suffisamment avancer l'outil pour atteindre le régime stationnaire ; cela se révèle d'autant plus long que l'outil est large, ce qui nous amène au calcul très intensif. L'influence de paramètres géométriques comme l'inclinaison de la dent ou l'angle d'étrave du système a été mise en évidence, elle est qualitativement conforme aux observations expérimentales. L'étude de l'influence des paramètres des modèles de comportement montre la prépondérance de ceux liés à l'état critique, i.e. aux propriétés mécaniques après de grandes déformations. Enfin nous avons validé qualitativement le code en termes de modes d'écoulements et de répartition des efforts pour des outils multi-dents. La comparaison quantitative des efforts reste à affiner en revenant sur le comportement choisi et son implémentation.

Le travail présenté dans cette thèse a été motivé par des problématiques posées par l'évaluation de contrats échangés sur le marché du gaz: les contrats de stockage et d'approvisionnement en gaz. Ceux-ci incorporent de l'optionalité et des contraintes, ce qui rend leur évaluation difficile dans un contexte de prix de matières premières aléatoires. L'évaluation de ces contrats mène à des problèmes de contrôle stochastique complexes: switching optimal ou contrôle impulsionnel et contrôle stochastique en grande dimension. La première partie de cette thèse est une revue relativement exhaustive de la littérature, mettant en perspective les différentes approches d'évaluation existantes. Dans une deuxième partie, nous considérons une méthode numérique de résolution de problèmes de contrôle impulsionnel basée sur leur représentation comme solution d'EDSRs à sauts contraints. Nous proposons une approximation à temps discret utilisant une pénalisation pour traiter la contrainte et donnons un taux de convergence de l'erreur introduite. Combinée avec des techniques Monte Carlo, cette méthode a été testée numériquement sur trois problèmes: gestion optimale de biomasse, évaluation d'options Swing et de contrats de stockage gaz. Dans une troisième partie, nous proposons une méthode pour l'évaluation d'options dont le payoff dépend de moyennes mobiles de prix sous-jacents. Elle utilise sur une approximation à dimension finie de la dynamique des processus de moyenne mobile, basée sur un développement en série de Laguerre tronquée. Les résultats numériques fournis incluent des exemples de contrats Swing gaziers à prix d'exercice indexés sur moyennes mobiles de prix pétroliers.

Cette thèse est consacrée à la conception et l'analyse théorique de différents composants en optique intégrée. Nous présentons un modèle de milieu effectif pour les cristaux photoniques (CPs) 2D qui rend compte des effets de surface, puis un algorithme pour la conception de réseaux antiréfléchissants grand-angle pour ces CPs. Ces réseaux permettent d'améliorer significativement la transmission à travers une lentille plate d'indice négatif. Nous proposons une nouvelle génération de circulateurs magnétooptiques compacts, fonctionnant dans un champ magnétique extérieur uniforme et constitués d'une cavité résonnante en anneaux circulaires couplée directement à des guides d'ondes standards. Nous généralisons la méthode multipolaire 2D aux matériaux gyrotropiques et la formulons sans « lattice sums » pour les structures périodiques. Enfin, nous décrivons en détail la méthode des éléments finis pour le calcul des modes propres des cavités 3D en anneaux circulaires et matériaux gyrotropiques.

In this thesis, we construct a new hydrodynamical algorithm able of handling general compressible reactive flow problems, based on a finite-volume method inspired by the original MUSCL scheme of van Leer (1979). The algorithm is of second-order in the smooth part of the flow and avoids dimensional splitting. It uses MPI to achieve parallelism, and includes an astrophysical equation of state and a nuclear reaction network. It proves to be robust to tests cases. In particular it reproduces quite well the reactive and non-reactive results obtained with two different numerical methods (Fryxell & al. 1989, Busegnies & al. 2007). Moreover the time-dependent results are in agreement with the corresponding steady state solution. This gives us confidence in applying it to an astrophysical situation which has never been studied, the propagation of a detonation in conditions relevant to superbursts. The algorithm is described in (Noel & al. 2007). In a firt step we obtain the detonation profiles in pure carbon and in a mixture of carbon and iron. In both cases we underline the large difference between the total reaction length and the length on which some species burn. This difference leads to enormous numerical difficulties because all the length scales cannot be resolved at the same time in a single simulation. We show that the carbon detonation might be studied in a partial resolution approach like the one of Gamezo & al. (1999). In a second step we construct a new reduced nuclear reaction network able to reproduce the energy production due to the photo-disintegrations of heavy elements, like ruthenium, which are thought to occur during superbursts in mixed H/He accreting systems. Using this new nuclear network we simulate detonations in mixture of carbon and ruthenium. An interesting feature is that, in this case, all the reaction lengths can be resolved in the same simulation. This makes the C/Ru detonations easier to study in future multi-dimensional simulations than the pure carbon ones (Noel & al. 2007b). Finally we perform some numerical experiments which show that our algorithm is able to deal with initially inhomogeneous medium, and that the multi-dimensional simulations are attainable even if they are quite computational time consuming. - B. Van Leer, J. Comp. Phys., 21, 101, 1979 - Fryxell, B.A., Muller, E., and Arnett, W.D., Technical report MPA 449, 1989 - Busegnies, Y., Francois, J. and Paulus, G., Shock Waves, 11, 2007 - Gamezo, V.N., Wheeler, J.C., Khokhlov, A.M., and Oran, E.S., ApJ, 512, 827, 1999 - Noël, C., Busegnies, Y., Papalexandris, M.V. & al., A&A, 470, 653, 2007 - Noël, C., Goriely, S., Busegnies, Y. & Papalexandris, M.V., submitted to A&A, 2007b / Un algorithme parallèle basé sur une méthode aux volumes finis inspirée du schéma MUSCL de Van Leer (1979) a été construit. Il a été développé sur base de la méthode de Lappas & al. (1999) qui permet de résoudre simultanément toutes les dimensions spatiales. Cette méthode se base sur la construction de surfaces appropriées dans l'espace-temps, le long desquelles les équations de bilan se découplent en équations plus simples à intégrer. Cet algorithme est actuellement le seul à éviter le "splitting" des dimensions spatiales. Dans les modèles conventionnels (PPM, FCT, etc.), l'intégration spatiale des équations est réalisée de manière unidimensionnelle pour chaque direction. Un réseau de réactions nucléaires ainsi qu'une équation d'état astrophysique ont été inclus dans l'algorithme et celui-ci a ensuite été soumis à une grande variété de cas tests réactifs et non réactifs. Il a été comparé à d'autres codes généralement utilisés en astrophysique (Fryxell & al. 1989, Fryxell & al. 2000, Busegnies & al. 2007) et il reproduit correctement leurs résultats. L'algorithme est décrit dans Noël & al. (2007). Sur base de cet algorithme, les premières simulations de détonation dans des conditions thermodynamiques représentatives des Superbursts ont été réalisées. Différentes compositions du milieu ont été envisagées (carbone pur, mélange de carbone et de fer, mélange de carbone et de cendres du processus rp). Dans la plupart des systèmes où des Superbursts ont été observés, la matière accrétée est un mélange d'hydrogène et d'hélium. Dans ce cas, des phases de combustion précédant le Superburst produisent des nucléides plus lourd que le fer (Schatz & al. 2003). Ces nucléides peuvent être photodésintégrés durant le Superburst. Pour prendre en compte ces réactions endothermiques de photodésintégration, nous avons construit un nouveau réseau réduit de réactions nucléaires qui a été incorporé dans l'algorithme hydrodynamique (Noël & al. 2007b). Ce réseau réduit reproduit globalement l'énergétique d'un réseau complet et a permis de faire la première simulation numérique de détonation dans des conditions caractéristiques de systèmes accréteurs d'un mélange hydrogène-hélium. Finallement quelques simulations multidimensionelles préliminaires ont éte réalisées. - Busegnies, Y., Francois, J. and Paulus, G., Shock Waves, 11, 2007 - Fryxell, B.A., Muller, E., and Arnett, W.D., Technical report MPA 449, 1989 - Fryxell, B.A., Olson, K., Ricker, P. & al., ApJS, 131, 273, 2000 - Lappas, T., Leonard, A. and Dimotakis, P.E., SIAM J. Sci. Comput., 20, 1064, 1999 - Noël, C., Busegnies, Y., Papalexandris, M.V. & al., A&A, 470, 653, 2007 - Noël, C., Goriely, S., Busegnies, Y. & Papalexandris, M.V., submitted to A&A, 2007b - Röpke, F. K. PhD thesis, Technischen Universitat Munchen, 2003 - Schatz, H., Bildsten, L., Cumming, A. and Ouellette, M., Nuclear Physics A, 718, 247, 2003 - Van Leer, B. Comp. Phys., 21, 101, 1979 - Weinberg, N.N. and Bildsten, L., ArXiv e-prints, 0706.3062, 2007

Les défauts d'aspect des ouvrants automobiles sont des perturbations géométriques de la surface extérieure dont les dimensions caractéristiques ont une profondeur de l'ordre de quelques dizaines de micromètres et une étendue de quelques dizaines de millimètres. Ils entrainent des irrégularités de réflexion de la lumière sur la surface de la carrosserie et influent fortement sur la qualité perçue du véhicule. Le terme " ouvrants automobiles " regroupe les portes, le capot moteur et le coffre à bagages, pièces propices à l'apparition des défauts d'aspect à cause de leurs caractéristiques géométriques tourmentées. Cette étude s'appuie sur des essais expérimentaux pour caractériser les défauts d'aspect et sur la simulation par éléments finis pour les prédire numériquement. Les matériaux étudiés sont des aciers d'épaisseur comprise entre 0,6 et 0,8 mm. Leur carac- térisation mécanique est réalisée à partir d'essais de traction uniaxiale et de cisaillement simple (monotone et Bauschinger), avec pour objectif d'identifier les paramètres matériau des lois de comportement élastoplastiques utilisées dans la simulation numérique. Un dispositif de laboratoire a été conçu pour reproduire à petite échelle le défaut d'aspect apparaissant lors de l'opération de tombage d'un coin de porte. Sa modularité permet de faire varier différents paramètres comme la hauteur et le rayon de tombage. La géométrie des éprou- vettes tombées est mesurée avec une machine à mesurer tridimensionnelle. Une nouvelle méthode de détection et de caractérisation expérimentale est proposée à partir d'un calcul de courbure de sections. Cette méthode, moins subjective que la méthode traditionnelle par déméritage, permet de caractériser les défauts d'aspect au travers d'une cartographie des profondeurs sur des sections de la surface étudiée par des plans sécants. On montre que la simulation numérique reproduit l'allure du défaut d'aspect de coin de porte à condition d'utiliser un maillage fin, mais minimise sa profondeur, quel que soit le type d'élément fini utilisé. La méthode a été validée sur un coin de porte industriel présentant un défaut d'aspect lors de l'opération de tombage. Le défaut est observé qualitativement et mesuré par digitalisation optique. On montre que l'allure du défaut est reproduite par la simulation numérique sous réserve d'un maillage raffiné et de la prise en compte de toute l'histoire de déformation. Enfin, les défauts d'aspect apparaissant autour d'une poignée de porte ont aussi été caractérisés par digitalisation optique. La simulation numérique de ce défaut récurrent sur les automobiles, sous les mêmes conditions que pour le coin de porte, permet de prédire la majorité des défauts d'aspect observés.

.

La compréhension de l'état turbulent d'un plasma de fusion, responsable du faible temps de confinement observé, constitue un enjeu fondamental vers la production d'énergie par cette voie. Pour les machines les plus performantes, les tokamaks, les conductivités thermiques ionique et électronique mesurées sont du même ordre de grandeur. Les sources potentielles de la turbulence sont les forts gradients de température, de densité,... présents au coeur d'un plasma de tokamak. Si les pertes de chaleur par le canal ionique sont relativement bien comprises, l'origine du fort transport de chaleur électronique est quant à elle largement inconnue. En plus des fluctuations de vitesses électrostatiques, il existe des fluctuations de vitesses magnétiques, auxquelles des particules rapides sont particulièrement sensibles. Expérimentalement, le temps de confinement peut s'exprimer en fonction de paramètres non adimensionnels. Ces lois d'échelle sont encore trop imprécises, néanmoins de fortes dépendances en fonction du rapport de la pression cinétique à la pression magnétique, β et du rayon de Larmor normalisé, ρ* sont prédites.<br /><br />La thèse proposée ici cherche à déterminer la pertinence d'un modèle fluide non linéaire, électromagnétique, tridimensionnel, basé sur une instabilité particulière pour décrire les pertes de chaleur par le canal électronique et de déterminer les dépendances du transport turbulent associé en fonction de paramètres adimensionnels, dont β et ρ*. L'instabilité choisie est une instabilité d'échange générée par le gradient de température électronique (Electron Temperature Gradient (ETG) driven turbulence en anglais). Ce modèle non linéaire est construit à partir des équations de Braginskii. Le code de simulation développé est global au sens qu'un flux de chaleur entrant est imposé, laissant les gradients libres d'évoluer.<br /><br />A partir des simulations non linéaires, nous avons pu mettre en évidence trois caractéristiques principales pour le modèle ETG fluide: le transport de chaleur turbulente est essentiellement électrostatique; les fluctuations de potentiel et de pression forment des structures radialement allongées; le niveau de transport observé est beaucoup plus faible que celui mesuré expérimentalement.<br /><br />L'étude de la dépendance du transport de chaleur en fonction du rapport de la pression cinétique à la pression magnétique a montré un faible impact de ce paramètre mettant ainsi en défaut la loi empirique d'Ohkawa. En revanche, il a été montré sans ambiguïté le rôle important du rayon de Larmor électronique normalisé dans le tranport de chaleur: le temps de confinement est inversement proportionnel à ce paramètre. Enfin, la faible dépendance du transport de chaleur turbulent en fonction du cisaillement magnétique et de l'inverse du rapport d'aspect a été mise en évidence.<br /><br />Bien que le niveau de transport observé dans les simulations soit plus faible que celui mesuré expérimentalement, nous avons tenté une confrontation directe avec un choc de Tore Supra. Ce tokamak est particulièrement bien désigné pour étudier les pertes de chaleur électronique. En conservant la plupart des paramètres d'un choc bien référencé de Tore Supra, la simulation non linéaire obtenue donne un seuil en gradient de température proche de la valeur expérimentale. Le niveau de transport observé est plus faible d'un facteur cinquante environ que le transport mesuré. Un paramètre important qui n'a pu être conservé est le rayon de Larmor normalisé.<br /><br />La limitation en ρ* devra être franchie afin de confirmer ces résultats. Enfin une rigoureuse confrontation avec des simulations girocinétiques permettra de disqualifier ou non l'instabilité ETG pour rendre compte des pertes de chaleur observées.<br /><br />Mots-clés: fusion thermonucléaire, tokamak, plasma, turbulence ETG, simulations numériques

Les écoulements diphasiques en milieu poreux sont des phénomènes complexes et qui concernent de nombreux problèmes industriels. EDF travaille sur la faisabilité et la sécurité d'un stockage en couche géologique profonde de déchets nucléaires. Dans ce domaine la simulation des écoulements diphasiques en milieu poreux est particulièrement importante dans au moins trois domaines : tout d'abord lors de la phase de ventilation des galeries du stockage qui pourrait désaturer la roche présente et ainsi en modifier ses propriétés de rétention, mais également lors de la phase de resaturation des matériaux et enfin lors de l'arrivée de l'eau sur les parties métalliques contenues dans le stockage qui entraînera alors des phénomènes de corrosion et un dégagement d'hydrogène. Dans ce contexte, EDF souhaite se doter de méthodes numériques performantes et robustes ne nécessitant pas de conditions restrictives sur la forme des mailles. Ce travail s'inscrivant dans cette problématique, est consacré dans un premier temps au développement du schéma volumes finis SUSHI (Scheme Using Stabilization and Hybrid Interfaces) dans le code de mécanique d'EDF, Code_Aster afin de modéliser les écoulements diphasique en milieu poreux. Ce schéma a été développé en 2D et en 3D. Parallèlement une nouvelle formulation qui permet de traiter de manière uniforme les écoulements en milieu saturé et insaturé pour des problèmes miscibles et immiscibles est proposée. Différentes études modélisant des difficultés liées aux problématiques du stockage de déchets radioactifs en couches géologiques profondes ont été traitées. On peut citer l'étude d'un bi-matériau qui met en avant le ré-équilibrage capillaire d'un matériau par un autre possédant des propriétés et des conditions initiales en saturation très hétérogènes. On citera également l'étude de l'injection d'hydrogène dans un milieu initialement saturé en eau pure qui est tirée du benchmark « Ecoulement diphasique » proposé par le GNR MOMAS. Cette étude avait pour objectif de mettre en évidence le bon traitement de l'apparition d'une phase dans un milieu saturé et donc la pertinence de notre nouvelle formulation à traiter d'une manière unifié un problème d'écoulement saturé et un problème d'écoulement insaturé / The two-phase flow in porous media is a complex phenomenon and which relate to many industrial problems. EDF works on the feasibility and the safety of a storage in deep geologic layer of nuclear waste. In this domain the simulation of the two-phase flow in porous media is particularly important in at least three domains : first of all during the phase of ventilation of the galleries of the storage which could desaturate the rock and so modify its properties, but also during the phase of resaturation of the materials and finally during the arrival of the water on the metal parts contained in the storage which will then involve phenomena of corrosion and a hydrogen release. In this context, EDF wishes to obtain robust numerical methods without restrictive condition on the mesh. This work is dedicated at first to the development of the finite volume scheme SUSHI (Scheme Using Stabilization and Hybrid Interfaces) in the code of mechanics of EDF, Code_Aster in order to simulate the two-phase flow in porous media. This scheme was developed in 2D and in 3D. At the same time a new formulation which allows to simulate in a uniform way the flows in saturated and unsaturated porous media for miscible and immiscible problems is proposed. Various studies simulating difficulties related to the problems of the storage of nuclear waste in deep geological layers were study. We can quote the study of a bi-material which advances the capillary rebalancing of a material by an other one possessing properties and initial very heterogeneous conditions in saturation. We will also quote the study of the injection of hydrogen in an porous media initially saturated in pure water which is proposed by the benchmark "two-phase Flow " proposed by the GNR MOMAS. This study had for objective to bring to light the good treatment of the appearance of a phase in a saturated porous media and thus the relevance of our new formulation to study with a way unified a problem of saturated flow and a problem of unsaturated flow

Depuis quinze ans, de multiples instruments de mesure ont été envoyés vers Mars et explorent aujourd'hui sa surface ou y sont en orbite. Les données, renvoyées par ces sondes ont fourni une grande quantité d'informations quant aux conditions régnant aujourd'hui sur Mars, mais aussi sur la géologie des terrains et sur la minéralogie associée. Les minéraux découverts sur Mars se sont formés dans des conditions tout à fait différentes de celles qui y règnent aujourd'hui.Une approche associant modélisations numériques et expériences est proposée pour permettre d'identifier l'influence de divers processus comme le rayonnement UVC, la température et la pression partielle de CO2 ou de SO2. Des expériences en boîte à gants sont réalisées à des pressions partielles de CO2 de 10-6 bar, 10-3 bar et 10-2 bar, en présence ou non de rayonnement UVC. Conjointement à ces expériences, des expériences en capsules en or sont menées pour des températures de 20°C et 50°C avec des pressions partielles de CO2 d'environ 10-3 bar, 10-2 bar et 1 bar ainsi qu'avec une pression partielle de SO2 d'environ 1 bar. Les résultats ont montré qu'il était possible de reproduire expérimentalement et par les codes de calculs, la minéralogie supposée résultant de l'altération d'un verre basaltique de composition proche de roches martiennes. L'observation de carbonates de calcium dans les expériences réalisées à des pressions partielles de CO2 supérieures à 10-3 bar est en accord avec la prédiction de ces mêmes minéraux par des simulations numériques réalisées avec le logiciel DISSOL. Néanmoins, l'abondance de ces carbonates dans les deux types de modélisations, ainsi que leur composition à dominance calcique, contrastent avec les observations des rares minéraux carbonatés magnésiens à la surface de Mars. Ce contraste pose les questions concernant la composition de l'atmosphère, la composition initiale du verre et la préservation des carbonates après leur précipitation. Bien que des argiles n'aient pas été observées dans les expériences à cause d'un taux de dissolution du verre trop faible, des montmorillonites ferrifères et magnésiennes sont prédites dans les simulations numériques quelle que soit la pression partielle de CO2. De la beidellite sodique est également prédite dans le cas où les conditions atmosphériques martiennes seraient réductrices quelle que soit la pression partielle de CO2. Contrairement aux observations de la surface de Mars, aucune nontronite ne précipite. Les sulfates sont absents des paragenèses minérales obtenues lors des expériences avec une atmosphère à CO2. Dans les modélisations, aucun minéral sulfaté n'est néoformé pour les pressions partielles de CO2 imposées (> 10-6 bar) et la teneur en soufre du verre expérimental de synthèse. Toutefois, la précipitation de gypse peut être observée dans les simulations si la teneur en soufre est 100 fois supérieure à celle du verre utilisé pour ce travail. L'absence de sulfates magnésiens, de même que de carbonates magnésiens suggère que la composition chimique choisie pour le verre soumis à l'altération ne correspond pas totalement celle de la roche qui a pu donner les différents assemblages minéralogiques observés à la surface de Mars. En revanche, des sulfates de calcium et, des sulfates de fer et d'aluminium sont observés pour des expériences en capsules avec une atmosphère à SO2. Le rôle des atmosphères est essentiel quant à l'abondance des minéraux carbonatés et des minéraux sulfatés. Les expériences montrent que les radiations UVC accélèrent la dissolution du verre et favorise l'oxydation partielle du soufre du verre en pyrite / Process of deterioration on the primitive surface of Mars: experimental and digital simulations and mineralogical implications

Introduction: Le but de l’étude était d’examiner l’effet des matériaux à empreintes sur la précision et la fiabilité des modèles d’études numériques. Méthodes: Vingt-cinq paires de modèles en plâtre ont été choisies au hasard parmi les dossiers de la clinique d’orthodontie de l’Université de Montréal. Une empreinte en alginate (Kromopan 100), une empreinte en substitut d’alginate (Alginot), et une empreinte en PVS (Aquasil) ont été prises de chaque arcade pour tous les patients. Les empreintes ont été envoyées chez Orthobyte pour la coulée des modèles en plâtre et la numérisation des modèles numériques. Les analyses de Bolton 6 et 12, leurs mesures constituantes, le surplomb vertical (overbite), le surplomb horizontal (overjet) et la longueur d’arcade ont été utilisés pour comparaisons. Résultats : La corrélation entre mesures répétées était de bonne à excellente pour les modèles en plâtre et pour les modèles numériques. La tendance voulait que les mesures répétées sur les modèles en plâtre furent plus fiables. Il existait des différences statistiquement significatives pour l’analyse de Bolton 12, pour la longueur d’arcade mandibulaire, et pour le chevauchement mandibulaire, ce pour tous les matériaux à empreintes. La tendance observée fut que les mesures sur les modèles en plâtre étaient plus petites pour l’analyse de Bolton 12 mais plus grandes pour la longueur d’arcade et pour le chevauchement mandibulaire. Malgré les différences statistiquement significatives trouvées, ces différences n’avaient aucune signification clinique. Conclusions : La précision et la fiabilité du logiciel pour l’analyse complète des modèles numériques sont cliniquement acceptables quand on les compare avec les résultats de l’analyse traditionnelle sur modèles en plâtre. / Introduction: The purpose of this study was to examine the effects of different impression materials on the accuracy and reliability of digital models. Methods: Models from 25 patients selected at random from the files of the Department of Orthodontics at the University of Montreal were used in this study. One alginate (Kromopan 100 alginate, Italy), 1 alginate alternative (Alginot, Kerr Dentistry, Orange, CA), and 1 PVS (Aquasil, Dentsply Caulk, Milford, DE) impression was taken of both arches of each patient and sent to Orthobyte (Othobyte Digital Technology inc., Calgary, AB) for fabrication of a plaster model and scanning for production of a digital model. The Bolton 6 and 12 analyses and their constituent measurements, overbite, overjet, and arch length were used for the comparison. Results: The repeatability of measurements using both the plaster and digital methods was good to excellent, with the plaster measurements tending to be more reliable. There were statistically significant differences in the Bolton 12 and mandibular arch length and spacing measurements for all impression materials, with the plaster models tending to give a smaller measurement for the Bolton 12 and a higher measurement for the mandibular arch length and spacing. Although statistically significant differences in some measurements were found for the reliability and validity of digital models, none was clinically significant. Conclusions: The accuracy and reliability of the software for comprehensive cast analysis is clinically acceptable and reproducible when compared with traditional plaster study model analysis.