La nature hétérogène et parfois karstifiée des milieux carbonatés a un impact direct sur nombre de problématiques environnementales, comme celle des aquifères côtiers affectés par une intrusion d'eau salée. Au sud de l'ile de Majorque (Baléares, Espagne), la plate-forme carbonatée de Llucmajor offre un site privilégié pour l'étude de telles intrusions notamment avec le site expérimental de Campos, qui rencontre l'intrusion à faible profondeur. La caractérisation des propriétés hétérogènes du réservoir passe par une description préalable des propriétés des réservoirs carbonatés et notamment celle de la porosité. Au regard de la large gamme d'échelles d'hétérogénéités allant du micromètre à la centaine de mètre, l'utilisation de méthodes multi-échelles est nécessaire. Ce travail de thèse a porté sur la caractérisation multi-échelle de la porosité, allant de la microporosité à la porosité karstique, et l'établissement de relations entre les changements d'échelles d'observation de ces propriétés. A très petite échelle, la microtomographie haute-résolution à rayons X extrait à l'échelle du micron la microporosité et la connectivité des réseaux microporeux. A plus basse fréquence (celle du mm), la tomographie est également utilisée afin d'étudier le réseau micro- et macroporeux à l'échelle des carottes, à laquelle sont réalisées nombre d'observations et de mesures, notamment pétrophysiques. A l'échelle du puits, les images de paroi permettent d'identifier les zones de dissolution allant des macropores centimétriques aux cavités karstiques plurimétriques. A plus grande échelle, la géophysique de surface permet d'étudier l'extension latérale de ces karsts et hétérogénéités. Enfin, la structure poreuse du site étudié sur cette vaste gamme d'échelles est intégrée à la géologie régionale de la plate-forme de Llucmajor. Ainsi, l'observatoire installé sur le site pour surveiller la propagation de l'eau salée dans l'aquifère a pu être replacé dans son contexte géologique et pétrophysique, afin d'affiner l'interprétation de ses mesures de variations de salinité au cours du temps.

Cette étude concerne la fatigue et le vieillissement des élastomères en environnements marin et thermique. Elle vise plus particulièrement à proposer des outils applicables à l'échelle industrielle permettant d'intégrer les effets du vieillissement sur l'évolution des propriétés mécaniques, en particulier sur les propriétés en fatigue. Une vaste campagne expérimentale a ainsi été lancée pour deux matériaux industriels et pour deux types d'environnements (marin et thermique). Le vieillissement est caractérisé à toutes les échelles, depuis l'étude des mécanismes de dégradation (causes) jusqu'aux évolutions du comportement mécanique (conséquences). Les méthodes classiques d'extrapolation de la durée de vie sont ensuite appliquées et confrontées à des données issues de structures vieillies en service. Le second aspect de cette étude concerne la proposition d'une méthode d'identification rapide des propriétés à la fatigue. Pour cela, une méthode spécifique, dédiée aux matériaux élastomères, est proposée. Elle repose sur l'utilisation de mesures d'auto- échauffement par thermographie infrarouge, sur la description des populations de défauts sous chargement cyclique par micro-tomographie aux rayons X et sur l'utilisation d'un modèle de terme source pour alimenter un critére d'endurance énergétique. Enfin, dans la dernière partie, nous présentons les outils manquants à l'extension de ces résultats aux structures industrielles permettant de prendre en compte les gradients de propriétés induits par le vieillissement.

