The modulational instability in deep water under the action of wind and dissipation

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Métodos analíticos, numéricos e experimentais para o cálculo de ondas de impacto em meios líquidos

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Dynamik orographischer Bannerwolken

Der Begriff "Bannerwolke" bezeichnet ein eindrucksvolles Phänomen aus dem Bereich der Gebirgsmeteorologie. Bannerwolken können gelegentlich im Hochgebirge im Bereich steiler Bergspitzen oder langgezogener Bergrücken, wie z.B. dem Matterhorn in den Schweizer Alpen oder dem Zugspitzgrat in den Bayrischen Alpen beobachtet werden. Der Begriff bezeichnet eine Banner- oder Fahnen-ähnliche Wolkenstruktur, welche an der windabgewandten Seite des Berges befestigt zu sein scheint, während die windzugewandte Seite vollkommen wolkenfrei ist. Bannerwolken fanden bislang, trotz ihres relativ häufigen Auftretens in der wissenschaftlichen Literatur kaum Beachtung. Entsprechend wenig ist über ihren Entstehungsmechanismus und insbesondere die relative Bedeutung dynamischer gegenüber thermodynamischer Prozesse bekannt. In der wissenschaftlichen Literatur wurden bislang 3 unterschiedliche Mechanismen postuliert, um die Entstehung von Bannerwolken zu erklären. Demnach entstehen Bannerwolken durch (a) den Bernoulli-Effekt, insbesondere durch die lokale adiabatische Kühlung hervorgerufen durch eine Druckabnahme entlang quasi-horizontal verlaufender, auf der windzugewandten Seite startender Trajektorien, (b) durch isobare Mischung bodennaher kälterer Luft mit wärmerer Luft aus höheren Schichten, oder (c) durch erzwungene Hebung im aufsteigenden Ast eines Leerotors. Ziel dieser Arbeit ist es, ein besseres physikalisches Verständnis für das Phänomen der Bannerwolke zu entwickeln. Das Hauptaugenmerk liegt auf dem dominierenden Entstehungsmechanismus, der relativen Bedeutung dynamischer und thermodynamischer Prozesse, sowie der Frage nach geeigneten meteorologischen Bedingungen. Zu diesem Zweck wurde ein neues Grobstruktursimulations (LES)-Modell entwickelt, welches geeignet ist turbulente, feuchte Strömungen in komplexem Terrain zu untersuchen. Das Modell baut auf einem bereits existierenden mesoskaligen (RANS) Modell auf. Im Rahmen dieser Arbeit wurde das neue Modell ausführlich gegen numerische Referenzlösungen und Windkanal-Daten verglichen. Die wesentlichen Ergebnisse werden diskutiert, um die Anwendbarkeit des Modells auf die vorliegende wissenschaftliche Fragestellung zu überprüfen und zu verdeutlichen. Die Strömung über eine idealisierte pyramidenförmige Bergspitze wurde für Froude-Zahlen Fr >> 1 sowohl auf Labor- als auch atmosphärischer Skala mit und ohne Berücksichtigung der Feuchtephysik untersucht. Die Simulationen zeigen, dass Bannerwolken ein primär dynamisches Phänomen darstellen. Sie entstehen im Lee steiler Bergspitzen durch dynamisch erzwungene Hebung. Die Simulationen bestätigen somit die Leerotor-Theorie. Aufgrund des stark asymmetrischen, Hindernis-induzierten Strömungsfeldes können Bannerwolken sogar im Falle horizontal homogener Anfangsbedingungen hinsichtlich Feuchte und Temperatur entstehen. Dies führte zu der neuen Erkenntnis, dass zusätzliche leeseitige Feuchtequellen, unterschiedliche Luftmassen in Luv und Lee, oder Strahlungseffekte keine notwendige Voraussetzung für die Entstehung einer Bannerwolke darstellen. Die Wahrscheinlichkeit der Bannerwolkenbildung steigt mit zunehmender Höhe und Steilheit des pyramidenförmigen Hindernisses und ist in erster Näherung unabhängig von dessen Orientierung zur Anströmung. Simulationen mit und ohne Berücksichtigung der Feuchtephysik machen deutlich, dass thermodynamische Prozesse (insbes. die Umsetzung latenter Wärme) für die Dynamik prototypischer (nicht-konvektiver) Bannerwolken zweitrangig ist. Die Verstärkung des aufsteigenden Astes im Lee und die resultierende Wolkenbildung, hervorgerufen durch die Freisetzung latenter Wärme, sind nahezu vernachlässigbar. Die Feuchtephysik induziert jedoch eine Dipol-ähnliche Struktur im Vertikalprofil der Brunt-Väisälä Frequenz, was zu einem moderaten Anstieg der leeseitigen Turbulenz führt. Es wird gezeigt, dass Gebirgswellen kein entscheidendes Ingredienz darstellen, um die Dynamik von Bannerwolken zu verstehen. Durch eine Verstärkung der Absinkbewegung im Lee, haben Gebirgswellen lediglich die Tendenz die horizontale Ausdehnung von Bannerwolken zu reduzieren. Bezüglich geeigneter meteorologischer Bedingungen zeigen die Simulationen, dass unter horizontal homogenen Anfangsbedingungen die äquivalentpotentielle Temperatur in der Anströmung mit der Höhe abnehmen muss. Es werden 3 notwendige und hinreichende Kriterien, basierend auf dynamischen und thermodynamischen Variablen vorgestellt, welche einen weiteren Einblick in geeignete meteorologische Bedingungen geben.

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.

