Multilevel Domain Decomposition Algorithms for Monolithic Fluid-Structure Interaction Problems with Application to Haemodynamics

Finite element techniques for solving the problem of fluid-structure interaction of an elastic solid material in a laminar incompressible viscous flow are described. The mathematical problem consists of the Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian formulation coupled with a non-linear structure model, considering the problem as one continuum. The coupling between the structure and the fluid is enforced inside a monolithic framework which computes simultaneously for the fluid and the structure unknowns within a unique solver. We used the well-known Crouzeix-Raviart finite element pair for discretization in space and the method of lines for discretization in time. A stability result using the Backward-Euler time-stepping scheme for both fluid and solid part and the finite element method for the space discretization has been proved. The resulting linear system has been solved by multilevel domain decomposition techniques. Our strategy is to solve several local subproblems over subdomain patches using the Schur-complement or GMRES smoother within a multigrid iterative solver. For validation and evaluation of the accuracy of the proposed methodology, we present corresponding results for a set of two FSI benchmark configurations which describe the self-induced elastic deformation of a beam attached to a cylinder in a laminar channel flow, allowing stationary as well as periodically oscillating deformations, and for a benchmark proposed by COMSOL multiphysics where a narrow vertical structure attached to the bottom wall of a channel bends under the force due to both viscous drag and pressure. Then, as an example of fluid-structure interaction in biomedical problems, we considered the academic numerical test which consists in simulating the pressure wave propagation through a straight compliant vessel. All the tests show the applicability and the numerical efficiency of our approach to both two-dimensional and three-dimensional problems.

Biomedical and industrial applications of atmospheric pressure non-equilibrium plasmas

This dissertation will be focused on the characterization of an atmospheric pressure plasma jet source with an application oriented diagnostic approach and the description of processes supported by this plasma source. The plasma source investigated is a single electrode plasma jet. Schlieren images, optical emission spectra, temperature and heat flux profiles are analyzed to deeply investigate the fluid dynamic, the chemical composition and the thermal output of the plasma generated with a nanosecond-pulsed high voltage generator. The maximum temperature measured is about 45 °C and values close to the room temperature are reached 10 mm down the source outlet, ensuring the possibility to use the plasma jet for the treatment of thermosensitive materials, such as, for example, biological substrate or polymers. Electrospinning of polymeric solution allows the production of nanofibrous non-woven mats and the plasma pre-treatment of the solutions leads to the realization of defect free nanofibers. The use of the plasma jet allows the electrospinnability of a non-spinnable poly(L-lactic acid) (PLLA) solution, suitable for the production of biological scaffold for the wound dressing.

Numerical analysis of the relation between interactions and structure in a molecular fluid

Coarse graining is a popular technique used in physics to speed up the computer simulation of molecular fluids. An essential part of this technique is a method that solves the inverse problem of determining the interaction potential or its parameters from the given structural data. Due to discrepancies between model and reality, the potential is not unique, such that stability of such method and its convergence to a meaningful solution are issues.rnrnIn this work, we investigate empirically whether coarse graining can be improved by applying the theory of inverse problems from applied mathematics. In particular, we use the singular value analysis to reveal the weak interaction parameters, that have a negligible influence on the structure of the fluid and which cause non-uniqueness of the solution. Further, we apply a regularizing Levenberg-Marquardt method, which is stable against the mentioned discrepancies. Then, we compare it to the existing physical methods - the Iterative Boltzmann Inversion and the Inverse Monte Carlo method, which are fast and well adapted to the problem, but sometimes have convergence problems.rnrnFrom analysis of the Iterative Boltzmann Inversion, we elaborate a meaningful approximation of the structure and use it to derive a modification of the Levenberg-Marquardt method. We engage the latter for reconstruction of the interaction parameters from experimental data for liquid argon and nitrogen. We show that the modified method is stable, convergent and fast. Further, the singular value analysis of the structure and its approximation allows to determine the crucial interaction parameters, that is, to simplify the modeling of interactions. Therefore, our results build a rigorous bridge between the inverse problem from physics and the powerful solution tools from mathematics. rn

Design and fluid dynamic analysis of a custom drip chamber in a medical disposable device

