Les milieux granulaires présentent une grande diversité de comportements physiques que l'on peut relier aux différents états "conventionnels" de la matière : solide, liquide ou gaz. Les travaux décrits dans ce manuscrit se concentrent sur les écoulements granulaires denses présentant de nombreuses analogies avec les liquides moléculaires. Les expériences présentées dans ce manuscrit sont réalisées en utilisant des jets granulaires s'écoulant d'un réservoir, dans l'air et sous l'effet de la gravité. Quelques propriétés de tels jets sont étudiées, notamment l'évolution du diamètre au cours de la propagation. L'évolution de cette grandeur permet une description des jets granulaires comme des fluides incompressibles dans une certaine gamme de diamètres de sortie D et de diamètres des particules d, résumés en une grandeur réduite : le rapport d'aspect D/d. Lorsque ce dernier est inférieur à une grandeur critique proche de 15, les jets granulaires cessent de se comporter comme des liquides incompressibles et sont alors décrits comme des gaz en expansion.La transition entre les deux régimes d'écoulements est expliquée par un modèle faisant intervenir l'équilibre entre la température granulaire et sa dissipation par l'inélasticité des collisions entre les billes constituant le milieu granulaire. En utilisant les mesures de la température granulaire accessibles à la paroi dans le réservoir, ce modèle permet d'obtenir les profils de température et de fraction volumique du milieu en fonction du rapport d'aspect ainsi que leurs valeurs moyennes, présentant un bon accord avec les mesures de la fraction volumique en sortie du réservoir.Dans le régime "liquide" des jets granulaires mis en évidence ici, une tension de surface effective est mesurée par l'étude d'une instabilité similaire à celle de Rayleigh-Plateau. En soumettant le jet à une oscillation verticale apparaît une modulation du diamètre du jet dont il est possible d'étudier l'évolution au cours du temps. La mesure des taux de croissance des modes instables excités permet d'établir la relation de dispersion et d'obtenir une tension de surface effective dans les jets granulaires. Il est montré que cette tension de surface effective décroit lorsque la taille des particules augmente, décroit lorsque la pression de l'air ambiant diminue et augmente avec la distance de propagation.Le mécanisme proposé pour expliquer l'origine de cette tension de surface effective est l'interaction entre le jet granulaire et l'air ambiant. Les profils de vitesse de l'air entraîné par le jet sont mesurés pour différentes diamètres de particules, différentes pression ambiantes et différentes distances de propagation. L'entraînement de l'air et les profils de vitesses sont utilisés pour définir certaines forces s'exerçant sur les billes et construire une tension de surface effective : l'ordre de grandeur de cette tension de surface est incompatible avec les mesures réalisées précédemment. Une mesure de la pression de l'air à l'intérieur du jet granulaire permet d'obtenir une force exercée par la perméation de l'air à travers les pores du jet dont l'ordre de grandeur est correct. / Granular media exhibit a broad range of behaviours that can be related to the common states of matter: solid, liquid and gaseous. The experiments described in this thesis focus on dense granular flows ad their analogies with molecular liquid flows.The work presented focuses on granular material flowing from a reservoir, in ambient air under the effect of gravity. Some properties of such jets are presented, among them the evolution of the diameter of the jet along propagation. This evolution allows one to describe granular jets as incompressible liquid jets for a certain range of parameters D (exit diameter) and d (beads diameter). When using a reduced variable D/d (the aspect ratio), a transition appears close to a critical value (approx. 15) below which the granular jet exhibit a gas-like behaviour.This transition is explained using a model supposing balance between the heat flux (from the granular temperature imposed at the wall of the reservoir) and its dissipation through inelastic collisions between the beads in the granular medium. This model allows one to compute the temperature and volume fraction profiles as well as their mean values which are in good agreement with the experimental measurements (volume fraction at the exit).In the "liquid" regime defined previously, an effective surface tension is measured by studying the unstable modes of a Rayleigh-Plateau like instability. This instability is obtained by imposing a vertical oscillation on the jet which induces a diametrical modulation whose evolution is measured. Measurements of the growth rates of these unstable modes leads to the dispersion relation $\omega (kr_{0})$ and to an effective granular surface tension. This surface tension shows a diminution when the beads' size increases, a diminution when the air pressure is reduced and a growth along propagation.The proposed mechanism explaining this effective surface tension is the interaction between the granular jet and the surrounding air. Air velocity profiles are measured for different beads diameters, different air pressures and different distances below the outlet. Using these velocity profiles, various forces exerted on the beads are calculated, leading to a surface tension whose order of magnitude is lower than the one measured through the Rayleigh-Plateau instability. The right order of magnitude can be obtained by considering the air permeation through the pores on the surface the jet and the force exerted along the surface of said pores.

