Imaging and modeling of pore scale processes in porous media using X-ray computed tomography and lattice Boltzmann solver

In der Erdöl– und Gasindustrie sind bildgebende Verfahren und Simulationen auf der Porenskala im Begriff Routineanwendungen zu werden. Ihr weiteres Potential lässt sich im Umweltbereich anwenden, wie z.B. für den Transport und Verbleib von Schadstoffen im Untergrund, die Speicherung von Kohlendioxid und dem natürlichen Abbau von Schadstoffen in Böden. Mit der Röntgen-Computertomografie (XCT) steht ein zerstörungsfreies 3D bildgebendes Verfahren zur Verfügung, das auch häufig für die Untersuchung der internen Struktur geologischer Proben herangezogen wird. Das erste Ziel dieser Dissertation war die Implementierung einer Bildverarbeitungstechnik, die die Strahlenaufhärtung der Röntgen-Computertomografie beseitigt und den Segmentierungsprozess dessen Daten vereinfacht. Das zweite Ziel dieser Arbeit untersuchte die kombinierten Effekte von Porenraumcharakteristika, Porentortuosität, sowie die Strömungssimulation und Transportmodellierung in Porenräumen mit der Gitter-Boltzmann-Methode. In einer zylindrischen geologischen Probe war die Position jeder Phase auf Grundlage der Beobachtung durch das Vorhandensein der Strahlenaufhärtung in den rekonstruierten Bildern, das eine radiale Funktion vom Probenrand zum Zentrum darstellt, extrahierbar und die unterschiedlichen Phasen ließen sich automatisch segmentieren. Weiterhin wurden Strahlungsaufhärtungeffekte von beliebig geformten Objekten durch einen Oberflächenanpassungsalgorithmus korrigiert. Die Methode der „least square support vector machine” (LSSVM) ist durch einen modularen Aufbau charakterisiert und ist sehr gut für die Erkennung und Klassifizierung von Mustern geeignet. Aus diesem Grund wurde die Methode der LSSVM als pixelbasierte Klassifikationsmethode implementiert. Dieser Algorithmus ist in der Lage komplexe geologische Proben korrekt zu klassifizieren, benötigt für den Fall aber längere Rechenzeiten, so dass mehrdimensionale Trainingsdatensätze verwendet werden müssen. Die Dynamik von den unmischbaren Phasen Luft und Wasser wird durch eine Kombination von Porenmorphologie und Gitter Boltzmann Methode für Drainage und Imbibition Prozessen in 3D Datensätzen von Böden, die durch synchrotron-basierte XCT gewonnen wurden, untersucht. Obwohl die Porenmorphologie eine einfache Methode ist Kugeln in den verfügbaren Porenraum einzupassen, kann sie dennoch die komplexe kapillare Hysterese als eine Funktion der Wassersättigung erklären. Eine Hysterese ist für den Kapillardruck und die hydraulische Leitfähigkeit beobachtet worden, welche durch die hauptsächlich verbundenen Porennetzwerke und der verfügbaren Porenraumgrößenverteilung verursacht sind. Die hydraulische Konduktivität ist eine Funktion des Wassersättigungslevels und wird mit einer makroskopischen Berechnung empirischer Modelle verglichen. Die Daten stimmen vor allem für hohe Wassersättigungen gut überein. Um die Gegenwart von Krankheitserregern im Grundwasser und Abwässern vorhersagen zu können, wurde in einem Bodenaggregat der Einfluss von Korngröße, Porengeometrie und Fluidflussgeschwindigkeit z.B. mit dem Mikroorganismus Escherichia coli studiert. Die asymmetrischen und langschweifigen Durchbruchskurven, besonders bei höheren Wassersättigungen, wurden durch dispersiven Transport aufgrund des verbundenen Porennetzwerks und durch die Heterogenität des Strömungsfeldes verursacht. Es wurde beobachtet, dass die biokolloidale Verweilzeit eine Funktion des Druckgradienten als auch der Kolloidgröße ist. Unsere Modellierungsergebnisse stimmen sehr gut mit den bereits veröffentlichten Daten überein.

