Three-dimensional finite element modelling of stack pollutant emissions

[EN]In this paper we propose a finite element method approach for modelling the air quality in a local scale over complex terrain. The area of interest is up to tens of kilometres and it includes pollutant sources. The proposed methodology involves the generation of an adaptive tetrahedral mesh, the computation of an ambient wind field, the inclusion of the plume rise effect in the wind field, and the simulation of transport and reaction of pollutants. The methodology is used to simulate a fictitious pollution episode in La Palma island (Canary Island, Spain)…

The meccano method for finite element and isogeometric analysis

Congresos y conferencias

Local scale finite element modelling of stack pollutant emissions

Congresos y conferencias

Calibration and validation of an air quality finite element model around an electric power plant in Gran Canaria island

[EN]This work presents the calibration and validation of an air quality finite element model applied to the surroundings of Jinamar electric power plant in Gran Canaria island (Spain). The model involves the generation of an adaptive tetrahedral mesh, the computation of an ambient wind field, the inclusion of the plume rise effect in the wind field, and the simulation of transport and reaction of pollutants. The main advantage of the model is the treatment of complex terrains that introduces an alternative to the standard implementation of current models. In addition, it improves the computational cost through the use of unstructured meshes...

Adaptive Finite Element Method with a local element parametrization nested meses

[EN]This work presents a novel approach to solve a two dimensional problem by using an adaptive finite element approach. The most common strategy to deal with nested adaptivity is to generate a mesh that represents the geometry and the input parameters correctly, and to refine this mesh locally to obtain the most accurate solution. As opposed to this approach, the authors propose a technique using independent meshes : geometry, input data and the unknowns. Each particular mesh is obtained by a local nested refinement of the same coarse mesh at the parametric space…

Ensemble wind forecasting based on the HARMONIE model and adaptive finite elements in complex orography

[EN]Ensemble forecasting is a methodology to deal with uncertainties in the numerical wind prediction. In this work we propose to apply ensemble methods to the adaptive wind forecasting model presented in. The wind field forecasting is based on a mass-consistent model and a log-linear wind profile using as input data the resulting forecast wind from Harmonie, a Non-Hydrostatic Dynamic model used experimentally at AEMET with promising results. The mass-consistent model parameters are estimated by using genetic algorithms. The mesh is generated using the meccano method and adapted to the geometry…

Computer-algebraische und analytische Methoden zur Berechnung von Vertexfunktionen im Standardmodell

Das Standardmodell der elektroschwachen Wechselwirkung hatin den vergangenen Jahrzehnten beachtliche Erfolge erzielt.Die Suche nach

Erkennungssequenzen und strukturelle Elemente als Determinanten der regulierten Spaltung von Membranproteinen

Eine funktionell und strukturell diverse Gruppe von Transmembranproteinen wie beispielsweise Mediatoren und deren Rezeptoren können proteolytisch gespalten werden. Dieser Prozess wird als Shedding bezeichnet. Kürzlich konnte die proteolytische Aktivität identifiziert werden, die für die Prozessierung von proTNFa verantwortlich ist. Sie wurde TACE (TNF Alpha Converting Enzyme) genannt. In Experimenten mit TACE-/- Fibroblasten konnte ich herausfinden, dass das durch PMA induzierte Shedding des IL-6Rs stark reduziert war. Eine basale hydroxamatsensitive Freisetzung des IL-6Rs konnte allerdings noch detektiert werden. Um Unterschiede im Shedding von IL-6R und proTNFa zu untersuchen, generierte ich chimäre Proteine aus diesen beiden Proteinen, bei denen die Spaltstellenregionen gegeneinander vertauscht worden waren. TNFa Chimären zeigten nur sehr geringes Shedding. Im Gegensatz dazu wurden IL-6R Chimären, die die proTNFa Spaltstelle enthielten spontan gespalten. Die PMA-Induzierbarkeit war verloren gegangen. Daraufhin wurden verschiedene Chimären des unspaltbaren Proteins gp130 und der Spaltstellenpeptide aus TNFa, TGFa und IL-6R generiert. Hierbei wurde ein kurzes membranproximales Peptid aus gp130 gegen die Spaltstellen ausgetauscht. Diese Peptide übertrugen sowohl spontane ale auch PMA-induzierte Spaltbarkeit auf gp130. Um die minimalen Bedingungen für Shedding zu untersuchen, setzte ich verkürzte IL-6R Spaltstellenpeptide in gp130 ein. Die resultierenden Chimären waren empfänglich für reguliertes Shedding. Überaschenderweise konnten auch spaltbare Chimären durch das Ersetzen der membrannahem Region von gp130 durch die entsprechende Region aus dem ebenfalls nicht spaltbaren LIFR generiert werden.

