Aeolian transport of sand

The airborne transport of particles on a granular surface by the saltation mechanism is studied through numerical simulation of particles dragged by turbulent air flow. We calculate the saturated flux q_s and show that its dependence on the wind strength u_* is consistent with several empirical relations obtained from experimental measurements. We propose and explain a new relation for fluxes close to the threshold velocity u_t, namely, q_s=a(u_*-u_t)^α with α≈2. We also obtain the distortion of the velocity profile of the wind due to the drag of the particles and find a novel dynamical scaling relation. We also obtain a new expression for the dependence of the height of the saltation layer as function of the strength of the wind.     

Scaling laws in aeolian sand transport

We report on wind tunnel measurements on saltating particles in a turbulent boundary layer and provide evidence that over an erodible bed the particle velocity in the saltation layer and the saltation length are almost invariant with the wind strength, whereas over a nonerodible bed these quantities vary significantly with the air friction speed. It results that the particle transport rate over an erodible bed does not exhibit a cubic dependence with the air friction speed, as predicted by Bagnold, but a quadratic one. This contrasts with saltation over a nonerodible bed where the cubic Bagnold scaling holds. Our findings emphasize the crucial role of the boundary conditions at the bed and may have important practical consequences for aeolian sand transport in a natural environment.     

Jump at the onset of saltation

We reveal that the transition in the saturated flux for aeolian saltation is generically discontinuous by explicitly simulating particle motion in turbulent flow. This is the first time that a jump in the saturated flux has been observed. The discontinuity is followed by a coexistence interval with two metastable solutions. The modification of the wind profile due to momentum exchange exhibits a maximum at high shear strength.     

An Analytical Approach to the Turbulence-Relaxed State in Tokamak Plasmas

Starting from the continuity, temperature, and motion equations of the trapped electron fluid in generaltokamak magnetic field with positive or reversed shear and the definition of Lagrangian invariant, dL / dt = ( t u. )L =0, where u is convective velocity, the trapped electron dynamics is considered in the following two assumptions: (i) theturbulence is low frequency electrostatic, and (ii) L is a functional only of the density n, temperature T, and magneticfield B, and the effect of perturbation potential φ is included in the convective velocity u, i.e., u is a functional of n,T, B, and φ. The Lagrangian invariant hidden in the trapped electron dynamics is strictly found: L= ln[(n/B)c1(T/B2/3)c2], where c1 and c2 are dimensionless changeable parameters and c1 ∝ c2. From this Lagrangian invariant thewhich, in the limit of large aspect ratio, reduce to n(r)q(r) = const. and T3/2(r)q(r) = const., respectively. The lattertwo scaling laws are compared with existent experimental results, being in good agreement.     

Correlation time and diffusion coefficient imaging: application to a granular flow system

A parametric method for spatially resolved measurements for velocity autocorrelation functions, R(u)(tau) = , expressed as a sum of exponentials, is presented. The method is applied to a granular flow system of 2-mm oil-filled spheres rotated in a half-filled horizontal cylinder, which is an Ornstein-Uhlenbeck process with velocity autocorrelation function R(u)(tau) = e(- ||tau ||/tau(c)), where tau(c) is the correlation time and D = tau(c) is the diffusion coefficient. The pulsed-field-gradient NMR method consists of applying three different gradient pulse sequences of varying motion sensitivity to distinguish the range of correlation times present for particle motion. Time-dependent apparent diffusion coefficients are measured for these three sequences and tau(c) and D are then calculated from the apparent diffusion coefficient images. For the cylinder rotation rate of 2.3 rad/s, the axial diffusion coefficient at the top center of the free surface was 5.5 x 10(-6) m(2)/s, the correlation time was 3 ms, and the velocity fluctuation or granular temperature was 1.8 x 10(-3) m(2)/s(2). This method is also applicable to study transport in systems involving turbulence and porous media flows.     

Crossover from intermittent to continuum dynamics for locally driven colloids

