No description available.

abstract: The application of novel visualization and modeling methods to the study of cardiovascular disease is vital to the development of innovative diagnostic techniques, including those that may aid in the early detection and prevention of cardiovascular disorders. This dissertation focuses on the application of particle image velocimetry (PIV) to the study of intracardiac hemodynamics. This is accomplished primarily though the use of ultrasound based PIV, which allows for in vivo visualization of intracardiac flow without the requirement for optical access, as is required with traditional camera-based PIV methods. The fundamentals of ultrasound PIV are introduced, including experimental methods for its implementation as well as a discussion on estimating and mitigating measurement error. Ultrasound PIV is then compared to optical PIV; this is a highly developed technique with proven accuracy; through rigorous examination it has become the “gold standard” of two-dimensional flow visualization. Results show good agreement between the two methods. Using a mechanical left heart model, a multi-plane ultrasound PIV technique is introduced and applied to quantify a complex, three-dimensional flow that is analogous to the left intraventricular flow. Changes in ventricular flow dynamics due to the rotational orientation of mechanical heart valves are studied; the results demonstrate the importance of multi-plane imaging techniques when trying to assess the strongly three-dimensional intraventricular flow. The potential use of ultrasound PIV as an early diagnosis technique is demonstrated through the development of a novel elasticity estimation technique. A finite element analysis routine is couple with an ensemble Kalman filter to allow for the estimation of material elasticity using forcing and displacement data derived from PIV. Results demonstrate that it is possible to estimate elasticity using forcing data derived from a PIV vector field, provided vector density is sufficient. / Dissertation/Thesis / Doctoral Dissertation Mechanical Engineering 2015

abstract: Cerebral aneurysms are pathological balloonings of blood vessels in the brain, commonly found in the arterial network at the base of the brain. Cerebral aneurysm rupture can lead to a dangerous medical condition, subarachnoid hemorrhage, that is associated with high rates of morbidity and mortality. Effective evaluation and management of cerebral aneurysms is therefore essential to public health. The goal of treating an aneurysm is to isolate the aneurysm from its surrounding circulation, thereby preventing further growth and rupture. Endovascular treatment for cerebral aneurysms has gained popularity over traditional surgical techniques due to its minimally invasive nature and shorter associated recovery time. The hemodynamic modifications that the treatment effects can promote thrombus formation within the aneurysm leading to eventual isolation. However, different treatment devices can effect very different hemodynamic outcomes in aneurysms with different geometries. Currently, cerebral aneurysm risk evaluation and treatment planning in clinical practice is largely based on geometric features of the aneurysm including the dome size, dome-to-neck ratio, and parent vessel geometry. Hemodynamics, on the other hand, although known to be deeply involved in cerebral aneurysm initiation and progression, are considered to a lesser degree. Previous work in the field of biofluid mechanics has demonstrated that geometry is a driving factor behind aneurysmal hemodynamics. The goal of this research is to develop a more combined geometric/hemodynamic basis for informing clinical decisions. Geometric main effects were analyzed to quantify contributions made by geometric factors that describe cerebral aneurysms (i.e., dome size, dome-to-neck ratio, and inflow angle) to clinically relevant hemodynamic responses (i.e., wall shear stress, root mean square velocity magnitude and cross-neck flow). Computational templates of idealized bifurcation and sidewall aneurysms were created to satisfy a two-level full factorial design, and examined using computational fluid dynamics. A subset of the computational bifurcation templates was also translated into physical models for experimental validation using particle image velocimetry. The effects of geometry on treatment were analyzed by virtually treating the aneurysm templates with endovascular devices. The statistical relationships between geometry, treatment, and flow that emerged have the potential to play a valuable role in clinical practice. / Dissertation/Thesis / Doctoral Dissertation Bioengineering 2016