Identifizierung und Charakterisierung cis regulatorischer Elemente des humanen Gens für Keratin 20 in vitro und im transgenen Mausmodell

ZusammenfassungKeratin 20 (K20) ist ein Intermediärfilament, das als Strukturelement in den Epithelien des Intestinaltrakts und den Merkelzellen der Haut exprimiert wird. Im Rahmen dieser Arbeit wurden verschiedene K20 Expressionsvektoren generiert, mit denen regulatorische Elemente des humanen Gens für K20 charakterisiert wurden. Analysiert wurde der Bereich von –4,8 kb bzw. –21,5 kb bis 12,9 kb vom Transkriptionsstart. Die Vektoren, welche die Sequenz von –21,5 kb bis 12,9 kb umfassten, konnten die Expression des EGFP Reportergens in HT-29 Zellen signifikant steigern. Zwei Enhancer-Elemente zwischen –21,5 und –18,4 kb bzw. –18,4 und –14,9 kb verstärkten die Expression der Reporterkonstrukte in vitro signifikant. Die Analyse der Vektoren in transgenen Mäusen zeigte, dass diese das Transgen gewebespezifisch exprimieren. Mit dem Bereich von –4,8 kb bis 12,9 kb ließ sich transgenspezifische mRNA Expression in intestinalen Geweben mittels RT-PCR nachweisen. Die Verwendung der Sequenz von –21,5 kb bis 21,8 kb steigerte die Expression in vivo nicht weiter.Für die gewebespezifische Expression des K20 Vektors reicht die 4,8 kb 5´upstream Sequenz aus, der Bereich bis –21,5 kb verstärkt die Expression in vitro, allerdings fehlen für die starke gewebespezifische Expression von Transgenen in vivo noch weitere Kontrollelemente.

Generalized soluble groups of finite co-central rank

Eine Gruppe G hat endlichen Prüferrang (bzw. Ko-zentralrang) kleiner gleich r, wenn für jede endlich erzeugte Gruppe H gilt: H (bzw. H modulo seinem Zentrum) ist r-erzeugbar. In der vorliegenden Arbeit werden, soweit möglich, die bekannten Sätze über Gruppen von endlichem Prüferrang (kurz X-Gruppen), auf die wesentlich größere Klasse der Gruppen mit endlichem Ko-zentralrang (kurz R-Gruppen) verallgemeinert.Für lokal nilpotente R-Gruppen, welche torsionsfrei oder p-Gruppen sind, wird gezeigt, dass die Zentrumsfaktorgruppe eine X-Gruppe sein muss. Es folgt, dass Hyperzentralität und lokale Nilpotenz für R-Gruppen identische Bediungungen sind. Analog hierzu sind R-Gruppen genau dann lokal auflösbar, wenn sie hyperabelsch sind. Zentral für die Strukturtheorie hyperabelscher R-Gruppen ist die Tatsache, dass solche Gruppen eine aufsteigende Normalreihe abelscher X-Gruppen besitzen. Es wird eine Sylowtheorie für periodische hyperabelsche R-Gruppen entwickelt. Für torsionsfreie hyperabelsche R-Gruppen wird deren Auflösbarkeit bewiesen. Des weiteren sind lokal endliche R-Gruppen fast hyperabelsch. Für R-Gruppen fallen sehr große Gruppenklassen mit den fast hyperabelschen Gruppen zusammen. Hierzu wird der Begriff der Sektionsüberdeckung eingeführt und gezeigt, dass R-Gruppen mit fast hyperabelscher Sektionsüberdeckung fast hyperabelsch sind.