Negli ultimi anni, parallelamente allo sviluppo di calcolatori elettronici sempre più performanti, la fluidodinamica computazionale è diventata uno strumento di notevole utilità nell’analisi dei flussi e nello sviluppo di dispositivi medici. Quando impiegate nello studio di flussi di fluidi fisiologici, come il sangue, il vantaggio principale delle analisi CFD è che permettono di caratterizzare il comportamento fluidodinamico senza dover eseguire test in-vivo/in-vitro, consentendo quindi notevoli vantaggi in termini di tempo, denaro e rischio derivante da applicazioni mediche. Inoltre, simulazioni CFD offrono una precisa e dettagliata descrizione di ogni parametro di interesse permettendo, già in fase di progettazione, di prevedere quali modifiche al layout garantiranno maggiori vantaggi in termini di funzionalità. Il presente lavoro di tesi si è posto l’obiettivo di valutare, tramite simulazioni CFD, le performances fluidodinamiche del comparto sangue “camera venosa” di un dispositivo medico monouso Bellco impiegato nella realizzazione di trattamenti di emodialisi. Dopo una panoramica del contesto, è presentata una breve descrizione della disfunzione renale e dei trattamenti sostitutivi. Notevole impegno è stato in seguito rivolto allo studio della letteratura scientifica in modo da definire un modello reologico per il fluido non-Newtoniano preso in considerazione e determinarne i parametri caratteristici. Il terzo capitolo presenta lo stato dell’arte delle apparecchiature Bellco, rivolgendosi con particolare attenzione al componente “cassette” del dispositivo monouso. L’analisi fluidodinamica del compartimento “camera venosa” della cassette, che sarà presa in considerazione nei capitoli quinto e sesto, si inserisce nell’ambito della riprogettazione del dispositivo attualmente in commercio: il quarto capitolo si incentra sul suo nuovo design, ponendo specifico interesse sul layout della camera venosa di nuova generazione. Per lo studio dei flussi che si sviluppano internamente ad essa ci si è avvalsi del modulo CFD del software COMSOL multiphysics® (versione 5.0); la definizione del modello implementato e della tipologia di studio effettuato sono presi in considerazione nel quinto capitolo. Le problematiche di maggior impatto nella realizzazione di un trattamento di emodialisi sono l’emolisi e la coagulazione del sangue. Nell'evenienza che si verifichino massivamente occorre infatti interrompere il trattamento con notevoli disagi per il paziente, per questo devono essere evitate. Nel sesto capitolo i risultati ottenuti sono stati esaminati rivolgendo particolare attenzione alla verifica dell’assenza di fenomeni che possano portare alle problematiche suddette.

Atlas-Based Segmentation of Brain Tumor Images Using a Markov Random Field-Based Tumor Growth Model and Non-Rigid Registration

We propose a new and clinically oriented approach to perform atlas-based segmentation of brain tumor images. A mesh-free method is used to model tumor-induced soft tissue deformations in a healthy brain atlas image with subsequent registration of the modified atlas to a pathologic patient image. The atlas is seeded with a tumor position prior and tumor growth simulating the tumor mass effect is performed with the aim of improving the registration accuracy in case of patients with space-occupying lesions. We perform tests on 2D axial slices of five different patient data sets and show that the approach gives good results for the segmentation of white matter, grey matter, cerebrospinal fluid and the tumor.

Quantitative PCR to diagnose Pneumocystis pneumonia in immunocompromised non-HIV patients