Nous étudions expérimentalement les instabilités hélicoïdales et chaotiques, des défauts d'extrusion de polymères fondus, sur deux polystyrènes et un polypropylène à l'aide de deux dispositifs expérimentaux: un rhéomètre capillaire et une extrudeuse monovis, munie d'une filière plate à parois transparentes, permettant d'effectuer des mesures de biréfringence d'écoulement et de vélocimétrie laser-Doppler. Dans ce dernier cas, on dispose d'une géométrie d'écoulement constituée d'un réservoir, suivi d'une contraction et d'un canal d'écoulement. En rhéométrie capillaire, nous avons étudié l'influence des conditions d'écoulement (diamètre, longueur, angle d'entrée de la filière, température...) sur le déclenchement de l'instabilité hélicoïdale. Les essais de biréfringence ont permis l'observation d'oscillations transverses d'isochromes dans le réservoir: lors de ce régime instable, l'écoulement principal a lieu de façon non symétrique dans le canal. On montre également que l'utilisation de "convergent doux" permet de repousser la déstabilisation de l'écoulement d'entrée à des débits sept fois plus importants que dans le cas d'une contraction brusque. On s'attache à mesurer, sur l'axe de symétrie du réservoir, des premières différences de contraintes normales et de taux de déformation dans les conditions critiques propres à l'initiation du défaut, par biréfringence d'écoulement et vélocimétrie laser-Doppler, respectivement. Des calculs viscoélastiques de l'écoulement par éléments finis ont également permis l'estimation de la répartition des contraintes élongationnelles autour de la contraction, dans les conditions critiques d'apparition de l'instabilité.

Les fluides pétroliers extraits de gisements ont des teneurs plus ou moins importantes en composés non souhaités tels que le dioxyde de carbone, l'hydrogène sulfuré ou autres produits soufrés. Ils doivent subir nombre de transformations et purifications. Le dimensionnement des unités correspondantes ne peut être bien optimisé que si l'on dispose de données fiables et précises d'équilibre de phases et de propriétés volumétriques et bien entendu d'une modélisation elle aussi fiable et précise. Le présent travail a été consacré en partie à des mesures de propriétés volumétriques sur trois mélanges binaires (éthane - hydrogène sulfuré, éthane - dioxyde de carbone et propane - hydrogène sulfuré). Ces mesures ont été réalisées à l'aide d'un appareillage, comportant un densimètre à tube vibrant (Anton Paar, modèle DMA 512), qui a été conçu et construit tout spécialement pour ce travail. L'étude menée sur la modélisation des propriétés volumétriques a permis de mettre en évidence l'inadéquation de l'usage classique des équations d'état cubiques pour représenter simultanément propriétés volumétriques et équilibres de phases. Parmi les équations d'état présentées, une attention particulière a toutefois été portée sur les équations d'état cubiques en raison de leur très grande utilisation dans le domaine pétrolier. Un moyen a été testé pour adapter les équations d'état cubiques à la représentation simultanée et satisfaisante des propriétés volumétriques et équilibres de phases : Il s'agit du couplage de l'équation d'état cubique de Soave-Redlich-Kwong avec une correction de volume par réseau de neurones.

