Academic literature on the topic 'Tomographic PIV'

This article deals with a pioneering application of tomographic particle image velocimetry (tomographic PIV) for the hydrodynamic and hydroacoustic analysis of a marine propeller. The hydrodynamic study was mainly focused on the topological analysis of the propeller wake characteristics in the near field based on the vorticity field and on the tilting and stretching terms of the vorticity transport equation. Hydroacoustic analysis concerned the use of tomographic PIV in combination with the Powell's acoustic analogy. Tomographic PIV proved to be a valid tool for the detailed quantitative reconstruction of the complex vortex topology in the propeller wake and provided an accurate description of the source terms of the Powell's analogy. In particular, it was shown that the tip vortex perturbation represents the dominant nonlinear contribution to the radiated far-field noise in non-cavitating flow conditions.

Further understanding the small scale coherent structures which occur in high Reynolds number turbulence would be of enormous benefit. Therefore, the aim of the current project was to make well resolved three-dimensional flow measurements of the mixing flow between counter rotating impellers, using Tomographic Particle Image Velocimetry (TPIV).TPIV software was developed, with a novel approach permitting a significant reduction in processing time, and a series of numerical accuracy studies contributing to the fundamental understanding of this new technique. Basic flow characterisation determined the local isotropy, homogeneity and expected Reynolds number scaling. A favourable comparison between planar PIV and TPIV increased confidence in the latter, which was used to assess the dynamics and topology of the dissipation scale structures. In support of previous investigations similar topology, strain rate alignment, scale-invariance, and clustering behaviours are demonstrated. Correlated high enstrophy and dissipation regions occur in the periphery of larger structures, resulting in intermittency. Geometry characterisation indicates a predominance of tube-like structures, which are observed to form from larger ribbon-like structures through unsteady breakdown and vortex roll-up. Significant correlation between intermittent fields of dissipation and enstrophy describe the fine scales effects. These relationships should pave the way for more accurate models, capable of relating small scales and large scales during the prediction of dynamically important quantities.

PONTIFÍCIA UNIVERSIDADE CATÓLICA DO RIO DE JANEIRO
COORDENAÇÃO DE APERFEIÇOAMENTO DO PESSOAL DE ENSINO SUPERIOR
CONSELHO NACIONAL DE DESENVOLVIMENTO CIENTÍFICO E TECNOLÓGICO
PROGRAMA DE EXCELENCIA ACADEMICA
PROGRAMA DE DOUTORADO SANDUÍCHE NO EXTERIOR
O estudo fundamental de estruturas turbulentas próximas a paredes é de significativa importância devido ao seu papel dominante em inúmeras aplicações de engenharia. Em escoamentos turbulentos, estruturas de diferentes escalas interagem entre si através de um fenômeno complexo caracterizado por um processo contínuo e auto-sustentável. Embora a presença de padrões de movimento coerente pareça ser responsável pela manutenção da turbulência, a evolução e a interação destas estruturas ainda não são completamente entendidas. Para investigar tais fenômenos, foram realizados estudos de escoamentos turbulentos próximos a paredes sobre uma placa plana em um túnel de vento e dentro de um canal de água quadrado. O presente trabalho empregou uma técnica tomográfica de velocimetria por imagem de partícula de alta taxa de medição de campos de velocidade tridimensionais. O processamento da técnica tomo-PIV foi otimizado a fim de melhorar o desempenho do algoritmo e a exatidão dos campos vetoriais medidos. As medições realizadas permitiram a obtenção de estatísticas condicionais das estruturas coerentes típicas da camada limite turbulenta, tais como regiões de baixa e de alta velocidades, varreduras (sweeps), ejeções (ejections) e vórtices. Os resultados quantitativos mostraram vórtices distribuídos por todo o escoamento, entretanto, concentrados na vizinhança das regiões de baixa velocidade e ejeções. Os resultados indicaram que os vórtices estavam indiretamente correlacionados com o processo de produção de turbulência, apesar de serem os principais responsáveis pela dissipação da energia turbulenta.
The fundamental study of turbulent structures near a wall is of significant importance due to its dominant role in innumerous engineering applications. In turbulent flows, eddies of different scales interact with each other in a complex phenomenon of a continuous self-sustaining process. Although the presence of patterns of coherent motion seem to be responsible for the maintenance of turbulence, the evolution and the interaction of these structures are still not completely understood. To investigate such phenomena, near-wall turbulent flows over a smooth at plate in a wind tunnel and inside a square water channel were conducted. The present work employed a high-repetition tomographic PIV technique to measure three dimensional, time-resolved velocity fields. The tomo-PIV processing was optimized in order to improve the algorithm performance and the vector fields accuracy. Conditional statistics of coherent structures, namely, lowand high-speed regions, sweeps, ejections and vortices, were computed. Quantitative results showed vortical structures spread all over the flow, nevertheless, concentrated close to low-speed regions and ejections. Results indicated that vortical structures were indirectly correlated with the process of turbulence production, although they were the main responsible for the turbulent energy dissipation.
L étude fondamentale des structures turbulentes proche d une paroi est d une significative importance en raison de son rôle dominant dans nombreuses applications d ingênierie. Dans les écoulements turbulents, des structures d échelles différentes interagissent dans un phénomène complexe caractérisé par un processus autonome et continu. Bien que la présence de modéles de mouvements cohérents semble entretenir la turbulence, l évolution et l interaction de ces structures ne sont pas encore complétement compris. Pour analyser ces phénomènes, l étude d écoulements turbulents à proximité de paróis a été effectuées sur une plaque plane dans une soufflerie et à l intérieur d un canal hydraulique de section carrée. Une technique tomographique de vélocimétrie par images de particules pour mesurer de champs de vitesse tridimensionnel a été employée. Le traitement des données a été optimisé afin d améliorer la précision des champs de vecteurs mesurés. Les mesures realisés ont permis d obtenir des statistiques conditionnelles des structures cohérentes, à savoir, les régions à faible et à grande vitesse, les sweeps, les ejections et les tourbillons. Les résultats quantitatives ont montré des structures tourbillonnaires distribués sur tout le ux, cependant, concentrée au voisinage des régions à faible vitesse et des ejections. Les rèsultats ont indiqué que les structures tourbillonnaires étaient indirectement corrélés avec le processus de production de la turbulence, bien qu ils ont été les principaux responsables de la dissipation d énergie turbulente.