Este trabalho representa o início de uma série de pesquisas que visam o estudo da turbulência e de sua relação com processos de mistura e trocas gasosas entre ar e água, através de experimentos com grades oscilantes. Seu objetivo foi o projeto e a construção de um tanque de grade oscilante, equipamento que gera turbulência com intensidade controlável. Após a construção do tanque, experimentos para medições de velocidade turbulenta foram feitos, com uso de uma técnica de velocimetria a laser. Uma grade de 9x9 barras foi usada na agitação do fluido, com uma amplitude de oscilação de 3 cm, para 4 diferentes freqüências de oscilação. Adquiriu-se 9600 imagens do movimento do fluido, em 6 regiões do tanque, para a obtenção dos campos de velocidade turbulenta, calculados através do software Visiflow e de um programa computacional desenvolvido neste trabalho. Os gráficos criados a partir dos campos possibilitaram a observação do decaimento espacial da turbulência e da região de sua produção. Os campos médios de velocidade mostraram-se bem inferiores aos campos instantâneos, indicando a existência de baixo escoamento médio do fluido. As condições de isotropia e a homogeneidade espacial da turbulência são mais aproximadas à medida que se afasta da grade. A intensidade turbulenta produzida é diretamente relacionada com a freqüência de oscilação. O número de imagens para a obtenção de uma média representativa da velocidade turbulenta neste equipamento parece ser dependente da freqüência de oscilação da grade. Isto aponta para a necessidade de estabelecer corretamente as condições para os cálculos estatísticos em escoamentos turbulentos / This work presents the project and construction of a tank with an oscillating grid, equipment which provides for the experimental studies of turbulence and its relation to mixing processes and gas transfer across fluid interface. Experiments were carried out with the use of digital particle image velocimetry technique, to investigate the properties of the produced turbulence. A grid made of 9x9 square bars was used to stir the water, with a stroke of 3 cm and 4 different oscillation frequencies. A number of 9600 images were acquired, in 6 regions of the tank, for the generation of the turbulent velocity fields through the software Visiflow and a computational program developed in this work. The results showed that the current equipment, with a new concept for the grid oscillation system, can be conveniently useful for studies in this research field. Average velocity fields appeared to be considerably smaller than the instantaneous velocity fields, which leads to the existence of nearly-stationary turbulence in the water volume. Nearly-isotropic turbulence and spatial homogeneity were approximate as the measurements distanced from the grid. The turbulent intensity was directly dependent on the oscillation frequency. The spatial decay of the turbulence and the region near the grid where it is produced could also be observed. The number of images necessary for the calculus of reliable root-mean-square turbulent velocities seems to be dependent on the oscillation frequency of the grid. It results in the necessity of establishing correct statistical analysis of turbulent flows

Red blood cells (RBCs) form approximately 40 to 45% of the human blood volume, and their behaviour and characteristics are the main determinant of blood properties, such as viscosity. RBCs are deformable species and stack together under low shear rate to form aggregates or rouleaux. Flowing RBCs migrate away from the wall leaving a cell-depleted layer known as the cell-free layer (CFL). This layer contributes to the blood viscosity and exchange between the RBCs and the target cells: a thinner CFL enhances the exchange process by reducing the diffusion distance. The formation of this CFL, however, is not yet completely understood. The goal of this study is to improve the understanding of the formation of the CFL in the micro-flow. This was accomplished by studying the effects of changing both the flow rate and the microchannel geometry on blood flow in microchannels. In this work, 10% hematocrit human blood suspensions were prepared in native plasma and flowed through poly-dimethylsiloxane (PDMS) microchannels of 100 μm x 34 μm cross-section. Investigation of the flowing cells was performed by using micro particle image velocimetry (μPIV) coupled with a high-speed camera. First, the high-speed camera images were processed with customized Matlab programs to detect and measure the CFL thickness and the RBC aggregates sizes. Second, the blood flow velocity profiles were measured using μPIV in order to determine the actual flow rate, the RBCs’ centerline velocity, and the shear rate. The results showed that the increase in both flow rate and shear rate significantly reduced the CFL thickness and RBC aggregates size. Comparison of the upstream and downstream measurements in the bifurcating microchannel showed that the change in microchannel geometry did not significantly influence CFL thickness and RBC aggregate size, while within the daughter branches, RBCs tended to flow close to the inner wall resulting in an undetectable CFL at the inner wall and in a larger CFL at the outer wall of the branch. These in vitro results quantitatively relate CFL thickness and RBC aggregate size at different shear rates. The findings are of immediate interest regarding the understanding of microcirculation and improved designs of microchips.