The utility of quantitative Pneumocystis jirovecii PCR in clinical routine for diagnosing Pneumocystis pneumonia (PCP) in immunocompromised non-HIV patients is unknown. We analysed bronchoalveolar lavage fluid with real-time quantitative P. jirovecii PCR in 71 cases with definitive PCP defined by positive immunofluorescence (IF) tests and in 171 randomly selected patients with acute lung disease. In those patients, possible PCP cases were identified by using a novel standardised PCP probability algorithm and chart review. PCR performance was compared with IF testing, clinical judgment and the PCP probability algorithm. Quantitative P. jirovecii PCR values >1,450 pathogens·mL(-1) had a positive predictive value of 98.0% (95% CI 89.6-100.0%) for diagnosing definitive PCP. PCR values of between 1 and 1,450 pathogens·mL(-1) were associated with both colonisation and infection; thus, a cut-off between the two conditions could not be identified and diagnosis of PCP in this setting relied on IF and clinical assessment. Clinical PCP could be ruled out in 99.3% of 153 patients with negative PCR results. Quantitative PCR is useful for diagnosing PCP and is complementary to IF. PCR values of >1,450 pathogens·mL(-1) allow reliable diagnosis, whereas negative PCR results virtually exclude PCP. Intermediate values require additional clinical assessment and IF testing. On the basis of our data and for economic and logistical limitations, we propose a clinical algorithm in which IF remains the preferred first test in most cases, followed by PCR in those patients with a negative IF and strong clinical suspicion for PCP.

Rapid diagnosis of experimental meningitis by bacterial heat production in cerebrospinal fluid

BACKGROUND: Calorimetry is a nonspecific technique which allows direct measurement of heat generated by biological processes in the living cell. We evaluated the potential of calorimetry for rapid detection of bacterial growth in cerebrospinal fluid (CSF) in a rat model of bacterial meningitis. METHODS: Infant rats were infected on postnatal day 11 by direct intracisternal injection with either Streptococcus pneumoniae, Neisseria meningitidis or Listeria monocytogenes. Control animals were injected with sterile saline or heat-inactivated S. pneumoniae. CSF was obtained at 18 hours after infection for quantitative cultures and heat flow measurement. For calorimetry, 10 microl and 1 microl CSF were inoculated in calorimetry ampoules containing 3 ml trypticase soy broth (TSB). RESULTS: The mean bacterial titer (+/- SD) in CSF was 1.5 +/- 0.6 x 108 for S. pneumoniae, 1.3 +/- 0.3 x 106 for N. meningitidis and 3.5 +/- 2.2 x 104 for L. monocytogenes. Calorimetric detection time was defined as the time until heat flow signal exceeded 10 microW. Heat signal was detected in 10-microl CSF samples from all infected animals with a mean (+/- SD) detection time of 1.5 +/- 0.2 hours for S. pneumoniae, 3.9 +/- 0.7 hours for N. meningitidis and 9.1 +/- 0.5 hours for L. monocytogenes. CSF samples from non-infected animals generated no increasing heat flow (<10 microW). The total heat was the highest in S. pneumoniae ranging from 6.7 to 7.5 Joules, followed by L. monocytogenes (5.6 to 6.1 Joules) and N. meningitidis (3.5 to 4.4 Joules). The lowest detectable bacterial titer by calorimetry was 2 cfu for S. pneumoniae, 4 cfu for N. meningitidis and 7 cfu for L. monocytogenes. CONCLUSION: By means of calorimetry, detection times of <4 hours for S. pneumoniae and N. meningitidis and <10 hours for Listeria monocytogenes using as little as 10 microl CSF were achieved. Calorimetry is a new diagnostic method allowing rapid and accurate diagnosis of bacterial meningitis from a small volume of CSF.

Non-intrusive pressure measurement in microchannels

A non-intrusive interferometric measurement technique has been successfully developed to measure fluid compressibility in both gas and liquid phases via refractive index (RI) changes. The technique, consisting of an unfocused laser beam impinging a glass channel, can be used to separate and quantify cell deflection, fluid flow rates, and pressure variations in microchannels. Currently in fields such as microfluidics, pressure and flow rate measurement devices are orders of magnitude larger than the channel cross-sections making direct pressure and fluid flow rate measurements impossible. Due to the non-intrusive nature of this technique, such measurements are now possible, opening the door for a myriad of new scientific research and experimentation. This technique, adapted from the concept of Micro Interferometric Backscatter Detection (MIBD), boasts the ability to provide comparable sensitivities in a variety of channel types and provides quantification capability not previously demonstrated in backscatter detection techniques. Measurement sensitivity depends heavily on experimental parameters such as beam impingement angle, fluid volume, photodetector sensitivity, and a channel’s dimensional tolerances. The current apparatus readily quantifies fluid RI changes of 10-5 refractive index units (RIU) corresponding to pressures of approximately 14 psi and 1 psi in water and air, respectively. MIBD reports detection capability as low as 10-9 RIU and the newly adapted technique has the potential to meet and exceed this limit providing quantification in the place of detection. Specific device sensitivities are discussed and suggestions are provided on how the technique may be refined to provide optimal quantification capabilities based on experimental conditions.

