A subpixel edge detector applied to aortic dissection detection

[EN] The aortic dissection is a disease that can cause a deadly situation, even with a correct treatment. It consists in a rupture of a layer of the aortic artery wall, causing a blood flow inside this rupture, called dissection. The aim of this paper is to contribute to its diagnosis, detecting the dissection edges inside the aorta. A subpixel accuracy edge detector based on the hypothesis of partial volume effect is used, where the intensity of an edge pixel is the sum of the contribution of each color weighted by its relative area inside the pixel. The method uses a floating window centred on the edge pixel and computes the edge features. The accuracy of our method is evaluated on synthetic images of different hickness and noise levels, obtaining an edge detection with a maximal mean error lower than 16 percent of a pixel.

Angiogenesis and angioregression gene expression analyses in swine corpus luteum

The corpus luteum (CL) lifespan is characterized by a rapid growth, differentiation and controlled regression of the luteal tissue, accompanied by an intense angiogenesis and angioregression. Indeed, the CL is one of the most highly vascularised tissue in the body with a proliferation rate of the endothelial cells 4- to 20-fold more intense than in some of the most malignant human tumours. This angiogenic process should be rigorously controlled to allow the repeated opportunities of fertilization. After a first period of rapid growth, the tissue becomes stably organized and prepares itself to switch to the phenotype required for its next apoptotic regression. In pregnant swine, the lifespan of the CLs must be extended to support embryonic and foetal development and vascularisation is necessary for the maintenance of luteal function. Among the molecules involved in the angiogenesis, Vascular Endothelial Growth Factor (VEGF) is the main regulator, promoting endothelial cells proliferation, differentiation and survival as well as vascular permeability and vessel lumen formation. During vascular invasion and apoptosis process, the remodelling of the extracellular matrix is essential for the correct evolution of the CL, particularly by the action of specific class of proteolytic enzymes known as matrix metalloproteinases (MMPs). Another important factor that plays a role in the processes of angiogenesis and angioregression during the CL formation and luteolysis is the isopeptide Endothelin-1 (ET-1), which is well-known to be a potent vasoconstrictor and mitogen for endothelial cells. The goal of the present thesis was to study the role and regulation of vascularisation in an adult vascular bed. For this purpose, using a precisely controlled in vivo model of swine CL development and regression, we determined the levels of expression of the members of VEGF system (VEGF total and specific isoforms; VEGF receptor-1, VEGFR-1; VEGF receptor-2, VEGFR-2) and ET- 1 system (ET-1; endothelin converting enzyme-1, ECE-1; endothelin receptor type A, ET-A) as well as the activity of the Ca++/Mg++-dependent endonucleases and gelatinases (MMP-2 and MMP-9). Three experiments were conducted to reach such objectives in CLs isolated from ovaries of cyclic, pregnant or fasted gilts. In the Experiment I, we evaluated the influence of acute fasting on VEGF production and VEGF, VEGFR-2, ET-1, ECE-1 and ET-A mRNA expressions in CLs collected on day 6 after ovulation (midluteal phase). The results indicated a down-regulation of VEGF, VEGFR-2, ET-1 and ECE-1 mRNA expression, although no change was observed for VEGF protein. Furthermore, we observed that fasting stimulated steroidogenesis by luteal cells. On the basis of the main effects of VEGF (stimulation of vessel growth and endothelial permeability) and ET-1 (stimulation of endothelial cell proliferation and vasoconstriction, as well as VEGF stimulation), we concluded that feed restriction possibly inhibited luteal vessel development. This could be, at least in part, compensated by a decrease of vasal tone due to a diminution of ET-1, thus ensuring an adequate blood flow and the production of steroids by the luteal cells. In the Experiment II, we investigated the relationship between VEGF, gelatinases and Ca++/Mg++-dependent endonucleases activities with the functional CL stage throughout the oestrous cycle and at pregnancy. The results demonstrated differential patterns of expression of those molecules in correspondence to the different phases of the oestrous cycle. Immediately after ovulation, VEGF mRNA/protein levels and MMP-9 activity are maximal. On days 5–14 after ovulation, VEGF expression and MMP-2 and -9 activities are at basal levels, while Ca++/Mg++-dependent endonuclease levels increased significantly in relation to day 1. Only at luteolysis (day 17), Ca++/Mg++-dependent endonuclease and MMP-2 spontaneous activity increased significantly. At pregnancy, high levels of MMP-9 and VEGF were observed. These results suggested that during the very early luteal phase, high MMPs activities coupled with high VEGF levels drive the tissue to an angiogenic phenotype, allowing CL growth under LH (Luteinising Hormone) stimulus, while during the late luteal phase, low VEGF and elevate MMPs levels may play a role in the apoptotic tissue and extracellular matrix remodelling during structural luteolysis. In the Experiment III, we described the expression patterns of all distinct VEGF isoforms throughout the oestrous cycle. Furthermore, the mRNA expression and protein levels of both VEGF receptors were also evaluated. Four novel VEGF isoforms (VEGF144, VEGF147, VEGF182, and VEGF164b) were found for the first time in swine and the seven identified isoforms presented four different patterns of expression. All isoforms showed their highest mRNA levels in newly formed CLs (day 1), followed by a decrease during mid-late luteal phase (days 10–17), except for VEGF182, VEGF188 and VEGF144 that showed a differential regulation during late luteal phase (day 14) or at luteolysis (day 17). VEGF protein levels paralleled the most expressed and secreted VEGF120 and VEGF164 isoforms. The VEGF receptors mRNAs showed a different pattern of expression in relation to their ligands, increasing between day 1 and 3 and gradually decreasing during the mid-late luteal phase. The differential regulation of some VEGF isoforms principally during the late luteal phase and luteolysis suggested a specific role of VEGF during tissue remodelling process that occurs either for CL maintenance in case of pregnancy or for noncapillary vessel development essential for tissue removal during structural luteolysis. In summary, our findings allow us to determine relationships among factors involved in the angiogenesis and angioregression mechanisms that take place during the formation and regression of the CL. Thus, CL provides a very interesting model for studying such factors in different fields of the basic research.