Introduction to the Maïdo Lidar Calibration Campaign dedicated to the validation of upper air meteorological parameters

The first operations at the new High-altitude Maïdo Observatory at La Réunion began in 2013. The Maïdo Lidar Calibration Campaign (MALICCA) was organized there in April 2013 and has focused on the validation of the thermodynamic parameters (temperature, water vapor, and wind) measured with many instruments including the new very large lidar for water vapor and temperature profiles. The aim of this publication consists of providing an overview of the different instruments deployed during this campaign and their status, some of the targeted scientific questions and associated instrumental issues. Some specific detailed studies for some individual techniques were addressed elsewhere. This study shows that temperature profiles were obtained from the ground to the mesopause (80 km) thanks to the lidar and regular meteorological balloon-borne sondes with an overlap range showing good agreement. Water vapor is also monitored from the ground to the mesopause by using the Raman lidar and microwave techniques. Both techniques need to be pushed to their limit to reduce the missing range in the lower stratosphere. Total columns obtained from global positioning system or spectrometers are valuable for checking the calibration and ensuring vertical continuity. The lidar can also provide the vertical cloud structure that is a valuable complementary piece of information when investigating the water vapor cycle. Finally, wind vertical profiles, which were obtained from sondes, are now also retrieved at Maïdo from the newly implemented microwave technique and the lidar. Stable calibrations as well as a small-scale dynamical structure are required to monitor the thermodynamic state of the middle atmosphere, ensure validation of satellite sensors, study the transport of water vapor in the vicinity of the tropical tropopause and study their link with cirrus clouds and cyclones and the impact of small-scale dynamics (gravity waves) and their link with the mean state of the mesosphere.

Some physical and chemical proceddes in fluids

The ability to reproduce reduced gravity conditions for long periods is one of the reasons why the orbiting laboratory is so attractive. In this paper several fluid dynamics problem areas are reviewed in which zero-gravity conditions are of great importance. Although emphasis is placed on space processing, there are some older problems also in which gravity masks the phenomcna, impeding a reasonably simple approach to the solution. Three problems are considered: Thermal convection under reduced gravity. The dumping effect ofsurface gravity waves at the outset of convection induced by surface tractions is discussed in particular. The existence of convection is of concern for some satellite thermal control techniques presently used, and for most of the proposed manufacturing processes. Whereas convection should be normally avoided, problems related to the containerless stirring ofa melt constitute an exception. Secondly, gravity and chemical reactions. Although chemical reactions are independent of gravity because of the small mass of the molecules and atoms involved, in many cases the reaction rate dcpends on the arrival of the species to the reaction zone. When the arrival process is buoyancy-controlled, the net specd of the reaction will be affected by the gravity. Thirdly, two-phase flows under reduced gravity provkle interesting problems from boiling heat transfer to degasslng of melts. This part of the paper deals only with the measurement of sound veiocity in a liquid containing bubbles. It is suggested that such measurements should be mude under reduced gravity to provide reliable residís.

