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.

Risk factors contributing to symptomatic miniplate removal: a retrospective study of 153 bilateral sagittal split osteotomy patients

Removal of miniplates is a controversial topic in oral and maxillofacial surgery. Originally, miniplates were designed to be removed on completion of bone healing. The introduction of low profile titanium miniplates has led to the routine removal of miniplates becoming comparatively rare in many parts of the world. Few studies have investigated the reasons for non-routine removal of miniplates and the factors that affect osteosynthesis after osteotomy in large numbers of patients. The aim of the present study was to investigate complications related to osteosynthesis after bilateral sagittal split osteotomy (BSSO) in a large number (n=153) of patients. In addition to the rates of removal, emphasis was placed on investigating the reasons and risk factors associated with symptomatic miniplate removal. The rate of plate removal per patient was 18.6%, the corresponding rate per plate being 18.2%. Reasons for plate removal included plate-related complications in 16 patients and subjective discomfort in 13 patients. Half of the plates were removed during the first postoperative year. Smoking was the only significant predictor for plate removal. Patients undergoing orthognathic surgery should be screened with regard to smoking and encouraged and assisted to cease smoking, at least perioperatively.

Interaction of pregnancy and autoimmune rheumatic disease

During pregnancy, the fetus represents a natural allograft that is not normally rejected. While the maternal immune system retains the ability to respond to foreign antigens, tolerance mechanisms are up-regulated to protect the fetus from immunologic attacks by the mother. The profound immunologic adaptations during and after pregnancy do influence maternal autoimmune rheumatic diseases in several ways. One is triggering the onset of a rheumatic disease in the post partum period, the other influencing disease activity of established rheumatic disease. The review will discuss the mechanisms of increased susceptibility of rheumatoid arthritis (RA) in the first year post partum with a specific emphasis on the role of fetal cells or antigens persisting in the maternal circulation (so called microchimerism). Furthermore, the different influences of pregnancy on established rheumatic diseases will be highlighted. A marked beneficial effect of pregnancy is observed on RA whereas several other rheumatic diseases as ankylosing spondylitis (AS) and systemic lupus erythematosus (SLE) show either no particular effect or an aggravation of symptoms during pregnancy. Differences emerging in regard to modulation of disease symptoms during pregnancy seem related to response to hormones, the type of cytokine profile and immune response prevailing as well as further downstream interactions of molecular pathways that are important in disease pathogenesis.

Generalized eMC implementation for Monte Carlo dose calculation of electron beams from different machine types

The electron Monte Carlo (eMC) dose calculation algorithm available in the Eclipse treatment planning system (Varian Medical Systems) is based on the macro MC method and uses a beam model applicable to Varian linear accelerators. This leads to limitations in accuracy if eMC is applied to non-Varian machines. In this work eMC is generalized to also allow accurate dose calculations for electron beams from Elekta and Siemens accelerators. First, changes made in the previous study to use eMC for low electron beam energies of Varian accelerators are applied. Then, a generalized beam model is developed using a main electron source and a main photon source representing electrons and photons from the scattering foil, respectively, an edge source of electrons, a transmission source of photons and a line source of electrons and photons representing the particles from the scrapers or inserts and head scatter radiation. Regarding the macro MC dose calculation algorithm, the transport code of the secondary particles is improved. The macro MC dose calculations are validated with corresponding dose calculations using EGSnrc in homogeneous and inhomogeneous phantoms. The validation of the generalized eMC is carried out by comparing calculated and measured dose distributions in water for Varian, Elekta and Siemens machines for a variety of beam energies, applicator sizes and SSDs. The comparisons are performed in units of cGy per MU. Overall, a general agreement between calculated and measured dose distributions for all machine types and all combinations of parameters investigated is found to be within 2% or 2 mm. The results of the dose comparisons suggest that the generalized eMC is now suitable to calculate dose distributions for Varian, Elekta and Siemens linear accelerators with sufficient accuracy in the range of the investigated combinations of beam energies, applicator sizes and SSDs.

Fastener-Based Computational Models of Cold-Formed Steel Shear Walls

