Free turbulent shear layer in a point vortex gas as a problem in nonequilibrium statistical mechanics

This paper attempts to unravel any relations that may exist between turbulent shear flows and statistical mechanics through a detailed numerical investigation in the simplest case where both can be well defined. The flow considered for the purpose is the two-dimensional (2D) temporal free shear layer with a velocity difference Delta U across it, statistically homogeneous in the streamwise direction (x) and evolving from a plane vortex sheet in the direction normal to it (y) in a periodic-in-x domain L x +/-infinity. Extensive computer simulations of the flow are carried out through appropriate initial-value problems for a ``vortex gas'' comprising N point vortices of the same strength (gamma = L Delta U/N) and sign. Such a vortex gas is known to provide weak solutions of the Euler equation. More than ten different initial-condition classes are investigated using simulations involving up to 32 000 vortices, with ensemble averages evaluated over up to 10(3) realizations and integration over 10(4)L/Delta U. The temporal evolution of such a system is found to exhibit three distinct regimes. In Regime I the evolution is strongly influenced by the initial condition, sometimes lasting a significant fraction of L/Delta U. Regime III is a long-time domain-dependent evolution towards a statistically stationary state, via ``violent'' and ``slow'' relaxations P.-H. Chavanis, Physica A 391, 3657 (2012)], over flow time scales of order 10(2) and 10(4)L/Delta U, respectively (for N = 400). The final state involves a single structure that stochastically samples the domain, possibly constituting a ``relative equilibrium.'' The vortex distribution within the structure follows a nonisotropic truncated form of the Lundgren-Pointin (L-P) equilibrium distribution (with negatively high temperatures; L-P parameter lambda close to -1). The central finding is that, in the intermediate Regime II, the spreading rate of the layer is universal over the wide range of cases considered here. The value (in terms of momentum thickness) is 0.0166 +/- 0.0002 times Delta U. Regime II, extensively studied in the turbulent shear flow literature as a self-similar ``equilibrium'' state, is, however, a part of the rapid nonequilibrium evolution of the vortex-gas system, which we term ``explosive'' as it lasts less than one L/Delta U. Regime II also exhibits significant values of N-independent two-vortex correlations, indicating that current kinetic theories that neglect correlations or consider them as O(1/N) cannot describe this regime. The evolution of the layer thickness in present simulations in Regimes I and II agree with the experimental observations of spatially evolving (3D Navier-Stokes) shear layers. Further, the vorticity-stream-function relations in Regime III are close to those computed in 2D Navier-Stokes temporal shear layers J. Sommeria, C. Staquet, and R. Robert, J. Fluid Mech. 233, 661 (1991)]. These findings suggest the dominance of what may be called the Kelvin-Biot-Savart mechanism in determining the growth of the free shear layer through large-scale momentum and vorticity dispersal.

Discussion on the geometric structure of strange attractor

We discuss the transversal heteroclinic cycle formed by hyperbolic periodic pointes of diffeomorphism on the differential manifold. We point out that every possible kind of transversal heteroclinic cycle has the Smalehorse property and the unstable manifolds of hyperbolic periodic points have the closure relation mutually. Therefore the strange attractor may be the closure of unstable manifolds of a countable number of hyperbolic periodic points. The Henon mapping is used as an example to show that the conclusion is reasonable.

The geometric structure and dynamic properties of lauwerier attractor

 Introduction The strange chaotic attractor (ACS) is an important subject in the nonlinear field. On the basis of the theory of transversal heteroclinic cycles, it is suggested that the strange attractor is the closure of the unstable manifolds of countable infinite hyperbolic periodic points. From this point of view some nonlinear phenomena are explained reasonably. 更多还原

Nonexistence of looping trajectories in hydromagnetic waves of finite amplitude. Breaking of waves in a cold collision-free plasma in a magnetic field. On stabilities of periodic waves in a cold collision-free plasma in a magnetic field.

