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.

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.

Investigation of dangling bonds in al2o3 passivated si surface using electron paramagnetic resonance (epr)

Nel presente lavoro di tesi magistrale sono stati depositati e caratterizzati sottili film di ossido di alluminio, Al2O3, (di spessore compreso tra 3-30 nm) su un substrato di FZ-Si drogato p. La deposizione è avvenuta mediante plasma ALD (Atomic Layer Depostion). La tecnica spettroscopica EPR (Electron Paramagnetic Resonance) è stata utilizzata per studiare l’interfaccia Si/Al2O3 con lo scopo di scoprire l’origine della formazione di densità di carica negativa Qf all’interfaccia: tale carica negativa induce una passivazione per effetto di campo ed è quindi la ragione per cui il dielettrico Al2O3 risulta essere un ottimo materiale passivante. Si è deciso di variare alcuni parametri, come lo spessore dello strato di Al2O3, lo spessore dello strato intermedio di ossido di silicio, depositato mediante ossidazione termica (dry thermal oxidation), e la superficie del substrato di silicio. Sono stati realizzati cinque differenti gruppi di campioni: per ciascuno di essi sono state impiegate varie tecniche di caratterizzazione, come la QSSPC (Quasi Steady State Photoconuctance) e la tecnica di spettroscopia ottica SE (spettroscopic ellipsometry). Per ogni gruppo sono stati riportati gli spettri EPR ottenuti ed i rispettivi fit, da cui è stato possibile risalire ai fattori giromagnetici di spin g, riportati in tabelle con le loro possibili attribuzioni. E’ stato dimostrato che la presenza di uno strato di ossido di silicio tra il substrato di silicio e lo strato di ossido di alluminio risulta essere fondamentale per la formazione di densità di carica negativa all’interfaccia: aumentando lo spessore dello strato di SiOx (nel range 1-30 nm) si assiste ad una diminuzione di carica negativa Qf. Analizzando gli spettri EPR, è stato possibile concludere che all’interfaccia Si/Al2O3 sono presenti difetti caratteristici dell’interfaccia Si/SiOx. Le nostre osservazioni, dunque, sono coerenti con la formazione di uno strato di ossido di silicio tra Si e Al2O3.

Aeroelastic stability of structures: flutter analysis using Computational Fluid Dynamics

Thanks to the increasing slenderness and lightness allowed by new construction techniques and materials, the effects of wind on structures became in the last decades a research field of great importance in Civil Engineering. Thanks to the advances in computers power, the numerical simulation of wind tunnel tests has became a valid complementary activity and an attractive alternative for the future. Due to its flexibility, during the last years, the computational approach gained importance with respect to the traditional experimental investigation. However, still today, the computational approach to fluid-structure interaction problems is not as widely adopted as it could be expected. The main reason for this lies in the difficulties encountered in the numerical simulation of the turbulent, unsteady flow conditions generally encountered around bluff bodies. This thesis aims at providing a guide to the numerical simulation of bridge deck aerodynamic and aeroelastic behaviour describing in detail the simulation strategies and setting guidelines useful for the interpretation of the results.

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.

Structural analysis of the major lightharvesting complex II by electron paramagnetic resonance (EPR) : comparison between monomers and trimers

Der Haupt-Lichtsammenkomplex II (LHCII) höherer Pflanzen ist das häufigsternMembranprotein der Welt und in die chloroplastidäre Thylakoidmembran integriert. DerrnLHCII kann als Modellsystem genutzt werden, um die Funktionsweise vonrnMembranproteinen besser zu verstehen, da 96 % seiner Struktur kristallografisch aufgelöstrnist und er in rekombinanter Form in vitro rückgefaltet werden kann. Hierbei entsteht einrnvoll funktionaler Protein-Pigment.Komplex, der nahezu identisch mit der in vivo Varianternist.rnElektronenparamagnetischen Resonanz (EPR) Spektroskopie ist eine hoch sensitive undrnideal geeignete Methode, um die Strukturdynamik von Proteinen zu untersuchen. Hierzurnist eine ortsspezifische Markierung mit Spinsonden notwendig, die kovalent an Cysteinernbinden. Möglich wird dies, indem sorgfältig ausgewählte Aminosäuren gegen Cysteinerngetauscht werden, ohne dass die Funktionsweise des LHCII beeinträchtigt wird.rnIm Rahmen dieser Arbeit wurden die Stabilität des verwendeten Spinmarkers und diernProbenqualität verbessert, indem alle Schritte der Probenpräparation untersucht wurden.rnMithilfe dieser Erkenntnisse konnte sowohl die Gefahr einer Proteinaggregation als auchrnein Verlust des EPR Signals deutlich vermindert werden. In Kombination mit derrngleichzeitigen Etablierung des Q-Band EPR können nun deutlich geringer konzentrierternProben zuverlässig vermessen werden. Darüber hinaus wurde eine reproduzierbarernMethode entwickelt, um heterogene Trimere herzustellen. Diese bestehen aus einemrndoppelt markierten Monomer und zwei unmarkierten Monomeren und erlauben es, diernkristallografisch unvollständig aufgelöste N-terminale Domäne im monomeren undrntrimeren Assemblierungsgrad zu untersuchen. Die Ergebnisse konnten einerseits diernVermutung bestätigen, dass diese Domäne im Vergleich zum starren Proteinkern sehrrnflexibel ist und andererseits, dass sie in Monomeren noch mobiler ist als in Trimeren.rnZudem wurde die lumenale Schleifenregion bei unterschiedlichen pH Werten undrnvariierender Pigmentzusammensetzung untersucht, da dieser Bereich sehr kontroversrndiskutiert wird. Die Messergebnisse offenbarten, dass diese Region starre und flexiblerernSektionen aufweist. Während der pH Wert keinen Einfluss auf die Konformation hatte,rnzeigte sich, dass die Abwesenheit von Neoxanthin zu einer Änderung der Konformationrnführt. Weiterführende Analysen der strukturellen Dynamik des LHCII in einerrnLipidmembran konnten hingegen nicht durchgeführt werden, da dies eine gerichteternInsertion des rückgefalteten Proteins in Liposomen erfordert, was trotz intensiverrnVersuche nicht zum Erfolg führte.

