Lagrangian Dynamics in High-Dimensional Point-Vortex Systems

We study the Lagrangian dynamics of systems of N point vortices and passive particles in a two-dimensional, doubly periodic domain. The probability distribution function of vortex velocity, p_N, has a slow-velocity Gaussian component and a significant high-velocity tail caused by close vortex pairs. In the limit for N -> oo, p_N tends to a Gaussian. However, the form of the single-vortex velocity causes very slow convergence with N; for N ~ 10^6 the non-Gaussian high-velocity tails still play a significant role. At finite N, the Gaussian component is well modeled by an Ornstein-Uhlenbeck (OU) stochastic process with variance sigma_N = sqrt{N ln N /2 pi}. Considering in detail the case N=100, we show that at short times the velocity autocorrelation is dominated by the Gaussian component and displays an exponential decay with a short Lagrangian decorrelation time. The close pairs have a long correlation time and cause nonergodicity over at least the time of the integration. Due to close vortex dipoles the absolute dispersion differs significantly from the OU prediction, and shows evidence of long-time anomalous dispersion. We discuss the mathematical form of a new stochastic model for the Lagrangian dynamics, consisting of an OU model combined with long-lived close same-sign vortices engaged in rapid rotation and long-lived close dipoles engaged in ballistic motion. From a dynamical-systems perspective this work indicates that systems of dimension O(100) can have behavior which is a combination of both low-dimensional behavior, i.e. close pairs, and extremely high-dimensional behavior described by traditional stochastic processes.Comment: 48 pages, 17 figures, to appear in Physics of Fluid

Ideabook: Libraries for Families

The IDEABOOK is a research-based framework to guide and broaden family engagement in libraries.The framework helps libraries move beyond thinking of family engagement as random, individual activities or programs, but rather as a system where library leadership, activities, and resources that are linked to goals. The framework represents a theory of change that begins with a set of elements—leadership, engagement, and support services—that build a pathway for meaningful family engagement beginning in the early childhood years and extending through young adulthood.This IDEABOOK was developed for anyone who works in a library setting—from library directors and children's and youth librarians, to volunteers and support staff—and shares many innovative ways that libraries support and guide families in children's learning and development

A systematic correlation between two-dimensional flow topology and the abstract statistics of turbulence

Velocity differences in the direct enstrophy cascade of two-dimensional turbulence are correlated with the underlying flow topology. The statistics of the transverse and longitudinal velocity differences are found to be governed by different structures. The wings of the transverse distribution are dominated by strong vortex centers, whereas, the tails of the longitudinal differences are dominated by saddles. Viewed in the framework of earlier theoretical work this result suggests that the transfer of enstrophy to smaller scales is accomplished in regions of the flow dominated by saddles.Comment: 4 pages, 4 figure

Influence of bottom topography on integral constraints in zonal flows with parameterized potential vorticity fluxes

An integral constraint for eddy fluxes of potential vorticity (PV), corresponding to global momentum conservation, is applied to two-layer zonal quasi-geostrophic channel flow. This constraint must be satisfied for any type of parameterization of eddy PV fluxes. Bottom topography strongly influence the integral constraint compared to a flat bottom channel. An analytical solution for the mean flow solution has been found by using asymptotic expansion in a small parameter which is the ratio of the Rossby radius to the meridional extent of the channel. Applying the integral constraint to this solution, one can find restrictions for eddy PV transfer coefficients which relate the eddy fluxes of PV to the mean flow. These restrictions strongly deviate from restrictions for the channel with flat bottom topography

GiRaFFE: An Open-Source General Relativistic Force-Free Electrodynamics Code

We present GiRaFFE, the first open-source general relativistic force-free electrodynamics (GRFFE) code for dynamical, numerical-relativity generated spacetimes. GiRaFFE adopts the strategy pioneered by McKinney and modified by Paschalidis and Shapiro to convert a GR magnetohydrodynamic (GRMHD) code into a GRFFE code. In short, GiRaFFE exists as a modification of IllinoisGRMHD, a user-friendly, open-source, dynamical-spacetime GRMHD code. Both GiRaFFE and IllinoisGRMHD leverage the Einstein Toolkit's highly-scalable infrastructure to make possible large-scale simulations of magnetized plasmas in strong, dynamical spacetimes on adaptive-mesh refinement (AMR) grids. We demonstrate that GiRaFFE passes a large suite of both flat and curved-spacetime code tests passed by a number of other state-of-the-art GRFFE codes, and is thus ready for production-scale simulations of GRFFE phenomena of key interest to relativistic astrophysics.Comment: 23 pages, 4 figures. Consistent with published versio

Random Projections For Large-Scale Regression

Fitting linear regression models can be computationally very expensive in large-scale data analysis tasks if the sample size and the number of variables are very large. Random projections are extensively used as a dimension reduction tool in machine learning and statistics. We discuss the applications of random projections in linear regression problems, developed to decrease computational costs, and give an overview of the theoretical guarantees of the generalization error. It can be shown that the combination of random projections with least squares regression leads to similar recovery as ridge regression and principal component regression. We also discuss possible improvements when averaging over multiple random projections, an approach that lends itself easily to parallel implementation.Comment: 13 pages, 3 Figure

Binary black hole merger in the extreme mass ratio limit

We discuss the transition from quasi-circular inspiral to plunge of a system of two nonrotating black holes of masses $m_1$ and $m_2$ in the extreme mass ratio limit $m_1m_2\ll (m_1+m_2)^2$. In the spirit of the Effective One Body (EOB) approach to the general relativistic dynamics of binary systems, the dynamics of the two black hole system is represented in terms of an effective particle of mass $\mu\equiv m_1m_2/(m_1+m_2)$ moving in a (quasi-)Schwarzschild background of mass $M\equiv m_1+m_2$ and submitted to an ${\cal O}(\mu)$ radiation reaction force defined by Pad\'e resumming high-order Post-Newtonian results. We then complete this approach by numerically computing, \`a la Regge-Wheeler-Zerilli, the gravitational radiation emitted by such a particle. Several tests of the numerical procedure are presented. We focus on gravitational waveforms and the related energy and angular momentum losses. We view this work as a contribution to the matching between analytical and numerical methods within an EOB-type framework.Comment: 14 pages, six figures. Revised version. To appear in the CQG special issue based around New Frontiers in Numerical Relativity conference, Golm (Germany), July 17-21 200

Forced Stratified Turbulence: Successive Transitions with Reynolds Number

Numerical simulations are made for forced turbulence at a sequence of increasing values of Reynolds number, R, keeping fixed a strongly stable, volume-mean density stratification. At smaller values of R, the turbulent velocity is mainly horizontal, and the momentum balance is approximately cyclostrophic and hydrostatic. This is a regime dominated by so-called pancake vortices, with only a weak excitation of internal gravity waves and large values of the local Richardson number, Ri, everywhere. At higher values of R there are successive transitions to (a) overturning motions with local reversals in the density stratification and small or negative values of Ri; (b) growth of a horizontally uniform vertical shear flow component; and (c) growth of a large-scale vertical flow component. Throughout these transitions, pancake vortices continue to dominate the large-scale part of the turbulence, and the gravity wave component remains weak except at small scales.Comment: 8 pages, 5 figures (submitted to Phys. Rev. E