Mechanosensitivity in the myenteric plexus of the guinea pig ileum

The enteric nervous system regulates autonomously from the central nervous system all the reflex pathways that control blood flow, motility, water and electrolyte transport and acid secretion. The ability of the gut to function in isolation is one of the most intriguing phenomenons in neurogastroenterology. This requires coding of sensory stimuli by cells in the gut wall. Enteric neurons are prominent candidates to relay mechanosensitivity. Surprisingly, the identity of mechanosensitive neurons in the enteric nervous system as well as the appropriate stimulus modality is unknown despite the evidence that enteric neurons respond to sustained distension. Objectives: The aim of our study was to record from mechanosensitive neurons using physiological stimulus modalities. Identification of sensory neurons is of central importance to understand sensory transmission under normal conditions and in gut diseases associated with sensorimotor dysfunctions, such as Irritable Bowel Syndrome. Only then it will be possible to identify novel targets that help to normalise sensory functions. Methods: We used guinea-pig ileum myenteric plexus preparations and recorded responses of all neurons in a given ganglion with a fast neuroimaging technique based on voltage sensitive dyes. To evoke a mechanical response we used two different kinds of stimuli: firstly we applied a local mechanical distortion of the ganglion surface with von Frey hair. Secondarily we mimic the ganglia deformation during physiological movements of myenteric ganglia in a freely contracting ileal preparation. We were able to reliably and reproducibly mimic this distortion by intraganglionic injections of small volumes of oxygenated and buffered Krebs solution using stimulus parameters that correspond to single contractions. We also performed in every ganglion tested, electrical stimulations to evoke fast excitatory postsynaptic potentials. Immunohistochemistry reactions were done with antibodies against Calbindin and NeuN, considered markers for sensory neurons. Results: Recordings were performed in 46 ganglia from 31 guinea pigs. In every ganglion tested we found from 1 to 21 (from 3% to 62%) responding cells with a median value of 7 (24% of the total number of neurons). The response consisted of an almost instantaneous spike discharge that showed adaptation. The median value of the action potential frequency in the responding neurons was 2.0 Hz, with a recording time of 1255 ms. The spike discharge lasted for 302 ± 231 ms and occurred only during the initial deformation phase. During sustained deformation no spike discharge was observed. The response was reproducible and was a direct activation of the enteric neurons since it remained after synaptic blockade with hexamethonium or ω-conotoxin and after long time perfusion with capsaicin. Muscle tone appears not to be required for activation of mechanosensory neurons. Mechanosensory neurons showed a response to mechanical stimulation related to the stimulus strength. All mechanosensory neurons received fast synaptic inputs. There was no correlation between mechanosensitivity and Calbindin-IR and NeuN-IR (44% of mechanosensitive neurones Calb-IR-/NeuN-IR-). Conclusions: We identified mechanosensitive neurons in the myenteric plexus of the guinea pig ileum which responded to brief deformation. These cells appear to be rapidly accommodating neurons which respond to dynamic change. All mechanosensitive neurons received fast synaptic input suggesting that their activity can be highly modulated by other neurons and hence there is a low stimulus fidelity which allows adjusting the gain in a sensory network. Mechanosensitivity appears to be a common feature of many enteric neurons belonging to different functional classes. This supports the existence of multifunctional enteric neurons which may fulfil sensory, integrative and motor functions.