Bio-inspired polymer-supported nitrido molybdenum(VI) complexes for ammonia synthesis − reactivity towards protons and electrons : EPR investigations on paramagnetic molybdenum(V) and copper(II) complexes

Die biologische Stickstofffixierung durch Molybdän-haltige Nitrogenasen sowie die Erforschung des zugrundeliegenden komplexen Mechanismus (N2-Aktivierung an Metall-Zentren, 6-fache Protonierung und Reduktion, N–N Bindungsspaltung unter Bildung von Ammoniak) ist von erheblichem Interesse. Insbesondere Molybdän-Komplexe wurden bereits erfolgreich als Modellverbindungen für die Untersuchung elementarer Einzelschritte der N2-Aktivierung eingesetzt. Durch die Verwendung von Triamidoamin-Liganden ist es Schrock et al. sogar gelungen mehrere Katalysezyklen zu durchlaufen und einen Mechanismus zu formulieren. Trotz der sterisch anspruchsvollen Substituenten in den Schrock-Komplexen ist die Umsatzrate dieses homogenen Katalysators, aufgrund Komplex-Deaktivierung infolge intermolekularer Reaktionen wie Dimerisierung und Disproportionierung, limitiert. In der vorliegenden Arbeit wurden einige dieser Herausforderungen angegangen und die aktiven Spezies auf einer Festphase immobilisiert, um intermolekulare Reaktionen durch räumliche Isolierung der Komplexe zu unterdrücken.rnEin Polymer-verankertes Analogon des Schrock Nitrido-Molybdän(VI)-Komplexes wurde auf einem neuen Reaktionsweg synthetisiert. Dieser beinhaltet nur einen einzigen Reaktionsschritt, um die funktionelle Gruppe „MoN“ einzuführen. Protonierung des immobilisierten Nitrido-Molybdän(VI)-Komplexes LMoVIN (L = Polymer-verankerter Triamidoamin-Ligand) mit 2,6-Lutidinium liefert den entsprechenden Imido-Molybdän(VI)-Komplex. Durch anschließende Ein-Elektronen-Reduktion mit Cobaltocen wird der Polymer-angebundene Imido-Molybdän(V)-Komplex erhalten, bewiesen durch EPR-Spektroskopie (g1,2,3 = 1.989, 1.929, 1.902). Durch die Immobilisierung und die effektive räumliche Separation der Reaktionszentren auf der Festphase werden bimolekulare Nebenreaktionen, die oft in homogenen Systemen auftreten, unterdrückt. Dies ermöglicht zum ersten Mal die Darstellung des Imido-Molybdän(V)-Intermediates des Schrock-Zyklus.rnEPR-Spektren des als Spin-Label eingeführten immobilisierten Nitrato-Kupfer(II)-Komplexes wurden unter verschiedenen Bedingungen (Lösungsmittel, Temperatur) aufgenommen, wobei sich eine starke Abhängigkeit zwischen der Zugänglichkeit und Reaktivität der immobilisierten Reaktionszentren und der Art des Lösungsmittels zeigte. Somit wurde die Reaktivität von LMoVIN gegenüber Protonen und Elektronen, welches zur Bildung von NH3 führt, unter Verwendung verschiedener Lösungsmittel untersucht und optimiert. Innerhalb des kugelförmigen Polymers verläuft die Protonierung und Reduktion von LMoVIN stufenweise. Aktive Zentren, die sich an der „äußeren Schale“ des Polymers befinden, sind gut zugänglich und reagieren schnell nach H+/e− Zugabe. Aktive Zentren im „Inneren des Polymers“ hingegen sind schlechter zugänglich und zeigen langsame diffusions-kontrollierte Reaktionen, wobei drei H+/e− Schritte gefolgt von einer Ligandenaustausch-Reaktion erforderlich sind, um NH3 freizusetzen: LMoVIN  LMoVNH  LMoIVNH2  LMoIIINH3 und anschließender Ligandenaustausch führt zur Freisetzung von NH3.rnIn einem weiteren Projekt wurde der Bis(ddpd)-Kupfer(II)-Komplex EPR-spektroskopisch in Hinblick auf Jahn−Teller-Verzerrung und -Dynamik untersucht. Dabei wurden die EPR-Spektren bei variabler Temperatur (70−293 K) aufgenommen. Im Festkörperspektrum bei T < 100 K erscheint der Kupfer(II)-Komplex als gestreckter Oktaeder, wohingegen das EPR-Spektrum bei höheren Temperaturen g-Werte aufzeigt, die einer pseudo-gestauchten oktaedrischen Kupfer(II)-Spezies zuzuordnen sind. Diese Tatsache wird einem intramolekularen dynamischen Jahn−Teller Phänomen zugeschrieben, welcher bei 100 K eingefroren wird.

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.

