Analysis of Spatially-Coupled Counter Braids

A counter braid (CB) is a novel counter architecture introduced by Lu et al. in 2007 for per-flow measurements on high-speed links. CBs achieve an asymptotic compression rate (under optimal decoding) that matches the entropy lower bound of the flow size distribution. Spatially-coupled CBs (SC-CBs) have recently been proposed. In this work, we further analyze single-layer CBs and SC-CBs using an equivalent bipartite graph representation of CBs. On this equivalent representation, we show that the potential and area thresholds are equal. We also show that the area under the extended belief propagation (BP) extrinsic information transfer curve (defined for the equivalent graph), computed for the expected residual CB graph when a peeling decoder equivalent to the BP decoder stops, is equal to zero precisely at the area threshold. This, combined with simulations and an asymptotic analysis of the Maxwell decoder, leads to the conjecture that the area threshold is in fact equal to the Maxwell decoding threshold and hence a lower bound on the maximum a posteriori (MAP) decoding threshold. Finally, we present some numerical results and give some insight into the apparent gap of the BP decoding threshold of SC-CBs to the conjectured lower bound on the MAP decoding threshold.Comment: To appear in the IEEE Information Theory Workshop, Jeju Island, Korea, October 201

The baroclinic secondary instability of the two-dimensional shear layer

The focus of this study is on the numerical investigation of two-dimensional, isovolume, high Reynolds and Froude numbers, variable-density mixing layers. Lagrangian simulations, of both the temporal and the spatial models, are performed. They reveal the breaking-up of the strained vorticity and density-gradient braids, connecting two neighboring primary structures. The secondary instability arises where the vorticity has been intensified by the baroclinic torque. A simplified model of the braid of the variable-density mixing layer, consisting of a strained vorticity and density-gradient filament, is analyzed. It is concluded that the physical mechanism responsible for the secondary instability is the forcing of the vorticity field by the baroclinic torque, itself sensitive to perturbations. This mechanism suggests a rapid route to turbulence for the variable-density mixing layer

Sensitivity of mixing layers to three-dimensional forcing

It is well known that turbulent mixing layers are dominated by large scale, fairly coherent structures, and that these structures are related to the stability characteristics of the flow. These facts have led researchers to attempt controlling such flows by selectively forcing certain unstable modes, which can in addition have the effect of suppressing other modes. Much of the work on controlling the mixing layer has relied on forcing 2-D instabilities. The results of forcing 3-D instabilities are addressed. The objectives of the work are twofold: to understand how a mixing layer responds to 3-D perturbations, and to test the validity of an amplitude expansion in predicting the mixing layer development. The amplitude expansion could be very useful in understanding and predicting the 3-D response of the flow to a variety of initial conditions

Asymptotic Analysis and Spatial Coupling of Counter Braids

A counter braid (CB) is a novel counter architecture introduced by Lu et al. in 2007 for per-flow measurements on high-speed links which can be decoded with low complexity using message passing (MP). CBs achieve an asymptotic compression rate (under optimal decoding) that matches the entropy lower bound of the flow size distribution. In this paper, we apply the concept of spatial coupling to CBs to improve the performance of the original CBs and analyze the performance of the resulting spatially-coupled CBs (SC-CBs). We introduce an equivalent bipartite graph representation of CBs with identical iteration- by-iteration finite-length and asymptotic performance. Based on this equivalent representation, we then analyze the asymptotic performance of single-layer CBs and SC-CBs under the MP decoding algorithm proposed by Lu et al.. In particular, we derive the potential threshold of the uncoupled system and show that it is equal to the area threshold. We also derive the Maxwell decoder for CBs and prove that the potential threshold is an upper bound on the Maxwell decoding threshold, which, in turn, is a lower bound on the maximum a posteriori (MAP) decoding threshold. We then show that the area under the extended MP extrinsic information transfer curve (defined for the equivalent graph), computed for the expected residual CB graph when a peeling decoder equivalent to the MP decoder stops, is equal to zero precisely at the area threshold. This, combined with the analysis of the Maxwell decoder and simulation results, leads us to the conjecture that the potential threshold is in fact equal to the Maxwell decoding threshold and hence a lower bound on the MAP decoding threshold. Interestingly, SC-CBs do not show the well-known phenomenon of threshold saturation of the MP decoding threshold to the potential threshold characteristic of spatially-coupled low-density parity-check codes and other coupled systems. However, SC-CBs yield better MP decoding thresholds than their uncoupled counterparts. Finally, we also consider SC-CBs as a compressed sensing scheme and show that low undersampling factors can be achieved

Contour Dynamics Methods