La qualité des images en TEP est directement liée à l'algorithme de reconstruction utilisé. Les méthodes de reconstruction statistiques peuvent modéliser, via la matrice système, les différents phénomènes physiques qui détériorent la qualité des images. Notre approche consiste d'une part, à modéliser les différents phénomènes géométriques et physiques, au moyen de l'outil de simulation Monte-Carlo, et d'autre part, à exploiter les données dans un format "mode-liste" présentant de nombreux avantages, afin de les reconstruire à l'aide d'un algorithme de reconstruction 3D. Cette méthode, appelée SIMALIR, a été appliquée et étudiée sur la caméra MOSAIC(TM) de Philips, système TEP dédié au petit animal qui a fait l'objet d'une modélisation à l'aide de la plateforme GATE. Les principales performances du modèle simulé et particulièrement la distribution spatiale ont été testées et validées, afin de pouvoir les exploiter dans la modélisation de la matrice système. Les résultats précliniques et sur des objets tests ont montré que la méthode SIMALIR permet d'améliorer considérablement la qualité des images reconstruites par rapport aux programmes de reconstruction classiques.

Nous avons développé une simulation de Monte Carlo rapide pour la tomographie à émission de positons (TEP) fondée sur le code SimSET modélisant la propagation des photons gammas dans le patient et le scanner basé sur un design en blocs. La validation de ce simulateur avec un code bien validé, GATE, et des données acquises sur un scanner TEP GE Discovery STE a montré qu'il modélise précisément les spectres en énergie (erreur inférieure à 4,6%), la résolution spatiale (6,1%), la fraction de diffusé (3,5%), la sensibilité aux coïncidences primaires (2,3%) et les taux de comptage (12,7%). Nous avons ensuite développé une correction de la diffusion incorporant l'énergie des photons détectés en mode liste en plus de la distribution spatiale des coïncidences diffusées. Notre approche est basée sur une reformulation de la fonction de vraisemblance contenant l'information en énergie donnant lieu à un algorithme de reconstruction EM contenant des termes de correction de la diffusion spatiaux et énergétiques. Nous avons aussi proposé une nouvelle méthode de normalisation du sinogramme diffusé utilisant l'information en énergie. Enfin, nous avons développé une méthode d'estimation des spectres primaires et diffusés détectés dans différents secteurs du scanner TEP. Nous avons évalué notre méthode spatio-spectrale de correction de la diffusion ainsi que la méthode spatiale traditionnelle dans des simulations de Monte Carlo réalistes. Ces résultats montrent que notre approche réduit les biais de quantification de 60% dans les régions froides dans les patients obèses, donnant lieu à des erreurs de quantification inférieures à 13% même dans les patients les plus larges en mode 3D (comparé à une erreur de 65% avec la méthode conventionnelle).

Nous présentons d'abord un système de tomographie optique fonctionnant en lumière blanche, avec un interféromètre de Linnik. La résolution tomographique obtenue est de l'ordre de 1,2 m dans l'air. Les capacités de l'instrument sont illustrées par des images tomographiques d'intérêt biomédical. D'autre part, nous présentons un polarimètre de Mueller à cristaux liquides nématiques optimisé pour la propagation des erreurs de mesures brutes aux images finales. Les mesures sont rapides (de l'ordre de la seconde) et précises (erreurs de l'ordre de 1%). Les images polarimétriques de coupes histologiques montrent l'intérêt de ce dispositif pour l'étude microscopique d'échantillons contenant des macromolécules fibreuses (collagène). Enfin, des images en degré de polarisation (DOP) de cols de l'utérus humains ex vivo viennent confirmer les possibilités de l'imagerie polarimétrique pour la détection de dysplasies (DOP plus élevé que dans les zones saines).