Cold-formed steel (CFS) combined with wood sheathing, such as oriented strand board (OSB), forms shear walls that can provide lateral resistance to seismic forces. The ability to accurately predict building deformations in damaged states under seismic excitations is a must for modern performance-based seismic design. However, few static or dynamic tests have been conducted on the non-linear behavior of CFS shear walls. Thus, the purpose of this research work is to provide and demonstrate a fastener-based computational model of CFS wall models that incorporates essential nonlinearities that may eventually lead to improvement of the current seismic design requirements. The approach is based on the understanding that complex interaction of the fasteners with the sheathing is an important factor in the non-linear behavior of the shear wall. The computational model consists of beam-column elements for the CFS framing and a rigid diaphragm for the sheathing. The framing and sheathing are connected with non-linear zero-length fastener elements to capture the OSB sheathing damage surrounding the fastener area. Employing computational programs such as OpenSees and MATLAB, 4 ft. x 9 ft., 8 ft. x 9 ft. and 12 ft. x 9 ft. shear wall models are created, and monotonic lateral forces are applied to the computer models. The output data are then compared and analyzed with the available results of physical testing. The results indicate that the OpenSees model can accurately capture the initial stiffness, strength and non-linear behavior of the shear walls.

Association between symptoms and quality of life in patients with schizophrenia: a pooled analysis of changes over time

Quality of life is an important outcome in the treatment of patients with schizophrenia. It has been suggested that patients' quality of life ratings (referred to as subjective quality of life, SQOL) might be too heavily influenced by symptomatology to be a valid independent outcome criterion. There has been only limited evidence on the association of symptom change and changes in SQOL over time. This study aimed to examine the association between changes in symptoms and in SQOL among patients with schizophrenia. A pooled data set was obtained from eight longitudinal studies that had used the Brief Psychiatric Rating Scale (BPRS) for measuring psychiatric symptoms and either the Lancashire Quality of Life Profile or the Manchester Short Assessment of Quality of Life for assessing SQOL. The sample comprised 886 patients with schizophrenia. After controlling for heterogeneity of findings across studies using linear mixed models, a reduction in psychiatric symptoms was associated with improvements in SQOL scores. In univariate analyses, changes in all BPRS subscales were associated with changes in SQOL scores. In a multivariate model, only associations between changes in the BPRS depression/anxiety and hostility subscales and changes in SQOL remained significant, with 5% and 0.5% of the variance in SQOL changes being attributable to changes in depression/anxiety and hostility respectively. All BPRS subscales together explained 8.5% of variance. The findings indicate that SQOL changes are influenced by symptom change, in particular in depression/anxiety. The level of influence is limited and may not compromise using SQOL as an independent outcome measure.

Factors influencing subjective quality of life in patients with schizophrenia and other mental disorders: a pooled analysis

Subjective quality of life (SQOL) is an important outcome in the treatment of patients with schizophrenia. However, there is only limited evidence on factors influencing SQOL, and little is known about whether the same factors influence SQOL in patients with schizophrenia and other mental disorders. This study aimed to identify the factors associated with SQOL and test whether these factors are equally important in schizophrenia and other disorders. For this we used a pooled data set obtained from 16 studies that had used either the Lancashire Quality of Life Profile or the Manchester Short Assessment of Quality of Life for assessing SQOL. The sample comprised 3936 patients with schizophrenia, mood disorders, and neurotic disorders. After controlling for confounding factors, within-subject clustering, and heterogeneity of findings across studies in linear mixed models, patients with schizophrenia had more favourable SQOL scores than those with mood and neurotic disorders. In all diagnostic groups, older patients, those in employment, and those with lower symptom scores had higher SQOL scores. Whilst the strength of the association between age and SQOL did not differ across diagnostic groups, symptom levels were more strongly associated with SQOL in neurotic than in mood disorders and schizophrenia. The association of employment and SQOL was stronger in mood and neurotic disorders than in schizophrenia. The findings may inform the use and interpretation of SQOL data for patients with schizophrenia.

Monte Carlo simulation of a dynamic MLC based on a multiple source model