We simulate a colloid with charge q(d) driven through a disordered assembly of interacting colloids with charge q and show that, for q(d) approximately q, the velocity-force relation is nonlinear and the velocity fluctuations of the driven particle are highly intermittent with a 1/f characteristic. When g(d) >q , the average velocity drops, the velocity-force relation becomes linear, and the velocity fluctuations are Gaussian. We discuss the results in terms of a crossover from strongly intermittent heterogeneous dynamics to continuum dynamics. We also make several predictions for the transient response in the different regimes.     

Energy spectrum of particles accelerated in relativistic collisionless shocks

We analytically study diffusive particle acceleration in relativistic, collisionless shocks. We find a simple relation between the spectral index s and the anisotropy of the momentum distribution along the shock front. Based on this relation, we obtain s=(3beta(u)-2beta(u)beta(2)(d)+beta(3)(d))/(beta(u)-beta(d)) for isotropic diffusion, where beta(u) (beta(d)) is the upstream (downstream) fluid velocity normalized to the speed of light. This result is in agreement with previous numerical determinations of s for all (beta(u),beta(d)), and yields s=38/9 in the ultrarelativistic limit. The spectrum-anisotropy connection is useful for testing numerical studies and constraining anisotropic diffusion results. It suggests that the spectrum is highly sensitive to the form of the diffusion function for particles traveling along the shock front.     

Nonequilibrium probabilistic dynamics of the logistic map at the edge of chaos

We consider nonequilibrium probabilistic dynamics in logisticlike maps x(t+1)=1-a|x(t)|(z), (z>1) at their chaos threshold: We first introduce many initial conditions within one among W>1 intervals partitioning the phase space and focus on the unique value q(sen)<1 for which the entropic form S(q) identical with (1- summation operator Wp(q)(i))/(q-1) linearly increases with time. We then verify that S(q(sen))(t)-S(q(sen))( infinity ) vanishes like t(-1/[q(rel)(W)-1]) [q(rel)(W)>1]. We finally exhibit a new finite-size scaling, q(rel)( infinity )-q(rel)(W) proportional, variant W(-|q(sen)|). This establishes quantitatively, for the first time, a long pursued relation between sensitivity to the initial conditions and relaxation, concepts which play central roles in nonextensive statistical mechanics.     

薛定谔方程对单个粒子在均匀场中运动的因果描述

把薛定谔方程当成扩展了的经典力学中的雅科毕-哈密顿方程,对单个粒子在均匀场U(x)～±x中的运动进行因果描述。严格求解薛定谔方程,得到了上述两种情况下具有量子力学能级分立特性的粒子的速度随空间位置变化的曲线u(x),这两条速度曲线u(x)都可以遵循对应原理退化到与经典力学的速度曲线U_cla(x)重合。     

Saltation transport on Mars

We present the first calculation of saltation transport and dune formation on Mars and compare it to real dunes. We find that the rate at which grains are entrained into saltation on Mars is 1 order of magnitude higher than on Earth. With this fundamental novel ingredient, we reproduce the size and different shapes of Mars dunes, and give an estimate for the wind velocity on Mars.     

Scaling of directed dynamical small-world networks with random responses