Cette thèse porte sur le développement de la PIV tomographique (tomo-PIV) pour la mesure d’écoulements turbulents (Elsinga et al 2006). Cette méthode se fonde sur la reconstruction tomographique d’une distribution volumique d’intensité de particules traceuses, à partir de projections enregistrées par des caméras. Les distributions volumiques sont corrélées, fournissant ainsi un champ de déplacement 3D.Les principales avancées de la recherche sur cette technique sont présentées ainsi que les points bloquants. Les efforts ont principalement été portés sur la reconstruction tomographique. La principale difficulté est le bruit tomographique (particules fantômes) qui croît exponentiellement lorsqu’une forte densité de traceur est requise pour obtenir une résolution spatiale fine de la mesure, particulièrement pour les écoulements turbulents.Afin de mieux appréhender ce bruit de reconstruction, nous avons étudié numériquement les facteurs expérimentaux nuisant à la qualité de la reconstruction. Des considérations géométriques ont permis de quantifier l’impact de «particules ajoutées», qui se trouvent dans le volume de l’union mais pas dans le volume de l’intersection entre la zone laser et les champs de vue des caméras. La diminution du ratio signal-à-bruit dans les images, due à la diffusion de Mie et l’astigmatisme des optiques, a pour principal effet la perte de vraies particules dans la reconstruction.Étudier les conditions optiques de la tomo-PIV nous a permis de proposer une approche alternative à la reconstruction tomographique classique, qui vise à reconstruire une particule presque sur un unique voxel, plutôt que comme un agrégat de voxels de taille étendue, en se fondant sur une représentation particulaire des images. Nous nommons cette méthode Reconstruction Volumique de Particules (PVR). Après avoir été incorporée à un algorithme de reconstruction (SMART), il est possible d’élargir la représentation particulaire de PVR, afin d’obtenir des blobs de 2/3 voxels de diamètre requis par les algorithmes de corrélation de 3D-PIV. Des simulations numériques sur un large spectre de conditions génératrices, ont montré qu’utiliser PVR-SMART permettait des gains de performance par rapport à un algorithme classique comme tomo-SMART (Atkinson 2009).L’aspect vélocimétrie par corrélation de la méthode a aussi été pris en compte avec une extension sur GPU à la 3D (FOLKI-3D) de l’algorithme FOLKI-PIV (Champagnat et al. 2011). Le déplacement y est cherché en minimisant itérativement une fonctionnelle, du type des moindres carrés, par déformation de volume. Les tests synthétiques confirment que la réponse fréquentielle d’espace est semblable à celle d’autres algorithmes classiques itératifs de déformation de volume. Les simulations numériques de reconstruction tomographique ont permis de caractériser la robustesse de l’algorithme au bruit spécifique de la tomographie. Nous avons montré que FOLKI-3D était plus robuste aux particules fantômes cohérentes que les algorithmes classiques de déformation volumique. De plus, des gains de performance ont été observés en utilisant des schémas d’ordre élevé pour différents types de bruit.L’application de PVR-SMART sur des données expérimentales a été effectuée sur un jet d’air turbulent. Différentes densités de particules ont été utilisées pour comparer les performance de PVR-SMART avec tomo-SMART sur la région proche buse du jet. Avec le pré-traitement d’image utilisé, nous avons montré que les champs de vitesse de PVR-SMART étaient près de 50 % moins bruités que ceux de tomo-SMART. L’analyse sur les champs de vitesse comporte l’étude de quantités statistiques, de peak-locking, de divergence, du tenseur des gradients ainsi que de structures cohérentes.Enfin, nous concluons avec une synthèse des résultats obtenus au cours de cette étude, en envisageant de nouvelles perspectives de recherche dans le contexte de la PIV tomographique
This research dissertation focuses on the developments of tomographic PIV (tomo-PIV) for the measurement of turbulent flows (Elsinga et al. 2006). It is based on the tomographic reconstruction of a volumic intensity distribution of tracer particles from projections recorded on cameras. The corresponding volumic distributions are correlated to obtain 3D displacement fields.The present work surveys the state of advancement of the research conducted on this technique and the main issues it has been confronted with so far. The main research focus was on tomographic reconstruction. Indeed, its main limitation is the appearance of ghost particles, ie reconstruction noise, which occurs when high tracer concentrations are required for high spatial resolution measurements.For a thorough understanding of tomographic noise, we carried out a numerical study of experimental factors impacting the quality of tomographic reconstruction. Geometric considerations quantified the impact of "added particles" lying in the Union volume but not in the Intersection volume, between the camera fields of view and the illumination area. This phenomenon was shown to create ghost particles. The decrease in signal-to-noise ratio in the image was investigated, considering Mie scattering and defocusing effects. Particle image defocusing mainly results in the loss of real particles in reconstruction. Mie scattering’s main impact is also the loss of real particles due to the polydisperse nature of the seeding.This study of imaging conditions for tomo-PIV led us to propose an alternative approach to classical tomographic reconstruction. It seeks to recover nearly single voxel particles rather than blobs of extended size using a particle-based representation of image data. We term this approach Particle Volume Reconstruction (PVR). PVR underlies a more physical, sparse volumic representation of point particles, which lives halfway between infinitely small particles, and voxel blobs commonly used in tomo-PIV. From that representation, it is possible to smooth it to 2 voxel diameter blobs for a 3D-PIV use of PVR incorporated in a SMART algorithm. Numerical simulations showed that PVR-SMART outperforms tomo-SMART (Atkinson et al. 2009) on a variety generating conditions and a variety of metrics on volume reconstruction and displacement estimation, especially in the case of seeding density greater than 0.06 ppp.We introduce a cross-correlation technique for 3D-PIV (FOLKI-3D) as an extension to 3D of the FOLKI-PIV algorithm (Champagnat et al. 2011). The displacement is searched as the minimizer of a sum of squared differences, solved iteratively by using volume deformation. Numerical tests confirmed that spatial frequency response is similar to that of standard iterative deformation algorithms. Numerical simulations of tomographic reconstruction characterized the robustness of the algorithm to specific tomographic noise. FOLKI-3D was found more robust to coherent ghosts than standard deformation algorithms, while gains in accuracy of the high-order deformation scheme were obtained for various signal noises.The application of PVR-SMART on experimental data was performed on a turbulent air jet. Several seeding density conditions were used to compare the performance of tomo-SMART and PVR-SMART on the near field region of the jet. With the given image pre-processing, PVR-SMART was found to yield velocity fields that are about 50 % less noisy than tomo-SMART. The velocity field comparison included velocity field statistical properties, peak-locking study, flow divergence analysis, velocity gradient tensor and coherent structures exploration.Finally, conclusions are drawn from the main results of this dissertation and lead to potential research perspectives of our work with respect to the future of tomographic PIV