Non-invasive method to predict viscosity and size using total internal reflection fluorescence microscopy (TIRFM) system

This study develops an automated analysis tool by combining total internal reflection fluorescence microscopy (TIRFM), an evanescent wave microscopic imaging technique to capture time-sequential images and the corresponding image processing Matlab code to identify movements of single individual particles. The developed code will enable us to examine two dimensional hindered tangential Brownian motion of nanoparticles with a sub-pixel resolution (nanoscale). The measured mean square displacements of nanoparticles are compared with theoretical predictions to estimate particle diameters and fluid viscosity using a nonlinear regression technique. These estimated values will be confirmed by the diameters and viscosities given by manufacturers to validate this analysis tool. Nano-particles used in these experiments are yellow-green polystyrene fluorescent nanospheres (200 nm, 500 nm and 1000 nm in diameter (nominal); 505 nm excitation and 515 nm emission wavelengths). Solutions used in this experiment are de-ionized (DI) water, 10% d-glucose and 10% glycerol. Mean square displacements obtained near the surface shows significant deviation from theoretical predictions which are attributed to DLVO forces in the region but it conforms to theoretical predictions after ~125 nm onwards. The proposed automation analysis tool will be powerfully employed in the bio-application fields needed for examination of single protein (DNA and/or vesicle) tracking, drug delivery, and cyto-toxicity unlike the traditional measurement techniques that require fixing the cells. Furthermore, this tool can be also usefully applied for the microfluidic areas of non-invasive thermometry, particle tracking velocimetry (PTV), and non-invasive viscometry.

DEVELOPMENT OF AN IONIC LIQUID FERROFLUID ELECTROSPRAY SOURCE AND MODE SHAPE STUDIES OF A FERROFLUID IN A NON-UNIFORM MAGNETIC FIELD

An electrospray source has been developed using a novel new fluid that is both magnetic and conductive. Unlike conventional electrospray sources that required microfabricated structures to support the fluid to be electrosprayed, this new electrospray fluid utilizes the Rosensweig instability to create the structures in the magnetic fluid when an external magnetic field was applied. Application of an external electric field caused these magnetic fluid structures to spray. These fluid based structures were found to spray at a lower onset voltage than was predicted for electrospray sources with solid structures of similar geometry. These fluid based structures were also found to be resilient to damage, unlike the solid structures found in traditional electrospray sources. Further, experimental studies of magnetic fluids in non-uniform magnetic fields were conducted. The modes of Rosensweig instabilities have been studied in-depth when created by uniform magnetic fields, but little to no studies have been performed on Rosensweig instabilities formed due to non-uniform magnetic fields. The measured spacing of the cone-like structures of ferrofluid, in a non-uniform magnetic field, were found to agree with a proposed theoretical model.

Effect of fluid resuscitation on mortality and organ function in experimental sepsis models

Introduction Several recent studies have shown that a positive fluid balance in critical illness is associated with worse outcome. We tested the effects of moderate vs. high-volume resuscitation strategies on mortality, systemic and regional blood flows, mitochondrial respiration, and organ function in two experimental sepsis models. Methods 48 pigs were randomized to continuous endotoxin infusion, fecal peritonitis, and a control group (n = 16 each), and each group further to two different basal rates of volume supply for 24 hours [moderate-volume (10 ml/kg/h, Ringer's lactate, n = 8); high-volume (15 + 5 ml/kg/h, Ringer's lactate and hydroxyethyl starch (HES), n = 8)], both supplemented by additional volume boli, as guided by urinary output, filling pressures, and responses in stroke volume. Systemic and regional hemodynamics were measured and tissue specimens taken for mitochondrial function assessment and histological analysis. Results Mortality in high-volume groups was 87% (peritonitis), 75% (endotoxemia), and 13% (controls). In moderate-volume groups mortality was 50% (peritonitis), 13% (endotoxemia) and 0% (controls). Both septic groups became hyperdynamic. While neither sepsis nor volume resuscitation strategy was associated with altered hepatic or muscle mitochondrial complex I- and II-dependent respiration, non-survivors had lower hepatic complex II-dependent respiratory control ratios (2.6 +/- 0.7, vs. 3.3 +/- 0.9 in survivors; P = 0.01). Histology revealed moderate damage in all organs, colloid plaques in lung tissue of high-volume groups, and severe kidney damage in endotoxin high-volume animals. Conclusions High-volume resuscitation including HES in experimental peritonitis and endotoxemia increased mortality despite better initial hemodynamic stability. This suggests that the strategy of early fluid management influences outcome in sepsis. The high mortality was not associated with reduced mitochondrial complex I- or II-dependent muscle and hepatic respiration.