Ricerca, classificazione e archiviazione dei dati di elettrofisiologia cellulare cardiaca umana per l'identificazione di modelli matematici

L’obiettivo di questo lavoro di tesi è stato quindi quello di ricercare e archiviare tutti i dati sperimentali di correnti ioniche umane ventricolari presenti in letteratura fino ad oggi, per costruire uno strumento di facile utilizzo per chiunque abbia la necessità di sviluppare o validare modelli matematici di potenziale d’azione. Partendo da una fase iniziale di ricerca vera e propria degli articoli in letteratura, utilizzando il motore di ricerca PubMed come strumento principale, sono stati estratti e archiviati tutti i dati di interesse, divisi per tipo di corrente, memorizzando le informazioni principali in un foglio di lavoro e salvando i dati sia come immagini che come vettori, per consentirne in futuro una rapida consultazione e un facile utilizzo.

Classificazione delle galassie: caratteristiche morfologiche, fotometriche e cinematiche

Questo elaborato scritto tratta la classificazione e la modellizzazione della fisica delle galassie. La prima parte illustra la classificazione morfologica delle galassie. Descrivo la classificazione di Hubble perché è la prima in ordine cronologico, la più semplice e la più importante. In seguito espongo le modifiche introdotte da Sandage e De Vaucouleurs. In particolare cerco di giustificare l’importanza delle classificazioni morfologiche come punto di partenza per una comprensione e una modellizzazione della fisica delle galassie. Nella seconda parte mi concentro sulle caratteristiche fotometriche e cinematiche interne relative a due particolari tipi di galassie: le Ellittiche e le Spirali. Approfondisco due argomenti in particolare: i profili di brillanza delle Ellittiche per quanto riguarda le caratteristiche fotometriche delle Ellittiche e le curve di rotazione delle Spirali, per le caratteristiche cinematiche delle Spirali. Questi due argomenti sono stati scelti perché il primo fornisce un modello analitico per descrivere la luminosità delle Ellittiche e il secondo permette di rappresentare, sempre tramite una modellizzazione, la cinematica delle Spirali, in particolar modo del loro disco.

Applicazione del metodo QDanet_PRO alla classificazione di dati omici

Il presente lavoro di tesi si pone nell'ambito dell'analisi dati attraverso un metodo (QDanet_PRO), elaborato dal Prof. Remondini in collaborazine coi Dott. Levi e Malagoli, basato sull'analisi discriminate a coppie e sulla Teoria dei Network, che ha come obiettivo la classificazione di dati contenuti in dataset dove il numero di campioni è molto ridotto rispetto al numero di variabili. Attraverso questo studio si vogliono identificare delle signature, ovvero un'insieme ridotto di variabili che siano in grado di classificare correttamente i campioni in base al comportamento delle variabili stesse. L'elaborazione dei diversi dataset avviene attraverso diverse fasi; si comincia con una un'analisi discriminante a coppie per identificare le performance di ogni coppia di variabili per poi passare alla ricerca delle coppie più performanti attraverso un processo che combina la Teoria dei Network con la Cross Validation. Una volta ottenuta la signature si conclude l'elaborazione con una validazione per avere un'analisi quantitativa del successo o meno del metodo.

Analisi spettrale in banda X e proprietà multi-frequenza della galassia ULIRG IRAS 13120-5453