Algorithmische Methoden zur Berechnung von Vierbeinfunktionen

Die vorliegende Arbeit beschäftigt sich mit derAutomatisierung von Berechnungen virtuellerStrahlungskorrekturen in perturbativen Quantenfeldtheorien.Die Berücksichtigung solcher Korrekturen aufMehrschleifen-Ebene in der Störungsreihenentwicklung istheute unabdingbar, um mit der wachsenden Präzisionexperimenteller Resultate Schritt zu halten. Im allgemeinen kinematischen Fall können heute nur dieEinschleifen-Korrekturen als theoretisch gelöst angesehenwerden -- für höhere Ordnungen liegen nur Teilergebnissevor. In Mainz sind in den letzten Jahren einige neuartigeMethoden zur Integration von Zweischleifen-Feynmandiagrammenentwickelt und im xloops-Paket in algorithmischer Formteilweise erfolgreich implementiert worden. Die verwendetenVerfahren sind eine Kombination exakter symbolischerRechenmethoden mit numerischen. DieZweischleifen-Vierbeinfunktionen stellen in diesem Rahmenein neues Kapitel dar, das durch seine große Anzahl vonfreien kinematischen Parametern einerseits leichtunüberschaubar wird und andererseits auch auf symbolischerEbene die bisherigen Anforderungen übersteigt. Sie sind ausexperimenteller Sicht aber für manche Streuprozesse vongroßem Interesse. In dieser Arbeit wurde, basierend auf einer Idee von DirkKreimer, ein Verfahren untersucht, welches die skalarenVierbeinfunktionen auf Zweischleifen-Niveau ganz ohneRandbedingungen an den Parameterraum zu integrierenversucht. Die Struktur der nach vier Residuenintegrationenauftretenden Terme konnte dabei weitgehend geklärt und dieKomplexität der auftretenden Ausdrücke soweit verkleinertwerden, dass sie von heutigen Rechnern darstellbar sind.Allerdings ist man noch nicht bei einer vollständigautomatisierten Implementierung angelangt. All dies ist dasThema von Kapitel 2. Die Weiterentwicklung von xloops über Zweibeinfunktionenhinaus erschien aus vielfältigen Gründen allerdings nichtmehr sinnvoll. Im Rahmen dieser Arbeit wurde daher einradikaler Bruch vollzogen und zusammen mit C. Bauer und A.Frink eine Programmbibliothek entworfen, die als Vehikel fürsymbolische Manipulationen dient und es uns ermöglicht,übliche symbolische Sprachen wie Maple durch C++ zuersetzen. Im dritten Kapitel wird auf die Gründeeingegangen, warum diese Umstellung sinnvoll ist, und dabeidie Bibliothek GiNaC vorgestellt. Im vierten Kapitel werdenDetails der Implementierung dann im Einzelnen vorgestelltund im fünften wird sie auf ihre Praxistauglichkeituntersucht. Anhang A bietet eine Übersicht über dieverwendeten Hilfsmittel komplexer Analysis und Anhang Bbeschreibt ein bewährtes numerisches Instrument.

Informationsstrukturen der protokollbasierten Medizin

