No Arabic abstract
We use direct statistical simulation (DSS) to find the low-order statistics of the well-known dynamical system, the Lorenz63 model. Instead of accumulating statistics from numerical simulation of the dynamical systems, we solve the equations of motion for the statistics themselves after closing them by making several different choices for the truncation. Fixed points of the statistics are obtained either by time evolving, or by iterative methods. Statistics so obtained are compared to those found by the traditional approach.
The dynamics of macroscopically homogeneous sheared suspensions of neutrally buoyant, non-Brownian spheres is investigated in the limit of vanishingly small Reynolds numbers using Stokesian dynamics. We show that the complex dynamics of sheared suspensions can be characterized as a chaotic motion in phase space and determine the dependence of the largest Lyapunov exponent on the volume fraction $phi$. The loss of memory at the microscopic level of individual particles is also shown in terms of the autocorrelation functions for the two transverse velocity components. Moreover, a negative correlation in the transverse particle velocities is seen to exist at the lower concentrations, an effect which we explain on the basis of the dynamics of two isolated spheres undergoing simple shear. In addition, we calculate the probability distribution function of the velocity fluctuations and observe, with increasing $phi$, a transition from exponential to Gaussian distributions. The simulations include a non-hydrodynamic repulsive interaction between the spheres which qualitatively models the effects of surface roughness and other irreversible effects, such as residual Brownian displacements, that become particularly important whenever pairs of spheres are nearly touching. We investigate the effects of such a non-hydrodynamic interparticle force on the scaling of the particle tracer diffusion coefficient $D$ for very dilute suspensions, and show that, when this force is very short-ranged, $D$ becomes proportional to $phi^2$ as $phi to 0$. In contrast, when the range of the non-hydrodynamic interaction is increased, we observe a crossover in the dependence of $D$ on $phi$, from $phi^2$ to $phi$ as $phi to 0$.
We propose a new look at the heat bath for two Brownian particles, in which the heat bath as a `system is both perturbed and sensed by the Brownian particles. Non-local thermal fluctuation give rise to bath-mediated static forces between the particles. Based on the general sum-rule of the linear response theory, we derive an explicit relation linking these forces to the friction kernel describing the particles dynamics. The relation is analytically confirmed in the case of two solvable models and could be experimentally challenged. Our results point out that the inclusion of the environment as a part of the whole system is important for micron- or nano-scale physics.
We study the evolution of a reactive field advected by a one-dimensional compressible velocity field and subject to an ignition-type nonlinearity. In the limit of small molecular diffusivity the problem can be described by a spatially discretized system, and this allows for an efficient numerical simulation. If the initial field profile is supported in a region of size l < lc one has quenching, i.e., flame extinction, where lc is a characteristic length-scale depending on the system parameters (reacting time, molecular diffusivity and velocity field). We derive an expression for lc in terms of these parameters and relate our results to those obtained by other authors for different flow settings.
Linking thermodynamic variables like temperature $T$ and the measure of chaos, the Lyapunov exponents $lambda$, is a question of fundamental importance in many-body systems. By using nonlinear fluid equations in one and three dimensions, we prove that in thermalised flows $lambda propto sqrt{T}$, in agreement with results from frustrated spin systems. This reveals an underlying universality and provides evidence for recent conjectures on the thermal scaling of $lambda$. We also reconcile seemingly disparate effects -- equilibration on one hand and pushing systems out-of-equilibrium on the other -- of many-body chaos by relating $lambda$ to $T$ through the dynamical structures of the flow.
We numerically study the two-dimensional, area preserving, web map. When the map is governed by ergodic behavior, it is, as expected, correctly described by Boltzmann-Gibbs statistics, based on the additive entropic functional $S_{BG}[p(x)] = -kint dx,p(x) ln p(x)$. In contrast, possible ergodicity breakdown and transitory sticky dynamical behavior drag the map into the realm of generalized $q$-statistics, based on the nonadditive entropic functional $S_q[p(x)]=kfrac{1-int dx,[p(x)]^q}{q-1}$ ($q in {cal R}; S_1=S_{BG}$). We statistically describe the system (probability distribution of the sum of successive iterates, sensitivity to the initial condition, and entropy production per unit time) for typical values of the parameter that controls the ergodicity of the map. For small (large) values of the external parameter $K$, we observe $q$-Gaussian distributions with $q=1.935dots$ (Gaussian distributions), like for the standard map. In contrast, for intermediate values of $K$, we observe a different scenario, due to the fractal structure of the trajectories embedded in the chaotic sea. Long-standing non-Gaussian distributions are characterized in terms of the kurtosis and the box-counting dimension of chaotic sea.