Comparison of two non-bronchoscopic methods for evaluating inflammation in patients with acute hypoxaemic respiratory failure

INTRODUCTION: The simple bedside method for sampling undiluted distal pulmonary edema fluid through a normal suction catheter (s-Cath) has been experimentally and clinically validated. However, there are no data comparing non-bronchoscopic bronchoalveolar lavage (mini-BAL) and s-Cath for assessing lung inflammation in acute hypoxaemic respiratory failure. We designed a prospective study in two groups of patients, those with acute lung injury (ALI)/acute respiratory distress syndrome (ARDS) and those with acute cardiogenic lung edema (ACLE), designed to investigate the clinical feasibility of these techniques and to evaluate inflammation in both groups using undiluted sampling obtained by s-Cath. To test the interchangeability of the two methods in the same patient for studying the inflammation response, we further compared mini-BAL and s-Cath for agreement of protein concentration and percentage of polymorphonuclear cells (PMNs). METHODS: Mini-BAL and s-Cath sampling was assessed in 30 mechanically ventilated patients, 21 with ALI/ARDS and 9 with ACLE. To analyse agreement between the two sampling techniques, we considered only simultaneously collected mini-BAL and s-Cath paired samples. The protein concentration and polymorphonuclear cell (PMN) count comparisons were performed using undiluted sampling. Bland-Altman plots were used for assessing the mean bias and the limits of agreement between the two sampling techniques; comparison between groups was performed by using the non-parametric Mann-Whitney-U test; continuous variables were compared by using the Student t-test, Wilcoxon signed rank test, analysis of variance or Student-Newman-Keuls test; and categorical variables were compared by using chi-square analysis or Fisher exact test. RESULTS: Using protein content and PMN percentage as parameters, we identified substantial variations between the two sampling techniques. When the protein concentration in the lung was high, the s-Cath was a more sensitive method; by contrast, as inflammation increased, both methods provided similar estimates of neutrophil percentages in the lung. The patients with ACLE showed an increased PMN count, suggesting that hydrostatic lung edema can be associated with a concomitant inflammatory process. CONCLUSIONS: There are significant differences between the s-Cath and mini-BAL sampling techniques, indicating that these procedures cannot be used interchangeably for studying the lung inflammatory response in patients with acute hypoxaemic lung injury.

Neurally adjusted ventilatory assist decreases ventilator-induced lung injury and non-pulmonary organ dysfunction in rabbits with acute lung injury