Auf der Sonne sind viele Phänomene zu sehen die mit der solaren magnetischen Aktivität zusammenhängen. Das dafür zuständige Magnetfeld wird durch einen Dynamo erzeugt, der sich vermutlich am Boden der Konvektionszone in der sogenannten Tachocline befindet. Angetrieben wird der Dynamo teils von der differenziellen Rotation, teils von den magnetischen Turbulenzen in der Konvektionszone. Die differentielle Rotation kann an der Sonnenoberfläche durch beobachten der Sonnenfleckbewegungen gemessen werden.<br>Um einen größeren Parameterraum zum Testen von Dynamotheorien zu erhalten, kann man diese Messungen auch auf andere Sterne ausdehnen. Das primäre Problem dabei ist, dass die Oberflächen von Sternen nicht direkt beobachtet werden können. Indirekt kann man dies jedoch mit Hilfe der Doppler-imaging Methode erreichen, die die Doppler-Verbreitung der Spektrallinien von schnell rotierenden Sternen benützt. Um jedoch ein Bild der Sternoberfläche zu erhalten, bedarf es vieler hochaufgelöster spektroskopischer Beobachtungen, die gleichmäßig über eine Sternrotation verteilt sein müssen. Für Sterne mit langen Rotationsperioden sind diese Beobachtungen nur schwierig durchzuführen. Das neue robotische Observatorium STELLA adressiert dieses Problem und bietet eine auf Dopplerimaging abgestimmte Ablaufplanung der Beobachtungen an. Dies wird solche Beobachtungen nicht nur leichter durchführbar machen, sondern auch effektiver gestalten.<br>Als Vorschau welche Ergebnisse mit STELLA erwartet werden können dient eine Studie an sieben Sternen die allesamt eine lange (zwischen sieben und 25 Tagen) Rotationsperiode haben. Alle Sterne zeigen differentielle Rotation, allerdings sind die Messfehler aufgrund der nicht zufriedenstellenden Datenqualität von gleicher Größenordnung wie die Ergebnisse, ein Problem das bei STELLA nicht auftreten wird. Um die Konsistenz der Ergebnisse zu prüfen wurde wenn möglich sowohl eine Kreuzkorrelationsanalyse als auch die sheared-image Methode angewandt. Vier von diesen sieben Sternen weisen eine differentielle Rotation in umgekehrter Richtung auf als auf der Sonne zu sehen ist. Die restlichen drei Sterne weisen schwache, aber in der Richtung sonnenähnliche differentielle Rotation auf.<br>Abschließend werden diese neuen Messungen mit bereits publizierten Werten kombiniert, und die so erhaltenen Daten auf Korrelationen zwischen differentieller Rotation, Rotationsperiode, Evolutionsstaus, Spektraltyp und Vorhandensein eines Doppelsterns überprüft. Alle Sterne zusammen zeigen eine signifikante Korrelation zwischen dem Betrag der differenziellen Rotation und der Rotationsperiode. Unterscheidet man zwischen den Richtungen der differentiellen Rotation, so bleibt nur eine Korrelation der Sterne mit antisolarem Verhalten. Darüberhinaus zeigt sich auch, dass Doppelsterne schwächer differentiell rotieren. / The sun shows a wide variety of magnetic-activity related phenomena. The magnetic field responsible for this is generated by a dynamo process which is believed to operate in the tachocline, which is located at the bottom of the convection zone. This dynamo is driven in part by differential rotation and in part by magnetic turbulences in the convection zone. The surface differential rotation, one key ingredient of dynamo theory, can be measured by tracing sunspot positions.<br>To extend the parameter space for dynamo theories, one can extend these measurements to other stars than the sun. The primary obstacle in this endeavor is the lack of resolved surface images on other stars. This can be overcome by the Doppler imaging technique, which uses the rotation-induced Doppler-broadening of spectral lines to compute the surface distribution of a physical parameter like temperature. To obtain the surface image of a star, high-resolution spectroscopic observations, evenly distributed over one stellar rotation period are needed. This turns out to be quite complicated for long period stars. The upcoming robotic observatory STELLA addresses this problem with a dedicated scheduling routine, which is tailored for Doppler imaging targets. This will make observations for Doppler imaging not only easier, but also more efficient.<br>As a preview of what can be done with STELLA, we present results of a Doppler imaging study of seven stars, all of which show evidence for differential rotation, but unfortunately the errors are of the same order of magnitude as the measurements due to unsatisfactory data quality, something that will not happen on STELLA. Both, cross-correlation analysis and the sheared image technique where used to double check the results if possible. For four of these stars, weak anti-solar differential rotation was found in a sense that the pole rotates faster than the equator, for the other three stars weak differential rotation in the same direction as on the sun was found.<br>Finally, these new measurements along with other published measurements of differential rotation using Doppler imaging, were analyzed for correlations with stellar evolution, binarity, and rotation period. The total sample of stars show a significant correlation with rotation period, but if separated into antisolar and solar type behavior, only the subsample showing anti-solar differential rotation shows this correlation. Additionally, there is evidence for binary stars showing less differential rotation as single stars, as is suggested by theory. All other parameter combinations fail to deliver any results due to the still small sample of stars available.

