We present a mechanism of amplification of phonons by phonons on the basis of nonlinear band-gap transmission phenomenon. As a concept the idea may be applied to the various number of systems; however we introduce the specific idea of creating an amplification scenario in the chain of coupled cantilever arrays. One chain is driven at the constant frequency located in the upper band of the ladder system, thus no wave enters the system. However the frequency is specifically chosen to be very close to the maximum value of the frequency corresponding to the dispersion relation of the system. An amplification scenario happens when a counter phase pulse of the same frequency with a small amplitude is introduced to the second chain. If both signals exceed a threshold amplitude for the band-gap transmission a large amplitude soliton enters the system—therefore we have an amplifier. Although the concept may be applied in a variety of contexts, all-optical or all-magnonic systems, we choose the system of coupled cantilever arrays and represent a clear example of the application of the presented conceptual idea. Logical operations is the other probable field, where such a mechanism could be used, which might significantly broaden the horizon of the considered applications of band-gap soliton dynamics.

The generalized $H\left(n\right)$ Hirsch index of order $n$ has been recently introduced and shown to interpolate between the degree and the $K$-core centrality in networks. We provide a detailed analytic characterization of the properties of sets of nodes having the same $H\left(n\right)$, within the annealed network approximation. The connection between the Hirsch indices and the degree is highlighted. Numerical tests in synthetic uncorrelated networks and real-world correlated ones validate the findings. We also test the use of the Hirsch index for the identification of influential spreaders in networks, finding that it is in general outperformed by the recently introduced nonbacktracking centrality.

We study a two-dimensional nonvibrating granular system as a model of a magnetorheological fluid. The system is composed of magnetic steel particles on a horizontal plane under a vertical sinusoidal magnetic field and a horizontal static magnetic field. When the amplitude of the horizontal field is zero, we find that the motion of the particles has characteristics similar to those of Brownian particles. A slowing down of the dynamics is observed as the particle concentration increases or the magnitude of the vertical magnetic field decreases. When the amplitude of the horizontal field is nonzero, the particles interact through effective dipolar interactions. Above a threshold in the amplitude of the horizontal field, particles form chains that become longer and more stable as time increases. For some conditions, at short time intervals, the average chain length as a function of time exhibits scaling behavior. The chain length distribution at a given time is a decreasing exponential function. The behavior of this granular system is consistent with theoretical and experimental results for magnetorheological fluids.

Pumping of electrolytes in microchannels can be achieved with the use of microelectrodes subjected to AC potentials. Experiments have shown an influence of Faradaic currents in the pumping performance, and theoretical studies for asymmetric electrolytes suggest that induced charges in the diffusion layer play an important role. In this work we consider the case of a diffusion layer induced by an array of electrodes subjected to a traveling wave potential and we include Faradaic currents. Previous theoretical studies considered the case of very small applied voltages, which allowed for two major simplifications: (i) Butler-Volmer (B-V) equation was linearized, and (ii) the presence of gradients in ion concentration was neglected. We extend previous results and used the full nonlinear B-V equation. A comparison with the linear limit shows that the flow rate in both cases coincides for voltages around and below $\approx 0.25$ V. For voltages larger than this, the nonlinear equations show that gradients in ion concentration appear and have an important influence, therefore, the predictions deviate from the linear model. We show that the electrical force in the diffusion layer can induce pumping either in the same or the opposite direction of the applied traveling-wave potential and it could be responsible for the reversal of the flow as observed in experiments.

The collision of a spherical grain with a granular bed is commonly parametrized by the splash function, which provides the velocity of the rebounding grain and the velocity distribution and number of ejected grains. Starting from elementary geometric considerations and physical principles, like momentum conservation and energy dissipation in inelastic pair collisions, we derive a rebound parametrization for the collision of a spherical grain with a granular bed. Combined with a recently proposed energy-splitting model [Ho et al., Phys. Rev. E 85, 052301 (2012)] that predicts how the impact energy is distributed among the bed grains, this yields a coarse-grained but complete characterization of the splash as a function of the impact velocity and the impactor-bed grain-size ratio. The predicted mean values of the rebound angle, total and vertical restitution, ejection speed, and number of ejected grains are in excellent agreement with experimental literature data and with our own discrete-element computer simulations. We extract a set of analytical asymptotic relations for shallow impact geometries, which can readily be used in coarse-grained analytical modeling or computer simulations of geophysical particle-laden flows.