This dissertation consists of three parts. In Part I, it is shown that looping trajectories cannot exist in finite amplitude stationary hydromagnetic waves propagating across a magnetic field in a quasi-neutral cold collision-free plasma. In Part II, time-dependent solutions in series expansion are presented for the magnetic piston problem, which describes waves propagating into a quasi-neutral cold collision-free plasma, ensuing from magnetic disturbances on the boundary of the plasma. The expansion is equivalent to Picard's successive approximations. It is then shown that orbit crossings of plasma particles occur on the boundary for strong disturbances and inside the plasma for weak disturbances. In Part III, the existence of periodic waves propagating at an arbitrary angle to the magnetic field in a plasma is demonstrated by Stokes expansions in amplitude. Then stability analysis is made for such periodic waves with respect to side-band frequency disturbances. It is shown that waves of slow mode are unstable whereas waves of fast mode are stable if the frequency is below the cutoff frequency. The cutoff frequency depends on the propagation angle. For longitudinal propagation the cutoff frequency is equal to one-fourth of the electron's gyrofrequency. For transverse propagation the cutoff frequency is so high that waves of all frequencies are stable.

Earthquakes of the Nepal Himalaya : towards a physical model of the seismic cycle

Home to hundreds of millions of souls and land of excessiveness, the Himalaya is also the locus of a unique seismicity whose scope and peculiarities still remain to this day somewhat mysterious. Having claimed the lives of kings, or turned ancient timeworn cities into heaps of rubbles and ruins, earthquakes eerily inhabit Nepalese folk tales with the fatalistic message that nothing lasts forever. From a scientific point of view as much as from a human perspective, solving the mysteries of Himalayan seismicity thus represents a challenge of prime importance. Documenting geodetic strain across the Nepal Himalaya with various GPS and leveling data, we show that unlike other subduction zones that exhibit a heterogeneous and patchy coupling pattern along strike, the last hundred kilometers of the Main Himalayan Thrust fault, or MHT, appear to be uniformly locked, devoid of any of the “creeping barriers” that traditionally ward off the propagation of large events. The approximately 20 mm/yr of reckoned convergence across the Himalaya matching previously established estimates of the secular deformation at the front of the arc, the slip accumulated at depth has to somehow elastically propagate all the way to the surface at some point. And yet, neither large events from the past nor currently recorded microseismicity nearly compensate for the massive moment deficit that quietly builds up under the giant mountains. Along with this large unbalanced moment deficit, the uncommonly homogeneous coupling pattern on the MHT raises the question of whether or not the locked portion of the MHT can rupture all at once in a giant earthquake. Univocally answering this question appears contingent on the still elusive estimate of the magnitude of the largest possible earthquake in the Himalaya, and requires tight constraints on local fault properties. What makes the Himalaya enigmatic also makes it the potential source of an incredible wealth of information, and we exploit some of the oddities of Himalayan seismicity in an effort to improve the understanding of earthquake physics and cipher out the properties of the MHT. Thanks to the Himalaya, the Indo-Gangetic plain is deluged each year under a tremendous amount of water during the annual summer monsoon that collects and bears down on the Indian plate enough to pull it away from the Eurasian plate slightly, temporarily relieving a small portion of the stress mounting on the MHT. As the rainwater evaporates in the dry winter season, the plate rebounds and tension is increased back on the fault. Interestingly, the mild waggle of stress induced by the monsoon rains is about the same size as that from solid-Earth tides which gently tug at the planets solid layers, but whereas changes in earthquake frequency correspond with the annually occurring monsoon, there is no such correlation with Earth tides, which oscillate back-and-forth twice a day. We therefore investigate the general response of the creeping and seismogenic parts of MHT to periodic stresses in order to link these observations to physical parameters. First, the response of the creeping part of the MHT is analyzed with a simple spring-and-slider system bearing rate-strengthening rheology, and we show that at the transition with the locked zone, where the friction becomes near velocity neutral, the response of the slip rate may be amplified at some periods, which values are analytically related to the physical parameters of the problem. Such predictions therefore hold the potential of constraining fault properties on the MHT, but still await observational counterparts to be applied, as nothing indicates that the variations of seismicity rate on the locked part of the MHT are the direct expressions of variations of the slip rate on its creeping part, and no variations of the slip rate have been singled out from the GPS measurements to this day. When shifting to the locked seismogenic part of the MHT, spring-and-slider models with rate-weakening rheology are insufficient to explain the contrasted responses of the seismicity to the periodic loads that tides and monsoon both place on the MHT. Instead, we resort to numerical simulations using the Boundary Integral CYCLes of Earthquakes algorithm and examine the response of a 2D finite fault embedded with a rate-weakening patch to harmonic stress perturbations of various periods. We show that such simulations are able to reproduce results consistent with a gradual amplification of sensitivity as the perturbing period get larger, up to a critical period corresponding to the characteristic time of evolution of the seismicity in response to a step-like perturbation of stress. This increase of sensitivity was not reproduced by simple 1D-spring-slider systems, probably because of the complexity of the nucleation process, reproduced only by 2D-fault models. When the nucleation zone is close to its critical unstable size, its growth becomes highly sensitive to any external perturbations and the timings of produced events may therefore find themselves highly affected. A fully analytical framework has yet to be developed and further work is needed to fully describe the behavior of the fault in terms of physical parameters, which will likely provide the keys to deduce constitutive properties of the MHT from seismological observations.