La croissance récente d'Internet lors deux dernières décennies a conduit à un besoin croissant de techniques permettant de mesurer la structure et la performance d'Internet. Les techniques de mesures de réseaux peuvent être classifiées en méthodes passives qui utilisent des données collectées au niveau des routeurs, et les méthodes actives, reposant sur l'injection active et l'observation de paquets-sondes. Les méthodes actives, qui sont la motivation principale de ce doctorat, sont particulièrement adaptées aux utilisateurs finaux, qui ne peuvent pas accéder aux données mesurées par les routeurs avec l'architecture actuelle d'Internet. Sur un autre plan, la théorie des réseaux se développe depuis un siècle, et de nombreux outils permettent de prédire la performance d'un système, en fonction de quelques paramètres clés. La théorie des files d'attentes émerge comme une solution particulièrement fructueuse, que ce soit pour les réseaux téléphoniques ou pour les réseaux filaires à commutation de paquet. Dans ce dernier cas, elle s'intéresse au mécanisme à l'échelle des paquets, et prédit des statistiques à ce niveau. À l'échelle des flots de paquets, la théorie des réseaux à partage de bande passante permet une abstraction de tout schéma d'allocation de bande passante, y compris le partage implicite résultant du protocole TCP. De nombreux travaux ont montré comment les résultats provenant de ces théories peuvent s'appliquer aux réseaux réels, et en particulier à Internet, et dans quels aspects le comportement de réseaux réels diffère des prédictions théoriques. Cependant, il y a eu peu de travaux établissant des liens entre le point de vue théorique d'un réseau et le problème pratique consistant à le mesurer. Le but de ce manuscrit est de bâtir quelques ponts entre le monde des méthodes de mesure par sondes actives et le monde de la théorie des réseaux. Nous adoptons l'approche des problèmes inverses, qui peuvent être vus en opposition aux problèmes directs. Un problème direct prédit l'évolution d'un système défini, en fonction des conditions initiales et d'une équation d'évolution connue. Un problème inverse observe une partie de la trajectoire d'un système défini, et cherche à estimer les conditions initiales ou paramètres pouvant conduire à cette trajectoire. Les données des méthodes de mesure par sondes actives sont les séries temporelles des pertes et délais des sondes, c'est-à-dire précisément une partie de la "trajectoire" d'un réseau. Ainsi, les méthodes de mesures par sondes actives peuvent être considérées comme des problèmes inverses pour une théorie des réseaux qui permettrait une prédiction exacte de l'évolution des réseaux. Nous montrons dans ce document comment les méthodes de mesures par sondes actives sont reliées aux problèmes inverses dans la théories des files d'attentes. Nous spécifions comment les contraintes de mesures peuvent être incluses dans les problèmes inverses, quels sont les observables, et détaillons les étapes successives pour un problème inverse dans la théorie des files d'attentes. Nous classifions les problèmes en trois catégories différentes, en fonction de la nature de leur résultat et de leur généralité, et donnons des exemples simples pour illustrer leurs différentes propriétés. Nous étudions en détail un problème inverse spécifique, où le réseau se comporte comme un réseau dit "de Kelly" avecK serveurs en tandem. Dans ce cas précis, nous calculons explicitement la distribution des délais de bout en bout des sondes, en fonction des capacités résiduelles des serveurs et de l'intensité des sondes. Nous montrons que l'ensemble des capacités résiduelles peut être estimé à partir du délai moyen des sondes pour K intensités de sondes différentes. Nous proposons une méthodes d'inversion alternative, à partir de la distribution des délais des sondes pour une seule intensité de sonde. Dans le cas à deux serveurs, nous donnons une caractérisation directe de l'estimateur du maximum de vraisemblance des capacités résiduelles. Dans le cas général, nous utilisons l'algorithme Espérance-Maximisation (E-M). Nous prouvons que dans le cas à deux serveurs, la suite des estimations de E-M converge vers une limite finie, qui est une solution de l'équation de vraisemblance. Nous proposons une formule explicite pour le calcul de l'itération quand K = 2 ou K = 3, et prouvons que la formule reste calculable quelque soit le nombre de serveurs. Nous évaluons ces techniques numériquement. À partir de simulations utilisant des traces d'un réseau réel, nous étudions indépendamment l'impact de chacune des hypothèses d'un réseau de Kelly sur les performances de l'estimateur, et proposons des facteurs de correction simples si besoin. Nous étendons l'exemple précédant au cas des réseaux en forme d'arbre. Les sondes sont multicast, envoyées depuis la racine et à destination des feuilles. À chaque noeud, elles attendent un temps aléatoire distribué de façon exponentielle. Nous montrons que ce modèle est relié au modèle des réseaux de Kelly sur une topologie d'arbre, avec du trafic transverse unicast et des sondes multicast, et calculons une formule explicite pour la vraisemblance des délais joints. Nous utilisons l'algorithme E-M pour calculer l'estimateur de vraisemblance du délai moyen à chaque noeud, et calculons une formule explicite pour la combinaison des étapes E et M. Des simulations numériques illustrent la convergence de l'estimateur et ses propriétés. Face à la complexité de l'algorithme, nous proposons une technique d'accélération de convergence, permettant ainsi de considérer des arbres beaucoup plus grands. Cette technique contient des aspects innovant dont l'intérêt peut dépasser le cadre de ces travaux. Finalement, nous explorons le cas des problèmes inverses dans la théorie des réseaux à partage de bande passante. À partir de deux exemples simples, nous montrons comment un sondeur peut mesurer le réseau en faisant varier le nombre de flots de sondes, et en mesurant le débit associé aux flots dans chaque cas. En particulier, si l'allocation de bande passante maximise une fonction d'utilité -équitable, l'ensemble des capacités des réseaux et leur nombre de connections associé peut être identifié de manière unique dans la plupart des cas. Nous proposons un algorithme pour effectuer cette inversion, avec des exemples illustrant ses propriétés numériques.