We report an experimental and computational study investigating the effects of laser preheat on the hydrodynamic behavior of a material layer. In particular, we find that perturbation of the surface of the layer results in a complex interaction, in which the bulk of the layer develops density, pressure, and temperature structure and in which the surface experiences instability-like behavior, including mode coupling. A uniform one-temperature preheat model is used to reproduce the experimentally observed behavior, and we find that this model can be used to capture the evolution of the layer, while also providing evidence of complexities in the preheat behavior. This result has important consequences for inertially confined fusion plasmas, which can be difficult to diagnose in detail, as well as for laser hydrodynamics experiments, which generally depend on assumptions about initial conditions in order to interpret their results.

The Weibel instability driven by two symmetric counterstreaming relativistic electron plasmas, also referred to as current-filamentation instability, is studied in a constant and uniform external magnetic field aligned with the plasma flows. Both the linear and nonlinear stages of the instability are investigated using analytical modeling and particle-in-cell simulations. While previous studies have already described the stabilizing effect of the magnetic field, we show here that the saturation stage is only weakly affected. The different mechanisms responsible for the saturation are discussed in detail in the relativistic cold fluid framework considering a single unstable mode. The application of an external field leads to a slight increase of the saturation level for large wavelengths, while it does not affect the small wavelengths. Multimode and temperature effects are then investigated. While at high temperature the saturation level is independent of the external magnetic field, at low but finite temperature the competition between different modes in the presence of an external magnetic field leads to a saturation level lower with respect to the unmagnetized case.

In this paper we provide a theoretical model describing the evolution of the charged-particle distribution function in a system with nonlinear wave-particle interactions. Considering a system with strong electrostatic waves propagating in an inhomogeneous magnetic field, we demonstrate that individual particle motion can be characterized by the probability of trapping into the resonance with the wave and by the efficiency of scattering at resonance. These characteristics, being derived for a particular plasma system, can be used to construct a kinetic equation (or generalized Fokker-Planck equation) modeling the long-term evolution of the particle distribution. In this equation, effects of charged-particle trapping and transport in phase space are simulated with a nonlocal operator. We demonstrate that solutions of the derived kinetic equations agree with results of test-particle tracing. The applicability of the proposed approach for the description of space and laboratory plasma systems is also discussed.

The lattice Boltzmann method is a simulation technique in computational fluid dynamics. In its standard formulation, it is restricted to regular computation grids, second-order spatial accuracy, and a unity Courant-Friedrichs-Lewy (CFL) number. This paper advances the standard lattice Boltzmann method by introducing a semi-Lagrangian streaming step. The proposed method allows significantly larger time steps, unstructured grids, and higher-order accurate representations of the solution to be used. The appealing properties of the approach are demonstrated in simulations of a two-dimensional Taylor-Green vortex, doubly periodic shear layers, and a three-dimensional Taylor-Green vortex.

Shear banding and stick-slip instabilities have been long observed in sheared granular materials. Yet, their microscopic underpinnings, interdependencies, and variability under different loading conditions have not been fully explored. Here we use a nonequilibrium thermodynamics model, the Shear Transformation Zone theory, to investigate the dynamics of strain localization and its connection to stability of sliding in sheared, dry, granular materials. We consider frictional and frictionless grains as well as the presence and absence of acoustic vibrations. Our results suggest that at low and intermediate strain rates, persistent shear bands develop only in the absence of vibrations. Vibrations tend to fluidize the granular network and delocalize slip at these rates. Stick-slip is observed only for frictional grains, and it is confined to the shear band. At high strain rates, stick-slip disappears and the different systems exhibit similar stress-slip response. Changing the vibration intensity, duration or time of application alters the system response and may cause long-lasting rheological changes. We analyze these observations in terms of possible transitions between rate strengthening and rate weakening response facilitated by a competition between shear-induced dilation and vibration-induced compaction. We discuss the implications of our results on dynamic triggering, quiescence, and strength evolution in gouge-filled fault zones.