Periodic nanovoid structures via femtosecond laser irradiation

We have observed periodically aligned nanovoid structures inside a conventional borosilicate glass induced by a single femtosecond (fs) laser beam for the first time, to our knowledge. The spherical voids of nanosized diameter were aligned spontaneously with a period along the propagation direction of the laser beam. The period, the number of voids, and the whole length of the aligned void structure were controlled by changing the laser power, the pulse number, and the position of the focal point.

Functional significance of stiffness in adaptation of multijoint arm movements to stable and unstable dynamics.

This study compared the mechanisms of adaptation to stable and unstable dynamics from the perspective of changes in joint mechanics. Subjects were instructed to make point to point movements in force fields generated by a robotic manipulandum which interacted with the arm in either a stable or an unstable manner. After subjects adjusted to the initial disturbing effects of the force fields they were able to produce normal straight movements to the target. In the case of the stable interaction, subjects modified the joint torques in order to appropriately compensate for the force field. No change in joint torque or endpoint force was required or observed in the case of the unstable interaction. After adaptation, the endpoint stiffness of the arm was measured by applying displacements to the hand in eight different directions midway through the movements. This was compared to the stiffness measured similarly during movements in a null force field. After adaptation, the endpoint stiffness under both the stable and unstable dynamics was modified relative to the null field. Adaptation to unstable dynamics was achieved by selective modification of endpoint stiffness in the direction of the instability. To investigate whether the change in endpoint stiffness could be accounted for by change in joint torque or endpoint force, we estimated the change in stiffness on each trial based on the change in joint torque relative to the null field. For stable dynamics the change in endpoint stiffness was accurately predicted. However, for unstable dynamics the change in endpoint stiffness could not be reproduced. In fact, the predicted endpoint stiffness was similar to that in the null force field. Thus, the change in endpoint stiffness seen after adaptation to stable dynamics was directly related to changes in net joint torque necessary to compensate for the dynamics in contrast to adaptation to unstable dynamics, where a selective change in endpoint stiffness occurred without any modification of net joint torque.

Some remarks about periodic feedback stabilization

Periodic feedback stabilization is a very natural solution to overcome the topological obstructions which may occur when one tries to asymptotically (locally) stabilize a (locally) controllable nonlinear system around an equilibrium point. The object of this paper is to give a simple geometric interpretation of this fact, to show that one obtains a weakened form of those obstructions when periodic feedback is used, and to illustrate the success of periodic feedback stabilization on a representative system which contains a drift.

Symmetry-extended counting rules for periodic frameworks.

A symmetry-adapted version of the Maxwell rule appropriate to periodic bar-and-joint frameworks is obtained, and is further extended to body-and-joint systems. The treatment deals with bodies and forces that are replicated in every unit cell, and uses the point group isomorphic to the factor group of the space group of the framework. Explicit expressions are found for the numbers and symmetries of detectable mechanisms and states of self-stress in terms of the numbers and symmetries of framework components. This approach allows detection and characterization of mechanisms and states of self-stress in microscopic and macroscopic materials and meta-materials. Illustrative examples are described. The notion of local isostaticity of periodic frameworks is extended to include point-group symmetry.

Theoretical and experimental investigations of flexural wave propagation in periodic grid structures designed with the idea of phononic crystals

The propagation characteristics of fiexural waves in periodic grid structures designed with the idea of phononic crystals are investigated by combining the Bloch theorem with the finite element method. This combined analysis yields phase constant surfaces, which predict the location and the extension of band gaps, as well as the directions and the regions of wave propagation at assigned frequencies. The predictions are validated by computation and experimental analysis of the harmonic responses of a finite structure with 11 × 11 unit cells. The fiexural wave is localized at the point of excitation in band gaps, while the directional behaviour occurs at particular frequencies in pass bands. These studies provide guidelines to designing periodic structures for vibration attenuation.

Dynamical Equations For The Period Vectors In A Periodic System Under Constant External Stress