Les inversions conjointes sont maintenant communément utilisées en Sciences de la Terre. Elles ont été développées afin d'améliorer notre connaissance et compréhension de la structure interne de la Terre puisqu'elles apportent de plus en plus de contraintes sur les paramètres inversés. Dans cette thèse, nous proposons un nouveau processus d'inversion conjointe qui prend en compte trois paramètres différents et qui mène à l'obtention d'un modèle lithosphérique tridimensionnel. Notre méthode utilise trois types de données gravimétriques et sismologiques qui présentent une bonne complémentarité : (1) Les fonctions récepteur P afin d'obtenir les variations de profondeur du Moho, (2) les délais de temps d'arrivées P des téléséismes pour retrouver les anomalies de vitesse dans la croûte et le manteau supérieur, (3) les anomalies gravimétriques qui donnent accès aux variations de densité à l'échelle lithosphérique. Dans le schéma d'inversion proposé, nous inversons d'abord les fonctions récepteur. Les variations de profondeur du Moho qui en résultent sont alors interpolées puis incorporées comme information a priori dans le processus d'inversion conjointe gravimétrie - tomographie télésismique. Que ce soit pour les fonctions récepteur ou pour la partie conjointe, notre méthode d'inversion se base sur un algorithme de recherche afin de trouver la structure qui minimise l'écart aux trois types de données. Nous avons opté pour l'algorithme de voisinage qui présente l'avantage de concentrer sa recherche dans les régions qui minimisent l'écart aux données sans toutefois abandonner les autres régions afin de limiter les risques de convergence vers un minimum local. Nous utilisons un modèle 3-D constitué de noeuds de vitesse afin de modéliser les délais de temps. Ces noeuds de vitesse sont liés aux noeuds de densité grâce à une relation linéaire entre vitesse et densité. Cette relation peut évoluer avec la profondeur de manière à simuler les effets dus aux variations de pression. Nous avons réalisé plusieurs tests synthétiques afin d'évaluer le comportement de notre nouvelle méthode. Ils ont permis de montrer que les variations du Moho ainsi que les anomalies de vitesse et densité sont bien retrouvées. Ils ont également montré l'importance de la prise en compte de ces trois types de données dans un même schéma d'inversion, notamment en ce qui concerne la faculté de discerner les variations d'interfaces des anomalies volumiques. Nous avons ensuite appliqué la méthode dans une région continentale en phase de rifting où les processus de déformation sont complexes et mal connus. Notre nouvelle méthode a apporté davantage de contraintes sur les paramètres inversés et a permis d'obtenir un modèle correctement résolu de la croûte et du manteau supérieur sous le dôme de Hangai (Mongolie Centrale). Une anomalie asthénosphérique sous le dôme est souvent invoquée pour expliquer sa topographie élevée mais c'est une anomalie (en vitesse et densité) crustale (croûte inférieure) qui s'exprime dans nos modèles. Ses effets sur la topographie et son origine (compositionnelle et/ou thermique) sont encore débattus et nécessiteraient des études plus poussées (paramétriques, VP =VS, ...).