In an early paper on the stability of fluid layers with uniform vorticity Rayleigh (1880) remarks: "... In such cases, the velocity curve is composed of portions of straight lines which meet each other at finite angles. This state of things may be supposed to be slightly disturbed by bending the surfaces of transition, and the determination of the subsequent motion depends upon that of the form of these surfaces. For co retains its constant value throughout each layer unchanged in the absence of friction, and by a well-known theorem the whole motion depends upon [omega]." We can now recognize this as essentially the principal of contour dynamics (CD), where [omega] is the uniform vorticity. The theorem referred to is the Biot-Savart law. Nearly a century later Zabusky et al (1979) presented numerical CD calculations of nonlinear vortex patch evolution. Subsequently, owing to its compact form conferring a deceptive simplicity, CD has become a widely used method for the investigation of two-dimensional rotational flow of an incompressible inviscid fluid. The aim of this article is to survey the development, technical details, and vortex-dynamic applications of the CD method in an effort to assess its impact on our understanding of the mechanics of rotational flow in two dimensions at ultrahigh Reynolds numbers. The study of the dynamics of two- and three-dimensional vortex mechanics by computational methods has been an active research area for more than two decades. Quite apart from many practical applications in the aerodynamics of separated flows, the theoretical and numerical study of vortices in incompressible fluids has been stimulated by the idea that turbulent fluid motion may be viewed as comprising ensembles of more or less coherent laminar vortex structures that interact via relatively simple dynamics and by the appeal of the vorticity equation, which does not contain the fluid pressure. Two-dimensional vortex interactions have been perceived as supposedly relevant to the origins of coherent structures observed experimentally in mixing layers, jets, and wakes, and for models of large-scale atmospheric and oceanic turbulence. Interest has often focused on the limit of infinite Reynolds number, where in the absence of boundaries, the inviscid Euler equations are assumed to properly describe the flow dynamics. The numerous surveys of progress in the study of vorticity and the use of numerical methods applied to vortex mechanics include articles by Saffman & Baker (1979) and Saffman (1981) on inviscid vortex interactions and Aref (1983) on two-dimensional flows. Numerical methods have been surveyed by Chorin (1980), and Leonard (1980, 1985). Caflisch (1988) describes work on the mathematical aspects of the subject. Zabusky (1981), Aref (1983), and Melander et al (1987b) discuss various aspects of CD. The review of Dritschel (1989) gives emphasis to numerical issues in CD and to recent computations with contour surgery. This article is confined to a discussion of vortices on a two-dimensional surface. We generally follow Saffman & Baker (1979) in matters of definition. In two dimensions a vortex sheet is a line of discontinuity in velocity while a vortex jump is a line of discontinuity in vorticity. We shall, however, use filament to denote a two-dimensional ribbon of vorticity surrounded by fluid with vorticity of different magnitude (which may be zero), rather than the more usual three-dimensional idea of a vortex tube. The ambiguity is unfortunate but is already in the literature. Additionally, a vortex patch is a finite, singly connected area of uniform vorticity while a vortex strip is an infinite strip of uniform vorticity with finite thickness, or equivalently, an infinite filament. Contour Dynamics will refer to the numerical solution of initial value problems for piecewise constant vorticity distributions by the Lagrangian method of calculating the evolution of the vorticity jumps. Such flows are often related to corresponding solutions of the Euler equations that are steady in some translating or rotating frame of reference. These solutions will be called vortex equilibria, and the numerical technique for computing their shapes based on CD is often referred to as contour statics. The mathematical foundation for the study of vorticity was laid primarily by the well-known investigations of Helmholtz, Kelvin, J. J. Thomson, Love, and others. In our century, efforts to produce numerical simulations of flows governed by the Euler equations have utilized a variety of Eulerian, Lagrangian, and hybrid methods. Among the former are the class of spectral methods that now comprise the prevailing tool for large-scale two- and three-dimensional calculations (see Hussaini & Zang 1987). The Lagrangian methods for two-dimensional flows have been predominantly vortex tracking techniques based on the Helmholtz vorticity laws. The first initial value calculations were those of Rosenhead (193l) and Westwater (1935) who attempted to calculate vortex sheet evolution by the motion of O(10) point vortices. Subsequent efforts by Moore (1974) (see also Moore 1983, 1985) and others to produce more refined computations for vortex sheets have failed for reasons related to the tendency for initially smooth vortex sheet data to produce singularities (Moore 1979). Discrete vortex methods used to study the nonlinear dynamics of vortex patches and layers have included the evolution of assemblies of point vortices by direct summation (e.g. Acton 1976) and the cloud in cell method (Roberts & Christiansen 1972, Christiansen & Zabusky 1973, Aref & Siggia 1980, 1981). For reviews see Leonard (1980) and Aref (1983). These techniques have often been criticized for their lack of accuracy and numerical convergence and because they may be subject to grid scale dispersion. However, many qualitative vortex phenomena observed in nature and in experiments, such as amalgamation events and others still under active investigation (e.g. filamentation) were first simulated numerically with discrete vortices. The contour dynamics approach is attractive because it appears to allow direct access, at least for small times, to the inviscid dynamics for vorticity distributions smoother than those of either point vortices or vortex sheets, while at the same time enabling the mapping of the two-dimensional Euler equations to a one-dimensional Lagrangian form. In Section 2 we discuss the formulation and numerical implementation of contour dynamics for the Euler equations in two dimensions. Section 3 is concerned with applications to isolated and multiple vortex systems and to vortex layers. An attempt is made to relate this work to calculations of the relevant vortex equilibria and to results obtained with other methods. Axisymmetric contour dynamics and the treatment of the multi-layer model of quasigeostrophic flows are described in Section 4 while Section 5 is devoted to a discussion of the tendency shown by vorticity jumps to undergo the strange and subtle phenomenon of filamentation