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.

The Equidistant Method - a novel hip joint simulation algorithm for detection of femoroacetabular impingement

A novel computerized algorithm for hip joint motion simulation and collision detection, called the Equidistant Method, has been developed. This was compared to three pre-existing methods having different properties regarding definition of the hip joint center and behavior after collision detection. It was proposed that the Equidistant Method would be most accurate in detecting the location and extent of femoroacetabular impingement.

Choosing Smoothness Parameters for Smoothing Splines by Minimizing and Estimate of Risk

Smoothing splines are a popular approach for non-parametric regression problems. We use periodic smoothing splines to fit a periodic signal plus noise model to data for which we assume there are underlying circadian patterns. In the smoothing spline methodology, choosing an appropriate smoothness parameter is an important step in practice. In this paper, we draw a connection between smoothing splines and REACT estimators that provides motivation for the creation of criteria for choosing the smoothness parameter. The new criteria are compared to three existing methods, namely cross-validation, generalized cross-validation, and generalization of maximum likelihood criteria, by a Monte Carlo simulation and by an application to the study of circadian patterns. For most of the situations presented in the simulations, including the practical example, the new criteria out-perform the three existing criteria.

Mortality from neurodegenerative disease and exposure to extremely low-frequency magnetic fields: 31 years of observations on Swiss railway employees

AIMS: The objective of the present study was to investigate the relationship between extremely low-frequency magnetic field (ELF-MF) exposure and mortality from several neurodegenerative conditions in Swiss railway employees. METHODS: We studied a cohort of 20,141 Swiss railway employees with 464,129 person-years of follow-up between 1972 and 2002. For each individual, cumulative exposure was calculated from on-site measurements and modelling of past exposure. We compared cause-specific mortality in highly exposed train drivers (mean exposure: 21 microT) with less exposed occupational groups (for example station masters: 1 microT). RESULTS: The hazard ratio for train drivers compared to station masters was 1.96 [95% confidence interval (CI) = 0.98-3.92] for senile dementia and 3.15 (95% CI = 0.90-11.04) for Alzheimer's disease. For every 10 microT years of cumulative exposure senile dementia mortality increased by 5.7% (95% CI = 1.3-10.4), Alzheimer's disease by 9.4% (95% CI = 2.7-16.4) and amyotrophic lateral sclerosis by 2.1% (95% CI = -6.8 to 11.7). There was no evidence for an increase in mortality from Parkinson's disease and multiple sclerosis. CONCLUSIONS: This study suggests a link between exposure to ELF-MF and Alzheimer's disease and indicates that ELF-MF might act in later stages of the disease process.

Leukaemia, brain tumours and exposure to extremely low frequency magnetic fields: cohort study of Swiss railway employees

AIMS: To investigate the relationship between extremely low frequency magnetic field (ELF-MF) exposure and mortality from leukaemia and brain tumour in a cohort of Swiss railway workers. METHODS: 20,141 Swiss railway employees with 464,129 person-years of follow-up between 1972 and 2002 were studied. Mortality rates for leukaemia and brain tumour of highly exposed train drivers (21 muT average annual exposure) were compared with medium and low exposed occupational groups (i.e. station masters with an average exposure of 1 muT). In addition, individual cumulative exposure was calculated from on-site measurements and modelling of past exposures. RESULTS: The hazard ratio (HR) for leukaemia mortality of train drivers was 1.43 (95% CI 0.74 to 2.77) compared with station masters. For myeloid leukaemia the HR of train drivers was 4.74 (95% CI 1.04 to 21.60) and for Hodgkin's disease 3.29 (95% CI 0.69 to 15.63). Lymphoid leukaemia, non-Hodgkin's disease and brain tumour mortality were not associated with magnetic field exposure. Concordant results were obtained from analyses based on individual cumulative exposure. CONCLUSIONS: Some evidence of an exposure-response association was found for myeloid leukaemia and Hodgkin's disease, but not for other haematopoietic and lymphatic malignancies and brain tumours.

Critical evaluation of diameter increment modelling in the Great Lakes region