L’interruzione dell’attività di formazione stellare nelle galassie attraverso l’azione di venti guidati dalla potenza dell’AGN e/o dall’attività di starburst è una fase prevista da tutti i modelli di coevoluzione tra i buchi neri super massicci e la galassia ospite. La recente scoperta di venti molecolari massivi in alcune ULIRG locali ha fornito una prova forte a favore dell’esistenza del meccanismo di feedback negativo nelle galassie. Ad oggi non è chiaro il ruolo che l'attività di AGN e di starburst hanno nella produzione dei venti: tuttavia, la maggior parte delle ULIRG in cui sono stati osservati i venti molecolari presentano elevati tassi e velocità di outflow, suggerendo che l'AGN giochi un ruolo principale. IRAS 13120-5453 rappresenta un caso particolare: la sua attività e dominata dalla formazione stellare, ma sembra anche ospitare un AGN oscurato al suo interno. Tuttavia, non presenta elevati tassi di outflow e velocità dei venti: non è quindi chiaro quale dei due fenomeni vi sia all'origine. In questo lavoro di tesi si è operata un'analisi multibanda di questa sorgente con lo scopo di studiarne le principali proprietà fisiche. Lo studio in banda X attraverso l'analisi dei dati provenienti dai satelliti XMM-Newton, Chandra e NuSTAR ha permesso di conoscere parametri importanti come il grado di oscuramento della sorgente e la potenza dell'AGN. Con l'analisi in banda IR, è stato possibile conoscere il contributo dell'AGN e della starburst alla luminosità IR e i principali parametri fisici di questa galassia. L’obiettivo di tale lavoro è quello di capire il ruolo svolto da questi due principali fenomeni e quale possa essere la connessione con i venti molecolari.

Two dimensional Lattice Boltzmann simulation of fluid flow through an idealized micro-structure of disordered media

This technical report discusses the application of Lattice Boltzmann Method (LBM) in the fluid flow simulation through porous filter-wall of disordered media. The diesel particulate filter (DPF) is an example of disordered media. DPF is developed as a cutting edge technology to reduce harmful particulate matter in the engine exhaust. Porous filter-wall of DPF traps these soot particles in the after-treatment of the exhaust gas. To examine the phenomena inside the DPF, researchers are looking forward to use the Lattice Boltzmann Method as a promising alternative simulation tool. The lattice Boltzmann method is comparatively a newer numerical scheme and can be used to simulate fluid flow for single-component single-phase, single-component multi-phase. It is also an excellent method for modelling flow through disordered media. The current work focuses on a single-phase fluid flow simulation inside the porous micro-structure using LBM. Firstly, the theory concerning the development of LBM is discussed. LBM evolution is always related to Lattice gas Cellular Automata (LGCA), but it is also shown that this method is a special discretized form of the continuous Boltzmann equation. Since all the simulations are conducted in two-dimensions, the equations developed are in reference with D2Q9 (two-dimensional 9-velocity) model. The artificially created porous micro-structure is used in this study. The flow simulations are conducted by considering air and CO2 gas as fluids. The numerical model used in this study is explained with a flowchart and the coding steps. The numerical code is constructed in MATLAB. Different types of boundary conditions and their importance is discussed separately. Also the equations specific to boundary conditions are derived. The pressure and velocity contours over the porous domain are studied and recorded. The results are compared with the published work. The permeability values obtained in this study can be fitted to the relation proposed by Nabovati [8], and the results are in excellent agreement within porosity range of 0.4 to 0.8.

LATTICE BOLTZMANN METHOD AND CELLULAR AUTOMATA SIMULATION OF PARTICLE MOTION AND DEPOSITION IN 2-D CASE

This technical report discusses the application of the Lattice Boltzmann Method (LBM) and Cellular Automata (CA) simulation in fluid flow and particle deposition. The current work focuses on incompressible flow simulation passing cylinders, in which we incorporate the LBM D2Q9 and CA techniques to simulate the fluid flow and particle loading respectively. For the LBM part, the theories of boundary conditions are studied and verified using the Poiseuille flow test. For the CA part, several models regarding simulation of particles are explained. And a new Digital Differential Analyzer (DDA) algorithm is introduced to simulate particle motion in the Boolean model. The numerical results are compared with a previous probability velocity model by Masselot [Masselot 2000], which shows a satisfactory result.

Carnot Cycles: Traditional and Stefan Boltzmann, and the First Planck Argument

The Stefan Boltzmann equation is obtained using a non-traditional Carnot Engine. In addition, the original Planck argument for radiation density is given.

Measured Stark widths of several spectral lines of Pb III