Background: Nowadays, one of the medical problem is rolling cell adhesion. Rolling cell adhesion is a complex process that requires the analysis of the challenging environment such as body fluid and is the process responsible for recruiting the cell to specific organs. In order to explore the rolling cell adhesion, mathematical model is proposed. Different image processing methods are created, such as optical flow - Lucas Kanade algorithm, and other type of methods related to mechanical fluid, namely PIV (Particle Image Velocimetry). Aim: The aim of this master thesis is the identification of challenges while using PIV in live cell images and propose the algorithm, which may analyze the rolling cell adhesion problem. Methods: In order to understand properly the rolling cell adhesion problem from biological site, literature review combined with the expert consultation is performed. According to gather information, mathematical model is proposed. Particle Image Velocimetry is explained according to literature review, where at the beginning the expert recommends some books as a primary research. As a result of this research, PIV challenges are identified and generally PIV idea is explained. Then two experiments are performed. The first experiment evaluates detection algorithms and the second one, analyses track algorithm vs. PIV. In order to evaluate the mentioned algorithms, some evaluation method are selected and some criteria are defined. Unfortunately the found methods are not perfect, therefore a new method related to performance evaluation using time series is proposed. Thesis result: The result of this thesis is a proposition of the algorithm, which can be used in the rolling cell adhesion. The algorithm is formed according to the detailed exploration of the rolling cell adhesion and analysis of the selected algorithms related to the image analysis during the theoretical research and experiments.

D'une part, les équations de Navier-Stokes permettent de décrire les écoulements fluides, la littérature est riche de méthodes numériques permettant la résolution de celle-ci. D'autre part, nous sommes capables de mesurer de manière non-intrusive différentes caractéristique d'un écoulement (champ de vitesse et pression, etc.). Dans le cadre de cette thèse, nous nous intéressons aux techniques d'assimilation de données qui combinent les modèles numériques avec les observations afin de déterminer une meilleure approximation du système. Cette thèse s'articule autour de l'assimilation de donnée variationnelle (4DVar) qui est plus précise par construction. Nous avons mené une première application sur la reconstruction de la hauteur et vitesse de la surface libre d'un fluide contenu dans un récipient rectangulaire à fond plat. L'écoulement est modélisé par les équations de shallow water et résolues numériquement. Les observations de l'évolution de la hauteur de la surface libre ont été prélevées par un capteur de profondeur (Kinect). Nous avons comparé les résultats de la reconstruction par 4DVar avec plusieurs version de la méthode d'assimilation hybride 4DEnVar. Enfin, nous avons appliqué la technique 4DVar à la reconstruction volumique de l'aval d'un sillage de cylindre à Reynolds 300. L'écoulement turbulent a été simulé par un code DNS parallèle Incompact3D. La reconstruction a été effectué en combinant tout d'abord des observations synthétiques en trois dimension, puis en combinant des observations de plans orthogonales en stéréo PIV. / In the one hand, flow dynamics are usually described by the NavierStokes equations and the literature provides a wide range of techniques to solve such equations. On the other hand, we can nowadays measure different characteristics of a flow (velocity, pressure, temperature etc...) with non-intrusive Particle Image Velocimetry techniques. Within this thesis, we take interest in the data assimilation techniques, that combine a dynamics model with measurements to determine a better approximation of the system. This thesis focus on the classic variational assimilation technique (4DVar) which ensures a high accuracy of the solution by construction. We carry out a first application of the 4DVar technique to reconstruct the characteristics (height and velocity field) of a uni directional wave at its free surface. The fluid evolution is simulated by the shallow water equations and solved numerically. We use a simple experimental setup envolving a depth sensor (Kinect sensor) to extract the free surface height. We compared the results of the 4DVar reconstruction with different versions of the hybrid data assimilation technique 4DEnVar. Finally, we apply the 4DVar technique to reconstruct the downstream of a three dimensional cylinder wake at Reynolds 300. The turbulent flow is simulated by the high-performance multi-threading DNS code Incompact3d. This dynamics model is first combined with synthetic three dimensional observations, then with real orthogonal-plane stereo PIV observations to reconstruct the full three dimensional flow.