Outcome respiratorio a lungo termine nei soggetti affetti da cardiopatia congenita sottoposti ad intervento cardiochirurgico

Introduction: In the last years cardiac surgery for congenital heart disease (CHD) reduced dramatically mortality modifying prognosis, but, at the same time, increased morbidity in this patient population. Respiratory and cardiovascular systems are strictly anatomically and functionally connected, so that alterations of pulmonary hemodynamic conditions modify respiratory function. While very short-term alterations of respiratory mechanics after surgery were investigated by many authors, not as much works focused on long-term changes. In these subjects rest respiratory function may be limited by several factor: CHD itself (fetal pulmonary perfusion influences vascular and alveolar development), extracorporeal circulation (CEC), thoracotomy and/or sternotomy, rib and sternal contusions, pleural adhesions and pleural fibrosis, secondary to surgical injury. Moreover inflammatory cascade, triggered by CEC, can cause endothelial damage and compromise gas exchange. Aims: The project was conceived to 1) determine severity of respiratory functional impairement in different CHD undergone to surgical correction/palliation; 2) identify the most and the least CHD involved by pulmonary impairement; 3) find a correlation between a specific hemodynamic condition and functional anomaly, and 4) between rest respiratory function and cardiopulmonary exercise test. Materials and methods: We studied 113 subjects with CHD undergone to surgery, and distinguished by group in accord to pulmonary blood flow (group 0: 28 pts with normal pulmonary flow; group 1: 22 pts with increased flow; group 2: 43 pts with decreased flow; group 3: 20 pts with total cavo-pulmonary anastomosis-TCPC) followed by the Pediatric Cardiology and Cardiac Surgery Unit, and we compare them to 37 age- and sex-matched healthy subjects. In Pediatric Pulmonology Unit all pts performed respiratory function tests (static and dynamic volumes, flow/volume curve, airway resistances-raw- and conductance-gaw-, lung diffusion of CO-DLCO- and DLCO/alveolar volume), and CHD pts the same day had cardiopulmonary test. They all were examined and had allergological tests, and respiratory medical history. Results: restrictive pattern (measured on total lung capacity-TLC- and vital capacity-VC) was in all CHD groups, and up to 45% in group 2 and 3. Comparing all groups, we found a significant difference in TLC between healthy and group 2 (p=0.001) and 3 (p=0.004), and in VC between group 2 and healthy (p=0.001) and group 1(p=0.034). Inspiratory capacity (IC) was decreased in group 2 related to healthy (p<0.001) and group 1 (p=0.037). We showed a direct correlation between TLC and VC with age at surgery (p=0.01) and inverse with number of surgical interventions (p=0.03). Reduced FEV1/FVC ratio, Gaw and increased Raw were mostly present in group 3. DLCO was impaired in all groups, but up to 80% in group 3 and 50% in group 2; when corrected for alveolar volume (DLCO/VA) reduction persisted in group 3 (20%), 2 (6.2%) and 0 (7.1%). Exercise test was impaired in all groups: VO2max and VE markedly reduced in all but especially in group 3, and VE/VCO2 slope, marker of ventilatory response to exercise, is increased (<36) in 62.5% of group 3, where other pts had anyway value>32. Comparing group 3 and 2, the most involved categories, we found difference in VO2max and VE/VCO2 slope (respectively p=0.02 and p<0.0001). We evidenced correlation between rest and exercise tests, especially in group 0 (between VO2max and FVC, FEV1, VC, IC; inverse relation between VE/VCO2slope and FVC, FEV1 and VC), but also in group 1 (VO2max and IC), group 2 (VO2max and FVC and FEV1); never in group 3. Discussion: According with literature, we found a frequent impairment of rest pulmonary function in all groups, but especially in group 2 and 3. Restrictive pattern was the most frequent alteration probably due to compromised pulmonary (vascular and alveolar) development secondary to hypoperfusion in fetal and pre-surgery (and pre-TCPC)life. Parenchymal fibrosis, pleural adhesions and thoracic deformities can add further limitation, as showed by the correlation between group 3 and number of surgical intervention. Exercise tests were limited, particularly in group 3 (complex anatomy and lost of chronotropic response), and we found correlations between rest and exercise tests in all but group 3. We speculate that in this patients hemodynamic exceeds respiratory contribution, though markedly decreased.