Simulations of forest stand dynamics in a modelling framework including Forest Vegetation Simulator (FVS) are diameter driven, thus the diameter or basal area increment model needs a special attention. This dissertation critically evaluates diameter or basal area increment models and modelling approaches in the context of the Great Lakes region of the United States and Canada. A set of related studies are presented that critically evaluate the sub-model for change in individual tree basal diameter used in the Forest Vegetation Simulator (FVS), a dominant forestry model in the Great Lakes region. Various historical implementations of the STEMS (Stand and Tree Evaluation and Modeling System) family of diameter increment models, including the current public release of the Lake States variant of FVS (LS-FVS), were tested for the 30 most common tree species using data from the Michigan Forest Inventory and Analysis (FIA) program. The results showed that current public release of the LS-FVS diameter increment model over-predicts 10-year diameter increment by 17% on average. Also the study affirms that a simple adjustment factor as a function of a single predictor, dbh (diameter at breast height) used in the past versions, provides an inadequate correction of model prediction bias. In order to re-engineer the basal diameter increment model, the historical, conceptual and philosophical differences among the individual tree increment model families and their modelling approaches were analyzed and discussed. Two underlying conceptual approaches toward diameter or basal area increment modelling have been often used: the potential-modifier (POTMOD) and composite (COMP) approaches, which are exemplified by the STEMS/TWIGS and Prognosis models, respectively. It is argued that both approaches essentially use a similar base function and neither is conceptually different from a biological perspective, even though they look different in their model forms. No matter what modelling approach is used, the base function is the foundation of an increment model. Two base functions – gamma and Box-Lucas – were identified as candidate base functions for forestry applications. The results of a comparative analysis of empirical fits showed that quality of fit is essentially similar, and both are sufficiently detailed and flexible for forestry applications. The choice of either base function in order to model diameter or basal area increment is dependent upon personal preference; however, the gamma base function may be preferred over the Box-Lucas, as it fits the periodic increment data in both a linear and nonlinear composite model form. Finally, the utility of site index as a predictor variable has been criticized, as it has been widely used in models for complex, mixed species forest stands though not well suited for this purpose. An alternative to site index in an increment model was explored, using site index and a combination of climate variables and Forest Ecosystem Classification (FEC) ecosites and data from the Province of Ontario, Canada. The results showed that a combination of climate and FEC ecosites variables can replace site index in the diameter increment model.

Erosion Damage Mapping: Assessing Current Soil Erosion Damage in Switzerland

Within the scope of a comprehensive assessment of the degree of soil erosion in Switzerland, common methods have been used in the past including test plot measurements, artificial rainfall simulation, and erosion modelling. In addition, mapping guidelines for all visible erosion features have been developed since the 1970s and are being successfully applied in many research and soil conservation projects. Erosion damage has been continuously mapped over a period of 9 years in a test region in the central Bernese plateau. In 2005, two additional study areas were added. The present paper assesses the data gathered and provides a comparison of the three study areas within a period of one year (from October 2005 to October 2006), focusing on the on-site impacts of soil erosion. During this period, about 11 erosive rainfall events occurred. Average soil loss rates mapped at each study site amounted to 0.7 t ha-1, 1.2 t ha-1 and 2.3 t ha-1, respectively. About one fourth of the total arable land showed visible erosion damage. Maximum soil losses of about 70 t ha-1 occurred on individual farm plots. Average soil erosion patterns are widely used to underline the severity of an erosion problem (e.g. impacts on water bodies). But since severe rainfall events, wheel tracks, headlands, and other “singularities” often cause high erosion rates, analysis of extreme erosion patterns such as maximum values led to a more differentiated understanding and appropriate conclusions for planning and design of soil protection measures. The study contains an assessment of soil erosion in Switzerland, emphasizing questions about extent, frequency and severity. At the same time, the effects of different types of land management are investigated in the field, aiming at the development of meaningful impact indicators of (un-)sustainable agriculture/soil erosion risk as well as the validation of erosion models. The results illustrate that conservation agriculture including no-till, strip tillage and in-mulch seeding plays an essential role in reducing soil loss as compared to conventional tillage.

Learning for sustainability - LforS: Simulation games - a creative tool for interactive learning

Description of simulation and training games as tool for awareness and capacity development in multi steakeholder processes

Calibrating Wireless Sensor Network Simulation Models with Real-World Experiments

This paper studies the energy-efficiency and service characteristics of a recently developed energy-efficient MAC protocol for wireless sensor networks in simulation and on a real sensor hardware testbed. This opportunity is seized to illustrate how simulation models can be verified by cross-comparing simulation results with real-world experiment results. The paper demonstrates that by careful calibration of simulation model parameters, the inevitable gap between simulation models and real-world conditions can be reduced. It concludes with guidelines for a methodology for model calibration and validation of sensor network simulation models.

The MIRELA framework: modeling and analyzing mixed reality applications using timed automata

Mixed Reality (MR) aims to link virtual entities with the real world and has many applications such as military and medical domains [JBL+00, NFB07]. In many MR systems and more precisely in augmented scenes, one needs the application to render the virtual part accurately at the right time. To achieve this, such systems acquire data related to the real world from a set of sensors before rendering virtual entities. A suitable system architecture should minimize the delays to keep the overall system delay (also called end-to-end latency) within the requirements for real-time performance. In this context, we propose a compositional modeling framework for MR software architectures in order to specify, simulate and validate formally the time constraints of such systems. Our approach is first based on a functional decomposition of such systems into generic components. The obtained elements as well as their typical interactions give rise to generic representations in terms of timed automata. A whole system is then obtained as a composition of such defined components. To write specifications, a textual language named MIRELA (MIxed REality LAnguage) is proposed along with the corresponding compilation tools. The generated output contains timed automata in UPPAAL format for simulation and verification of time constraints. These automata may also be used to generate source code skeletons for an implementation on a MR platform. The approach is illustrated first on a small example. A realistic case study is also developed. It is modeled by several timed automata synchronizing through channels and including a large number of time constraints. Both systems have been simulated in UPPAAL and checked against the required behavioral properties.