The purpose of this paper is to derive the dynamical equations for the period vectors of a periodic system under constant external stress. The explicit starting point is Newton’s second law applied to halves of the system. Later statistics over indistinguishable translated states and forces associated with transport of momentum are applied to the resulting dynamical equations. In the final expressions, the period vectors are driven by the imbalance between internal and external stresses. The internal stress is shown to have both full interaction and kinetic-energy terms.

Drifting periodic structures in a degenerate-planar bent-rod nematic liquid crystal beyond the dielectric inversion frequency

We report on the electric-field-generated effects in the nematic phase of a twin mesogen formed of bent-core and calamitic units, aligned homeotropically in the initial ground state and examined beyond the dielectric inversion point. The bend-Freedericksz (BF) state occurring at the primary bifurcation and containing a network of umbilics is metastable; we focus here on the degenerate planar (DP) configuration that establishes itself at the expense of the BF state in the course of an anchoring transition. In the DP regime, normal rolls, broad domains, and chevrons (both defect-mediated and defect-free types) form at various linear defect-sites, in different regions of the frequency-voltage plane. A significant novel aspect common to all these patterned states is the sustained propagative instability, which does not seem explicable on the basis of known driving mechanisms.

Dynamics of nonlinearly coupled magnetic-field-aligned electromagnetic electron-cyclotron waves near the zero-group-dispersion point in magnetized plasmas

The nonlinear coupling between two magnetic-field-aligned electromagnetic electron-cyclotron (EMEC) waves in plasmas is considered. Evaluating the ponderomotive coupling between the EMEC waves and quasistationary plasma density perturbations, a pair of coupled nonlinear Schrodinger equations (CNLSEs) is obtained. The CNLSEs are then used to investigate the occurrence of modulational instability in magnetized plasmas. Waves in the vicinity of the zero-group-dispersion point are considered, so that the group dispersion terms may either bear the same or different signs. It is found that a stable EMEC wave can be destabilized due to its nonlinear interactions with an unstable one, while a pair of unstable EMEC waves yields an increased instability growth rate. Individually stable waves remain stable while interacting with one another. Stationary nonlinear solutions of the coupled equations are presented. The relevance of our investigation to nonlinear phenomena in space plasmas is discussed. (c) 2005 American Institute of Physics.

Localizationlike effect in two-dimensional alternate quantum walks with periodic coin operations

Exploiting multidimensional quantum walks as feasible platforms for quantum computation and quantum simulation attracts constantly growing attention from a broad experimental physics community. Here, we propose a two-dimensional quantum walk scheme with a single-qubit coin that presents, in the considered regimes, a strong localizationlike effect on the walker. The result could provide new possible directions for the implementation of quantum algorithms or from the point of view of quantum simulation. We characterize the localizationlike effect in terms of the parameters of a step-dependent qubit operation that acts on the coin space after any standard coin operation, showing that a proper choice can guarantee a nonnegligible probability of finding the walker in the origin even for large times. We finally discuss the robustness to imperfections, a qualitative relation with coherences behavior, and possible experimental realizations of this model with the current state-of-the-art settings.

Monte Carlo study of the XY-model on quasi-periodic lattices

Monte Carlo Simulations were carried out using a nearest neighbour ferromagnetic XYmodel, on both 2-D and 3-D quasi-periodic lattices. In the case of 2-D, both the unfrustrated and frustrated XV-model were studied. For the unfrustrated 2-D XV-model, we have examined the magnetization, specific heat, linear susceptibility, helicity modulus and the derivative of the helicity modulus with respect to inverse temperature. The behaviour of all these quatities point to a Kosterlitz-Thouless transition occuring in temperature range Te == (1.0 -1.05) JlkB and with critical exponents that are consistent with previous results (obtained for crystalline lattices) . However, in the frustrated case, analysis of the spin glass susceptibility and EdwardsAnderson order parameter, in addition to the magnetization, specific heat and linear susceptibility, support a spin glass transition. In the case where the 'thin' rhombus is fully frustrated, a freezing transition occurs at Tf == 0.137 JlkB , which contradicts previous work suggesting the critical dimension of spin glasses to be de > 2 . In the 3-D systems, examination of the magnetization, specific heat and linear susceptibility reveal a conventional second order phase transition. Through a cumulant analysis and finite size scaling, a critical temperature of Te == (2.292 ± 0.003) JI kB and critical exponents of 0:' == 0.03 ± 0.03, f3 == 0.30 ± 0.01 and I == 1.31 ± 0.02 have been obtained.