Extreme deviations and applications

Stretched exponential probability density functions (pdf), having the form of the exponential of minus a fractional power of the argument, are commonly found in turbulence and other areas. They can arise because of an underlying random multiplicative process. For this, a theory of extreme deviations is developed, devoted to the far tail of the pdf of the sum $X$ of a finite number $n$ of independent random variables with a common pdf $e^{-f(x)}$. The function $f(x)$ is chosen (i) such that the pdf is normalized and (ii) with a strong convexity condition that $f''(x)>0$ and that $x^2f''(x)\to +\infty$ for $|x|\to\infty$. Additional technical conditions ensure the control of the variations of $f''(x)$. The tail behavior of the sum comes then mostly from individual variables in the sum all close to $X/n$ and the tail of the pdf is $\sim e^{-nf(X/n)}$. This theory is then applied to products of independent random variables, such that their logarithms are in the above class, yielding usually stretched exponential tails. An application to fragmentation is developed and compared to data from fault gouges. The pdf by mass is obtained as a weighted superposition of stretched exponentials, reflecting the coexistence of different fragmentation generations. For sizes near and above the peak size, the pdf is approximately log-normal, while it is a power law for the smaller fragments, with an exponent which is a decreasing function of the peak fragment size. The anomalous relaxation of glasses can also be rationalized using our result together with a simple multiplicative model of local atom configurations. Finally, we indicate the possible relevance to the distribution of small-scale velocity increments in turbulent flow.Comment: 26 pages, 1 figure ps (now available), addition and discussion of mathematical references; appeared in J. Phys. I France 7, 1155-1171 (1997

Non-unique factorization of polynomials over residue class rings of the integers

We investigate non-unique factorization of polynomials in Z_{p^n}[x] into irreducibles. As a Noetherian ring whose zero-divisors are contained in the Jacobson radical, Z_{p^n}[x] is atomic. We reduce the question of factoring arbitrary non-zero polynomials into irreducibles to the problem of factoring monic polynomials into monic irreducibles. The multiplicative monoid of monic polynomials of Z_{p^n}[x] is a direct sum of monoids corresponding to irreducible polynomials in Z_p[x], and we show that each of these monoids has infinite elasticity. Moreover, for every positive integer m, there exists in each of these monoids a product of 2 irreducibles that can also be represented as a product of m irreducibles.Comment: 11 page

Experimental evidence of accelerated energy transfer in turbulence

We investigate the vorticity dynamics in a turbulent vortex using scattering of acoustic waves. Two ultrasonic beams are adjusted to probe simultaneously two spatial scales in a given volume of the flow, thus allowing a dual channel recording of the dynamics of coherent vorticity structures. Our results show that this allows to measure the average energy transfer time between different spatial length scales, and that such transfer goes faster at smaller scales.Comment: 5 pages, 5 figure

An order (n) algorithm for the dynamics simulation of robotic systems

The formulation of an Order (n) algorithm for DISCOS (Dynamics Interaction Simulation of Controls and Structures), which is an industry-standard software package for simulation and analysis of flexible multibody systems is presented. For systems involving many bodies, the new Order (n) version of DISCOS is much faster than the current version. Results of the experimental validation of the dynamics software are also presented. The experiment is carried out on a seven-joint robot arm at NASA's Goddard Space Flight Center. The algorithm used in the current version of DISCOS requires the inverse of a matrix whose dimension is equal to the number of constraints in the system. Generally, the number of constraints in a system is roughly proportional to the number of bodies in the system, and matrix inversion requires O(p exp 3) operations, where p is the dimension of the matrix. The current version of DISCOS is therefore considered an Order (n exp 3) algorithm. In contrast, the Order (n) algorithm requires inversion of matrices which are small, and the number of matrices to be inverted increases only linearly with the number of bodies. The newly-developed Order (n) DISCOS is currently capable of handling chain and tree topologies as well as multiple closed loops. Continuing development will extend the capability of the software to deal with typical robotics applications such as put-and-place, multi-arm hand-off and surface sliding

Universal decay of scalar turbulence

The asymptotic decay of passive scalar fields is solved analytically for the Kraichnan model, where the velocity has a short correlation time. At long times, two universality classes are found, both characterized by a distribution of the scalar -- generally non-Gaussian -- with global self-similar evolution in time. Analogous behavior is found numerically with a more realistic flow resulting from an inverse energy cascade.Comment: 4 pages, 3 Postscript figures, submitted to PR

Effect of helicity and rotation on the free decay of turbulent flows

The self-similar decay of energy in a turbulent flow is studied in direct numerical simulations with and without rotation. Two initial conditions are considered: one non-helical (mirror-symmetric), and one with maximal helicity. The results show that, while in the absence of rotation the energy in the helical and non-helical cases decays with the same rate, in rotating flows the helicity content has a major impact on the decay rate. These differences are associated with differences in the energy and helicity cascades when rotation is present. Properties of the structures that arise in the flow at late times in each time are also discussed.Comment: 4 pages, 4 figure

Entire solutions of hydrodynamical equations with exponential dissipation

We consider a modification of the three-dimensional Navier--Stokes equations and other hydrodynamical evolution equations with space-periodic initial conditions in which the usual Laplacian of the dissipation operator is replaced by an operator whose Fourier symbol grows exponentially as \ue ^{|k|/\kd} at high wavenumbers $|k|$. Using estimates in suitable classes of analytic functions, we show that the solutions with initially finite energy become immediately entire in the space variables and that the Fourier coefficients decay faster than \ue ^{-C(k/\kd) \ln (|k|/\kd)} for any $C<1/(2\ln 2)$. The same result holds for the one-dimensional Burgers equation with exponential dissipation but can be improved: heuristic arguments and very precise simulations, analyzed by the method of asymptotic extrapolation of van der Hoeven, indicate that the leading-order asymptotics is precisely of the above form with $C= C_\star =1/\ln2$. The same behavior with a universal constant $C_\star$ is conjectured for the Navier--Stokes equations with exponential dissipation in any space dimension. This universality prevents the strong growth of intermittency in the far dissipation range which is obtained for ordinary Navier--Stokes turbulence. Possible applications to improved spectral simulations are briefly discussed.Comment: 29 pages, 3 figures, Comm. Math. Phys., in pres

Helical rotating turbulence. Part II. Intermittency, scale invariance and structures

We study the intermittency properties of the energy and helicity cascades in two 1536^3 direct numerical simulations of helical rotating turbulence. Symmetric and anti-symmetric velocity increments are examined, as well as probability density functions of the velocity field and of the helicity density. It is found that the direct cascade of energy to small scales is scale invariant and non-intermittent, whereas the direct cascade of helicity is highly intermittent. Furthermore, the study of structure functions of different orders allows us to identify a recovery of isotropy of strong events at very small scales in the flow. Finally, we observe the juxtaposition in space of strong laminar and persistent helical columns next to time-varying vortex tangles, the former being associated with the self-similarity of energy and the latter with the intermittency of helicity.Comment: 11 pages, 10 figure