Flow pattern analysis for magnetic resonance velocity imaging

Blood flow in the heart is highly complex. Although blood flow patterns have been investigated by both computational modelling and invasive/non-invasive imaging techniques, their evolution and intrinsic connection with cardiovascular disease has yet to be explored. Magnetic resonance (MR) velocity imaging provides a comprehensive distribution of multi-directional in vivo flow distribution so that detailed quantitative analysis of flow patterns is now possible. However, direct visualisation or quantification of vector fields is of little clinical use, especially for inter-subject or serial comparison of changes in flow patterns due to the progression of the disease or in response to therapeutic measures. In order to achieve a comprehensive and integrated description of flow in health and disease, it is necessary to characterise and model both normal and abnormal flows and their effects. To accommodate the diversity of flow patterns in relation to morphological and functional changes, we have described in this thesis an approach of detecting salient topological features prior to analytical assessment of dynamical indices of the flow patterns. To improve the accuracy of quantitative analysis of the evolution of topological flow features, it is essential to restore the original flow fields so that critical points associated with salient flow features can be more reliably detected. We propose a novel framework for the restoration, abstraction, extraction and tracking of flow features such that their dynamic indices can be accurately tracked and quantified. The restoration method is formulated as a constrained optimisation problem to remove the effects of noise and to improve the consistency of the MR velocity data. A computational scheme is derived from the First Order Lagrangian Method for solving the optimisation problem. After restoration, flow abstraction is applied to partition the entire flow field into clusters, each of which is represented by a local linear expansion of its velocity components. This process not only greatly reduces the amount of data required to encode the velocity distribution but also permits an analytical representation of the flow field from which critical points associated with salient flow features can be accurately extracted. After the critical points are extracted, phase portrait theory can be applied to separate them into attracting/repelling focuses, attracting/repelling nodes, planar vortex, or saddle. In this thesis, we have focused on vortical flow features formed in diastole. To track the movement of the vortices within a cardiac cycle, a tracking algorithm based on relaxation labelling is employed. The constraints and parameters used in the tracking algorithm are designed using the characteristics of the vortices. The proposed framework is validated with both simulated and in vivo data acquired from patients with sequential MR examination following myocardial infarction. The main contribution of the thesis is in the new vector field restoration and flow feature abstraction method proposed. They allow the accurate tracking and quantification of dynamic indices associated with salient features so that inter- and intra-subject comparisons can be more easily made. This provides further insight into the evolution of blood flow patterns and permits the establishment of links between blood flow patterns and localised genesis and progression of cardiovascular disease.Open acces

Illustrative Flow Visualization of 4D PC-MRI Blood Flow and CFD Data

