Scaling and the prediction of energy spectra in decaying hydrodynamic turbulence

Few rigorous results are derived for fully developed turbulence. By applying the scaling properties of the Navier-Stokes equation we have derived a relation for the energy spectrum valid for unforced or decaying isotropic turbulence. We find the existence of a scaling function $\psi$. The energy spectrum can at any time by a suitable rescaling be mapped onto this function. This indicates that the initial (primordial) energy spectrum is in principle retained in the energy spectrum observed at any later time, and the principle of permanence of large eddies is derived. The result can be seen as a restoration of the determinism of the Navier-Stokes equation in the mean. We compare our results with a windtunnel experiment and find good agreement.Comment: 4 pages, 1 figur

The two-dimensional magnetohydrodynamic Kelvin-Helmholtz instability: compressibility and large-scale coalescence effects

The Kelvin-Helmholtz (KH) instability occurring in a single shear flow configuration that is embedded in a uniform flow-aligned magnetic field, is revisited by means of high resolution two-dimensional (2D) magnetohydrodynamic (MHD) simulations. First, the calculations extend previous studies of magnetized shear flows to a higher compressibility regime. The nonlinear evolution of an isolated KH billow emerging from the fastest growing linear mode for a convective sonic Mach number $M_{cs}=0.7$ layer is in many respects similar to its less compressible counterpart (Mach $M_{cs}=0.5$). In particular, the disruptive regime where locally amplified, initially weak magnetic fields, control the nonlinear saturation process is found for Alfv\'en Mach numbers 4\simlt M_A \simlt 30. The most notable difference between $M_{cs}=0.7$ versus $M_{cs}=0.5$ layers is that higher density contrasts and fast magnetosonic shocklet structures are observed. Second, the use of adaptive mesh refinement allows to parametrically explore much larger computational domains, including up to 22 wavelengths of the linearly dominant mode. A strong process of large-scale coalescence is found, whatever the magnetic field regime. It proceeds through continuous pairing/merging events between adjacent vortices up to the point where the final large-scale vortical structure reaches the domain dimensions. This pairing/merging process is attributed to the growth of subharmonic modes and is mainly controlled by relative phase differences between them. These grid-adaptive simulations demonstrate that even in very weak magnetic field regimes ($M_A \simeq 30$), the large-scale KH coalescence process can trigger tearing-type reconnection events previously identified in cospatial current-vortex sheets.Comment: Published in Physics of Plasmas, figures absent due to file sizes, full version at http://www.phys.uu.nl/~toth/ (follow Publications

On small time asymptotics for rough differential equations driven by fractional Brownian motions

We survey existing results concerning the study in small times of the density of the solution of a rough differential equation driven by fractional Brownian motions. We also slightly improve existing results and discuss some possible applications to mathematical finance.Comment: This is a survey paper, submitted to proceedings in the memory of Peter Laurenc

Geodetic, teleseismic, and strong motion constraints on slip from recent southern Peru subduction zone earthquakes

We use seismic and geodetic data both jointly and separately to constrain coseismic slip from the 12 November 1996 M_w 7.7 and 23 June 2001 M_w 8.5 southern Peru subduction zone earthquakes, as well as two large aftershocks following the 2001 earthquake on 26 June and 7 July 2001. We use all available data in our inversions: GPS, interferometric synthetic aperture radar (InSAR) from the ERS-1, ERS-2, JERS, and RADARSAT-1 satellites, and seismic data from teleseismic and strong motion stations. Our two-dimensional slip models derived from only teleseismic body waves from South American subduction zone earthquakes with M_w > 7.5 do not reliably predict available geodetic data. In particular, we find significant differences in the distribution of slip for the 2001 earthquake from models that use only seismic (teleseismic and two strong motion stations) or geodetic (InSAR and GPS) data. The differences might be related to postseismic deformation or, more likely, the different sensitivities of the teleseismic and geodetic data to coseismic rupture properties. The earthquakes studied here follow the pattern of earthquake directivity along the coast of western South America, north of 5°S, earthquakes rupture to the north; south of about 12°S, directivity is southerly; and in between, earthquakes are bilateral. The predicted deformation at the Arequipa GPS station from the seismic-only slip model for the 7 July 2001 aftershock is not consistent with significant preseismic motion

Adaptive density estimation for stationary processes

We propose an algorithm to estimate the common density $s$ of a stationary process $X_1,...,X_n$. We suppose that the process is either $\beta$ or $\tau$-mixing. We provide a model selection procedure based on a generalization of Mallows' $C_p$ and we prove oracle inequalities for the selected estimator under a few prior assumptions on the collection of models and on the mixing coefficients. We prove that our estimator is adaptive over a class of Besov spaces, namely, we prove that it achieves the same rates of convergence as in the i.i.d framework

Adaptive estimation in circular functional linear models

We consider the problem of estimating the slope parameter in circular functional linear regression, where scalar responses Y1,...,Yn are modeled in dependence of 1-periodic, second order stationary random functions X1,...,Xn. We consider an orthogonal series estimator of the slope function, by replacing the first m theoretical coefficients of its development in the trigonometric basis by adequate estimators. Wepropose a model selection procedure for m in a set of admissible values, by defining a contrast function minimized by our estimator and a theoretical penalty function; this first step assumes the degree of ill posedness to be known. Then we generalize the procedure to a random set of admissible m's and a random penalty function. The resulting estimator is completely data driven and reaches automatically what is known to be the optimal minimax rate of convergence, in term of a general weighted L2-risk. This means that we provide adaptive estimators of both the slope function and its derivatives

The Crab pulsar light curve in the soft gamma ray range: FIGARO II results

The FIGARO II experiment (a large area, balloon borne, crystal scintillator detector working from 0.15 to 4.3 MeV) observed the Crab pulsar on 1990 Jul. 9 for about seven hours. The study of the pulse profile confirms some structures detected with a low significance during the shorter observation of 1986, and adds new important elements to the picture. In particular, between the two main peaks, two secondary peaks appear centered at phase values 0.1 and 0.3, in the energy range 0.38 to 0.49 MeV; in the same energy range, a spectral feature at 0.44 MeV, interpreted as a redshifted positron annihilation line, was observed during the same balloon flight in the phase interval including the second main peak and the neighboring secondary peak. If the phase interval considered is extended to include also the other secondary peak, the significance of the spectral line appears to increase