This study elucidates the origin of the multiplicity of stable oscillatory flows detected by time integration in tall rectangular cavities heated from the side. By using continuation techniques for periodic orbits, it is shown that initially unstable branches, arising at Hopf bifurcations of the basic steady flow, become stable after crossing Neimark-Sacker points. There are no saddle-node or pitchfork bifurcations of periodic orbits, which could have been alternative mechanisms of stabilization. According to the symmetries of the system, the orbits are either fixed cycles, which retain at any time the center symmetry of the steady flow, or symmetric cycles involving a time shift in the global invariance of the orbit. The bifurcation points along the branches of periodic flows are determined. By using time integrations, with unstable periodic solutions as initial conditions, we determine which of the bifurcations at the limits of the intervals of stable periodic orbits are sub- or supercritical.

We describe the physical mechanism responsible for electrohydrodynamic (EHD) instability of a fluid layer with collinear conductivity gradient and electric field. In particular, we resolve the ambiguity in literature regarding the cause for switching between stationary and oscillatory modes of EHD instability. Using linear stability analysis, we show that a small perturbation in conductivity field perturbs the local electric field and also induces a perturbation charge. The coupling of base-state electric field with the perturbation charge leads to a force which causes overstability. Whereas, the coupling of base-state free charge with perturbation electric field leads to a force which causes EHD instability via a stationary mode. The proposed mechanism correctly explains the existence of stationary and oscillatory modes for varying conductivity gradients and wave number of disturbances, depending upon the relative magnitude of these two forces.

An idea for evaluating transition probabilities in chemical reaction systems is proposed, which is efficient for repeated calculations with various rate constants. The idea is based on duality relations; instead of direct time evolutions of the original reaction system, the dual process is dealt with. Usually, if one changes rate constants of the original reaction system, the direct time evolutions should be performed again, using the new rate constants. On the other hands, only one solution of an extended dual process can be reused to calculate the transition probabilities for various rate constant cases. The idea is demonstrated in a parameter estimation problem for the Lotka-Volterra system.

We have intensively examined the dynamic behavior of flame front instability in a lean swirling premixed flame generated by a change in gravitational orientation [H. Gotoda, T. Miyano, and I. G. Shepherd, Phys. Rev. E 81, 026211 (2010)] from the viewpoints of complex networks, symbolic dynamics, and statistical complexity. Here, we considered the permutation entropy in combination with the surrogate data method, the permutation spectrum test, and the multiscale complexity-entropy causality plane incorporating a scale-dependent approach, none of which have been considered in the study of flame front instabilities. Our results clearly show the possible presence of chaos in flame front dynamics induced by the coupling of swirl-buoyancy interaction in inverted gravity. The flame front dynamics also possesses a scale-free structure, which is reasonably shown by the probability distribution of the degree in $\epsilon$-recurrence networks.

We study the dynamical response of a diatomic periodic chain of rotors coupled by springs, whose unit cell breaks spatial inversion symmetry. In the continuum description, we derive a nonlinear field theory which admits topological kinks and antikinks as nonlinear excitations but where a topological boundary term breaks the symmetry between the two and energetically favors the kink configuration. Using a cobweb plot, we develop a fixed-point analysis for the kink motion and demonstrate that kinks propagate without the Peierls-Nabarro potential energy barrier typically associated with lattice models. Using continuum elasticity theory, we trace the absence of the Peierls-Nabarro barrier for the kink motion to the topological boundary term which ensures that only the kink configuration, and not the antikink, costs zero potential energy. Further, we study the eigenmodes around the kink and antikink configurations using a tangent stiffness matrix approach appropriate for prestressed structures to explicitly show how the usual energy degeneracy between the two no longer holds. We show how the kink-antikink asymmetry also manifests in the way these nonlinear excitations interact with impurities introduced in the chain as disorder in the spring stiffness. Finally, we discuss the effect of impurities in the (bond) spring length and build prototypes based on simple linkages that verify our predictions.