[pt] Neste trabalho foi desenvolvido um sistema para medição de campos instantâneos de velocidade em regiões extensas do escoamento de fluidos. A técnica utilizada foi a velocimetria por imagem de partículas. Nesta técnica, as medidas de velocidade são obtidas a partir do registro da imagem de partículas traçadoras previamente distribuídas no fluido e iluminadas externamente por um plano de luz pulsada. As imagens capturadas são digitalizadas e processadas com algoritmos numéricos especialmente desenvolvidos para este fim. Estes algoritmos empregam técnicas de correlação cruzada de imagens ou de autocorrelação de imagens. Testes do programa desenvolvido foram realizados com imagens de escoamentos geradas em computador. Um sistema experimental foi desenvolvido para capturar imagens de diversos escoamentos teste de baixa velocidade. Os resultados obtidos foram satisfatórios. / [en] The present work describes the development and testo f a whole-field velocimetry system for measuring transient velocity fields in extensive flow regions. The technique employed is known as Particle Image Velocimetry. In this technique, double-exposure images of tracer particles distributed in the fluid and externally illuminated are registered in film or electronic camera. The images are digitized and processed by specially developed computer algorithms based on cross-correlation and autocorrelation techniques. The computer program developed was tested against computer generated flow images. A experimental setup was constructed to capture and analyze several real flow images. The tests conducted with several low-velocity flows yielded satisfactory results.

In this study, experimental data was collected in order to create a test case that can be used to validate computational fluid dynamics (CFD) simulations and the individual models used therein for gas turbine combustor applications. In many cases, the CFD results of gas turbine combustors do not correlate well with experimental results. For this reason, there is a requirement to test the simulation method used before CFD can successfully be used for combustor design. This test case encompasses all the features of a gas turbine combustor such as a swirler, primary, secondary and dilution holes as well as cooling rings. Experiments were performed on the same combustor geometry for both non-reacting and reacting flows. The non-reacting flow experiments consisted of stereoscopic particle image velocimetry (PIV) measurements performed at various planes in the three zones of the combustor. Data was collected on planes, both in line with the holes and in between the holes of each zone. For the reacting experiments, the temperatures on the outlet plane were measured using a thermocouple rake, thus a temperature contour plot on the outlet plane was produced. Further, the combustor can was modified with passive inserts, which were tested to determine their influence on the outlet temperature distribution during reacting runs. In this set-up, the outlet velocity profiles were also measured using a Pitot tube during both non-reacting and reacting flows. In addition to the outlet temperature distribution and velocity profiles, images of the flame patterns were captured, which showed the positions of flame tongues, fluctuating flames and steady flames. Carbon burn patterns on the walls of the combustor liner were also captured. From the data collected during the reacting runs, the pattern factor, profile factor, overall pressure loss and pressure loss factor were calculated. The non-reacting experiments performed using the PIV, produced three-dimensional velocity vector fields throughout the combustor. These experiments were performed at various flow rates, which gave an indication of which features of the combustor flow were affected by the flow rate. When comparing the individual PIV images alongside one another, the temporal nature of the combustor flow was also evident. The reacting experiments revealed a hot region of exhaust gas around the outer edge of the exhaust while there was a cooler region in the centre of the outlet flow. The PIV flowfield results revealed the reason for then hot outer ring-like region was due to the path the hot gasses would take. The hot combustor gas from the primary zone diverges outwards in the secondary zone then is further forced to the outside by the dilution recirculation zone. The hot flow then leaves the combustor along the wall while the cooler air from the jets leaves the combustor in the centre. The experiments performed produced a large variety of data that can be used to validate a number of aspects of combustor simulation using CFD. The non-reacting experimental data can be used to validate the turbulence models used and to evaluate how well the flow features were modelled or captured during the non-reacting stage of the combustor simulation process. The typical flow features such as jet penetration depths and the position and size of the recirculation regions are provided for effective comparison. The thermal results presented on the outlet plane of the combustor can be used for comparison with CFD results once combustion is modelled. Copyright / Dissertation (MEng)--University of Pretoria, 2010. / Mechanical and Aeronautical Engineering / unrestricted