Die protokollbasierte Medizin stellt einen interdisziplinären Brennpunkt der Informatik dar. Als besonderer Ausschnitt der medizinischen Teilgebiete erlaubt sie die relativ formale Spezifikation von Prozessen in den drei Bereichen der Prävention, Diagnose und Therapie.Letzterer wurde immer besonders fokussiert und gilt seit jeher im Rahmen klinischer Studien als Projektionsfläche für informationstechnologische Konzepte. Die Euphorie der frühen Jahre ernüchtert sich jedoch bei jeder Bilanz. Nur sehr wenige der unzählbaren Projekte haben ihre Routine in der alltäglichen Praxis gefunden. Die meisten Vorhaben sind an der Illusion der vollständigen Berechenbarkeit medizinischer Arbeitsabläufe gescheitert. Die traditionelle Sichtweise der klinischen Praxis beruht auf einer blockorientierten Vorstellung des Therapieausführungsprozesses. Sie entsteht durch seine Zerlegung in einzelne Therapiezweige, welche aus vordefinierten Blöcken zusammengesetzt sind. Diese können sequentiell oder parallel ausgeführt werden und sind selbst zusammengesetzt aus jeweils einer Menge von Elementen,welche die Aktivitäten der untersten Ebene darstellen. Das blockorientierte Aufbaumodell wird ergänzt durch ein regelorientiertes Ablaufmodell. Ein komplexes Regelwerk bestimmt Bedingungen für die zeitlichen und logischen Abhängigkeiten der Blöcke, deren Anordnung durch den Ausführungsprozeß gebildet wird. Die Modellierung der Therapieausführung steht zunächst vor der grundsätzlichen Frage, inwieweit die traditionelle Sichtweise für eine interne Repräsentation geeignet ist. Das übergeordnete Ziel besteht in der Integration der unterschiedlichen Ebenen der Therapiespezifikation. Dazu gehört nicht nur die strukturelle Komponente, sondern vorallem die Ablaufkomponente. Ein geeignetes Regelmodell ist erforderlich, welches den spezifischen Bedürfnissen der Therapieüberwachung gerecht wird. Die zentrale Aufgabe besteht darin, diese unterschiedlichen Ebenen zusammenzuführen. Eine sinnvolle Alternative zur traditionellen Sichtweise liefert das zustandsorientierte Modell des Therapieausführungsprozesses. Das zustandsorientierte Modell beruht auf der Sichtweise, daß der gesamte Therapieausführungsprozeß letztendlich eine lineare Folge von Zuständen beschreibt, wobei jeder Zustandsübergang durch ein Ereignis eingeleitet wird, an bestimmte Bedingungen geknüpft ist und bestimmte Aktionen auslösen kann. Die Parallelität des blockorientierten Modells tritt in den Hintergrund, denn die Menge der durchzuführenden Maßnahmen sind lediglich Eigenschaften der Zustände und keine strukturellen Elemente der Ablaufspezifikation. Zu jedem Zeitpunkt ist genau ein Zustand aktiv, und er repräsentiert eine von endlich vielen klinischen Situationen, mit all ihren spezifischen Aktivitäten und Ausführungsregeln. Die Vorteile des zustandsorientierten Modells liegen in der Integration. Die Grundstruktur verbindet die statische Darstellung der möglichen Phasenanordnungen mit der dynamischen Ausführung aktiver Regeln. Die ursprünglichen Inhalte des blockorientierten Modells werden als gewöhnliche Eigenschaften der Zustände reproduziert und stellen damit nur einen Spezialfall der zustandsbezogenen Sicht dar.Weitere Möglichkeiten für die Anreicherung der Zustände mit zusätzlichen Details sind denkbar wie sinnvoll. Die Grundstruktur bleibt bei jeder Erweiterung jedoch die gleiche. Es ergibt sich ein wiederverwendbares Grundgerüst,ein gemeinsamer Nenner für die Erfüllung der Überwachungsaufgabe.

Computer simulations of charged systems in partially periodic geometries