Le comportement mécanique du polyéthylène est déterminé en traction et en fluage à l'aide du dispositif VidéoTraction. Les résultats les plus inattendus concernent la variation de volume : après une expansion élastique, on constate une compétition entre un phénomène de compaction et un effet de dilatation. D'autres techniques et une recherche des erreurs du dispositif d'extensométrie ont confirmé la tendance à la dilatation et la faible contribution de la compaction. Les propriétés de la matrice polymère et celles des microcavités ont été distinguées. La matière est le siège de deux processus volumiques opposés se compensant presque exactement : la perte de cristallinité et l'orientation des chaînes amorphes. Dès lors, la cavitation est prédominante. Un modèle prédisant l'intensité de la cavitation est développé à partir des données microstructurales. Enfin, le comportement mécanique global est décomposé en une partie intrinsèque (matrice polymère) et une partie volumique (cavitation).

L’étude des écoulements compressibles, notamment supersoniques, traversant les chambres de combustion de certains moteurs ramjet et scramjet, requiert la prise en compte de différents dispositifs complexes pour l’amélioration du processus de combustion et en particulier de leur stabilisation. La connaissance des interactions entre la turbulence, les effets de compressibilité, et les interaction fluide-solide dans ce type d’écoulement reste imparfaite. Ce travail de thèse est dédié à l’amélioration de notre compréhension de ce type d’écoulement dans un certain nombre de configurations canoniques par la biais de la simulation numérique directe. L’ensemble des simulations conduites s’appuie sur l’emploi d’un outil de simulation numérique haute fidélité : CREAMS (Compressible REActive Multi-species Solver) développé à l’Institut Pprime. Ce code de calcul met en oeuvre des schémas numériques d’ordre élevé : schéma Runge–Kutta d’ordre 3 pour l’intégration temporelle combiné à un schéma WENO d’ordre 7 et centré d’ordre 8 pour la discrétisation spatiale. Dans un premier temps, nous présentons une nouvelle méthode de frontières immergées pour le calcul d’écoulement d’un fluide visqueux compressible dans des géométries irrégulières. La méthode développée dans le cadre de cette thèse est basée sur la combinaison de l’approche appelée "Direct forcing" et celle de "Ghost-Point-Forcing". L’originalité de cette méthode réside dans sa capacité à simuler des écoulements subsoniques et supersoniques à différents nombres de Reynolds. L’examen de précision de cette méthode a permis d’établir un ordre supérieur à deux et sa robustesse est éprouvée par l’étude d’un bon nombre de cas tests. Dans un second temps, une configuration canonique idéalisée d’interaction choc-turbulence est étudiée pour mettre en lumière les mécanismes physiques fondamentaux caractéristiques du phénomène d’interaction entre une turbulence homogène isotrope et une onde de choc droite. Cette étude est complétée par une étude d’interaction choc-mélange scalaire pour étudier l’impact du choc normal sur le processus du mélange. Ce travail permet de mettre en place une base de données de résultats susceptibles d’être confrontés ultérieurement à des calculs basés sur l’emploi de modèles de turbulence. Enfin, nous nous sommes intéressés à l’effet des propriétés de transport moléculaire, en particulier celles de la viscosité volumique, sur le développement d’une couche de mélange impactée par un choc oblique. Les simulations réalisées dans cette configuration ont permis d’étudier la validité de l’hypothèse de Stokes consistant à négliger l’effet de la viscosité volumique. / The study of compressible flows, especially supersonic, passing through the combustion chambers of ramjet and scramjet engines, requires the consideration of various complex devices for improving the combustion processand in particular its stabilization. Indeed, the knowledge of the interactions between turbulence, compressibility effects, and fluid-solid interactions in this type of flow still remains imperfect. This thesis is dedicated to improving our understanding of this type of flow in a number of canonical flow configurations through direct numerical simulation. All the simulations that have been conducted are based on the use of a high-precision numerical simulations tool, called CREAMS (Compressible REActive Multi-species Solver), developed at the Pprime Institute. This computational solver makes use of high precision numerical schemes: a 3rd order Runge–Kutta scheme for time integration combined with a 7th order WENO and 8th order centered scheme for spatial discretization. In a first step of this study, we present a new immersed boundary method for calculating the flow of compressible viscous fluids in irregular geometries. The method developed in this thesis is based on the combination of the "Directforcing" approach with the "Ghost-Point-Forcing" strategy. The originality of this method lies in its ability to simulate subsonic and supersonic flows at different Reynolds numbers. The accuracy of this method is found to be slightly larger than second order and its robustness is investigated by considering a large set of benchmarks. Ina second step, an idealized canonical configuration of shock-turbulence interaction is studied to highlight the fundamental physical mechanisms that are characteristic of the interaction between an isotropic homogeneous turbulence and a normal shock wave. This analysis is complemented by a scalar shock-mixing interaction study to investigate the impact of normal shock on the mixing process properties. Through this work, a database is made available. It can be used to assess and improve turbulence models. Finally, we investigated the effect of molecular transport properties, more specifically the volume viscosity, on the development of a mixture layer impacted by an oblique shock. The simulations performed in this configuration allow to scrutinize the validity of the Stokes hypothesis that is based on the neglection of the volume viscosity.