Detailed knowledge of the characteristics of the radiation field shaped by a multileaf collimator (MLC) is essential in intensity modulated radiotherapy (IMRT). A previously developed multiple source model (MSM) for a 6 MV beam was extended to a 15 MV beam and supplemented with an accurate model of an 80-leaf dynamic MLC. Using the supplemented MSM and the MC code GEANT, lateral dose distributions were calculated in a water phantom and a portal water phantom. A field which is normally used for the validation of the step and shoot technique and a field from a realistic IMRT treatment plan delivered with dynamic MLC are investigated. To assess possible spectral changes caused by the modulation of beam intensity by an MLC, the energy spectra in five portal planes were calculated for moving slits of different widths. The extension of the MSM to 15 MV was validated by analysing energy fluences, depth doses and dose profiles. In addition, the MC-calculated primary energy spectrum was verified with an energy spectrum which was reconstructed from transmission measurements. MC-calculated dose profiles using the MSM for the step and shoot case and for the dynamic MLC case are in very good agreement with the measured data from film dosimetry. The investigation of a 13 cm wide field shows an increase in mean photon energy of up to 16% for the 0.25 cm slit compared to the open beam for 6 MV and of up to 6% for 15 MV, respectively. In conclusion, the MSM supplemented with the dynamic MLC has proven to be a powerful tool for investigational and benchmarking purposes or even for dose calculations in IMRT.

Simple beam models for Monte Carlo photon beam dose calculations in radiotherapy

Monte Carlo (code GEANT) produced 6 and 15 MV phase space (PS) data were used to define several simple photon beam models. For creating the PS data the energy of starting electrons hitting the target was tuned to get correct depth dose data compared to measurements. The modeling process used the full PS information within the geometrical boundaries of the beam including all scattered radiation of the accelerator head. Scattered radiation outside the boundaries was neglected. Photons and electrons were assumed to be radiated from point sources. Four different models were investigated which involved different ways to determine the energies and locations of beam particles in the output plane. Depth dose curves, profiles, and relative output factors were calculated with these models for six field sizes from 5x5 to 40x40cm2 and compared to measurements. Model 1 uses a photon energy spectrum independent of location in the PS plane and a constant photon fluence in this plane. Model 2 takes into account the spatial particle fluence distribution in the PS plane. A constant fluence is used again in model 3, but the photon energy spectrum depends upon the off axis position. Model 4, finally uses the spatial particle fluence distribution and off axis dependent photon energy spectra in the PS plane. Depth dose curves and profiles for field sizes up to 10x10cm2 were not model sensitive. Good agreement between measured and calculated depth dose curves and profiles for all field sizes was reached for model 4. However, increasing deviations were found for increasing field sizes for models 1-3. Large deviations resulted for the profiles of models 2 and 3. This is due to the fact that these models overestimate and underestimate the energy fluence at large off axis distances. Relative output factors consistent with measurements resulted only for model 4.

Paleoenvironmental reconstructions from cave sediments of the Moravian Karst, Czech Republic / by Pavel Sroubek.

This thesis presents a paleoclimatic/paleoenvironmental study conducted on clastic cave sediments of the Moravian Karst, Czech Republic. The study is based on environmental magnetic techniques, yet a wide range of other scientific methods was used to obtain a clearer picture of the Quaternary climate. My thesis also presents an overview of the significance of cave deposits for paleoclimatic reconstructions, explains basic environmental magnetic techniques and offers background information on the study area – a famous karst region in Central Europe with a rich history. In Kulna Cave magnetic susceptibility variations and in particular variations in pedogenic susceptibility yield a detailed record of the palaeoenvironmental conditions during the Last Glacial Stage. The Kulna long-term climatic trends agree with the deep-sea SPECMAP record, while the short-term oscillations correlate with rapid changes in the North Atlantic sea surface temperatures. Kulna Cave sediments reflect the intensity of pedogenesis controlled by short-term warmer events and precipitation over the mid-continent and provide a link between continental European climate and sea surface temperatures in the North Atlantic during the Last Glacial Stage. Given the number of independent climate proxies determined from the entrance facies of the cave and their high resolution, Kulna is an extremely important site for studying Late Pleistocene climate. In the interior of Spiralka Cave, a five meter high section of fine grained sediments deposited during floods yields information on the climatic and environmental conditions of the last millenium. In the upper 1.5 meters of this profile, mineral magnetic and other non-magnetic data indicate that susceptibility variations are controlled by the concentration of magnetite and its magnetic grain size. Comparison of our susceptibility record to the instrumental record of winter temperature anomalies shows a remarkable correlation. This correlation is explained by coupling of the flooding events, cultivation of land and pedogenetic processes in the cave catchment area. A combination of mineral magnetic and geochemical proxies yields a detail picture of the rapidly evolving climate of the near past and tracks both natural and human induced environmental changes taking place in the broader region.