Les systèmes de climatisation et de ventilation sont souvent composés de configurations type jets impactant, sur leur partie terminale. Ainsi, les flux d’air soufflés viennent impacter des obstacles munis de fentes, de différentes formes, afin d’améliorer le mélange. Les conditions de confinement et de soufflage provoquent parfois un inconfort au niveau acoustique. Les nuisances acoustiques générées sont dues à un phénomène de boucles de rétroaction se traduisant par l’apparition des sons auto-entretenus. La production du son par un écoulement fluide en champ libre ou en interaction avec une structure a fait l’objet de nombreuses études. Dans le cas d’un champ acoustique et pour un écoulement à faible nombre de Mach la résolution du corollaire énergétique de Howe permet d’évaluer la puissance acoustique générée ou absorbée par les interactions entre le champ acoustique et l’écoulement. Le calcul de cette puissance nécessite la connaissance de trois paramètres : la vorticité, la vitesse et la vitesse acoustique par des méthodes analytiques ou en utilisant des données expérimentales. Expérimentalement, la mesure du champ cinématique, pour en déduire la vorticité, nécessite une technique de mesure tridimensionnelle. Pour cela une plate-forme expérimentale, utilisant de la vélocimétrie Laser, a été développée, et équipée pour générer les écoulements d’un jet plan. Les champs cinématiques de ces écoulements ont été mesurés en utilisant la technique PIV, avec un dispositif de PIV stéréoscopique. Les champs cinématiques de trente plans parallèles ont été mesurés afin d’étudier les champs de vitesses correspondants. Deux méthodes de reconstruction ont été appliquées à ces plans : la POD et la moyenne de phase. Le volume est obtenu par une interpolation des plans reconstruits donnant accès aux trois composantes de la vitesse. Pour valider ces méthodes de reconstruction en 3D à faible coût, elles étaient comparées à des mesures expérimentales réalisées par le même dispositif expérimental, dans les mêmes conditions, par la PIV tomographique donnant accès aux champs cinématiques tridimensionnels
Air conditioning and ventilation systems are often composed of jets having a configuration of an impinging jet, on their end part. Thus, the blown airflows impact slotted obstacles of different shapes to improve mixing. The conditions of confinement and blowing sometimes cause acoustic incompatibility. The acoustic noises generated are due to a phenomenon of feedback loops resulting of the appearance of self-sustained sounds. The production of sound by a free flow or in interaction with a structure has been the subject of many studies. In the case of an acoustic field and for a flow of low Mach number, Howe's energetic correlation is used to evaluate the acoustic power generated or absorbed by the interactions between the acoustic field and the flow. The calculation of this power requires the knowledge of three parameters : vorticity, velocity and acoustic velocity by analytical methods or by using experimental data. Experimentally, the measurement of the kinematic field, to deduce the vorticity, requires a three-dimensional measurement technique. For this purpose, an experimental platform, using laser velocimetry, has been developed and equipped to generate flows of a plane jet. The kinematic fields of these flows were measured using the PIV technique, with a stereoscopic PIV device. The kinematic fields of thirty parallel planes were measured to study the corresponding velocity fields. Two reconstruction methods have been applied to these plans : the POD and the phase average. The volume is obtained by an interpolation of the reconstructed planes giving access to the three components of the velocity. To validate these low-cost 3D reconstruction methods, they were compared to experimental measurements made by the same experimental setup, under the same conditions, by using the tomographic PIV giving access to the three-dimensional kinematic fields