Les écoulements autour de structures en mouvement font l'objet de plusieurs travaux numériques et expérimentaux. L'objectif de ce travail de thèse consiste à montrer la pertinence de la combinaison de la pénalisation volumique avec la méthode de lattice Boltzmann (LBM), dans l'étude du mouvement d'obstacles mobiles dans un écoulement, et de leur interaction avec celui-ci. La LBM,qui est simple et précise à mettre en œuvre, a prouvé ces dernières années son efficacité dans le domaine de la mécanique des fluides. Par ailleurs, la méthode de pénalisation volumique consiste à introduire un terme de pénalisation dans l'équation que l'on souhaite résoudre, afin de prendre en compte l'influence de l'obstacle sur le domaine fluide. Comme cette équation est résolue sur l'ensemble du domaine composé du fluide et du solide, les conditions aux limites à l'interface fluide-solide sont appliquées de façon naturelle. Il semble donc aisé de combiner cette technique avec la méthode de lattice Boltzmann. Nous avons dans un premier temps rappelé les notions de base et les principales caractéristiques de la méthode de lattice Boltzmann. On a présenté quelques exemples d'applications sur des cas tests, que nous avons programmés. Ensuite, une étude bibliographique faisant état des différentes approches qui utilisent la LBM dans l'étude des problèmes d'interaction fluide structure (IFS) a été réalisée. Puis, la combinaison de la pénalisation volumique avec la LBM a été testée avec succès sur l'équation de Burgers monodimensionnelle. La validation s'est portée en premier lieu, sur un écoulement autour d'un solide fixe, puis sur un écoulement autour d'une structure dont le mouvement est imposé, et finalement sur un problème d'IFS de type masse-ressort. La méthode développée a été ensuite testée sur les équations de Navier-Stokes, en considérant un fluide incompressible et une structure rigide. La validation s'est portée tout d'abord sur un écoulement autour d'obstacles immobiles (carré et cylindre), puis autour d'un cylindre mobile en oscillations forcées et libres. Enfin, une dernière application a été portée sur un écoulement entre deux plaques mobiles dans un canal. Nous avons montré que pour tous les cas étudiés, l'approche développée donne de bons résultats, elle reproduit avec précision les résultats de référence. / Flows around moving bodies are the subject of several numerical and experimental studies. The work presented in this document deals with the implementation of a volume penalization technique in a lattice Boltzmann model (LBM), in order to compute flows around moving obstacles. The LBM, which is accurate and easy to implement, has been successfully applied in fluid mechanics during the last decades. It was thus chosen in the present work, for flow computation. Furthermore, the volume penalization technique consists in introducing a volume penalization term into the equation that needs to be solved, in order to take into account the influence of the obstacle on the fluid domain. Since this equation is solved on both fluid and solid domains, the boundary conditions at the fluid-solid interface are naturally applied. Hence this technique seems easy to implement in a lattice Boltzmann framework. In the first chapter, the foundations and the main features of the lattice Boltzmann method are recalled, and several test cases that we simulated are presented. The second chapter deals with a literature review of the techniques developed for the simulation of fluid structure interaction problems in combination with the LBM. In the third chapter, the volume penalization method combined with the LBM was first applied to the one dimensional Burgers equation, considering motionless and moving obstacles (forced motion, and coupling between the fluid force calculated with the penalized Burgers equation and the motion of the obstacle). The combination of the volume penalization approach and the LBM was then employed to solve the incompressible NavierStokes equations, for cases of flows past motionless obstacles (flows over a square obstacle, and past a circular cylinder), and past an oscillating cylinder (where forced and free oscillations of the cylinder were simulated). Finally, this method was also applied to a symmetric Couette flow. For all these simulated cases, a good agreement with numerical results obtained with other techniques, and with results found in literature, was obtained.