L'un des aspects les plus critiques dans les planifications dosimétriques des traitements par radiothérapie est la délinéation des limites de la tumeur. Cette délinéation se fait généralement sur les images anatomiques de tomodensitométrie (TDM). Mais récemment, il est recommandé de faire cette délinéation pour les cancers broncho-pulmonaires non à petites cellules (CBNPC) sur les images fonctionnelles de Tomographie par Émission de Positon (TEP) pour prendre en compte les caractéristiques biologiques de la cible. Jusqu'à ce jour, aucune technique de segmentation ne s'est révélée satisfaisante pour les images TEP en application clinique. Une solution pour ce problème est proposée dans cette étude. Méthodes : Les optimisations de notre méthode ont consisté principalement à faire l'ajustement des seuils directement à partir des corps des patients au lieu de le faire à partir du fantôme. Résultats : Pour les lésions de grands axes supérieurs à 20 mm, notre technique de segmentation a montré une bonne estimation des mesures histologiques (la moyenne de différence de diamètre entre données mesurées et déterminées avec notre technique = +1,5 ± 8,4 %) et une estimation acceptable des mesures TDM. Pour les lésions de grands axes inférieurs ou égaux à 20 mm, cette méthode a montré un écart avec les mesures dérivées des données histologiques ou bien des données TDM. Conclusion : Cette nouvelle méthode d'ajustement montre une bonne précision pour la délimitation des lésions de grands axes compris entre 2 et 4,5 cm. Néanmoins, elle n'évalue pas correctement les lésions les plus petites, cela peut être dû à l'effet du volume partiel