Les mécanismes d’évolution spatio-temporelle des structures turbulentes instationnaires tridimensionnelles, et en particulier ceux rencontrés aux plus grandes échelles, sont à l’origine de phénomènes d’instabilité qui conduisent très souvent à une diminution de la performance des systèmes énergétiques. C’est le cas des variations cycle-à-cycle dans le moteur à combustion interne. Malgré les progrès substantiels réalisés par la simulation numérique en mécanique des fluides, les approches expérimentales demeurent essentielles pour l’analyse et la compréhension des phénomènes physiques ayant lieu. Dans ce travail de thèse, deux types de vélocimétrie par image de particules (PIV) ont été appliqués et adaptés au banc moteur optique du laboratoire Coria pour étudier l’écoulement en fonction de six conditions de fonctionnement du moteur. La PIV Haute Cadence 2D2C a permis d’abord d’obtenir un suivi temporel de l’écoulement dans le cylindre durant un même cycle moteur ainsi qu’identifier ces variations cycliques. La PIV Tomographique 3D3C a permis ensuite d’étendre les données mesurées vers l’espace tridimensionnel. La Tomo-PIV fait intervenir 4 caméras en position angulaire visualisant un environnement de géométrie complexe, confinée, ayant un accès optique restreint et introduisant des déformations optiques importantes. Cela a nécessité une attention particulière vis-à-vis du processus de calibration 3D des modèles de caméras. Des analyses conditionnées 2D et 3D de l’écoulement sont effectuées en se basant principalement sur la décomposition propre orthogonale (POD) permettant de séparer les différentes échelles de structure et le critère Γ permettant l’identification des centres des tourbillons
The unsteady evolution of three-dimensional large scale flow structures can often lead to a decrease in the performance of energetic systems. This is the case of cycle-to-cycle variations occurring in the internal combustion engine. Despite the substantial advancement made by numerical simulations in fluid mechanics, experimental measurements remain a requirement to validate any numerical model of a physical process. In this thesis, two types of particle image velocimetry (PIV) were applied and adapted to the optical engine test bench of the Coria laboratory in order to study the in-cylinder flow with respect to six operating conditions. First, the Time-Resolved PIV (2D2C) allowed obtaining a temporal tracking of the in-cylinder flow and identifying cyclic variabilities. Then tomographic PIV (3D3C) allowed extending the measured data to the three-dimensional domain. The Tomo-PIV setup consisted of 4 cameras in angular positioning, visualizing a confined environment with restricted optical access and important optical deformations. This required a particular attention regarding the 3D calibration process of camera models. 2D and 3D conditional analyses of the flow were performed using the proper orthogonal decomposition (POD) allowing to separate the different scales of flow structures and the Γ criterion allowing the identification of vortices centres

Die 3D Particle Tracking Velocimetry (3D PTV) ist eine Methode zur bildbasierten Bestimmung von Geschwindigkeitsfeldern in Gas- oder Flüssigkeitsströmungen. Dazu wird die Strömung mit Partikeln markiert und durch ein Mehrkamerasystem beobachtet. Das Ergebnis der Datenauswertung sind 3D Trajektorien einer großen Anzahl von Partikeln, die zur statistischen Analyse der Strömung genutzt werden können. In der vorliegenden Arbeit werden verschiedene neu entwickelte Modelle gezeigt, die das Einsatzspektrum vergrößern und die Leistungsfähigkeit der 3D PTV erhöhen. Wesentliche Neuerungen sind der Einsatz eines Spiegelsystems zur Generierung eines virtuellen Kamerasystems, die Modellierung von komplex parametrisierten Trennflächen der Mehrmedienphotogrammetrie, eine wahrscheinlichkeitsbasierte Trackingmethode sowie eine neuartige Methode zur tomographischen Rekonstruktion von Rastervolumendaten. Die neuen Modelle sind an drei realen Experimentieranlagen und mit synthetischen Daten getestet worden. Durch den Einsatz eines Strahlteilers vor dem Objektiv einer einzelnen Kamera und vier Umlenkspiegeln, positioniert im weiteren Strahlengang, werden vier virtuelle Kameras generiert. Diese Methode zeichnet sich vor allem durch die Wirtschaftlichkeit als auch durch die nicht notwendige Synchronisation aus. Vor allem für die Anwendung im Hochgeschwindigkeitsbereich sind diese beiden Faktoren entscheidend. Bei der Beobachtung von Phänomenen in Wasser kommt es an den Trennflächen verschiedener Medien zur optischen Brechung. Diese muss für die weitere Auswertung zwingend modelliert werden. Für komplexe Trennflächen sind einfache Ansätze über zusätzliche Korrekturterme nicht praktikabel. Der entwickelte Ansatz basiert auf der mehrfachen Brechung jedes einzelnen Bildstrahls. Dazu müssen die Trennflächenparameter und die Kameraorientierungen im selben Koordinatensystem bekannt sein. Zumeist wird die Mehrbildzuordnung von Partikeln durch die Verwendung von Kernlinien realisiert. Auf Grund von instabilen Kameraorientierungen oder bei einer sehr hohen Partikeldichte sind diese geometrischen Eigenschaften nicht mehr ausreichend, um die Mehrbildzuordnung zu lösen. Unter der Ausnutzung weiterer geometrischer, radiometrischer und physikalischer Eigenschaften kann die Bestimmung der 3D Trajektorien dennoch durchgeführt werden. Dabei werden durch die Analyse verschiedener Merkmale diejenigen ausgewählt, welche sich für die spatio-temporale Zuordnung eignen. Die 3D PTV beruht auf der Diskretisierung der Partikelabbildungen im Bildraum und der anschließenden Objektkoordinatenbestimmung. Eine rasterbasierte Betrachtungsweise stellt die tomographische Rekonstruktion des Volumens dar. Hierbei wird die Intensitätsverteilung wird im Volumen rekonstruiert. Die Bewegungsinformationen werden im Anschluss aus den Veränderungen aufeinander folgender 3D-Bilder bestimmt. Durch dieses Verfahren können Strömungen mit einer höheren Partikeldichte im Volumen analysiert werden. Das entwickelte Verfahren basiert auf der schichtweisen Entzerrung und Zusammensetzung der Kamerabilder. Die entwickelten Modelle und Ansätze sind an verschiedenen Versuchsanlagen erprobt worden. Diese unterschieden sich stark in der Größe (0,5 dm³ – 20 dm³ – 130 m³) und den vorherrschenden Strömungsgeschwindigkeiten (0,3 m/s – 7 m/s – 0,5 m/s)
3D Particle Tracking Velocimetry (3D PTV) is an image based method for flow field determination. It is based on seeding a flow with tracer particles and recording the flow with a multi camera system. The results are 3D trajectories of a large number of particles for a statistical analysis of the flow. The thesis shows different novel models to increase the spectrum of applications and to optimize efficiency of 3D PTV. Central aspects are the use of the mirror system to generate a virtual multi camera system, the modelling of complex interfaces of multimedia photogrammetry, a probability based tracking method and a novel method for tomographic reconstruction of volume raster data. The improved models are tested in three real testing facilities and with synthetic data. Using a beam splitter in front of the camera lens and deflecting mirrors arranged in the optical path, a four headed virtual camera system can be generated. This method is characterised by its economic efficiency and by the fact that a synchronisation is not necessary. These facts are important especially when using high speed cameras. When observing phenomena in water, there will be refraction at the different interfaces. This has to be taken into account and modelled for each application. Approaches which use correction terms are not suitable to handle complex optical interfaces. The developed approach is based on a multiple refraction ray tracing with known interface parameters and camera orientations. Mostly the multi image matching of particles is performed using epipolar geometry. Caused by the not stable camera orientation or a very high particle density this geometric properties are not sufficient to solve the ambiguities. Using further geometrical radiometrical and physical properties of particles, the determination of the 3D trajectories can be performed. After the analysis of different properties those of them are chosen which are suitable for spatio-temporal matching. 3D PTV bases on the discretisation of particle images in image space and the following object coordinate determination. A raster based approach is the tomographic reconstruction of the volume. Here the light intensity distribution in the volume will be reconstructed. Afterwards the flow information is determined from the differences in successive 3D images. Using tomographic reconstruction techniques a higher particle density can be analysed. The developed approach bases on a slice by slice rectification of the camera images and on a following assembly of the volume. The developed models and approaches are tested at different testing facilities. These differ in size (0.5 dm³ – 20 dm³ – 130 m³) and flow velocities (0.3 m/s – 7 m/s – 0.5 m/s)