A dynamical model of small-world networks, with directed links which describe various correlations in social and natural phenomena, is presented. Random responses of sites to the input message are introduced to simulate real systems. The interplay of these ingredients results in the collective dynamical evolution of a spinlike variable S(t) of the whole network. The global average spreading length (s) and average spreading time (s) are found to scale as p(-alpha)ln(N with different exponents. Meanwhile, S(t) behaves in a duple scaling form for N>N(*): S approximately f(p(-beta)q(gamma)t), where p and q are rewiring and external parameters, alpha, beta, and gamma are scaling exponents, and f(t) is a universal function. Possible applications of the model are discussed.     

Multifractal PDF analysis for intermittent systems

The formula for probability density functions (PDFs) has been extended to include PDF for energy dissipation rates in addition to other PDFs such as for velocity fluctuations, velocity derivatives, fluid particle accelerations, energy transfer rates, etc., and it is shown that the formula actually explains various PDFs extracted from direct numerical simulations and experiments performed in a wind tunnel. It is also shown that the formula with appropriate zooming increment corresponding to experimental situation gives a new route to obtain the scaling exponents of velocity structure function, including intermittency exponent, out of PDFs of velocity fluctuations.     

A Mechanical Model of Brownian Motion for One Massive Particle Including Slow Light Particles

We provide a connection between Brownian motion and a classical mechanical system. Precisely, we consider a system of one massive particle interacting with an ideal gas, evolved according to non-random mechanical principles, via interaction potentials, without any assumption requiring that the initial velocities of the environmental particles should be restricted to be “fast enough”. We prove the convergence of the (position, velocity)-process of the massive particle under a certain scaling limit, such that the mass of the environmental particles converges to 0 while the density and the velocities of them go to infinity, and give the precise expression of the limiting process, a diffusion process.     

Visualization of saltating sand particle movement near a flat ground surface

Movement of natural sand particles (d=200−300 μm) in a simulated atmospheric boundary layer was visualized using a digital high-speed camera. The consecutive particle images recorded at 2000 fps (frame per second) enabled us to observe the particle transport in detail, especially near the flat sand surface. Various modes of sand saltation were identified. The transverse motion of particles, often ignored in previous studies, was also visualized. In addition, instantaneous velocity fields of saltating particles were obtained using a particle tracking velocimetry (PTV), and statistical analysis of saltating particle trajectories was performed. The qualitative and quantitative results of the present study will be useful for understanding the basic physics of transport of saltating sand particles.     

Relating high-energy Lepton-Hadron, proton-nucleus, and nucleus-nucleus collisions through geometric scaling

A characteristic feature of small-x lepton-proton data from HERA is geometric scaling: the fact that in the region of small Bjorken variable x, x less, similar 0.01, all data can be described by a single variable Q(2)/Q(2)(s,p)(x), with all x dependence encoded in the so-called saturation momentum Q(s,p)(x). Here, we observe that the same scaling ansatz accounts for nuclear photoabsorption cross sections and favors the nuclear dependence Q(2)(s,A) proportional, variant A(alpha)Q(2)(s,p), alpha approximately 4/9. We then make the empirical finding that the same A dependence accounts for the centrality evolution of the multiplicities measured in Au+Au collisions at RHIC. It also allows one to parametrize the high-p(t) particle suppression in d+Au collisions at forward rapidities. If these geometric scaling properties have a common dynamical origin, then this A dependence of Q(2)(s,A) should emerge as a consequence of the underlying dynamical model.     

An Analytical APproach to the Turbulence—Relaxed State in Tokamak Plasmas

Starting from the continuity, temperature, and motion equations of the trapped electron fluid in general tokamak magnetic field with positive or reversed shear and the definition of Lagrangian invariant, dL/dt≡(\partial_t+ u•▽)L=0, where u is convective velocity, the trapped electron dynamics is considered in the following two assumptions: (i) the turbulence is low frequency electrostatic, and (ii) L is a functional only of the density n, temperature T, and magnetic field B, and the effect of perturbation potential φ is included in the convective velocity u, i.e., u is a functional of n, T, B, and φ. The Lagrangian invariant hidden in the trapped electron dynamics is strictly found: L=ln[(n/B)^c₁(T/B^2/3)^c₂], where c₁ and c₂ are dimensionless changeable parameters and c₁∝c₂. From this Lagrangian invariant the turbulent particle and electron thermal transport scaling laws are derived: 〈n>_ψq(ψ)=const. and 〈T^3/2>_ψq(ψ)=const., which, in the limit of large aspect ratio, reduce to n(r)q(r)=const. and T^3/2(r)q(r)=const., respectively. The latter two scaling laws are compared with existent experimental results, being in good agreement.     

Large-eddy simulation of formation of three-dimensional aeolian sand ripples in a turbulent field

With the method of large-eddy simulation, the equation of spherule motion and the method of immersed boundary condition, numerical simulations of three-dimensional turbulent aeolian motion and the formation of sand ripples under three-dimensional turbulent wind and the mutual actions of saltation and creeping motion were carried out. The resulting sand ripples have the form that is flat on the upwind side and steep on the leeward, which is identical to the sand ripples in nature. We also realized the self-restoration process of three-dimensional sand ripples, which shows the correctness of the method of numerical simulation and the models of saltation and creeping. Finally, We analyzed the influence of sand ripples on the three-dimensional turbulent wind field, and found that due to the appearance and development of sand ripples, in the normal direction of ground there exists stronger energy exchange, and moreover, there is close correspondence between the forms of sand ripples and the vorticity close to the ground surface. Supported by the Key Project of National Natural Science Foundation of China (Grant No. 10532040)     

100～500 GPa 4.2Ni2.45Fe0.35CoW合金的冲击压缩性和物态方程

 用二级轻气炮驱动飞片技术及对称碰撞技术，测量了两相合金4.2Ni2.45Fe0.35CoW（以下简称93W）的雨贡纽线。其压力范围为100~500 GPa。实测的冲击波速度D和粒子速度u可用直线关系式D=4.008+1.277u (km/s)描述。实验结果与用混合物雨贡纽线叠加原理的计算结果符合甚好。文中还给出了由实验数据计算得到的物态方程数据。     

A scalar Euclidean theory of gravitation: motion in a static spherically symmetric field

The motion of a particle in a static, spherically symmetric gravitational field is investigated in Euclidean space. The gravitational effects are described as due to a scalar field: To every point in space there is assigned a refractive index deciding the velocity of light in that point. The motion of light in the vacuum is described by the equation of classical optics. An equation of motion for material test particles is then derived by employing the usual Lagrangian formalism. The motion of the planets around the sun is explained, in particular the perihelion motion of Mercury. The present theory fully explains the four classical tests of general relativity in a mathematically far simpler way, and it can be equivalent to the Schwarzschild solution. It is also found that the effect of gravitation depends on the velocity of the particle, becoming repulsive for radial velocities larger thanc/ (c is the velocity of light). This seemingly odd result can also be obtained from the equations of general relativity, as was shown by Cavalleri and Spinelli.     

Embedding of Particle Waves in a Schwarzschild Metric Background

The special and general relativity theories are used to demonstrate that the velocity of an unradiative particle in a Schwarzschild metric background, and in an electrostatic field, is the group velocity of a wave that we call a particle wave, which is a monochromatic solution of a standard equation of wave motion and possesses the following properties. It generalizes the de Broglie wave. The rays of a particle wave are the possible particle trajectories, and the motion equation of a particle can be obtained from the ray equation. The standing particle wave equation generalizes the Schrödinger equation of wave amplitudes. The particle wave motion equation generalizes the Klein–Gordon equation; this result enables us to analyze the essence of the particle wave frequency. The equation of the eikonal of a particle wave generalizes the Hamilton–Jacobi equation; this result enables us to deduce the general expression for the linear momentum. The Heisenberg uncertainty relation expresses the diffraction of the particle wave, and the uncertainty relation connecting the particle instant of presence and energy results from the fact that the group velocity of the particle wave is the particle velocity. A single classical particle may be considered as constituted of geometrical particle wave; reciprocally, a geometrical particle wave may be considered as constituted of classical particles. The expression for a particle wave and the motion equation of the particle wave remain valid when the particle mass is zero. In that case, the particle is a photon, the particle wave is a component a classical electromagnetic wave that is embedded in a Schwarzschild metric background, and the motion equation of the wave particle is the motion equation of an electromagnetic wave in a Schwarzschild metric background. It follows that a particle wave possesses the same physical reality as a classical electromagnetic wave. This last result and the fact that the particle velocity is the group velocity of its wave are in accordance with the opinions of de Broglie and of Schrödinger. We extend these results to the particle subjected to any static field of forces in any gravitational metric background. Therefore we have achieved a synthesis of undulatory mechanics, classical electromagnetism, and gravitation for the case where the field of forces and the gravitational metric background are static, and this synthesis is based only on special and general relativity.