Fluidodinamica biomedica computazionale del flusso sanguigno in presenza di affezioni strumentali e/o patologiche

La tesi di Dottorato studia il flusso sanguigno tramite un codice agli elementi finiti (COMSOL Multiphysics). Nell’arteria è presente un catetere Doppler (in posizione concentrica o decentrata rispetto all’asse di simmetria) o di stenosi di varia forma ed estensione. Le arterie sono solidi cilindrici rigidi, elastici o iperelastici. Le arterie hanno diametri di 6 mm, 5 mm, 4 mm e 2 mm. Il flusso ematico è in regime laminare stazionario e transitorio, ed il sangue è un fluido non-Newtoniano di Casson, modificato secondo la formulazione di Gonzales & Moraga. Le analisi numeriche sono realizzate in domini tridimensionali e bidimensionali, in quest’ultimo caso analizzando l’interazione fluido-strutturale. Nei casi tridimensionali, le arterie (simulazioni fluidodinamiche) sono infinitamente rigide: ricavato il campo di pressione si procede quindi all’analisi strutturale, per determinare le variazioni di sezione e la permanenza del disturbo sul flusso. La portata sanguigna è determinata nei casi tridimensionali con catetere individuando tre valori (massimo, minimo e medio); mentre per i casi 2D e tridimensionali con arterie stenotiche la legge di pressione riproduce l’impulso ematico. La mesh è triangolare (2D) o tetraedrica (3D), infittita alla parete ed a valle dell’ostacolo, per catturare le ricircolazioni. Alla tesi sono allegate due appendici, che studiano con codici CFD la trasmissione del calore in microcanali e l’ evaporazione di gocce d’acqua in sistemi non confinati. La fluidodinamica nei microcanali è analoga all’emodinamica nei capillari. Il metodo Euleriano-Lagrangiano (simulazioni dell’evaporazione) schematizza la natura mista del sangue. La parte inerente ai microcanali analizza il transitorio a seguito dell’applicazione di un flusso termico variabile nel tempo, variando velocità in ingresso e dimensioni del microcanale. L’indagine sull’evaporazione di gocce è un’analisi parametrica in 3D, che esamina il peso del singolo parametro (temperatura esterna, diametro iniziale, umidità relativa, velocità iniziale, coefficiente di diffusione) per individuare quello che influenza maggiormente il fenomeno.

Inverse problems in classical and quantum physics

The subject of this thesis is in the area of Applied Mathematics known as Inverse Problems. Inverse problems are those where a set of measured data is analysed in order to get as much information as possible on a model which is assumed to represent a system in the real world. We study two inverse problems in the fields of classical and quantum physics: QCD condensates from tau-decay data and the inverse conductivity problem. Despite a concentrated effort by physicists extending over many years, an understanding of QCD from first principles continues to be elusive. Fortunately, data continues to appear which provide a rather direct probe of the inner workings of the strong interactions. We use a functional method which allows us to extract within rather general assumptions phenomenological parameters of QCD (the condensates) from a comparison of the time-like experimental data with asymptotic space-like results from theory. The price to be paid for the generality of assumptions is relatively large errors in the values of the extracted parameters. Although we do not claim that our method is superior to other approaches, we hope that our results lend additional confidence to the numerical results obtained with the help of methods based on QCD sum rules. EIT is a technology developed to image the electrical conductivity distribution of a conductive medium. The technique works by performing simultaneous measurements of direct or alternating electric currents and voltages on the boundary of an object. These are the data used by an image reconstruction algorithm to determine the electrical conductivity distribution within the object. In this thesis, two approaches of EIT image reconstruction are proposed. The first is based on reformulating the inverse problem in terms of integral equations. This method uses only a single set of measurements for the reconstruction. The second approach is an algorithm based on linearisation which uses more then one set of measurements. A promising result is that one can qualitatively reconstruct the conductivity inside the cross-section of a human chest. Even though the human volunteer is neither two-dimensional nor circular, such reconstructions can be useful in medical applications: monitoring for lung problems such as accumulating fluid or a collapsed lung and noninvasive monitoring of heart function and blood flow.

Entwicklungen zur Absolutquantifizierung der Lungenperfusion mittels kontrastmittelverstärkter Magnetresonanztomographie