This work presents algorithms for the calculation of the electrostatic interaction in partially periodic systems. The framework for these algorithms is provided by the simulation package ESPResSo, of which the author was one of the main developers. The prominent features of the program are listed and the internal structure is described. In the following, algorithms for the calculation of the Coulomb sum in three dimensionally periodic systems are described. These methods are the foundations for the algorithms for partially periodic systems presented in this work. Starting from the MMM2D method for systems with one non-periodic coordinate, the ELC method for these systems is developed. This method consists of a correction term which allows to use methods for three dimensional periodicity also for the case of two periodic coordinates. The computation time of this correction term is neglible for large numbers of particles. The performance of MMM2D and ELC are demonstrated by results from the implementations contained in ESPResSo. It is also discussed, how different dielectric constants inside and outside of the simulation box can be realized. For systems with one periodic coordinate, the MMM1D method is derived from the MMM2D method. This method is applied to the problem of the attraction of like-charged rods in the presence of counterions, and results of the strong coupling theory for the equilibrium distance of the rods at infinite counterion-coupling are checked against results from computer simulations. The degree of agreement between the simulations at finite coupling and the theory can be characterized by a single parameter gamma_RB. In the special case of T=0, one finds under certain circumstances flat configurations, in which all charges are located in the rod-rod plane. The energetically optimal configuration and its stability are determined analytically, which depends on only one parameter gamma_z, similar to gamma_RB. These findings are in good agreement with results from computer simulations.

Finite-strain analysis in orthogneiss of the Gran Paradiso massif, Western Alps, Italy

A finite-strain study in the Gran Paradiso massif of the Italian Western Alps has been carried out to elucidate whether ductile strain shows a relationship to nappe contacts and to shed light on the nature of the subhorizontal foliation typical of the gneiss nappes in the Alps. The Rf/_ and Fry methods used on feldspar porphyroclasts from 143 augengneiss and 11 conglomerate samples of the Gran Paradiso unit (upper tectonic unit of the Gran Paradiso massif), as well as, 9 augengneiss (Erfaulet granite) and 3 quartzite conglomerate samples from the underlying Erfaulet unit (lower unit of the Gran Paradiso massif), and 1 sample from mica schist. Microstructures and thermobarometric data show that feldspar ductility at temperatures >~450°C occurred only during high-pressure metamorphism, when the rocks were underplated beneath the overriding Adriatic plate. Therefore, the finite-strain data can be related to high-pressure metamorphism in the Alpine subduction zone. The augen gneiss was heterogeneously deformed and axial ratios of the strain ellipse in XZ sections range from 2.1 to 69.8. The long axes of the finite-strain ellipsoids trend W/WNW and the short axes are subvertical associated with a subhorizontal foliation. The strain magnitudes do not increase towards the nappe contacts. Geochemical work shows that the accumulation of finite strain was not associated with any significant volume strain. Hence, the data indicate flattening strain type in the Gran Paradiso unit and constrictional strain type in the Erfaulet unit and prove deviations from simple shear. In addition, electron microprobe work was undertaken to determine if the analysed fabrics formed during high-P metamorphism. The chemistry of phengites in the studied samples suggests that deformation and final structural juxtaposition of the Gran Paradiso unit against the Erfaulet took place during high-pressure metamorphism. On the other hand, nappe stacking occurred early during subduction probably by brittle imbrication and that ductile strain was superimposed on and modified the nappe structure during high-pressure underplating in the Alpine subduction zone. The accumulation of ductile strain during underplating was not by simple shear and involved a component of vertical shortening, which caused the subhorizontal foliation in the Gran Paradiso massif. It is concluded that this foliation formed during thrusting of the nappes onto each other suggesting that nappe stacking was associated with vertical shortening. The primary evidence for this interpretation is an attenuated metamorphic section with high-pressure metamorphic rocks of the Gran Paradiso unit juxtaposed against the Erfaulet unit. Therefore, the exhumation during high-pressure metamorphism in the Alpine subduction zone involved a component of vertical shortening, which is responsible for the subhorizontal foliation within the nappes.