Die 3D Particle Tracking Velocimetry (3D PTV) ist eine Methode zur bildbasierten Bestimmung von Geschwindigkeitsfeldern in Gas- oder Flüssigkeitsströmungen. Dazu wird die Strömung mit Partikeln markiert und durch ein Mehrkamerasystem beobachtet. Das Ergebnis der Datenauswertung sind 3D Trajektorien einer großen Anzahl von Partikeln, die zur statistischen Analyse der Strömung genutzt werden können. In der vorliegenden Arbeit werden verschiedene neu entwickelte Modelle gezeigt, die das Einsatzspektrum vergrößern und die Leistungsfähigkeit der 3D PTV erhöhen. Wesentliche Neuerungen sind der Einsatz eines Spiegelsystems zur Generierung eines virtuellen Kamerasystems, die Modellierung von komplex parametrisierten Trennflächen der Mehrmedienphotogrammetrie, eine wahrscheinlichkeitsbasierte Trackingmethode sowie eine neuartige Methode zur tomographischen Rekonstruktion von Rastervolumendaten. Die neuen Modelle sind an drei realen Experimentieranlagen und mit synthetischen Daten getestet worden. Durch den Einsatz eines Strahlteilers vor dem Objektiv einer einzelnen Kamera und vier Umlenkspiegeln, positioniert im weiteren Strahlengang, werden vier virtuelle Kameras generiert. Diese Methode zeichnet sich vor allem durch die Wirtschaftlichkeit als auch durch die nicht notwendige Synchronisation aus. Vor allem für die Anwendung im Hochgeschwindigkeitsbereich sind diese beiden Faktoren entscheidend. Bei der Beobachtung von Phänomenen in Wasser kommt es an den Trennflächen verschiedener Medien zur optischen Brechung. Diese muss für die weitere Auswertung zwingend modelliert werden. Für komplexe Trennflächen sind einfache Ansätze über zusätzliche Korrekturterme nicht praktikabel. Der entwickelte Ansatz basiert auf der mehrfachen Brechung jedes einzelnen Bildstrahls. Dazu müssen die Trennflächenparameter und die Kameraorientierungen im selben Koordinatensystem bekannt sein. Zumeist wird die Mehrbildzuordnung von Partikeln durch die Verwendung von Kernlinien realisiert. Auf Grund von instabilen Kameraorientierungen oder bei einer sehr hohen Partikeldichte sind diese geometrischen Eigenschaften nicht mehr ausreichend, um die Mehrbildzuordnung zu lösen. Unter der Ausnutzung weiterer geometrischer, radiometrischer und physikalischer Eigenschaften kann die Bestimmung der 3D Trajektorien dennoch durchgeführt werden. Dabei werden durch die Analyse verschiedener Merkmale diejenigen ausgewählt, welche sich für die spatio-temporale Zuordnung eignen. Die 3D PTV beruht auf der Diskretisierung der Partikelabbildungen im Bildraum und der anschließenden Objektkoordinatenbestimmung. Eine rasterbasierte Betrachtungsweise stellt die tomographische Rekonstruktion des Volumens dar. Hierbei wird die Intensitätsverteilung wird im Volumen rekonstruiert. Die Bewegungsinformationen werden im Anschluss aus den Veränderungen aufeinander folgender 3D-Bilder bestimmt. Durch dieses Verfahren können Strömungen mit einer höheren Partikeldichte im Volumen analysiert werden. Das entwickelte Verfahren basiert auf der schichtweisen Entzerrung und Zusammensetzung der Kamerabilder. Die entwickelten Modelle und Ansätze sind an verschiedenen Versuchsanlagen erprobt worden. Diese unterschieden sich stark in der Größe (0,5 dm³ – 20 dm³ – 130 m³) und den vorherrschenden Strömungsgeschwindigkeiten (0,3 m/s – 7 m/s – 0,5 m/s).
3D Particle Tracking Velocimetry (3D PTV) is an image based method for flow field determination. It is based on seeding a flow with tracer particles and recording the flow with a multi camera system. The results are 3D trajectories of a large number of particles for a statistical analysis of the flow. The thesis shows different novel models to increase the spectrum of applications and to optimize efficiency of 3D PTV. Central aspects are the use of the mirror system to generate a virtual multi camera system, the modelling of complex interfaces of multimedia photogrammetry, a probability based tracking method and a novel method for tomographic reconstruction of volume raster data. The improved models are tested in three real testing facilities and with synthetic data. Using a beam splitter in front of the camera lens and deflecting mirrors arranged in the optical path, a four headed virtual camera system can be generated. This method is characterised by its economic efficiency and by the fact that a synchronisation is not necessary. These facts are important especially when using high speed cameras. When observing phenomena in water, there will be refraction at the different interfaces. This has to be taken into account and modelled for each application. Approaches which use correction terms are not suitable to handle complex optical interfaces. The developed approach is based on a multiple refraction ray tracing with known interface parameters and camera orientations. Mostly the multi image matching of particles is performed using epipolar geometry. Caused by the not stable camera orientation or a very high particle density this geometric properties are not sufficient to solve the ambiguities. Using further geometrical radiometrical and physical properties of particles, the determination of the 3D trajectories can be performed. After the analysis of different properties those of them are chosen which are suitable for spatio-temporal matching. 3D PTV bases on the discretisation of particle images in image space and the following object coordinate determination. A raster based approach is the tomographic reconstruction of the volume. Here the light intensity distribution in the volume will be reconstructed. Afterwards the flow information is determined from the differences in successive 3D images. Using tomographic reconstruction techniques a higher particle density can be analysed. The developed approach bases on a slice by slice rectification of the camera images and on a following assembly of the volume. The developed models and approaches are tested at different testing facilities. These differ in size (0.5 dm³ – 20 dm³ – 130 m³) and flow velocities (0.3 m/s – 7 m/s – 0.5 m/s).