Die regionale Bestimmung der Durchblutung (Perfusion) ermöglicht differenzierte Aussagen über den Gesundheitszustand und die Funktionalität der Lunge. Durch neue Messverfahren ermöglicht die Magnetresonanztomographie (MRT) eine nicht-invasive und strahlungsfreie Untersuchung der Perfusion. Obwohl die Machbarkeit qualitativer MRT-Durchblutungsmessungen bereits gezeigt wurde, fehlt bisher eine validierte quantitative Methode. Ziel dieser Arbeit war eine Optimierung der bestehenden Messprotokolle und mathematischen Modelle zur Absolutquantifizierung der Lungenperfusion mit Magnetresonanztomographie. Weiterhin sollte die Methodik durch Vergleich mit einem etablierten Referenzverfahren validiert werden. Durch Simulationen und Phantommessungen konnten optimale MRT-Messparameter und ein standardisiertes Protokoll festgelegt werden. Des Weiteren wurde eine verallgemeinerte Bestimmung der Kontrastmittelkonzentration aus den gemessenen Signalintensitäten vorgestellt, diskutiert und durch Probandenmessungen validiert. Auf der Basis dieser Entwicklungen wurde die MRT-Durchblutungsmessung der Lunge tierexperimentell mit der Positronenemissionstomographie (PET) intraindividuell verglichen und validiert. Die Ergebnisse zeigten nur kleine Abweichungen und eine statistisch hochsignifikante, stark lineare Korrelation. Zusammenfassend war es durch die Entwicklungen der vorgestellten Arbeit möglich, die kontrastmittelgestützte MRT-Durchblutungsmessung der Lunge zu optimieren und erstmals zu validieren.

Lesion detection and classification in breast cancer: evaluation of approaches based on morphological features, tracer kinetic modelling and semi-quantitative parameters in MR functional imaging (DCE-MRI)

The diagnosis, grading and classification of tumours has benefited considerably from the development of DCE-MRI which is now essential to the adequate clinical management of many tumour types due to its capability in detecting active angiogenesis. Several strategies have been proposed for DCE-MRI evaluation. Visual inspection of contrast agent concentration curves vs time is a very simple yet operator dependent procedure, therefore more objective approaches have been developed in order to facilitate comparison between studies. In so called model free approaches, descriptive or heuristic information extracted from time series raw data have been used for tissue classification. The main issue concerning these schemes is that they have not a direct interpretation in terms of physiological properties of the tissues. On the other hand, model based investigations typically involve compartmental tracer kinetic modelling and pixel-by-pixel estimation of kinetic parameters via non-linear regression applied on region of interests opportunely selected by the physician. This approach has the advantage to provide parameters directly related to the pathophysiological properties of the tissue such as vessel permeability, local regional blood flow, extraction fraction, concentration gradient between plasma and extravascular-extracellular space. Anyway, nonlinear modelling is computational demanding and the accuracy of the estimates can be affected by the signal-to-noise ratio and by the initial solutions. The principal aim of this thesis is investigate the use of semi-quantitative and quantitative parameters for segmentation and classification of breast lesion. The objectives can be subdivided as follow: describe the principal techniques to evaluate time intensity curve in DCE-MRI with focus on kinetic model proposed in literature; to evaluate the influence in parametrization choice for a classic bi-compartmental kinetic models; to evaluate the performance of a method for simultaneous tracer kinetic modelling and pixel classification; to evaluate performance of machine learning techniques training for segmentation and classification of breast lesion.

CT perfusion image processing: analysis of liver tumors

Perfusion CT imaging of the liver has potential to improve evaluation of tumour angiogenesis. Quantitative parameters can be obtained applying mathematical models to Time Attenuation Curve (TAC). However, there are still some difficulties for an accurate quantification of perfusion parameters due, for example, to algorithms employed, to mathematical model, to patient’s weight and cardiac output and to the acquisition system. In this thesis, new parameters and alternative methodologies about liver perfusion CT are presented in order to investigate the cause of variability of this technique. Firstly analysis were made to assess the variability related to the mathematical model used to compute arterial Blood Flow (BFa) values. Results were obtained implementing algorithms based on “ maximum slope method” and “Dual input one compartment model” . Statistical analysis on simulated data demonstrated that the two methods are not interchangeable. Anyway slope method is always applicable in clinical context. Then variability related to TAC processing in the application of slope method is analyzed. Results compared with manual selection allow to identify the best automatic algorithm to compute BFa. The consistency of a Standardized Perfusion Index (SPV) was evaluated and a simplified calibration procedure was proposed. At the end the quantitative value of perfusion map was analyzed. ROI approach and map approach provide related values of BFa and this means that pixel by pixel algorithm give reliable quantitative results. Also in pixel by pixel approach slope method give better results. In conclusion the development of new automatic algorithms for a consistent computation of BFa and the analysis and definition of simplified technique to compute SPV parameter, represent an improvement in the field of liver perfusion CT analysis.