Une étude expérimentale du développement de jets turbulents à masse volumique variable est menée. Les comportements axiaux et radiaux des propriétés statistiques moyennes et fluctuantes obtenues sur les champs de vitesse et de concentration sont analysés. Ces résultats apportent une meilleure compréhension des phénomènes de mélange et de l’effet de la variation de densité dans le but d’apporter des données complémentaires pour le développement des calculs numériques et accroître leurs précisions. Les cas étudiés dans cette thèse sont deux jets d’hélium marqués par de la vapeur d’acétone pour des nombres de Reynolds Re = 7000 et 11000 permettant d’explorer l’effet de fortes variations de masse volumique (S = 0,39 et 0,41, respectivement) par rapport à l’air ambiant. Le cas d’un jet d’air (Re = 16000) également porteur de vapeur d’acétone est utilisé comme cas de référence comparable au cas d’un contaminant passif avec S = 1,17. L’étude porte sur la région de proche sortie de jet, jusqu’à une distance de 40 fois le diamètre Dj de sortie du jet. La configuration adoptée est celle d’un jet de tube rond, libre, axisymétrique, vertical et ascendant dont le diamètre intérieur du tube est Dj = 3,5mm se développant dans l’air ambiant. La particularité de cette étude est la mise en œuvre d’un système de mesure par diagnostic optique qui permet un couplage spatial et temporel des mesures de vitesse et de concentration sur une région de l’écoulement de quelques cm2 et non sur un seul point. Pour ce faire, le couplage des mesures PIV, pour la mesure du champ de vitesse, et PLIF basée sur la fluorescence de la vapeur d’acétone, pour la mesure du champ de concentration, a été étudié, conçu et testé. Il nous a permis d’obtenir une base de données, encore trop rares aujourd’hui, sur l’évolution des flux turbulents croisés de vitesse et concentration. Ces données ont mis en évidence une évolution plus rapide du jet le plus léger. Cependant, les coefficients de corrélation semblent identiques lorsque l’on atteint la zone autosimilaire du jet. Une approche basée sur les probabilités de densité conjointes vitesse-concentration a permis de mettre en évidence des différences dans la région extérieure des jets, où l’intermittence de frontière donne son empreinte sur les propriétés du mélange. / An experimental study of the development of turbulent jets with variable density is presented. The axial and radial behaviours of the mean and fluctuating statistical properties obtained on the velocity and concentration fields are analysed. These results provide a better understanding of mixing phenomena and the effect of density variation in order to provide complementary data for the development of numerical calculations and to increase their precision. The cases studied in this thesis are two helium jets carrying acetone vapor for Reynolds numbers Re = 7000 and 11000 to explore the effect of large density variations (S = 0.39 and 0.41, respectively) relative to the ambient air. The case of an air jet (Re = 16000) also carrying acetone vapor is used as a reference case comparable to the case of a passive contaminant with S = 1.17. The study focuses on the region of near jet exit, up to a distance of 40 times the jet outlet diameter Dj. The adopted configuration is a round jet tube, free, axisymmetric, vertical and ascending whose internal diameter of the tube is Dj = 3.5 mm developing in the ambient air. The particularity of this study is the implementation of an optical diagnostic measurement system that allows a spatial and temporal coupling of speed and concentration measurements over a region of the flow of a few cm2. To do this, the coupling of the PIV measurements for the measurement of the velocity field and the PLIF based on the fluorescence of the acetone vapor for the measurement of the concentration field has been studied, designed and tested. It allowed us to obtain a database, still too rare today, on the evolution of the turbulent flow of speed and concentration. These data have shown a faster evolution of the lightest jet. However, the correlation coefficients appear to be identical when the self-similar zone of the jet is reached. An approach based on joint speed-concentration density probabilities has revealed differences in the outer region of the jets, where the intermittent boundary gives its imprint on the properties of the mixture.

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