Acute fluid ingestion in the treatment of orthostatic intolerance - important implications for daily practice

Rapid water ingestion improves orthostatic intolerance (OI) in multiple system atrophy (MSA) and postural tachycardia syndrome (PoTS). We compared haemodynamic changes after water and clear soup intake, the latter being a common treatment strategy for OI in daily practice.

Formation of caries-like lesions in vitro on the root surfaces of human teeth in solutions simulating plaque fluid

Lesion formation on root surfaces of human posterior teeth was studied in acetate/lactate buffers with a background electrolyte composition based on plaque fluid analyses. Lesion depth after 28 days at 37 degrees C was measured in relation to: the presence or absence of cementum; the concentration of undissociated buffer; the presence or absence of magnesium ions at plaque fluid concentration. Each factor was evaluated at several values of -log(ion activity product for hydroxyapatite): pI(HA). Solutions were formulated to minimize variation in pH, which varied by < or =0.03 for a given comparison (individual pI(HA)) and by 0.42-0.82 over the range of pI(HA) within experiments. Lesions on surfaces from which cementum had been ground were significantly deeper than on intact surfaces, but this is considered to be due to subsurface mechanical damage and not to a solubility difference. Neither the concentration of undissociated buffer nor the presence of magnesium ions significantly affected lesion depth. Lesion depth was strongly influenced by the correlated variations in pI(HA) and pH. At pI(HA) 54 and 55, only extremely shallow lesions formed. From pI(HA) 56, lesion depth increased with increasing pI(HA). The results confirm that the solubility of the mineral of root tissues is higher than that of hydroxyapatite, but indicate that it is probably lower than suggested by Hoppenbrouwers et al. [Arch Oral Biol 1987;32:319-322]. For calcium concentrations of 3-12 mM, the critical pH for root tissue mineral was calculated as 5.22-5.66 assuming solubility equivalent to pI(HA) 54 and 5.08-5.51 assuming pI(HA) 55.

Is detection of bacterial DNA in ascitic fluid of clinical relevance?

In patients with cirrhosis, bacterial DNA has been found in ascites reflecting bacterial translocation. However, the clinical relevance of this finding is ill-defined especially compared with the standard diagnostics for detection of spontaneous bacterial peritonitis (SBP). Furthermore, other DNA tests have not been sufficiently evaluated.

Lung fluid movements in hypoxia

Pulmonary edema is a problem of major clinical importance resulting from a persistent imbalance between forces that drive water into the airspace of the lung and the biological mechanisms for its removal. Here, we will first review the fundamental mechanisms implicated in the regulation of lung fluid homeostasis, namely, the Starling forces and the respiratory transepithelial sodium transport. Second, we will discuss the contribution of hypoxia to the perturbation of this fine balance and the role of such perturbations in the development of high-altitude pulmonary edema, a disease characterized by a very high morbidity and mortality. Finally, we will review possible interventions aimed to maintain/restore lung fluid homeostasis and their importance for the prevention/treatment of pulmonary edema.

Petrous bone pneumatization is a risk factor for cerebrospinal fluid fistula following vestibular schwannoma surgery

For the prevention of postoperative CSF fistula a better understanding of origins and risk factors is necessary.

Attenuation of cerebrospinal fluid inflammation by the nonbacteriolytic antibiotic daptomycin versus that by ceftriaxone in experimental pneumococcal meningitis

Antibiotic-induced bacteriolysis exacerbates inflammation and brain damage in bacterial meningitis. Here the quality and temporal kinetics of cerebrospinal fluid (CSF) inflammation were assessed in an infant rat pneumococcal meningitis model for the nonbacteriolytic antibiotic daptomycin versus ceftriaxone. Daptomycin led to lower CSF concentrations of interleukin 1beta (IL-1beta), IL-10, IL-18, monocyte chemoattractant protein 1 (MCP-1), and macrophage inflammatory protein 1 alpha (MIP-1alpha) (P < 0.05). In experimental pneumococcal meningitis, daptomycin treatment resulted in more rapid bacterial killing, lower CSF inflammation, and less brain damage than ceftriaxone treatment.