La neoangiogenesi nel reimpianto autologo di tessuto ovarico in pazienti a rischio di fallimento ovarico precoce: ottimizzazione dei metodi di crioconservazione di corticle ovarica

Introduzione L’efficacia dei chemio/radioterapici ha aumentato notevolmente l’aspettativa di vita delle pazienti oncologiche, tuttavia, questi trattamenti possono compromettere la funzionalità ovarica. La crioconservazione di tessuto ovarico, con il successivo reimpianto, ha lo scopo di preservare la fertilità delle pazienti a rischio di fallimento ovarico precoce. Scopo dello studio Definire la migliore procedura di crioconservazione e reimpianto in grado di ottenere la neovascolarizzazione del tessuto reimpiantato nel minor tempo possibile al fine di diminuire la perdita follicolare causata dall’ischemia durante la procedura. Materiali e metodi Per ciascuna paziente (3) le biopsie ovariche, sono state prelevate laparoscopicamente e crioconservate secondo il protocollo di congelamento lento/scongelamento rapido. Campioni di corticale ovarica sono stati processati per l’analisi istologica, ultrastrutturale, immuistochimica e confocale per valutare la preservazione morfologiaca del tessuto. Le fettine di corticale ovarica sono state scongelate e reimpiantate ortotopicamente (2), nelle ovaia e in due tasche peritoneali, o eterotopicamente (1), in due tasche create nel sottocute sovrapubico. Risultati Le analisi di microscopia hanno mostrato il mantenimento di una discreta morfologia dello stroma, e dei vasi criopreservati e un lieve ma non significativo danneggiamento dei follicoli scongelati. Tutte le pazienti hanno mostrato la ripresa della funzionalità endocrina rispettivamente dopo 2/4 mesi dal reimpianto. Il color-doppler, inoltre ha rivelato un significativo aumento della vascolarizzazione ovarica rispetto alla quasi totale assenza di vascolarizzazione prima del reimpianto, quando le pazienti mostravano una conclamata menopausa. Conclusioni Lo studio ha confermato la ripresa della vascolarizzazione dell’ovaio in seguito a reimpianto avascolare di fettine di corticale, senza l’impiego di fattori esogeni o meccanici aggiuntivi, in tempi concordanti con i dati della letteratura. I risultati sono incoraggianti e l’avanzare degli studi e della ricerca potranno contribuire allo sviluppo di nuove metodologie di reimpianto che possano avere un successo clinico ed una sicurezza superiori a quelle finora ottenute.

Fluid-structure interaction by a coupled lattice Boltzmann-finite element approach

In this thesis, a strategy to model the behavior of fluids and their interaction with deformable bodies is proposed. The fluid domain is modeled by using the lattice Boltzmann method, thus analyzing the fluid dynamics by a mesoscopic point of view. It has been proved that the solution provided by this method is equivalent to solve the Navier-Stokes equations for an incompressible flow with a second-order accuracy. Slender elastic structures idealized through beam finite elements are used. Large displacements are accounted for by using the corotational formulation. Structural dynamics is computed by using the Time Discontinuous Galerkin method. Therefore, two different solution procedures are used, one for the fluid domain and the other for the structural part, respectively. These two solvers need to communicate and to transfer each other several information, i.e. stresses, velocities, displacements. In order to guarantee a continuous, effective, and mutual exchange of information, a coupling strategy, consisting of three different algorithms, has been developed and numerically tested. In particular, the effectiveness of the three algorithms is shown in terms of interface energy artificially produced by the approximate fulfilling of compatibility and equilibrium conditions at the fluid-structure interface. The proposed coupled approach is used in order to solve different fluid-structure interaction problems, i.e. cantilever beams immersed in a viscous fluid, the impact of the hull of the ship on the marine free-surface, blood flow in a deformable vessels, and even flapping wings simulating the take-off of a butterfly. The good results achieved in each application highlight the effectiveness of the proposed methodology and of the C++ developed software to successfully approach several two-dimensional fluid-structure interaction problems.

Studio delle applicazioni metodologiche di analisi d'immagini in tc perfusionali