La tomodensitométrie avec faisceau conique (CBCT) est actuellement utilisée en radiothérapie externe pour visualiser le patient dans la salle de traitement. Le mouvement respiratoire des patients y est encore difficilement pris en compte et des avancées sur le sujet pourraient améliorer la précision des traitements. En ce sens, l’obtention d’une séquence imageant les mouvements dans la région d’intérêt serait souhaitable. Ce mémoire présente le développement d’un algorithme de reconstruction 4D pour CBCT qui tente de répondre à certains besoins cliniques, soit l’obtention d’une qualité d’image suffisante, une facilité de mise en place clinique et une rapidité de calcul. L’algorithme 4D développé se base sur l’algorithme itératif convexe avec sous-ensembles ordonnés et régularisation de variation totale. Cette méthode a été choisie pour sa rapidité d’exécution, procurée par l’utilisation de sous-ensembles et la parallélisation des calculs sur carte graphique, et pour sa capacité à diminuer les artéfacts de rayons, communs en imagerie 4D, procurée par la régularisation de variation totale. La méthode développée pour obtenir une image 4D à partir d’acquisitions CBCT standards a fait appel à l’algorithme Amsterdam Shroud pour déduire le mouvement respiratoire de l’ensemble de projections CBCT. Elle a été validée sur un fantôme numérique et sur des acquisitions cliniques. Les résultats obtenus démontrent le potentiel de l’algorithme, puisqu’une image de résolution spatiale et temporelle satisfaisante a été reconstruite en moins de 5 minutes. Un tel temps de calcul se compare avantageusement à d’autres méthodes disponibles et ouvre la porte à une visualisation rapide du mouvement respiratoire en salle de traitement. / Cone beam computed tomography (CBCT) is currently used to visualize patients directly in the treatment room. However, the respiratory movement is still hardly taken into account and new developments could improve the precision of treatment. To this end, obtaining a film imaging movements in the region of interest would be beneficial. This master’s thesis presents the development of a reconstruction algorithm for 4D CBCT which seeks to respond to particular clinical needs, namely sufficient image quality, clinical implementation simplicity and high computational speed. The developed 4D algorithm is based on the ordered subsets convex iterative algorithm combined with the total variation minimization regularization technique. This method was chosen for its fast execution time, enabled by the use of subsets and the parallelization on GPU, and for its capability to reduce streaking artifacts, common on 4D imaging, enabled by the total variation regularization. The method developed to reconstruct a 4D image from standard CBCT scans employed the Amsterdam Shroud algorithm to deduce respiratory movement of a CBCT projections’ set. Its validation was performed on a numerical phantom and on clinical datasets. Results demonstrate the potential of the algorithm, since an image with sufficient spatial and temporal resolution was reconstructed in less than 5 minutes. Such computational times can be compared favorably with other available methods and could allow for online applications.