Space-time Measurements of Oceanic Sea States

Stereo video techniques are effective for estimating the space–time wave dynamics over an area of the ocean. Indeed, a stereo camera view allows retrieval of both spatial and temporal data whose statistical content is richer than that of time series data retrieved from point wave probes. We present an application of the Wave Acquisition Stereo System (WASS) for the analysis of offshore video measurements of gravity waves in the Northern Adriatic Sea and near the southern seashore of the Crimean peninsula, in the Black Sea. We use classical epipolar techniques to reconstruct the sea surface from the stereo pairs sequentially in time, viz. a sequence of spatial snapshots. We also present a variational approach that exploits the entire data image set providing a global space–time imaging of the sea surface, viz. simultaneous reconstruction of several spatial snapshots of the surface in order to guarantee continuity of the sea surface both in space and time. Analysis of the WASS measurements show that the sea surface can be accurately estimated in space and time together, yielding associated directional spectra and wave statistics at a point in time that agrees well with probabilistic models. In particular, WASS stereo imaging is able to capture typical features of the wave surface, especially the crest-to-trough asymmetry due to second order nonlinearities, and the observed shape of large waves are fairly described by theoretical models based on the theory of quasi-determinism (Boccotti, 2000). Further, we investigate space–time extremes of the observed stationary sea states, viz. the largest surface wave heights expected over a given area during the sea state duration. The WASS analysis provides the first experimental proof that a space–time extreme is generally larger than that observed in time via point measurements, in agreement with the predictions based on stochastic theories for global maxima of Gaussian fields.

AMS C-14 ages of coretops collected on RRS Charles Darwin cruise CD159, July 2004, N. Atlantic, for the NERC RAPID programme

Material and data were collected at 41 sites in the subpolar North Atlantic Ocean between Scotland and Newfoundland, during the RRS CharlesDarwin CD159 cruise in July 2004 (McCave, 2005). Sites were selected to reflect the major inputs of water that becomes the North Atlantic Deep Water (NADW); the Iceland-Scotland Overflow Water (ISOW), the Denmark Strait Overflow Water (DSOW) and the Labrador Sea Water (LSW). Areas cored were the south Iceland Rise, SE Greenland slope/rise and Eirik Drift, and the Labrador margin. A total of 29 box cores, 19 piston cores, 6 kasten cores, 9 short gravity cores and 20 CTD casts as well as 28 surface water samples were collected during the cruise. Here we present sediment core-top sample ages. The cores were sampled at 1 or 0.5 cm intervals and we used the top 1 or 2 cm, depending on availability of foraminifera in the samples. Sediment samples were disaggregated on an end-over-end wheel, wet sieved at >63 um, and dry sieved to 63-150 and >150 um. Accelerator Mass Spectrometer (AMS) radiocarbon dating was done for each core top based on between 900-1600 monospecific planktonic foraminifera (Globigerina bulloides or Neogloboquadrina pachyderma (sinistral)). All dates were of modern or late Holocene age except site RAPID-08-5B (9806 ± 38 uncorrected 14C years BP) and site RAPID-14-10B (11543 ± 40 uncorrected 14C years BP). The >150 um fraction was split until approximately 300 foraminifera remained and counted for number of lithic grains, benthic foraminifera, planktonic foraminifera and foraminifera fragments. In all but the shallowest sample (Greenland rise, 761m water depth) benthic foraminifera constituted less than 2% of the total >150 um fraction of the sample.