Utilizing Gaussian-Schell model beams to mitigate atmospheric turbulence in free space optical communications / by Kyle R. Drexler.

Turbulence affects traditional free space optical communication by causing speckle to appear in the received beam profile. This occurs due to changes in the refractive index of the atmosphere that are caused by fluctuations in temperature and pressure, resulting in an inhomogeneous medium. The Gaussian-Schell model of partial coherence has been suggested as a means of mitigating these atmospheric inhomogeneities on the transmission side. This dissertation analyzed the Gaussian-Schell model of partial coherence by verifying the Gaussian-Schell model in the far-field, investigated the number of independent phase control screens necessary to approach the ideal Gaussian-Schell model, and showed experimentally that the Gaussian-Schell model of partial coherence is achievable in the far-field using a liquid crystal spatial light modulator. A method for optimizing the statistical properties of the Gaussian-Schell model was developed to maximize the coherence of the field while ensuring that it does not exhibit the same statistics as a fully coherent source. Finally a technique to estimate the minimum spatial resolution necessary in a spatial light modulator was developed to effectively propagate the Gaussian-Schell model through a range of atmospheric turbulence strengths. This work showed that regardless of turbulence strength or receiver aperture, transmitting the Gaussian-Schell model of partial coherence instead of a fully coherent source will yield a reduction in the intensity fluctuations of the received field. By measuring the variance of the intensity fluctuations and the received mean, it is shown through the scintillation index that using the Gaussian-Schell model of partial coherence is a simple and straight forward method to mitigate atmospheric turbulence instead of traditional adaptive optics in free space optical communications.

Age modeling of young non-varved lake sediments: methods and limits. Examples from two lakes in Central Chile

High-resolution and highly precise age models for recent lake sediments (last 100–150 years) are essential for quantitative paleoclimate research. These are particularly important for sedimentological and geochemical proxies, where transfer functions cannot be established and calibration must be based upon the relation of sedimentary records to instrumental data. High-precision dating for the calibration period is most critical as it determines directly the quality of the calibration statistics. Here, as an example, we compare radionuclide age models obtained on two high-elevation glacial lakes in the Central Chilean Andes (Laguna Negra: 33°38′S/70°08′W, 2,680 m a.s.l. and Laguna El Ocho: 34°02′S/70°19′W, 3,250 m a.s.l.). We show the different numerical models that produce accurate age-depth chronologies based on 210Pb profiles, and we explain how to obtain reduced age-error bars at the bottom part of the profiles, i.e., typically around the end of the 19th century. In order to constrain the age models, we propose a method with five steps: (i) sampling at irregularly-spaced intervals for 226Ra, 210Pb and 137Cs depending on the stratigraphy and microfacies, (ii) a systematic comparison of numerical models for the calculation of 210Pb-based age models: constant flux constant sedimentation (CFCS), constant initial concentration (CIC), constant rate of supply (CRS) and sediment isotope tomography (SIT), (iii) numerical constraining of the CRS and SIT models with the 137Cs chronomarker of AD 1964 and, (iv) step-wise cross-validation with independent diagnostic environmental stratigraphic markers of known age (e.g., volcanic ash layer, historical flood and earthquakes). In both examples, we also use airborne pollutants such as spheroidal carbonaceous particles (reflecting the history of fossil fuel emissions), excess atmospheric Cu deposition (reflecting the production history of a large local Cu mine), and turbidites related to historical earthquakes. Our results show that the SIT model constrained with the 137Cs AD 1964 peak performs best over the entire chronological profile (last 100–150 years) and yields the smallest standard deviations for the sediment ages. Such precision is critical for the calibration statistics, and ultimately, for the quality of the quantitative paleoclimate reconstruction. The systematic comparison of CRS and SIT models also helps to validate the robustness of the chronologies in different sections of the profile. Although surprisingly poorly known and under-explored in paleolimnological research, the SIT model has a great potential in paleoclimatological reconstructions based on lake sediments