Hydroelectricity plays an important role to balance the stability of grid network.  In order to improve the stability of presently high loaded grids, hydropower plants are being operated over a wide range of operations and experiencing frequent start-stop, load rejection, and load acceptance. The turbines need to sustain sudden change in their operating condition to balance the grid frequency. Francis turbines have been widely used because of their wider operating range and higher stability in operation during rapid load variation. This has resulted in severe damage to the turbines as they are not normally designed to operate under such transient conditions. Several low and high frequency pressure fluctuations prevail during transients operating conditions. Generally, wall pressure measurements are performed which may not provide sufficient information to investigate the flow instabilities related to these fluctuations. Thus, the main objective of the present work is to simplify and perform optical measurements in a turbine during transient operating conditions to investigate the flow field. The measurements have been performed at the Water Power Laboratory using a high head model Francis turbine. The turbine is a 1:5.1 scale down model of a prototype operating at the Tokke Power Plant, Norway. The model runner diameter, net head, and discharge at the best efficiency point (BEP) were 0.349 m, 12 m, and 0.2 m3 s-1, respectively. A total ten pressure sensors were mounted at different locations namely, turbine inlet, vaneless space, and draft tube. The data were acquired at a sampling rate of 5 kHz. The instruments and sensors have been calibrated according to guidelines available in IEC standards. The determined total uncertainty in the measurement of hydraulic efficiency was ±0.15% at BEP condition. The velocity measurements in the draft tube cone were performed using a 2D PIV system and the images were sampled at a rate of 40 Hz.      Steady state measurements were carried out considering the realistic design and off-design operating conditions of the prototype turbine. Therefore, the angular speed of the runner was maintained constant for all steady state conditions during the measurements. The maximum hydraulic efficiency (92.4%) was observed at nED = 0.18, QED = 0.15, and a = 9.8º, which is named BEP. It is observed that the turbine experiences significant pressure fluctuations at the vaneless space, runner, and the draft tube. The fluctuations due to rotor-stator interaction (RSI) were observed to be most dominating at high load condition, however, fluctuations due to the rotating vortex rope (RVR) at part load (PL) condition. Two different modes (synchronous and asynchronous) modes of vortex rope are observed at PL condition of the turbine. An asymmetry in the flow leaving the runner was detected at both design and off-design conditions, with a stronger effect during off-design operating condition. Numerical simulations of the model turbine were carried out at PL operating condition. The simulations were performed using two turbulence models, standard k-ε and SST k-ω, with high-resolution advection scheme. The numerical pressure values obtained with both standard k-ε model and SST k-ω showed a small difference with the experimental values. The amplitudes of numerical pressure values were higher (~2.8%) in the vaneless space and lower (~5.0%) in the draft tube than the experimental values. The frequencies of the RSI and RVR were well captured in the turbine but the amplitudes were overestimated.   During load rejection from BEP to PL, the plunging mode of the vortex rope was observed to appear first in the system than that of the rotating mode. Whereas during the load acceptance from PL to BEP, both the modes were observed to disappear simultaneously from the system. In the velocity data, the axial velocity only contributed to the development of the plunging mode and radial velocity to the rotating mode. The region of low velocity, stagnation point, flow separation, recirculation, oscillating flow and high axial velocity gradients were well captured in the system during the transients. The induced high-velocity gradients during the load acceptance from BEP to HL was observed to develop a vortex core in the draft tube. During startup and shutdown, the guide vanes angular position was moved from one to another steady state condition to achieve the minimum load condition of the turbine. At this condition, the generator of the turbine was magnetized at the synchronous speed during startup and shutdown, respectively. The frequency of wave propagation was observed to vary with the runner angular speed during startup and complete shutdown of the turbine. Comparatively high-pressure fluctuations in the draft tube were observed during the guide vane movement from the high discharge conditions. Some unsteady phenomena such as the formation of dead velocity zone, backward flow, and flow oscillations were observed during startup and shutdown of the turbine.   The current work has been also used to continue a series of workshops, i.e., Francis-99. The first workshop was held on December 2014 with the cooperation of LTU and NTNU. The measurements performed in this work were used for the second workshop which was held on December 2016. The investigations presented in this thesis will be further explored in the third workshop scheduled for December 2018.