Análise dos microssismos registrados no arquipélago de São Pedro e São Paulo (ASPSP) e suas relações com variáveis oceanográficas

Microseisms are continuous vibrations pervasively recorded in the mili Hertz to 1 Hz frequency range. These vibrations are mostly composed of Rayleigh waves and are strongest in the 0.04 to 1 Hz frequency band. Their precise source mechanisms are still a matter of debate but it is agreed that they are related to atmospheric perturbations and ocean gravity waves. The Saint Peter Saint Paul Archipelago (SPSPA) is located in the equatorial region of the Atlantic Ocean about 1,100 km distant from the Brazilian northeastern coast. The SPSPA is composed by a set of several small rocky formations with a total area of approximately 17,000 m². Due to its remote distance from the continent and the lack of cultural noise, this location is a unique location for measuring microseismic noise and to investigate its relation with some climate and oceanographic variables. In the SPSPA we have recorded both primary microseisms (PM) at 0.04 – 0.12 Hz and the secondary microseisms (SM) at 0.12 – 0.4 Hz during 10 months in 2012 and 2013. Our analysis indicates a good correlation between the microseismic noise in the region and a seasonal dependency. In particular, the winter in the northern hemisphere. We have also shown that most of the PM is generated in the SPSPA itself. The SM source location depends with the seasonal climatic and oceanographic variables in the northern hemisphere

Estudos de distúrbios ionosféricos propagantes de larga escala (LSTID’S) na região equatorial durante grandes tempestades magnéticas: estudo de casos

The Ionospheric Disturbances – TIDs – are irregularities on the ionospheric plasma propagating in speeds in the order of tens to a few hundreds of meters per second. This present study detected and characterized the TIDs of LSTIDs (Large Scale Travelling Ionospheric Disturbance) type at low latitudes during intense geomagnetic storms and its propagation over the Brazilian sector. This work also shows as being the first to report systematically propagation of gravity waves over Natal. For this purpose, we used ionospheric records obtained from type of digisonde CADI (Canadiam Advanced Digital Ionosonde) located in Natal and the type DSP (Digisonde Portable Souder) located in Cachoeira Paulista, Fortaleza and São Luis, whereupon we used a dataset of 12 years collected by INPE (National Institute of Space Research). In this study, both calm days, that preceded the storms, and the geomagnetically disturbed days were related during the years 2000 and 1012, which cover a period of maximum and minimum solar activity. And it is presented the variations that happened in the electron density from region F of the ionosphere over the Brazilian sector, especially near the Equator (Natal, Fortaleza and São Luis), caused by ionospheric disturbances in the equatorial region during intense geomagnetic storms, because, as we know of the literature in this area, this phenomenon contributes positively to the emergence of LSTIDs in the auroral region, which may move to the equatorial region where a few cases have been documented and studied systematically. From the observation of signatures if TIDs in ionogram records, a study of the morphology of these events was performed and compared with the main characteristics of the wave of this phenomenon during great magnetic storms, i.e., DST <(-200 nT) and KP > 6. Thus, we obtained the main characteristics of TIDs over our region, i.e., period, vertical wavelength, phase and propagation speed, as well as the delay of these disturbances compared to the beginning of the magnetic storms to the Brazilian Sector.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission): an extended white paper

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high res- olution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM’s main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude bet- ter than COBE FIRAS. The data obtained will allow us to precisely measure the absolute sky brightness and polarization of all the components of the sky emission in the observed frequency range, separating the primordial and extragalactic components cleanly from the galactic and zodiacal light emissions. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM, which include: (1) the ultimate galaxy cluster survey using the Sunyaev-Zeldovich (SZ) e↵ect, detecting approximately 106 clusters extending to large redshift, including a char- acterization of the gas temperature of the brightest ones (through the relativistic corrections to the classic SZ template) as well as a peculiar velocity survey using the kinetic SZ e↵ect that comprises our entire Hubble volume; (2) a detailed characterization of the properties and evolution of dusty galaxies, where the most of the star formation in the universe took place, the faintest population of which constitute the di↵use CIB (Cosmic Infrared Background); (3) a characterization of the B modes from primordial gravity waves generated during inflation and from gravitational lensing, as well as the ultimate search for primordial non-Gaussianity using CMB polarization, which is less contaminated by foregrounds on small scales than thetemperature anisotropies; (4) a search for distortions from a perfect blackbody spectrum, which include some nearly certain signals and others that are more speculative but more informative; and (5) a study of the role of the magnetic field in star formation and its inter- action with other components of the interstellar medium of our Galaxy. These are but a few of the highlights presented here along with a description of the proposed instrument.