Die Aufklärung des kooperativen Materialtransports während des Sinterns von Metallen zielt auf eine verbesserte Vorhersage und theoretische Beschreibung des Sinterprozesses ab. Dem kooperativen Materialtransport wird in der Literatur ein signifikanter Beitrag zur Gesamtschwindung zugeschrieben. Lose oder gepresste Pulverschüttungen weisen unter Wärmebehandlung eine Volumenabnahme auf. Neben der Zentrumsannäherung und dem Wachstum der Sinterkontakte rührt diese Schwindung von spannungsbedingten Relativbewegungen der Partikel in Form von Translationen und insbesondere von Rotationen her. Die Sinterprozesse wurden bereits seit den 1940er Jahren für Zwei-Teilchen-Modelle, sowie für 1- und 2-dimensionale Versuchsanordnungen ausführlichen Analysen unterzogen. Die gewonnenen Erkenntnisse lassen sich aber nur bedingt auf 3-dimensionale Proben übertragen. Berechnungen der Gesamtschwindung basieren auf dem Zwei-Teilchen-Modell auf und berücksichtigen den kooperativen Materialtransport nur in unzureichendem Maße, so dass eine Diskrepanz zwischen berechneter und an realen Sinterkörpern beobachteter Schwindung auftritt. In dieser Arbeit wurde der kooperative Materialtransport im Sinteranfangsstadium (bis 1050°C) am Modellsystem Kupfer mittels Synchrotron-Computertomographie an 3-dimensionalen Proben untersucht. Die Tomographie bietet als erste Methode die Möglichkeit, Partikelbewegungen in 3-dimensionalen Proben über mehrere Sinterschritte hinweg zerstörungsfrei zu untersuchen und quantitative Daten zu erhalten. Auf diese Weise konnten die Bewegungen eines jeden Partikels in-situ über den Sinterverlauf verfolgt werden. Zur Bestimmung der Orientierungsänderung der sphärischen Partikel war es notwendig, die Partikel mittels Focused Ion Beam mit Markierungen zu versehen. Die Detektion und Auswertung der Partikelbewegung erfolgte über eine eigens entwickelte Analysesoftware basierend auf photogrammetrischen Methoden. Diese Software erlaubt es, Aussagen über Rotationsgeschwindigkeit und Rotationswinkel der Partikel sowie über die Dichte, Zentrumsannäherung, Änderung der Koordinationszahl und Gesamtgröße zu treffen. Neben der bereits in der Literatur dargestellten Rotation der Partikel in Bezug auf dessen Kontaktpartner konnte zusätzlich eine Eigenrotation um den eigenen Schwerpunkt festgestellt werden. Es wurde ein neues Modell entwickelt, das neben Diffusionsmechanismen auch Korngrenzengleiten zulässt und damit Eigenrotationen ermöglicht. / The investigation of the cooperative material transport of metals during sintering aims to improve the theoretical descriptions of the sintering process. In literature the cooperative material transport is ascribed to have a significant contribution to the total shrinkage. Loose or pressed powder packed beds show a volume shrinkage during heat treatment. This shrinkage is caused by centre approach and the growth of sintering necks but also results of particle movements like translations and particularly rotations. The sintering process was investigated in detail since the 1940’s on the basis of the two particle model as well as on 1- and 2- dimensional test arrangements. The obtained data can be transferred to 3-dimensional specimens only to a limited extent. Calculations of the total shrinkage are based on the two particle model and do hardly consider the cooperative material transport. That is why there is a difference between calculated and measured shrinkage in real sinter compacts. In this study the cooperative material transport was investigated in the initial sintering stage (up to 1050°C) in the model system copper. These measurements were conducted on 3-dimensional samples using synchrotron computed tomography. Tomography offers for the first time the opportunity to investigate particle movements in 3-dimensional samples non-destructively during several sintering steps to obtain quantitative data. Using this method it was possible to track the particle movements during the sintering process in-situ. To determine orientation rearrangements of spherical particles it was necessary to mark each particle by focused ion beam with a borehole. The detection and analysis of particle movements was carried out by custom image analysis software using photogrammetric methods. This software is able to determine the rotation speed and angle of particles, the local density, the centre approach, the change of the coordination number and the total size of a specimen. Beside the in literature described particle rotation in relation to contact partners an additional intrinsic rotation around the centre of mass of a particle was discovered. A new model was developed which permits not only diffusion mechanisms but also grain boundary sliding and therefore allows intrinsic rotations.