Diskontinuierliche Galerkin-Verfahren für die operationelle Wettervorhersage

In dieser Arbeit wird ein neuer Dynamikkern entwickelt und in das bestehendernnumerische Wettervorhersagesystem COSMO integriert. Für die räumlichernDiskretisierung werden diskontinuierliche Galerkin-Verfahren (DG-Verfahren)rnverwendet, für die zeitliche Runge-Kutta-Verfahren. Hierdurch ist ein Verfahrenrnhoher Ordnung einfach zu realisieren und es sind lokale Erhaltungseigenschaftenrnder prognostischen Variablen gegeben. Der hier entwickelte Dynamikkern verwendetrngeländefolgende Koordinaten in Erhaltungsform für die Orographiemodellierung undrnkoppelt das DG-Verfahren mit einem Kessler-Schema für warmen Niederschlag. Dabeirnwird die Fallgeschwindigkeit des Regens, nicht wie üblich implizit imrnKessler-Schema diskretisiert, sondern explizit im Dynamikkern. Hierdurch sindrndie Zeitschritte der Parametrisierung für die Phasenumwandlung des Wassers undrnfür die Dynamik vollständig entkoppelt, wodurch auch sehr große Zeitschritte fürrndie Parametrisierung verwendet werden können. Die Kopplung ist sowohl fürrnOperatoraufteilung, als auch für Prozessaufteilung realisiert.rnrnAnhand idealisierter Testfälle werden die Konvergenz und die globalenrnErhaltungseigenschaften des neu entwickelten Dynamikkerns validiert. Die Massernwird bis auf Maschinengenauigkeit global erhalten. Mittels Bergüberströmungenrnwird die Orographiemodellierung validiert. Die verwendete Kombination ausrnDG-Verfahren und geländefolgenden Koordinaten ermöglicht die Behandlung vonrnsteileren Bergen, als dies mit dem auf Finite-Differenzenverfahren-basierendenrnDynamikkern von COSMO möglich ist. Es wird gezeigt, wann die vollernTensorproduktbasis und wann die Minimalbasis vorteilhaft ist. Die Größe desrnEinflusses auf das Simulationsergebnis der Verfahrensordnung, desrnParametrisierungszeitschritts und der Aufteilungsstrategie wirdrnuntersucht. Zuletzt wird gezeigt dass bei gleichem Zeitschritt die DG-Verfahrenrnaufgrund der besseren Skalierbarkeit in der Laufzeit konkurrenzfähig zurnFinite-Differenzenverfahren sind.

Molekulardynamik-Simulation zur Struktur von SiO 2-Oberflächen mit adsorbiertem Wasser

In der vorliegenden Arbeit werdenMolekulardynamik-Simulationen zur Untersuchung derstatischen Eigenschaften von amorphenSiliziumdioxidoberflächen (Siliziumdioxid) durchgeführt. Da das von van Beest, Kramer und van Santen vorgeschlagene,sogenannte BKS-Potential für Bulksysteme optimiert wurde und an Oberflächen deutlichandere Ladungsverteilungenauftreten als im Bulk, ist die Anwendbarkeit diesesPotentials für Oberflächensystemefraglich. Aus diesem Grund haben wir untersucht, inwieweitsich die Oberflächeneigenschaften von Systemen, die mit Hilfe des BKS-Potentials äquilibriertwurden, durch ein Nachrelaxieren mit einer ab-initio-Simulation (Car-Parrinello-Methode)ändern. Mit Hilfe der Kombination aus BKS- und Car-Parrinello-Methode (CPMD)konnten wir feststellen, daß sich die Systeme aufgrund des Nachrelaxierens in z-Richtungweiter ausdehnen. Desweiteren zeigte sich insbesondere bei kleinen Ringen (kommen nur ander Oberfläche vor), daß es deutliche Abweichungen in den Geometrien (Atomabstände,Winkel usw.) zwischen der reinen BKS- und der kombinierten BKS-CPMD-Methode gibt. Anhand vonCPMD-Simulationen konnten wir zeigen, daß es durch die Wechselwirkung eines Wassermolekülsmit einem 2er-Ring zum Aufbrechen dieser Ringstruktur und zur Bildung von zweiSilanolgruppen (SiOH) kommt. Desweiteren stellten wir fest, daß es sich hierbei um eineexotherme Reaktion (Energiedifferenz 1.6 eV) handelt, für die eineEnergiebarriere von 1.1 eV überwunden werden muß. Ferner ergab sich, daß die an der Bildung des2er-Ringes beteiligten, stark deformierten Tetraeder nach dem Aufbrechen dieserRingstruktur eine nahezu ideale Tetraederform annehmen.

Novel simulation methods for Coulomb and hydrodynamic interactions

This thesis presents new methods to simulate systems with hydrodynamic and electrostatic interactions. Part 1 is devoted to computer simulations of Brownian particles with hydrodynamic interactions. The main influence of the solvent on the dynamics of Brownian particles is that it mediates hydrodynamic interactions. In the method, this is simulated by numerical solution of the Navier--Stokes equation on a lattice. To this end, the Lattice--Boltzmann method is used, namely its D3Q19 version. This model is capable to simulate compressible flow. It gives us the advantage to treat dense systems, in particular away from thermal equilibrium. The Lattice--Boltzmann equation is coupled to the particles via a friction force. In addition to this force, acting on {it point} particles, we construct another coupling force, which comes from the pressure tensor. The coupling is purely local, i.~e. the algorithm scales linearly with the total number of particles. In order to be able to map the physical properties of the Lattice--Boltzmann fluid onto a Molecular Dynamics (MD) fluid, the case of an almost incompressible flow is considered. The Fluctuation--Dissipation theorem for the hybrid coupling is analyzed, and a geometric interpretation of the friction coefficient in terms of a Stokes radius is given. Part 2 is devoted to the simulation of charged particles. We present a novel method for obtaining Coulomb interactions as the potential of mean force between charges which are dynamically coupled to a local electromagnetic field. This algorithm scales linearly, too. We focus on the Molecular Dynamics version of the method and show that it is intimately related to the Car--Parrinello approach, while being equivalent to solving Maxwell's equations with freely adjustable speed of light. The Lagrangian formulation of the coupled particles--fields system is derived. The quasi--Hamiltonian dynamics of the system is studied in great detail. For implementation on the computer, the equations of motion are discretized with respect to both space and time. The discretization of the electromagnetic fields on a lattice, as well as the interpolation of the particle charges on the lattice is given. The algorithm is as local as possible: Only nearest neighbors sites of the lattice are interacting with a charged particle. Unphysical self--energies arise as a result of the lattice interpolation of charges, and are corrected by a subtraction scheme based on the exact lattice Green's function. The method allows easy parallelization using standard domain decomposition. Some benchmarking results of the algorithm are presented and discussed.