Le travail mené dans le cadre de cette thèse porte sur l’étude expérimentale et numérique de l’écoulement au sein d’un éjecteur supersonique. Le régime d’écoulement qui s’installe dans ces appareils est très complexe du fait des phénomènes physiques qui les caractérisent comme la turbulence et les ondes de choc. Les méthodes expérimentales utilisées sont la mesure de la pression le long de l’axe de l’éjecteur `a l’aide d’une sonde développée à cet effet, la visualisation de l’écoulement par tomographie laser et la mesure de vitesse par PIV. Les simulations numériques sont réalisées à l’aide du code Ansys-Fluent en 2D axisymétrique et en 3D. Dans un premier temps, une étude de sensibilité du modèle numérique portant sur les paramètres de simulations et les modèles de turbulence est menée sur l’éjecteur fonctionnant sans flux induit. La validation des simulations repose sur une comparaison des résultats numériques avec des mesures de vitesse par PIV. Un modèle 3D s’est avéré incontournable pour l’étude de l’écoulement dans l’éjecteur avec flux induit à cause de sa géométrie complexe. Les outils expérimentaux et numériques développés permettent d’analyser finement l’interaction des flux moteur et induit, en particulier les processus de recompression par chocs obliques et de mélange. Une tentative de modélisation par LES des instabilités de l’écoulement détectées expérimentalement est également abordée
The work reported in this thesis relates to the experimental and numerical studies of the flow within a supersonic ejector. The flow pattern which occurs in these apparatuses is very complex because of the flow phenomena encountered like flow turbulence and shock waves. The experimental methods used are the measurement of the pressure along the axis of the ejector using a specific probe developed for this purpose, the flow visualization by laser tomography and the velocity measurement by PIV. The numerical simulations are carried out using the Ansys-Fluent code with 2D axisymmetric and 3D models. First, a study of sensitivity to the numerical parameters of simulation and to the turbulence models is carried out on the ejector operating without induced flow. The validation of the simulations is achieved by a comparison between the numerical results and velocity measurements by PIV. A 3D model is necessary for the simulation of the flow in the ejector operating with induced flow because of the complex ejector geometry. The experimental techniques and the numericalmodels developed make it possible to analyze the interaction of the primary and secondary flows, in particular the process of recompression by oblique shocks and the mixing process. An attempt at modeling by LES simulation the flow instabilities detected during experiments is also approached

Bones are multifunctional passive organs of movement that supports soft tissue and directly attached muscles. They also protect internal organs and are a reserve of calcium, phosphorus and magnesium. Each bone is covered with periosteum, and the adjacent bone surfaces are covered by articular cartilage. Histologically, the bone is an organ composed of many different tissues. The main component is bone tissue (cortical and spongy) composed of a set of bone cells and intercellular substance (mineral and organic), it also contains fat, hematopoietic (bone marrow) and cartilaginous tissue. Bones are a tissue that even in adult life retains the ability to change shape and structure depending on changes in their mechanical and hormonal environment, as well as self-renewal and repair capabilities. This process is called bone turnover. The basic processes of bone turnover are: • bone modeling (incessantly changes in bone shape during individual growth) following resorption and tissue formation at various locations (e.g. bone marrow formation) to increase mass and skeletal morphology. This process occurs in the bones of growing individuals and stops after reaching puberty • bone remodeling (processes involve in maintaining bone tissue by resorbing and replacing old bone tissue with new tissue in the same place, e.g. repairing micro fractures). It is a process involving the removal and internal remodeling of existing bone and is responsible for maintaining tissue mass and architecture of mature bones. Bone turnover is regulated by two types of transformation: • osteoclastogenesis, i.e. formation of cells responsible for bone resorption • osteoblastogenesis, i.e. formation of cells responsible for bone formation (bone matrix synthesis and mineralization) Bone maturity can be defined as the completion of basic structural development and mineralization leading to maximum mass and optimal mechanical strength. The highest rate of increase in pig bone mass is observed in the first twelve weeks after birth. This period of growth is considered crucial for optimizing the growth of the skeleton of pigs, because the degree of bone mineralization in later life stages (adulthood) depends largely on the amount of bone minerals accumulated in the early stages of their growth. The development of the technique allows to determine the condition of the skeletal system (or individual bones) in living animals by methods used in human medicine, or after their slaughter. For in vivo determination of bone properties, Abstract 10 double energy X-ray absorptiometry or computed tomography scanning techniques are used. Both methods allow the quantification of mineral content and bone mineral density. The most important property from a practical point of view is the bone’s bending strength, which is directly determined by the maximum bending force. The most important factors affecting bone strength are: • age (growth period), • gender and the associated hormonal balance, • genotype and modification of genes responsible for bone growth • chemical composition of the body (protein and fat content, and the proportion between these components), • physical activity and related bone load, • nutritional factors: – protein intake influencing synthesis of organic matrix of bone, – content of minerals in the feed (CA, P, Zn, Ca/P, Mg, Mn, Na, Cl, K, Cu ratio) influencing synthesis of the inorganic matrix of bone, – mineral/protein ratio in the diet (Ca/protein, P/protein, Zn/protein) – feed energy concentration, – energy source (content of saturated fatty acids - SFA, content of polyun saturated fatty acids - PUFA, in particular ALA, EPA, DPA, DHA), – feed additives, in particular: enzymes (e.g. phytase releasing of minerals bounded in phytin complexes), probiotics and prebiotics (e.g. inulin improving the function of the digestive tract by increasing absorption of nutrients), – vitamin content that regulate metabolism and biochemical changes occurring in bone tissue (e.g. vitamin D3, B6, C and K). This study was based on the results of research experiments from available literature, and studies on growing pigs carried out at the Kielanowski Institute of Animal Physiology and Nutrition, Polish Academy of Sciences. The tests were performed in total on 300 pigs of Duroc, Pietrain, Puławska breeds, line 990 and hybrids (Great White × Duroc, Great White × Landrace), PIC pigs, slaughtered at different body weight during the growth period from 15 to 130 kg. Bones for biomechanical tests were collected after slaughter from each pig. Their length, mass and volume were determined. Based on these measurements, the specific weight (density, g/cm3) was calculated. Then each bone was cut in the middle of the shaft and the outer and inner diameters were measured both horizontally and vertically. Based on these measurements, the following indicators were calculated: • cortical thickness, • cortical surface, • cortical index. Abstract 11 Bone strength was tested by a three-point bending test. The obtained data enabled the determination of: • bending force (the magnitude of the maximum force at which disintegration and disruption of bone structure occurs), • strength (the amount of maximum force needed to break/crack of bone), • stiffness (quotient of the force acting on the bone and the amount of displacement occurring under the influence of this force). Investigation of changes in physical and biomechanical features of bones during growth was performed on pigs of the synthetic 990 line growing from 15 to 130 kg body weight. The animals were slaughtered successively at a body weight of 15, 30, 40, 50, 70, 90, 110 and 130 kg. After slaughter, the following bones were separated from the right half-carcass: humerus, 3rd and 4th metatarsal bone, femur, tibia and fibula as well as 3rd and 4th metatarsal bone. The features of bones were determined using methods described in the methodology. Describing bone growth with the Gompertz equation, it was found that the earliest slowdown of bone growth curve was observed for metacarpal and metatarsal bones. This means that these bones matured the most quickly. The established data also indicate that the rib is the slowest maturing bone. The femur, humerus, tibia and fibula were between the values of these features for the metatarsal, metacarpal and rib bones. The rate of increase in bone mass and length differed significantly between the examined bones, but in all cases it was lower (coefficient b <1) than the growth rate of the whole body of the animal. The fastest growth rate was estimated for the rib mass (coefficient b = 0.93). Among the long bones, the humerus (coefficient b = 0.81) was characterized by the fastest rate of weight gain, however femur the smallest (coefficient b = 0.71). The lowest rate of bone mass increase was observed in the foot bones, with the metacarpal bones having a slightly higher value of coefficient b than the metatarsal bones (0.67 vs 0.62). The third bone had a lower growth rate than the fourth bone, regardless of whether they were metatarsal or metacarpal. The value of the bending force increased as the animals grew. Regardless of the growth point tested, the highest values were observed for the humerus, tibia and femur, smaller for the metatarsal and metacarpal bone, and the lowest for the fibula and rib. The rate of change in the value of this indicator increased at a similar rate as the body weight changes of the animals in the case of the fibula and the fourth metacarpal bone (b value = 0.98), and more slowly in the case of the metatarsal bone, the third metacarpal bone, and the tibia bone (values of the b ratio 0.81–0.85), and the slowest femur, humerus and rib (value of b = 0.60–0.66). Bone stiffness increased as animals grew. Regardless of the growth point tested, the highest values were observed for the humerus, tibia and femur, smaller for the metatarsal and metacarpal bone, and the lowest for the fibula and rib. Abstract 12 The rate of change in the value of this indicator changed at a faster rate than the increase in weight of pigs in the case of metacarpal and metatarsal bones (coefficient b = 1.01–1.22), slightly slower in the case of fibula (coefficient b = 0.92), definitely slower in the case of the tibia (b = 0.73), ribs (b = 0.66), femur (b = 0.59) and humerus (b = 0.50). Bone strength increased as animals grew. Regardless of the growth point tested, bone strength was as follows femur > tibia > humerus > 4 metacarpal> 3 metacarpal> 3 metatarsal > 4 metatarsal > rib> fibula. The rate of increase in strength of all examined bones was greater than the rate of weight gain of pigs (value of the coefficient b = 2.04–3.26). As the animals grew, the bone density increased. However, the growth rate of this indicator for the majority of bones was slower than the rate of weight gain (the value of the coefficient b ranged from 0.37 – humerus to 0.84 – fibula). The exception was the rib, whose density increased at a similar pace increasing the body weight of animals (value of the coefficient b = 0.97). The study on the influence of the breed and the feeding intensity on bone characteristics (physical and biomechanical) was performed on pigs of the breeds Duroc, Pietrain, and synthetic 990 during a growth period of 15 to 70 kg body weight. Animals were fed ad libitum or dosed system. After slaughter at a body weight of 70 kg, three bones were taken from the right half-carcass: femur, three metatarsal, and three metacarpal and subjected to the determinations described in the methodology. The weight of bones of animals fed aa libitum was significantly lower than in pigs fed restrictively All bones of Duroc breed were significantly heavier and longer than Pietrain and 990 pig bones. The average values of bending force for the examined bones took the following order: III metatarsal bone (63.5 kg) <III metacarpal bone (77.9 kg) <femur (271.5 kg). The feeding system and breed of pigs had no significant effect on the value of this indicator. The average values of the bones strength took the following order: III metatarsal bone (92.6 kg) <III metacarpal (107.2 kg) <femur (353.1 kg). Feeding intensity and breed of animals had no significant effect on the value of this feature of the bones tested. The average bone density took the following order: femur (1.23 g/cm3) <III metatarsal bone (1.26 g/cm3) <III metacarpal bone (1.34 g / cm3). The density of bones of animals fed aa libitum was higher (P<0.01) than in animals fed with a dosing system. The density of examined bones within the breeds took the following order: Pietrain race> line 990> Duroc race. The differences between the “extreme” breeds were: 7.2% (III metatarsal bone), 8.3% (III metacarpal bone), 8.4% (femur). Abstract 13 The average bone stiffness took the following order: III metatarsal bone (35.1 kg/mm) <III metacarpus (41.5 kg/mm) <femur (60.5 kg/mm). This indicator did not differ between the groups of pigs fed at different intensity, except for the metacarpal bone, which was more stiffer in pigs fed aa libitum (P<0.05). The femur of animals fed ad libitum showed a tendency (P<0.09) to be more stiffer and a force of 4.5 kg required for its displacement by 1 mm. Breed differences in stiffness were found for the femur (P <0.05) and III metacarpal bone (P <0.05). For femur, the highest value of this indicator was found in Pietrain pigs (64.5 kg/mm), lower in pigs of 990 line (61.6 kg/mm) and the lowest in Duroc pigs (55.3 kg/mm). In turn, the 3rd metacarpal bone of Duroc and Pietrain pigs had similar stiffness (39.0 and 40.0 kg/mm respectively) and was smaller than that of line 990 pigs (45.4 kg/mm). The thickness of the cortical bone layer took the following order: III metatarsal bone (2.25 mm) <III metacarpal bone (2.41 mm) <femur (5.12 mm). The feeding system did not affect this indicator. Breed differences (P <0.05) for this trait were found only for the femur bone: Duroc (5.42 mm)> line 990 (5.13 mm)> Pietrain (4.81 mm). The cross sectional area of the examined bones was arranged in the following order: III metatarsal bone (84 mm2) <III metacarpal bone (90 mm2) <femur (286 mm2). The feeding system had no effect on the value of this bone trait, with the exception of the femur, which in animals fed the dosing system was 4.7% higher (P<0.05) than in pigs fed ad libitum. Breed differences (P<0.01) in the coross sectional area were found only in femur and III metatarsal bone. The value of this indicator was the highest in Duroc pigs, lower in 990 animals and the lowest in Pietrain pigs. The cortical index of individual bones was in the following order: III metatarsal bone (31.86) <III metacarpal bone (33.86) <femur (44.75). However, its value did not significantly depend on the intensity of feeding or the breed of pigs.

The isothermal flow fields of injectors have undergone several computational and experimental investigations using point and planar measurement techniques,. Within the swirl induced vortex breakdown region it is only LES that has been able to predict fully the presence of a three dimensional helical vortex structure for a particular injector, and in certain conditions (no central fuel jet), a precessing vortex core. These structures can be elucidated from point and planar measurements and favorable comparisons of velocity statistics between experiment and LES predictions strengthen these findings. However, volumetric, 3-component measurement of velocity data has not been widely available to provide conclusive evidence of the exact three dimensional nature of the vortex structures that exist. An experimental setup utilizing time resolved tomographic PIV on a water flow rig is described in this paper. This is used to provide as high-quality aerodynamic study as possible of a single stream radially-fed air swirl gaseous fuel injector. The level of accuracy of the tomographic PIV technique is demonstrated by calculating the divergence of the velocity field as well as validating the results against a comprehensive 2 and 3 component standard PIV velocity database and other measurement techniques and predictions. Structure identification methods have been employed to visualise and understand the complex flow topology within the near field of the injector. The change in topology with and without the stabilising central jet is further investigated and agrees with findings of planar PIV results. While the velocity error associated with the tomo-PIV results is higher than the planar results the data agree well within the identified uncertainty bounds and are complimentary in understanding the flow field structure. Thus a full volumetric aerodynamic survey is available for the first time on this isothermal flow case.

The onset and progression of vortices generated by the hull of a surface combatant has been investigated for drift angles of 0, 10 and 20 degrees by means of tomographic particle image velocimetry. Vortices are seen to form on the sonar dome and fore-body of the ship model, as well as on the bilge keels, and on the after-body. Two major streamwise vortex systems were identified: the vortices generated by the sonar dome and the vortices generated on the fore-body of the hull. Leading edge and crossflow vortices form on the bilge keels, and on the windward side, these vortices merge. Vortices are also seen to form on the afterbody of the hull at increased drift angles. Frequency analysis of the sonar dome vortex at drift angle 20 degrees reveals a dominant frequency of approximately 3 Hz.