IMEX finite volume methods for the shallow water equations

Die Flachwassergleichungen (SWE) sind ein hyperbolisches System von Bilanzgleichungen, die adäquate Approximationen an groß-skalige Strömungen der Ozeane, Flüsse und der Atmosphäre liefern. Dabei werden Masse und Impuls erhalten. Wir unterscheiden zwei charakteristische Geschwindigkeiten: die Advektionsgeschwindigkeit, d.h. die Geschwindigkeit des Massentransports, und die Geschwindigkeit von Schwerewellen, d.h. die Geschwindigkeit der Oberflächenwellen, die Energie und Impuls tragen. Die Froude-Zahl ist eine Kennzahl und ist durch das Verhältnis der Referenzadvektionsgeschwindigkeit zu der Referenzgeschwindigkeit der Schwerewellen gegeben. Für die oben genannten Anwendungen ist sie typischerweise sehr klein, z.B. 0.01. Zeit-explizite Finite-Volume-Verfahren werden am öftersten zur numerischen Berechnung hyperbolischer Bilanzgleichungen benutzt. Daher muss die CFL-Stabilitätsbedingung eingehalten werden und das Zeitinkrement ist ungefähr proportional zu der Froude-Zahl. Deswegen entsteht bei kleinen Froude-Zahlen, etwa kleiner als 0.2, ein hoher Rechenaufwand. Ferner sind die numerischen Lösungen dissipativ. Es ist allgemein bekannt, dass die Lösungen der SWE gegen die Lösungen der Seegleichungen/ Froude-Zahl Null SWE für Froude-Zahl gegen Null konvergieren, falls adäquate Bedingungen erfüllt sind. In diesem Grenzwertprozess ändern die Gleichungen ihren Typ von hyperbolisch zu hyperbolisch.-elliptisch. Ferner kann bei kleinen Froude-Zahlen die Konvergenzordnung sinken oder das numerische Verfahren zusammenbrechen. Insbesondere wurde bei zeit-expliziten Verfahren falsches asymptotisches Verhalten (bzgl. der Froude-Zahl) beobachtet, das diese Effekte verursachen könnte.Ozeanographische und atmosphärische Strömungen sind typischerweise kleine Störungen eines unterliegenden Equilibriumzustandes. Wir möchten, dass numerische Verfahren für Bilanzgleichungen gewisse Equilibriumzustände exakt erhalten, sonst können künstliche Strömungen vom Verfahren erzeugt werden. Daher ist die Quelltermapproximation essentiell. Numerische Verfahren die Equilibriumzustände erhalten heißen ausbalanciert.rnrnIn der vorliegenden Arbeit spalten wir die SWE in einen steifen, linearen und einen nicht-steifen Teil, um die starke Einschränkung der Zeitschritte durch die CFL-Bedingung zu umgehen. Der steife Teil wird implizit und der nicht-steife explizit approximiert. Dazu verwenden wir IMEX (implicit-explicit) Runge-Kutta und IMEX Mehrschritt-Zeitdiskretisierungen. Die Raumdiskretisierung erfolgt mittels der Finite-Volumen-Methode. Der steife Teil wird mit Hilfe von finiter Differenzen oder au eine acht mehrdimensional Art und Weise approximniert. Zur mehrdimensionalen Approximation verwenden wir approximative Evolutionsoperatoren, die alle unendlich viele Informationsausbreitungsrichtungen berücksichtigen. Die expliziten Terme werden mit gewöhnlichen numerischen Flüssen approximiert. Daher erhalten wir eine Stabilitätsbedingung analog zu einer rein advektiven Strömung, d.h. das Zeitinkrement vergrößert um den Faktor Kehrwert der Froude-Zahl. Die in dieser Arbeit hergeleiteten Verfahren sind asymptotisch erhaltend und ausbalanciert. Die asymptotischer Erhaltung stellt sicher, dass numerische Lösung das "korrekte" asymptotische Verhalten bezüglich kleiner Froude-Zahlen besitzt. Wir präsentieren Verfahren erster und zweiter Ordnung. Numerische Resultate bestätigen die Konvergenzordnung, so wie Stabilität, Ausbalanciertheit und die asymptotische Erhaltung. Insbesondere beobachten wir bei machen Verfahren, dass die Konvergenzordnung fast unabhängig von der Froude-Zahl ist.

Hydraulic evaluation of the gyro electric wave energy converter

Motivation Thanks for a scholarship offered by ALma Mater Studiorum I could stay in Denmark for six months during which I could do physical tests on the device Gyro PTO at the Departmet of Civil Engineering of Aalborg University. Aim The goal of my thesis is an hydraulic evaluation of the device: Gyro PTO, a gyroscopic device for conversion of mechanical energy in ocean surface waves to electrical energy. The principle of the system is the application of the gyroscopic moment of flywheels equipped on a swing float excited by waves. The laboratory activities were carried out by: Morten Kramer, Jan Olsen, Irene Guaraldi, Morten Thøtt, Nikolaj Holk. The main purpose of the tests was to investigate the power absorption performance in irregular waves, but testing also included performance measures in regular waves and simple tests to get knowledge about characteristics of the device, which could facilitate the possibility of performing numerical simulations and optimizations. Methodology To generate the waves and measure the performance of the device a workstation was created in the laboratory. The workstation consist of four computers in each of wich there was a different program. Programs have been used : Awasys6, LabView, Wave lab, Motive optitrack, Matlab, Autocad Main Results Thanks to the obtained data with the tank testing was possible to make the process of wave analisys. We obtained significant wave height and period through a script Matlab and then the values of power produced, and energy efficiency of the device for two types of waves: regular and irregular. We also got results as: physical size, weight, inertia moments, hydrostatics, eigen periods, mooring stiffness, friction, hydrodynamic coefficients etc. We obtained significant parameters related to the prototype in the laboratory after which we scale up the results obtained for two future applications: one in Nissun Brending and in the North Sea. Conclusions The main conclusion on the testing is that more focus should be put into ensuring a stable and positive power output in a variety of wave conditions. In the irregular waves the power production was negative and therefore it does not make sense to scale up the results directly. The average measured capture width in the regular waves was 0.21 m. As the device width is 0.63 m this corresponds to a capture width ratio of: 0.21/0.63 * 100 = 33 %. Let’s assume that it is possible to get the device to produce as well in irregular waves under any wave conditions, and lets further assume that the yearly absorbed energy can be converted into electricity at a PTO-efficiency of 90 %. Under all those assumptions the results in table are found, i.e. a Nissum Bredning would produce 0.87 MWh/year and a North Sea device 85 MWh/year.