La Tomografia Computerizzata perfusionale (TCp) è una tecnica che permette di effettuare studi funzionali di natura emodinamica dei tessuti ispezionati, misurandone la perfusione attraverso l’analisi temporale della diffusione del mezzo di contrasto. Ricerche recenti indicano come alcuni parametri (ad esempio, il blood flow) della perfusione ematica dell’organo interessato rappresentino un valido strumento per la valutazione dell’efficacia delle terapie anti-angiogenetiche. Valutazione che, potendo avvenire già in fase precoce, senza attendere i più tardivi cambiamenti dimensionali, consente un eventuale adeguamento, altrettanto precoce, della terapia in caso di risultati parziali o di insuccesso. Tuttavia diversi problemi, tra cui la difficoltà ad ottenere misure riproducibili, dovuta alla estrema variabilità delle cause di errore, costituiscono un ostacolo alla standardizzazione e, conseguentemente, alla trasposizione in clinica della TCp. Il lavoro di tesi ha avuto come obiettivo quello di identificare le principali sorgenti di errore in tutto il processo di preparazione dell’esame, generazione e archiviazione dei dati TCp al fine di fornire a IRST-IRCCS un protocollo perfusionale che consente di tenere sotto controllo, innanzitutto, tutti i parametri relativi alla somministrazione del mezzo di contrasto e alla fase di acquisizione delle immagini. Successivamente, è stato creato un catalogo di esami perfusionali, facilmente consultabile sia dal personale medico sia da quello ingegneristico. Infine, è stato fornito il supporto, sia ai medici sia agli ingegneri, per la messa a punto di un metodo di validazione visiva dei risultati ottenuti. Contestualmente, uno dei principali risultati del lavoro di Tesi è stato quello di creare, all’interno di IRST-IRCCS, una figura che funga da “collettore” fra i radiologi e gli ingegneri dell’Università di Bologna. Proprio in questo modo è stato possibile standardizzare le procedure di esecuzione degli esami perfusionali e implementare, attraverso l’interfacciamento tra il personale medico ed i ricercatori dell’Università, i protocolli per la gestione ed il trattamento dei dati acquisiti.

Aspetti Ultrasonografici in Corso di Sviluppo Placentare e Adattamento Cardiovascolare nella Cavalla Gravida: ecografia bidimensionale, ecocardiografia ed ecocontrastografia (CEUS)

In corso di gravidanza normale avvengono modificazioni emodinamiche centrali e periferiche volte a garantire le crescenti richieste nutritive dell'unità feto-placentare. L’ecografia con mezzo di contrasto (CEUS-Contrast Enhanced Ultrasonography) a base di microbolle offre una nuova opportunità di monitorare e quantificare la perfusione utero-placentare in condizioni normali e patologiche. L’ecocardiografia è stata ampiamente usata in medicina umana per valutare l’adattamento morfo-funzionale cardiaco materno durante la gravidanza. Gli scopi di questo lavoro prospettico sono stati di applicare, per la prima volta nella specie equina, un mezzo di contrasto di II generazione (Sonovue®), al fine quantificare la perfusione utero-placentare in corso di gravidanza normale, valutandone gli effetti sul benessere materno-fetale e di descrivere le modificazioni nei parametri ecocardiografici morfometrici e funzionali cardiaci, in particolare relativi alla funzione del ventricolo sinistro nel corso di una gravidanza fisiologica. Due fattrici sane di razza Trottatore sono state monitorate ecograficamente in maniera seriale durante l’intero corso della gravidanza, tramite esame bidimensionale, ecocontrastografia dell'unità utero-placentare, flussimetria Doppler delle arterie uterine, ecocardiografia materna in modalità bidimensionale, M-mode, Doppler e Tissue Doppler Imaging. I neonati sono stati clinicamente monitorati e gli invogli fetali esaminati. Il pattern di microperfusione utero-placentare è valutabile quali-quantitativamente tramite la CEUS e dimostra un’aumento del flusso a livello di microvascolarizzazione uterina con l'avanzare della gravidanza; non è stata rilevata la presenza di microbolle a livello di strutture fetali nè effetti dannosi sul benessere materno-fetale. In questo studio sono state osservate delle modificazioni cardiache materne in corso di gravidanza fisiologica, relative all'aumento della FC, del CO ed in particolare all'aumento delle dimensioni dell'atrio sinistro ed a modificazioni nelle onde di velocità di flusso e tissutali di riempimento del ventricolo sinistro.

Numerical modelling of fluid-structure interaction with application in hemodynamics