We performed extensive simulations of the random-bond Ising model confined between walls where competitive surface fields act. By properly taking the thermodynamic limit we unambiguously determined wetting transition points of the system, as extrapolation of localization-delocalization transitions of the interface between domains of different orientation driven by the respective fields. The finite-size scaling theory for wetting with short-range fields [E. V. Albano and K. Binder, Phys. Rev. Lett. 109, 036101 (2012)] establishes that the average magnetization of the sample, with critical exponent $\beta$, is the proper order parameter for the study of wetting. While the hyperscaling relationship given by $\gamma +2\beta ={\nu }_{\parallel }+{\nu }_{\perp }$ requires $\beta =1/2$ ($\gamma =4, {\nu }_{\parallel }=3$, and ${\nu }_{\perp }=2$), the thermodynamic scaling establishes that ${\mathrm{\Delta }}_s=\gamma +\beta$, which in contrast requires $\beta =0$ (${\mathrm{\Delta }}_s=4$), where $\gamma , {\nu }_{\parallel }, {\nu }_{\perp }$, and ${\mathrm{\Delta }}_s$ are the critical exponents of the susceptibility, the correlation lengths parallel and perpendicular to the interface, and the gap exponent, respectively. So, we formulate a finite-size scaling theory for wetting without hyperscaling and perform numerical simulations that provide strong evidence of hyperscaling violation (i.e., $\beta =0$) and a direct measurement of the susceptibility critical exponent $\gamma /{\nu }_{\perp }=2.0\pm 0.2$, in agreement with theoretical results for the strong fluctuation regime of wetting transitions with quenched noise.

Viscoelasticity has been described since the time of Maxwell as an interpolation of purely viscous and purely elastic response, but its microscopic atomic-level mechanism in solids has remained elusive. We studied three model disordered solids: a random lattice, the bond-depleted fcc lattice, and the fcc lattice with vacancies. Within the harmonic approximation for central-force lattices, we applied sum rules for viscoelastic response derived on the basis of nonaffine atomic motions. The latter motions are a direct result of local structural disorder, and in particular, of the lack of inversion symmetry in disordered lattices. By defining a suitable quantitative and general atomic-level measure of nonaffinity and inversion symmetry, we show that the viscoelastic responses of all three systems collapse onto a master curve upon normalizing by the overall strength of inversion-symmetry breaking in each system. Close to the isostatic point for central-force lattices, power-law creep $G\left(t\right)\sim t^{-1/2}$ emerges as a consequence of the interplay between soft vibrational modes and nonaffine dynamics, and various analytical scalings, supported by numerical calculations, are predicted by the theory.

Modeling flow in porous media is usually focused on the governing equations for mass and momentum transport, which yield the velocity and pressure at the pore or Darcy scales. However, in many applications, it is important to determine the work (or power) needed to induce flow in porous media, and this can be achieved when the mechanical energy equation is taken into account. At the macroscopic scale, this equation may be postulated to be the result of the inner product of Darcy's law and the seepage velocity. However, near the porous medium boundaries, this postulate seems questionable due to the spatial variations of the effective properties (velocity, permeability, porosity, etc.). In this work we derive the macroscopic mechanical energy equation using the method of volume averaging for the simple case of incompressible single-phase flow in porous media. Our analysis shows that the result of averaging the pore-scale version of the mechanical energy equation at the Darcy scale is not, in general, the expected product of Darcy's law and the seepage velocity. As a matter of fact, this result is only applicable in the bulk region of the porous medium and, in the derivation of this result, the properties of the permeability tensor are determinant. Furthermore, near the porous medium boundaries, a more novel version of the mechanical energy equation is obtained, which incorporates additional terms that take into account the rapid variations of structural properties taking place in this particular portion of the system. This analysis can be applied to multiphase and compressible flows in porous media and in many other multiscale systems.

We present a molecular dynamics simulation of shock waves propagating in dense deuterium with the electron force field method [J. T. Su and W. A. Goddard, Phys. Rev. Lett. 99, 185003 (2007)], which explicitly takes the excitation of electrons into consideration. Nonequilibrium features associated with the excitation of electrons are systematically investigated. We show that chemical bonds in ${\mathrm{D}}_2$ molecules lead to a more complicated shock wave structure near the shock front, compared with the results of classical molecular dynamics simulation. Charge separation can bring about accumulation of net charges on large scales, instead of the formation of a localized dipole layer, which might cause extra energy for the shock wave to propagate. In addition, the simulations also display that molecular dissociation at the shock front is the major factor that accounts for the “bump” structure in the principal Hugoniot. These results could help to build a more realistic picture of shock wave propagation in fuel materials commonly used in the inertial confinement fusion.