The Cauchy-Lagrangian method for numerical analysis of Euler flow

A novel semi-Lagrangian method is introduced to solve numerically the Euler equation for ideal incompressible flow in arbitrary space dimension. It exploits the time-analyticity of fluid particle trajectories and requires, in principle, only limited spatial smoothness of the initial data. Efficient generation of high-order time-Taylor coefficients is made possible by a recurrence relation that follows from the Cauchy invariants formulation of the Euler equation (Zheligovsky & Frisch, J. Fluid Mech. 2014, 749, 404-430). Truncated time-Taylor series of very high order allow the use of time steps vastly exceeding the Courant-Friedrichs-Lewy limit, without compromising the accuracy of the solution. Tests performed on the two-dimensional Euler equation indicate that the Cauchy-Lagrangian method is more - and occasionally much more - efficient and less prone to instability than Eulerian Runge-Kutta methods, and less prone to rapid growth of rounding errors than the high-order Eulerian time-Taylor algorithm. We also develop tools of analysis adapted to the Cauchy-Lagrangian method, such as the monitoring of the radius of convergence of the time-Taylor series. Certain other fluid equations can be handled similarly.Comment: 30 pp., 13 figures, 45 references. Minor revision. In press in Journal of Scientific Computin

Photoassociation of a cold atom-molecule pair: long-range quadrupole-quadrupole interactions

The general formalism of the multipolar expansion of electrostatic interactions is applied to the calculation the potential energy between an excited atom (without fine structure) and a ground state diatomic molecule at large separations. Both partners exhibit a permanent quadrupole moment, so that their mutual quadrupole-quadrupole long-range interaction is attractive enough to bind trimers. Numerical results are given for an excited Cs(6P) atom and a ground state Cs2 molecule. The prospects for achieving photoassociation of a cold atom/dimer pair is thus discussed and found promising. The formalism can be easily generalized to the long-range interaction between molecules to investigate the formation of cold tetramers.Comment: 5 figure

Dynamo effect in parity-invariant flow with large and moderate separation of scales

It is shown that non-helical (more precisely, parity-invariant) flows capable of sustaining a large-scale dynamo by the negative magnetic eddy diffusivity effect are quite common. This conclusion is based on numerical examination of a large number of randomly selected flows. Few outliers with strongly negative eddy diffusivities are also found, and they are interpreted in terms of the closeness of the control parameter to a critical value for generation of a small-scale magnetic field. Furthermore, it is shown that, for parity-invariant flows, a moderate separation of scales between the basic flow and the magnetic field often significantly reduces the critical magnetic Reynolds number for the onset of dynamo action.Comment: 44 pages,11 figures, significantly revised versio

The Monge-Ampere equation: various forms and numerical methods

We present three novel forms of the Monge-Ampere equation, which is used, e.g., in image processing and in reconstruction of mass transportation in the primordial Universe. The central role in this paper is played by our Fourier integral form, for which we establish positivity and sharp bound properties of the kernels. This is the basis for the development of a new method for solving numerically the space-periodic Monge-Ampere problem in an odd-dimensional space. Convergence is illustrated for a test problem of cosmological type, in which a Gaussian distribution of matter is assumed in each localised object, and the right-hand side of the Monge-Ampere equation is a sum of such distributions.Comment: 24 pages, 2 tables, 5 figures, 32 references. Submitted to J. Computational Physics. Times of runs added, multiple improvements of the manuscript implemented

Persistence of small-scale anisotropy of magnetic turbulence as observed in the solar wind

The anisotropy of magnetophydrodynamic turbulence is investigated by using solar wind data from the Helios 2 spacecraft. We investigate the behaviour of the complete high-order moment tensors of magnetic field increments and we compare the usual longitudinal structure functions which have both isotropic and anisotropic contributions, to the fully anisotropic contribution. Scaling exponents have been extracted by an interpolation scaling function. Unlike the usual turbulence in fluid flows, small-scale magnetic fluctuations remain anisotropic. We discuss the radial dependence of both anisotropy and intermittency and their relationship.Comment: 7 pages, 2 figures, in press on Europhys. Let

Hybridization-related correction to the jellium model for fullerenes

We introduce a new type of correction for a more accurate description of fullerenes within the spherically symmetric jellium model. This correction represents a pseudopotential which originates from the comparison between an accurate ab initio calculation and the jellium model calculation. It is shown that such a correction to the jellium model allows one to account, at least partly, for the sp2-hybridization of carbon atomic orbitals. Therefore, it may be considered as a more physically meaningful correction as compared with a structureless square-well pseudopotential which has been widely used earlier.Comment: 16 pages, 10 figure

Passive scalar turbulence in high dimensions

Exploiting a Lagrangian strategy we present a numerical study for both perturbative and nonperturbative regions of the Kraichnan advection model. The major result is the numerical assessment of the first-order $1/d$-expansion by M. Chertkov, G. Falkovich, I. Kolokolov and V. Lebedev ({\it Phys. Rev. E}, {\bf 52}, 4924 (1995)) for the fourth-order scalar structure function in the limit of high dimensions $d$'s. %Two values of the velocity scaling exponent $\xi$ have been considered: %$\xi=0.8$ and $\xi=0.6$. In the first case, the perturbative regime %takes place at $d\sim 30$, while in the second at $d\sim 25$, %in agreement with the fact that the relevant small parameter %of the theory is $\propto 1/(d (2-\xi))$. In addition to the perturbative results, the behavior of the anomaly for the sixth-order structure functions {\it vs} the velocity scaling exponent, $\xi$, is investigated and the resulting behavior discussed.Comment: 4 pages, Latex, 4 figure

A note on the extension of the polar decomposition for the multidimensional Burgers equation

It is shown that the generalizations to more than one space dimension of the pole decomposition for the Burgers equation with finite viscosity and no force are of the form u = -2 viscosity grad log P, where the P's are explicitly known algebraic (or trigonometric) polynomials in the space variables with polynomial (or exponential) dependence on time. Such solutions have polar singularities on complex algebraic varieties.Comment: 3 pages; minor formatting and typos corrected. Submitted to Phys. Rev. E (Rapid Comm.

Statistical conservation laws in turbulent transport

We address the statistical theory of fields that are transported by a turbulent velocity field, both in forced and in unforced (decaying) experiments. We propose that with very few provisos on the transporting velocity field, correlation functions of the transported field in the forced case are dominated by statistically preserved structures. In decaying experiments (without forcing the transported fields) we identify infinitely many statistical constants of the motion, which are obtained by projecting the decaying correlation functions on the statistically preserved functions. We exemplify these ideas and provide numerical evidence using a simple model of turbulent transport. This example is chosen for its lack of Lagrangian structure, to stress the generality of the ideas