Characterization of microfluidic systems with Doppler optical coherence tomography

Doppler Optical Coherence Tomography (DOCT) is a biomedical imaging technique that allows simultaneous structural imaging and flow monitoring inside biological tissues and materials with spatial resolution in the micrometer scale. It has recently been applied to the characterization of microfluidic systems. Structural and flow imaging of novel microfluidics platforms for cytotoxicologic applications were obtained with a real-time, Near Infrared Spectral Domain DOCT system. Characteristics such as flow homogeneity in the chamber, which is one of the most important parameters for cell culture, are investigated. OCT and DOCT images were used to monitor flow inside a specific platform that is based on microchannel division for a better flow homogeneity. In particular, the evolution of flow profile at the transition between the microchannel structure and the chamber is studied.

Persistence of Artemisia steppe in the Tangra Yumco Basin, west-central Tibet, China: despite or in consequence of Holocene lake-level changes?

The closed Tangra Yumco Basin underwent the strongest Quaternary lake-level changes so far recorded on the Tibetan Plateau. It was hitherto unknown what effect this had on local Holocene vegetation development. A 3.6-m sediment core from a recessional lake terrace at 4,700 m a.s.l., 160 m above the present lake level of Tangra Yumco, was studied to reconstruct Holocene flooding phases (sedimentology and ostracod analyses), vegetation dynamics and human influence (palynology, charcoal and coprophilous fungi analyses). Peat at the base of the profile proves lake level was below 4,700 m a.s.l. during the Pleistocene/Holocene transition. A deep-lake phase started after 11 cal ka BP, but the ostracod record indicates the level was not higher than similar to 4,720 m a.s.l. (180 m above present) and decreased gradually after the early Holocene maximum. Additional sediment ages from the basin suggest recession of Tangra Yumco from the coring site after 2.6 cal ka BP, with a shallow local lake persisting at the site until similar to 1 cal ka BP. The final peat formation indicates drier conditions thereafter. Persistence of Artemisia steppe during the Holocene lake high-stand resembles palynological records from west Tibet that indicate early Holocene aridity, in spite of high lake levels that may have resulted from meltwater input. Yet pollen assemblages indicate humidity closer to that of present potential forest areas near Lhasa, with 500-600 mm annual precipitation. Thus, the early mid-Holocene humidity was sufficient to sustain at least juniper forest, but Artemisia dominance persisted as a consequence of a combination of environmental disturbances such as (1) strong early Holocene climate fluctuations, (2) inundation of habitats suitable for forest, (3) extensive water surfaces that served as barriers to terrestrial diaspore transport from refuge areas, (4) strong erosion that denuded the non-flooded upper slopes and (5) increasing human influence since the late glacial.

The cyclotron laboratory and the RFQ accelerator in Bern

Two proton accelerators have been recently put in operation in Bern: an 18 MeV cyclotron and a 2 MeV RFQ linac. The commercial IBA 18/18 cyclotron, equipped with a specifically conceived 6 m long external beam line ending in a separate bunker, will provide beams for routine 18-F and other PET radioisotope production as well as for novel detector, radiation biophysics, radioprotection, radiochemistry and radiopharmacy developments. The accelerator is embedded into a complex building hosting two physics laboratories and four Good Manufacturing Practice (GMP) laboratories. This project is the result of a successful collaboration between the Inselspital, the University of Bern and private investors, aiming at the constitution of a combined medical and research centre able to provide the most cutting-edge technologies in medical imaging and cancer radiation therapy. The cyclotron is complemented by the RFQ with the primary goals of elemental analysis via Particle Induced Gamma Emission (PIGE), and the detection of potentially dangerous materials with high nitrogen content using the Gamma-Resonant Nuclear Absorption (GRNA) technique. In this context, beam instrumentation devices have been developed, in particular an innovative beam profile monitor based on doped silica fibres and a setup for emittance measurements using the pepper-pot technique. On this basis, the establishment of a proton therapy centre on the campus of the Inselspital is in the phase of advanced study.