Simulationsgestützte Untersuchung der Spurtreue von Routenzügen

Für den innerbetrieblichen Transport und besonders für die Materialversorgung in der Produktion werden zunehmend Routenzüge eingesetzt, die aus einem Schleppfahrzeug und bis zu fünf Anhängern bestehen. Um gute Manövrierbarkeit und einen geringen Platzbedarf zu gewährleisten, sollten Routenzuganhänger möglichst spurtreu sein. In diesem Beitrag wird die Spurtreue exemplarisch für zwei am Markt verfügbare Anhänger für Routenzüge mit zwei ungelenkten und vier gelenkten Rädern untersucht. Hierfür wird zunächst ein Gütekriterium definiert, das die maximale Abweichung von der Spurtreue quantitativ erfasst. Außerdem werden Testszenarien vorgeschlagen, um die verschiedenen Fahrwerks- und Lenkkonzepte vergleichen zu können. Mit Hilfe eines entwickelten analytischen Modells werden die Abweichungen in der Spurtreue bei stationärer Kreisfahrt für die zwei gewählten Konzepte berechnet und dargestellt. Zusätzlich wird eine Mehrkörper-Simulation, die eine tiefere physikalische Modellierung und die Untersuchung komplexerer Fahrmanöver erlaubt, durchgeführt und mit den analytischen Ergebnissen verglichen. Es kann gezeigt werden, dass neben der Lenkkinematik weitere Parameter wie das Schräglaufverhalten der Reifen Einfluss auf die Spurtreue haben.

Design, Simulation und Erprobung einer aktiven Sitzaufhängung für Nutzfahrzeuge

Nutzfahrzeuge müssen oft durch sehr unebenes Gelände gefahren werden. In diesem Fall wird der Fahrer starken Vibrationen ausgesetzt, die von der Fahrzeugkarosserie durch die Sitzaufhängung auf ihn wirken. Um diese Schwingungen zu verringern, werden die Sitzaufhängungen in der Regel mit Feder-Dämpfer-Systemen ausgerüstet. Jedoch erreichen die passiven Systeme vor allem bei niederfrequenten Schwingungen ihre physikalischen Grenzen. Eine wesentliche Verbesserung des Sitzkomforts kann unter solchen Anregungsbedingungen nur mit einer aktiven Sitzaufhängung erreicht werden. In diesem Beitrag wird ein neuartiges aktives System für die Sitzaufhängung auf Basis von elektrorheologischen Flüssigkeiten vorgestellt. Außerdem werden die theoretischen Grundlagen für die Modellierung der beschriebenen aktiven Sitzaufhängung dargestellt. Anschließend werden die Simulationsergebnisse mit den Messergebnissen unter realen Betriebsbedingungen verglichen. Die Repräsentation der Ergebnisse mit Hilfe der im Bereich der Sitztechnik weit verbreiteten SEAT-Werten (Seat Effective Amplitude Transmissibility) zeigt das Potenzial des entwickelten Systems zur aktiven Reduktion der Schwingungsbelastung des Fahrers und ermöglicht seine objektive Bewertung.

A Routing Protocol Based on Energy and Link Quality for Internet of Things Applications

The Internet of Things (IoT) is attracting considerable attention from the universities, industries, citizens and governments for applications, such as healthcare, environmental monitoring and smart buildings. IoT enables network connectivity between smart devices at all times, everywhere, and about everything. In this context, Wireless Sensor Networks (WSNs) play an important role in increasing the ubiquity of networks with smart devices that are low-cost and easy to deploy. However, sensor nodes are restricted in terms of energy, processing and memory. Additionally, low-power radios are very sensitive to noise, interference and multipath distortions. In this context, this article proposes a routing protocol based on Routing by Energy and Link quality (REL) for IoT applications. To increase reliability and energy-efficiency, REL selects routes on the basis of a proposed end-to-end link quality estimator mechanism, residual energy and hop count. Furthermore, REL proposes an event-driven mechanism to provide load balancing and avoid the premature energy depletion of nodes/networks. Performance evaluations were carried out using simulation and testbed experiments to show the impact and benefits of REL in small and large-scale networks. The results show that REL increases the network lifetime and services availability, as well as the quality of service of IoT applications. It also provides an even distribution of scarce network resources and reduces the packet loss rate, compared with the performance of well-known protocols.