Das zentrale Thema dieser Dissertation ist die Anwendung illustrativer Methoden auf zwei bisher ungelöste Probleme der Strömungsvisualisierung. Das Ziel der Strömungsvisualisierung ist die Bereitstellung von Software, die Experten beim Auswerten ihrer Strömungsdaten und damit beim Erkenntnisgewinn unterstützt. Bei der illustrativen Visualisierung handelt es sich um einen Zweig der Visualisierung, der sich an der künstlerischen Arbeit von Illustratoren orientiert. Letztere sind darauf spezialisiert komplizierte Zusammenhänge verständlich und ansprechend zu vermitteln. Die angewendeten Techniken werden in der illustrativen Visualisierung auf reale Daten übertragen, um die Effektivität der Darstellung zu erhöhen. Das erste Problem, das im Rahmen dieser Dissertation bearbeitet wurde, ist die eingeschränkte Verständlichkeit von komplexen Stromflächen. Selbstverdeckungen oder Aufrollungen behindern die Form- und Strömungswahrnehmung und machen diese Flächen gerade in interessanten Strömungssituationen wenig nützlich. Auf Basis von handgezeichneten Strömungsdarstellungen haben wir ein Flächenrendering entwickelt, das Silhouetten, nicht-photorealistische Beleuchtung und illustrative Stromlinien verwendet. Interaktive Flächenschnitte erlauben die Exploration der Flächen und der Strömungen, die sie repräsentieren. Angewendet auf verschiedene Stromflächen ließ sich zeigen, dass die Methoden die Verständlichkeit erhöhen, v.a. in Bereichen komplexer Strömung mit Aufwicklungen oder Singularitäten. Das zweite Problem ist die Strömungsanalyse des Blutes aus 4D PC-MRI-Daten. An diese relativ neue Datenmodalität werden hohe Erwartungen für die Erforschung und Behandlung kardiovaskulärer Krankheiten geknüpft, da sie erstmals ein dreidimensionales, zeitlich aufgelöstes Abbild der Hämodynamik liefert. Bisher werden 4D PC-MRI-Daten meist mit Werkzeugen der klassischen Strömungsvisualisierung verarbeitet. Diese werden den besonderen Ansprüchen der medizinischen Anwender jedoch nicht gerecht, die in kurzer Zeit eine übersichtliche Darstellung der relevanten Strömungsaspekte erhalten möchten. Wir haben ein Werkzeug zur visuellen Analyse der Blutströmung entwickelt, welches eine einfache Detektion von markanten Strömungsmustern erlaubt, wie z.B. Jets, Wirbel oder Bereiche mit hoher Blutverweildauer. Die Grundidee ist hierbei aus vorberechneten Integrallinien mit Hilfe speziell definierter Linienprädikate die relevanten, d.h. am gefragten Strömungsmuster, beteiligten Linien ausgewählt werden. Um eine intuitive Darstellung der Resultate zu erreichen, haben wir uns von Blutflußillustrationen inspirieren lassen und präsentieren eine abstrakte Linienbündel- und Wirbeldarstellung. Die Linienprädikatmethode sowie die abstrakte Darstellung der Strömungsmuster wurden an 4D PC-MRI-Daten von gesunden und pathologischen Aorten- und Herzdaten erfolgreich getestet. Auch die Evaluierung durch Experten zeigt die Nützlichkeit der Methode und ihr Potential für den Einsatz in der Forschung und der Klinik.This thesis’ central theme is the use of illustrative methods to solve flow visualization problems. The goal of flow visualization is to provide users with software tools supporting them analyzing and extracting knowledge from their fluid dynamics data. This fluid dynamics data is produced in large amounts by simulations or measurements to answer diverse questions in application fields like engineering or medicine. This thesis deals with two unsolved problems in flow visualization and tackles them with methods of illustrative visualization. The latter is a subbranch of visualization whose methods are inspired by the art work of professional illustrators. They are specialized in the comprehensible and esthetic representation of complex knowledge. With illustrative visualization, their techniques are applied to real data to enhance their representation. The first problem dealt with in this thesis is the limited shape and flow perception of complex stream surfaces. Self-occlusion and wrap-ups hinder their effective use in the most interesting flow situations. On the basis of hand-drawn flow illustrations, a surface rendering method was designed that uses silhouettes, non-photorealistic shading, and illustrative surface stream lines. Additionally, geometrical and flow-based surface cuts allow the user an interactive exploration of the surface and the flow it represents. By applying this illustrative technique to various stream surfaces and collecting expert feedback, we could show that the comprehensibility of the stream surfaces was enhanced – especially in complex areas with surface wrap-ups and singularities. The second problem tackled in this thesis is the analysis of blood flow from 4D PC-MRI data. From this rather young data modality, medical experts expect many advances in the research of cardiovascular diseases because it delivers a three-dimensional and time-resolved image of the hemodynamics. However, 4D PC-MRI data are mainly processed with standard flow visualizaton tools, which do not fulfill the requirements of medical users. They need a quick and easy-to-understand display of the relevant blood flow aspects. We developed a tool for the visual analysis of blood flow that allows a fast detection of distinctive flow patterns, such as high-velocity jets, vortices, or areas with high residence times. The basic idea is to precalculate integral lines and use specifically designed line predicates to select and display only lines involved in the pattern of interest. Traditional blood flow illustrations inspired us to an abstract and comprehensible depiction of the resulting line bundles and vortices. The line predicate method and the illustrative flow pattern representation were successfully tested with 4D PC-MRI data of healthy and pathological aortae and hearts. Also, the feedback of several medical experts confirmed the usefulness of our methods and their capabilities for a future application in the clinical research and routine

NASA/ASEE Summer Faculty Fellowship Program

The contractor's report contains all sixteen final reports prepared by the participants in the 1989 Summer Faculty Fellowship Program. Reports describe research projects on a number of different topics. Interface software, metal corrosion, rocket triggering lightning, automatic drawing, 60-Hertz power, carotid-cardiac baroreflex, acoustic fields, robotics, AI, CAD/CAE, cryogenics, titanium, and flow measurement are discussed

Automated feature extraction in oceanographic visualization