Triggering of volcanic activity by large earthquakes

Statistical analyses of temporal relationships between large earthquakes and volcanic eruptions suggest seismic waves may trigger eruptions even over great (>1000 km) distances, although the causative mechanism is not well constrained. In this study the relationship between large earthquakes and subtle changes in volcanic activity was investigated in order to gain greater insight into the relationship between dynamic stresses propagated by surface waves and volcanic response. Daily measurements from the Ozone Monitoring Instrument (OMI), onboard the Aura satellite, provide constraints on volcanic sulfur-dioxide (SO2) emission rates as a measure of subtle changes in activity. Time series of SO2 emission rates were produced from OMI data for thirteen persistently active volcanoes from 1 October 2004 to 30 September 2010. In order to quantify the affect of earthquakes at teleseismic distances, we modeled surface-wave amplitudes from the source mechanisms of moment magnitude (Mw) ≥7 earthquakes, and calculated the Peak Dynamic Stress (PDS). We assessed the influence of earthquakes on volcanic activity in two ways: 1) by identifying increases in the SO2 time series data and looking for causative earthquakes and 2) by examining the average emission rate before and after each earthquake. In the first, the SO2 time series for each volcano was used to calculate a baseline threshold for comparison with post-earthquake emission. Next, we generated a catalog of responses based on sustained SO2 emission increases above this baseline. Delay times between each SO2 response and each prior earthquake were analyzed using both the actual earthquake catalog, and a randomly generated catalog of earthquakes. This process was repeated for each volcano. Despite varying multiple parameters, this analysis did not demonstrate a clear relationship between earthquake-generated PDS and SO2 emission. However, the second analysis, which was based on the occurrence of large earthquakes indicated a response at most volcanoes. Using the PDS calculations as a filtering criterion for the earthquake catalog, the SO2 mass for each volcano was analyzed in 28-day windows centered on the earthquake origin time. If the average SO2 mass after the earthquake was greater than an arbitrary percentage of pre-earthquake mass, we identified the volcano as having a response to the event. This window analysis provided insight on what type of volcanic activity is more susceptible to triggering by dynamic stress. The volcanoes with very open systems included in this study, Ambrym, Gaua, Villarrica, Erta Ale and, Turrialba, showed a clear response to dynamic stress while the volcanoes with more closed systems, Merapi, Semeru, Fuego, Pacaya, and Bagana, showed no response.

A MEMS Light Modulator Based on Diffractive Nanohole Gratings

We present the design, fabrication, and testing of a microelectromechanical systems (MEMS) light modulator based on pixels patterned with periodic nanohole arrays. Flexure-suspended silicon pixels are patterned with a two dimensional array of 150 nm diameter nanoholes using nanoimprint lithography. A top glass plate assembled above the pixel array is used to provide a counter electrode for electrostatic actuation. The nanohole pattern is designed so that normally-incident light is coupled into an in-plane grating resonance, resulting in an optical stop-band at a desired wavelength. When the pixel is switched into contact with the top plate, the pixel becomes highly reflective. A 3:1 contrast ratio at the resonant wavelength is demonstrated for gratings patterned on bulk Si substrates. The switching time is 0.08 ms and the switching voltage is less than 15V.

Decade Wind and Wave Data for NW Iberia, Cabo Sillerio

Seafloor sediment mobilization on the inner Northwest Iberian continental shelf is caused largely by ocean surface waves. The temporal and spatial variability in the wave height, wave period, and wave direction has a profound effect on local sediment mobilization, leading to distinct sediment mobilization scenarios. Six grain-size specific sediment mobilization scenarios, representing seasonal average and storm conditions, were simulated with a physics-based numerical model. Model inputs included meteorological and oceanographic data in conjunction with seafloor grain-size and the shelf bathymetric data. The results show distinct seasonal variations, most importantly in wave height, leading to sediment mobilization, specifically on the inner shelf shallower than 30 m water depth where up to 49% of the shelf area is mobilized. Medium to severe storm events are modeled to mobilize up to 89% of the shelf area above 150 m water depth. The frequency of each of these seasonal and storm-related sediment mobilization scenarios is addressed using a decade of meteorological and oceanographic data. The temporal and spatial patterns of the modeled sediment mobilization scenarios are discussed in the context of existing geological and environmental processes and conditions to assist scientific, industrial and environmental efforts that are directly affected by sediment mobilization. Examples, where sediment mobilization plays a vital role, include seafloor nutrient advection, recurrent arrival of oil from oil-spill-laden seafloor sediment, and bottom trawling impacts.

Abb. 3 Change of altitude in the EGIG-West-East-profile

The members of the two International Glaciological Greenland Expeditions (EGIG) in 1959 and 1968 determined the altitudes in the West-East-profile by levelling. The comparison of these two measurements gives evidence of the change of altitudes and of the surface-waves. The ways of measurement and of analysis are described here. The graphic representation of the result is discussed here in connexion with the possible causes.