Accurate series solutions for gravity-driven Stokes waves

In the past, high order series expansion techniques have been used to study the nonlinear equations that govern the form of periodic Stokes waves moving steadily on the surface of an inviscid fluid. In the present study, two such series solutions are recomputed using exact arithmetic, eliminating any loss of accuracy due to accumulation of round-off error, allowing a much greater number of terms to be found with confidence. It is shown that higher order behaviour of series generated by the solution casts doubt over arguments that rely on estimating the series’ radius of convergence. Further, the exact nature of the series is used to shed light on the unusual nature of convergence of higher order Pade approximants near the highest wave. Finally, it is concluded that, provided exact values are used in the series, these Pade approximants prove very effective in successfully predicting three turning points in both the dispersion relation and the total energy.

Short-Term Nonlinear Response of Tension Leg Platform in Random Sea Waves

Most of the existing mathematical models for analyzing the dynamic response of TLP are based on explicit or implicit assumptions that motions (translations and rotations) are small magnitude. However, when TLP works in severe adverse conditions, the a priori assumption on small displacements may be inadequate. In such situation, the motions should be regarded as finite magnitude. This paper will study stochastic nonlinear dynamic responses of TLP with finite displacements in random waves. The nonlinearities considered are: large amplitude motions, coupling the six degrees-of-freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. The nonlinear dynamic responses are calculated by using numerical integration procedure in the time domain. After the time histories of the dynamic responses are obtained, we carry out cycle counting of the stress histories of the tethers with rain-flow counting method to get the stress range distribution.

Short and long waves over a muddy seabed

The available experimental results have shown that in time-periodic motion the rheology of fluid mud displays complex viscoelastic behaviour. Based on the measured rheology of fluid mud from two field sites, we study the interaction of water waves and fluid mud by a two-layered model in which the water above is assumed to be inviscid and the mud below is viscoelastic. As the fluid-mud layer in shallow seas is usually much thinner than the water layer above, the sharp contrast of scales enables an approximate analytical theory for the interaction between fluid mud and small-amplitude waves with a narrow frequency band. It is shown that at the leading order and within a short distance of a few wavelengths, wave pressure from above forces mud motion below. Over a Much longer distance, waves are modified by the accumulative dissipation in mud. At the next order, infragravity waves owing to convective inertia (or radiation stresses) are affected indirectly by mud motion through the slow modulation of the short waves. Quantitative predictions are made for mud samples of several concentrations and from two different field sites.

Mono-energetic ion beam acceleration in solitary waves during relativistic transparency using high-contrast circularly polarized short-pulse laser and nanoscale targets

In recent experiments at the Trident laser facility, quasi-monoenergetic ion beams have been obtained from the interaction of an ultraintense, circularly polarized laser with a diamond-like carbon target of nm-scale thickness under conditions of ultrahigh laser pulse contrast. Kinetic simulations of this experiment under realistic laser and plasma conditions show that relativistic transparency occurs before significant radiation pressure acceleration and that the main ion acceleration occurs after the onset of relativistic transparency. Associated with this transition are a period of intense ion acceleration and the generation of a new class of ion solitons that naturally give rise to quasi-monoenergetic ion beams. An analytic theory has been derived for the properties of these solitons that reproduces the behavior observed in kinetic simulations and the experiments. © 2011 American Institute of Physics.