La présente thèse s’inscrit dans le domaine de la physique médicale et, plus précisément, de l’imagerie médicale tridimensionnelle (3D) et de la dosimétrie 3D pour la radiothérapie. L’objectif global du travail était de concevoir et évaluer un nouvel algorithme de reconstruction itératif rapide pour la tomodensitométrie (TDM) à faisceau conique, une modalité consistant à créer des images 3D des densités du sujet imagé à partir de mesures d’atténuation partielle d’un faisceau de radiation incidente. Cet algorithme a été implanté sur matériel graphique (GPU), une plate-forme de calcul hautement parallèle, menant à la conception de stratégies d’optimisation originales. En premier lieu, un nouvel algorithme itératif statistique régularisé, dénommé OSC-TV, a été conçu et implanté sur GPU. Il a été évalué sur des ensembles de projections synthétiques et cliniques de TDM à rayons X à faisceau conique. L’algorithme proposé a démontré une qualité d’image supérieure à celle de méthodes semblables pour des acquisitions basse-dose, ainsi que des temps de reconstruction compatibles avec les activités cliniques. L’impact principal de ce travail est la capacité d’offrir au patient une réduction de dose de radiation ionisante de deux à quatre fois par rapport aux protocoles d’acquisition usuels. En second lieu, cet algorithme a été testé sur des données expérimentales en tomographie optique à faisceau conique, donnant lieu à l’une des premières études de ce genre. La résolution spatiale des images 3D résultantes a été améliorée et le bruit a été réduit. L’on a aussi démontré l’importance de considérer le spectre de la source lumineuse afin d’assurer la justesse de l’estimation des densités. Le principal impact de l’étude est la démonstration de la supériorité de la reconstruction itérative pour des données affectées par les aberrations propres à la tomographie optique à faisceau conique, résultant potentiellement en l’amélioration de la dosimétrie 3D par gel radiochromique en radiothérapie. En troisième lieu, différentes approches de gestion de la matrice-système de type exact à rayons fins ont été évaluées pour la TDM à faisceau conique. Le pré-calcul et le stockage complet de la matrice-système dans la mémoire vive du GPU s’est montré comme l’approche la plus rapide, mais la moins flexible en termes de géométries représentables, en raison de la taille limitée de la mémoire vive. Le traçage de rayons à la volée est apparu très flexible, offrant aussi des temps de reconstruction raisonnables. En somme, les trois études ont permis de mettre en place et d’évaluer la méthode de reconstruction proposée pour deux modalités de tomographie, ainsi que de comparer différentes façons de gérer la matrice-système. / This thesis relates to the field of medical physics, in particular, three-dimensional (3D) imaging and 3D dosimetry for radiotherapy. The global purpose of the work was to design and evaluate a new fast iterative reconstruction algorithm for cone beam computed tomography (CT), an imaging technique used to create 3D maps of subject densities based on measurements of partial attenuation of a radiation beam. This algorithm was implemented for graphics processing units (GPU), a highly parallel computing platform, resulting in original optimization strategies. First, a new iterative regularized statistical method, dubbed OSC-TV, was designed and implemented for the GPU. It was evaluated on synthetic and clinical X ray cone beam CT data. The proposed algorithm yielded improved image quality in comparison with similar methods for low-dose acquisitions, as well as reconstruction times compatible with the clinical workflow. The main impact of this work is the capacity to reduce ionizing radiation dose to the patient by a factor of two to four, when compared to standard imaging protocols. Second, this algorithm was evaluated on experimental data from a cone beam optical tomography device, yielding one of the first studies of this kind. The spatial resolution of the resulting 3D images was improved, while the noise was reduced. The spectral properties of the light source were shown to be a key factor to take into consideration to ensure accurate density quantification. The main impact of the study was the demonstration of the superiority of iterative reconstruction for data affected by aberrations proper to cone beam optical tomography, resulting in a potential to improve 3D radiochromic gel dosimetry in radiotherapy. Third, different methods to handle an exact thin-ray system matrix were evaluated for the cone beam CT geometry. Using a GPU implementation, a fully pre-computed and stored system matrix yielded the fastest reconstructions, while being less flexible in terms of possible CT geometries, due to limited GPU memory capacity. On-the-fly ray-tracing was shown to be most flexible, while still yielding reasonable reconstruction times. Overall, the three studies resulted in the design and evaluation of the proposed reconstruction method for two tomographic modalities, as well as a comparison of the system matrix handling methods.