The Stark full widths at half of the maximal line intensity (FWHM, ω) have been measured for 25 spectrallines of PbIII (15 measured for the first time) arising from the 5d106s8s, 5d106s7p, 5d106s5f and 5d106s5g electronic configurations, in a lead plasma produced by ablation with a Nd:YAG laser. The optical emission spectroscopy from a laser-induced plasma generated by a 10 640 Å radiation, with an irradiance of 2 × 1010 W cm− 2 on a lead target (99.99% purity) in an atmosphere of argon was analysed in the wavelength interval between 2000 and 7000 Å. The broadening parameters were obtained with the target placed in argon atmosphere at 6 Torr and 400 ns after each laser light pulse, which provides appropriate measurement conditions. A Boltzmann plot was used to obtain the plasma temperature (21,400 K) and published values of the Starkwidths in Pb I, Pb II and PbIII to obtain the electron number density (7 × 1016 cm− 3); with these values, the plasma composition was determined by means of the Saha equation. Local Thermodynamic Equilibrium (LTE) conditions and plasma homogeneity has been checked. Special attention was dedicated to the possible self-absorption of the different transitions. Comparison of the new results with recent available data is also presented.

Conductance interaction identification by means of Boltzmann distribution and mutual information analysis in conductance-based neuron models

Conductance interaction identification by means of Boltzmann distribution and mutual information analysis in conductance-based neuron models.

Modulation of Sarco/Endoplasmic Reticulum Ca2+-ATPase 2 (SERCA2) function by acetylation following the treatment with histone deacetylase inhibitor suberoylanilide hydroxamic acid (SAHA)

Background: Acetylation and deacetylation at specific lysine (K) residues is mediated by histone acetylases (HATs) and deacetylases (HDACs), respectively. HATs and HDACs act on both histone and non-histone proteins, regulating various processes, including cardiac impulse propagation. Aim of the present work was to establish whether the function of the Ca2+ ATPase SERCA2, one of the major players in Ca2+ reuptake during excitation-contraction coupling in cardiac myocytes (CMs), could be modulated by direct K acetylation. Materials and methods: HL-1 atrial mouse cells (donated by Prof. Claycomb), zebrafish and Streptozotocin-induced diabetic rat CMs were treated with the pan-inhibitor of class I and II HDACs suberanilohydroxamic acid (SAHA) for 1.5 hour. Evaluation of SERCA2 acetylation was analyzed by co-immunoprecipitation. SERCA2 activity was measured on microsomes by pyruvate/NADH coupled reaction assay. SERCA2 mutants were obtained after cloning wild-type and mutated sequences into the pCDNA3 vector and transfected into HEK cells. Ca2+ transients in CMs (loading with Fluo3-AM, field stimulation, 0.5 Hz) and in transfected HEK cells (loading with FLUO-4, caffeine pulse) were recorded. Results: Co-Immunoprecipitation experiments performed on HL-1 cells demonstrated a significant increase in the acetylation of SERCA2 after SAHA-treatment (2.5 µM, n=3). This was associated with an increase in SERCA2 activity in microsomes obtained from HL-1 cells, after SAHA exposure (n=5). Accordingly, SAHA-treatment significantly shortened the Ca2+ reuptake time of adult zebrafish CMs. Further, SAHA 2.5 nM restored to control values the recovery time of Ca2+ transients decay in diabetic rat CMs. HDAC inhibition also improved contraction parameters, such as fraction of shortening, and increased pump activity in microsomes isolated from diabetic CMs (n=4). Notably, the K464, identified by bioinformatic tools as the most probable acetylation site on human SERCA2a, was mutated into Glutamine (Q) or Arginine (R) mimicking acetylation and deacetylation respectively. Measurements of Ca2+ transients in HEK cells revealed that the substitution of K464 with R significantly delayed the transient recovery time, thus indicating that deacetylation has a negative impact on SERCA2 function. Conclusions: Our results indicate that SERCA2 function can be improved by pro-acetylation interventions and that this mechanism of regulation is conserved among species. Therefore, the present work provides the basis to open the search for novel pharmacological tools able to specifically improve SERCA2 activity in diseases where its expression and/or function is impaired, such as diabetic cardiomyopathy.