OBJECTIVE: To determine if neurally adjusted ventilatory assist (NAVA) that delivers pressure in proportion to diaphragm electrical activity is as protective to acutely injured lungs (ALI) and non-pulmonary organs as volume controlled (VC), low tidal volume (Vt), high positive end-expiratory pressure (PEEP) ventilation. DESIGN: Prospective, randomized, laboratory animal study. SUBJECTS: Twenty-seven male New Zealand white rabbits. INTERVENTIONS: Anesthetized rabbits with hydrochloric acid-induced ALI were randomized (n = 9 per group) to 5.5 h NAVA (non-paralyzed), VC (paralyzed; Vt 6-ml/kg), or VC (paralyzed; Vt 15-ml/kg). PEEP was adjusted to hemodynamic goals in NAVA and VC6-ml/kg, and was 1 cmH2O in VC15-ml/kg. MEASUREMENTS AND MAIN RESULTS: PaO2/FiO2; lung wet-to-dry ratio; lung histology; interleukin-8 (IL-8) concentrations in broncho-alveolar-lavage (BAL) fluid, plasma, and non-pulmonary organs; plasminogen activator inhibitor type-1 and tissue factor in BAL fluid and plasma; non-pulmonary organ apoptosis rate; creatinine clearance; echocardiography. PEEP was similar in NAVA and VC6-ml/kg. During NAVA, Vt was lower (3.1 +/- 0.9 ml/kg), whereas PaO2/ FiO2, respiratory rate, and PaCO2 were higher compared to VC6-ml/kg (p<0.05 for all). Variables assessing ventilator-induced lung injury (VILI), IL-8 levels, non-pulmonary organ apoptosis rate, and kidney as well as cardiac performance were similar in NAVA compared to VC6-ml/kg. VILI and non-pulmonary organ dysfunction was attenuated in both groups compared to VC15-ml/kg. CONCLUSIONS: In anesthetized rabbits with early experimental ALI, NAVA is as effective as VC6-ml/kg in preventing VILI, in attenuating excessive systemic and remote organ inflammation, and in preserving cardiac and kidney function.

Free abdominal fluid without obvious solid organ injury upon CT imaging: an actual problem or simply over-diagnosing?

Whereas a non-operative approach for hemodynamically stable patients with free intraabdominal fluid in the presence of solid organ injury is generally accepted, the presence of free fluid in the abdomen without evidence of solid organ injury not only presents a challenge for the treating emergency physician but also for the surgeon in charge. Despite recent advances in imaging modalities, with multi-detector computed tomography (CT) (with or without contrast agent) usually the imaging method of choice, diagnosis and interpretation of the results remains difficult. While some studies conclude that CT is highly accurate and relatively specific at diagnosing mesenteric and hollow viscus injury, others studies deem CT to be unreliable. These differences may in part be due to the experience and the interpretation of the radiologist and/or the treating physician or surgeon.A search of the literature has made it apparent that there is no straightforward answer to the question what to do with patients with free intraabdominal fluid on CT scanning but without signs of solid organ injury. In hemodynamically unstable patients, free intraabdominal fluid in the absence of solid organ injury usually mandates immediate surgical intervention. For patients with blunt abdominal trauma and more than just a trace of free intraabdominal fluid or for patients with signs of peritonitis, the threshold for a surgical exploration - preferably by a laparoscopic approach - should be low. Based on the available information, we aim to provide the reader with an overview of the current literature with specific emphasis on diagnostic and therapeutic approaches to this problem and suggest a possible algorithm, which might help with the adequate treatment of such patients.

Flow disturbances in stent-related coronary evaginations: a computational fluid-dynamic simulation study

Aims: Angiographic ectasias and aneurysms in stented segments have been associated with late stent thrombosis. Using optical coherence tomography (OCT), some stented segments show coronary evaginations reminiscent of ectasias. The purpose of this study was to explore, using computational fluid-dynamic (CFD) simulations, whether OCT-detected coronary evaginations can induce local changes in blood flow. Methods and results: OCT-detected evaginations are defined as outward bulges in the luminal vessel contour between struts, with the depth of the bulge exceeding the actual strut thickness. Evaginations can be characterised cross ectionally by depth and along the stented segment by total length. Assuming an ellipsoid shape, we modelled 3-D evaginations with different sizes by varying the depth from 0.2-1.0 mm, and the length from 1-9 mm. For the flow simulation we used average flow velocity data from non-diseased coronary arteries. The change in flow with varying evagination sizes was assessed using a particle tracing test where the particle transit time within the segment with evagination was compared with that of a control vessel. The presence of the evagination caused a delayed particle transit time which increased with the evagination size. The change in flow consisted locally of recirculation within the evagination, as well as flow deceleration due to a larger lumen - seen as a deflection of flow towards the evagination. Conclusions: CFD simulation of 3-D evaginations and blood flow suggests that evaginations affect flow locally, with a flow disturbance that increases with increasing evagination size.