The thesis deals with numerical algorithms for fluid-structure interaction problems with application in blood flow modelling. It starts with a short introduction on the mathematical description of incompressible viscous flow with non-Newtonian viscosity and a moving linear viscoelastic structure. The mathematical model consists of the generalized Navier-Stokes equation used for the description of fluid flow and the generalized string model for structure movement. The arbitrary Lagrangian-Eulerian approach is used in order to take into account moving computational domain. A part of the thesis is devoted to the discussion on the non-Newtonian behaviour of shear-thinning fluids, which is in our case blood, and derivation of two non-Newtonian models frequently used in the blood flow modelling. Further we give a brief overview on recent fluid-structure interaction schemes with discussion about the difficulties arising in numerical modelling of blood flow. Our main contribution lies in numerical and experimental study of a new loosely-coupled partitioned scheme called the kinematic splitting fluid-structure interaction algorithm. We present stability analysis for a coupled problem of non-Newtonian shear-dependent fluids in moving domains with viscoelastic boundaries. Here, we assume both, the nonlinearity in convective as well is diffusive term. We analyse the convergence of proposed numerical scheme for a simplified fluid model of the Oseen type. Moreover, we present series of experiments including numerical error analysis, comparison of hemodynamic parameters for the Newtonian and non-Newtonian fluids and comparison of several physiologically relevant computational geometries in terms of wall displacement and wall shear stress. Numerical analysis and extensive experimental study for several standard geometries confirm reliability and accuracy of the proposed kinematic splitting scheme in order to approximate fluid-structure interaction problems.

Vergleich der retinalen und dermalen mikrovaskulären endothelialen Funktion bei Gesunden, Insulinresistenten und Typ-2-Diabetikern

Die vorliegende klinische Studie hatte zum Ziel, die mikrovaskuläre Endothelfunktion retinaler und dermaler Gefäße von Insulinresistenten und Typ 2- Diabetikern ohne Zeichen einer diabetischen Retinopathie mit einer gesunden insulinsensitiven nicht-diabetischen Kontrollgruppe hinsichtlich früher morphologischer und funktioneller Veränderungen zu vergleichen.rnrnMethode:rnEs wurden 54 Patienten ohne Nachweis einer diabetischen Retinopathie eingeschlossen und in 3 Gruppen entsprechend ihren metabolischen Ergebnissen eingeteilt: 1.) Gruppe K (Kontrollgruppe) setzte sich aus gesunden, nicht-diabetischen, insulin-sensitiven rn(HOMA ≤ 2) Probanden mit einem BMI ≤ 28 kg/m2 zusammen; 2.) Gruppe IR bestand aus den nicht-diabetischen, insulin-resistenten (HOMAs > 2), übergewichtigen Patienten mit einem BMI > 28 kg/m2 und 3.) Gruppe DM war definiert als Patienten mit einem manifesten Typ 2-Diabetes mellitus.rnrnDie mikrovaskuläre Funktion der Retina wurde mittels eines Laserdoppler-Verfahrens (Heidelberg Retina Flowmeter) untersucht und hierbei der retinale Blutfluss und das Verhältnis der Gefäßwand zum Lumen (WLR, wall-to-lumen-ratio) basal und nach Flickerlicht-Stimulation (10 Hz, Photo Stimulator 750) gemessen. Letzterer gilt als Marker für vaskuläre Schädigung. rnZusätzlich wurde die dermale Mikrozirkulation (Blutfluss, O2-Sättigung) als weiterer Faktor der mikrovaskulären Endothelfunktion in den 3 Studiengruppen untersucht und miteinander verglichen.rnErgebnisse:rnEs zeigte sich kein signifikanter Unterschied des retinalen Blutflusses zwischen den 3 Gruppen weder basal noch nach Flickerlicht-Stimulation. Es zeigte sich keine Korrelation zwischen der mikrovaskulären Funktion der Haut und der Retina. rnDie arterielle WLR zeigte nur geringe Unterscheide zwischen den 3 Gruppen.rnrnMit zunehmendem Grad der Insulinresistenz wurde jedoch eine Reduktion des basalen als auch des flickerlicht-stimulierten retinalen Blutflusses deutlich, dabei zeigte sich unerwarteter Weise eine Abnahme der WLR.rnrnDer (prä-ischämische) muskuläre Blutfluss war in der IR-Gruppe signifikant geringer als in der K-Gruppe. Auch war die postischämische dermale O2-Sättigung in der DM und IR-Gruppe signifikant niedriger im Vergleich zur K-Gruppe. Jedoch war die postischämische hyperämische dermale Reaktion in der IR und DM-Gruppe nur geringgradig weniger als in der K-Gruppe. rnrnSchlussfolgerung:rnEine Korrelation zwischen der Entwicklung der Insulinresistenz und retinaler sowie dermaler mikrovaskulärer endothelialer Funktion wurde bei der Studie deutlich. Mithilfe des neuen Verfahrens der Laser Scanner Flowmeter zur Messung der retinalen Endothelfunktion lassen sich sehr frühe morphologische Veränderungen des mikrovaskulären Blutflusses erfassen. rnDie fehlende Korrelation zwischen retinaler und dermaler mikrovaskulärer Funktion als auch die geringen Unterschiede der WLR sollte Gegenstand weiterer Studien seinrn