Thesis (S.M.)--Massachusetts Institute of Technology, Dept. of Ocean Engineering; and, (S.M.)--Massachusetts Institute of Technology, Dept. of Electrical Engineering and Computer Science, 2004.Includes bibliographical references (leaves 141-147).The ocean is characterized by a multitude of powerful, sporadic biophysical dynamical events; scientific research has reached the stage that their interpretation and prediction is now becoming possible. Ocean prediction, analogous to atmospheric weather prediction but combining biological, chemical and physical features is able to help us understand the complex coupled physics, biology and acoustics of the ocean. Applications of the prediction of the ocean environment include exploitation and management of marine resources, pollution control such as planning of maritime and naval operations. Given the vastness of ocean, it is essential for effective ocean prediction to employ adaptive sampling to best utilize the available sensor resources in order to minimize the forecast error. It is important to concentrate measurements to the regions where one can witness features of physical or biological significance in progress. Thus automated feature extraction in oceanographic visualization can facilitate adaptive sampling by presenting the physically relevant features directly to the operation planners. Moreover it could be used to help automate adaptive sampling. Vortices (eddies and gyres) and upwelling, two typical and important features of the ocean, are studied.(cont.) A variety of feature extraction methods are presented, and those more pertinent to this study are implemented, including derived field generation and attribute set extraction. Detection results are evaluated in terms of accuracy, computational efficiency, clarity and usability. Vortices, a very important flow feature is the primary focus of this study. Several point-based and set-based vortex detection methods are reviewed. A set-based vortex core detection method based on geometric properties of vortices is applied to both classical vortex models and real ocean models. The direction spanning property, which is a geometric property, guides the detection of all the vortex core candidates, and the conjugate pair eigenvalue method is responsible for filtering out the false positives from the candidate set. Results show the new method to be analytically accurate and practically feasible, and superior to traditional point-based vortex detection methods. Detection methods of streamlines are also discussed. Using the novel cross method or the winding angle method, closed streamlines around vortex cores can be detected.(cont.) Therefore, the whole vortex area, i.e., the combination of vortex core and surrounding streamlines, is detected. Accuracy and feasibility are achieved through automated vortex detection requiring no human inspection. The detection of another ocean feature, upwelling, is also discussed.by Da Guo.S.M

Numerical Investigation of Liquid Film Dynamics and Atomisation in Jet Engine Fuel Injectors

Today’s aerospace industry continues to exploit liquid hydrocarbon fossil fuels. Motivated by operational considerations, continued availability and cost, this is likely to be the case for many years, despite the obvious environmental concerns. The interplay of liquid atomisation, spray vaporisation and the combustion process are intricately linked. However, the physical process of fuel injection and its atomisation into tiny droplets prior to combustion remains poorly understood. Because atomisation governs the size of the fuel droplets, and therefore their subsequent evaporation rate, adjusting the injection sequence is of paramount importance and will have far-reaching repercussions on many aspects of the combustion process, for example pollutant formation. In the context of jet engines, kerosene is usually injected in its liquid form via an airblast-type fuel injector. A coflowing high-speed airstream destabilises the liquid fuel, which is thus sprayed into fine droplets into the combustion chamber. The prediction of this phenomenon for various operating conditions relevant to the aeronautical industry requires a deeper understanding of the mechanisms involved in the interaction of the two fluids. A key element in predicting the complex behaviour of spray formation and evolution in jet engines is accurate modelling of fuel atomisation. Atomisation represents one of the key challenges that remains to be undertaken to make predictive computational simulations possible. However, the inherent multi-physics and multi-scale nature of this process limits numerical investigations. Thanks to the steady progress in computer power and Computational Fluid Dynamics (CFD) methods, computational modelling of injection systems emerges as a promising tool that can drive the design of future devices. This research project sets out to investigate the atomisation process in detail, in particular in providing physical insight into the fundamental physics of the phenomenon, in conjunction with an analysis on wetting behaviours and liquid droplet tracking. High-fidelity numerical simulations are performed using a novel in-house state-of-the-art multiphase flow modelling capability, RCLSFoam. The performance of the numerical scheme is demonstrated on typical two-dimensional and three-dimensional benchmark test cases relevant to both multiphase flow modelling and atomisation, and validated against other computational methods. An informed and systematic qualitative assessment of the topological variations of the phase interface during primary atomisation of a liquid film is made through dynamical analysis, while investigating an extensive domain of operating conditions at ambient and aero-engine injection conditions relevant to industry. This analysis demonstrated the influence of shear-driven instabilities on the atomisation process. The shear stress and difference in inertia between liquid and gas are observed to play a significant role in the atomisation process. In addition, the key physical mechanisms and their competing effects have been mapped out in order to predict the evolution of the process according to the operating conditions of the injection system. The proposed cartography gathers four different atomisation mechanisms. In particular, for sufficiently high liquid injection speeds, three-dimensional wave modes were observed to co-exist (the “3-D wave mode” regime). For very low liquid flow rates, accumulated liquid at the atomising edge undergoes deformation by which droplets are generated (the “accumulation” regime). For an increasing gas injection speed and a fixed liquid velocity, the effects of surface tension were observed to result in the generation of streamwise ligaments only, which tend to pair up (the “ligament-merging” regime). Finally, “vortex action” is another observed mechanism by which the liquid film is fragmented. Overall, this research project culminated in (i) the study of dynamic wetting behaviours, with the implementation and validation against experimental data of the Kistler dynamic contact model; and (ii) the demonstration of an algorithm for droplet capture and subsequent post-processing analysis of the droplet characteristics.Rolls Royce plc. and EPSR