Fiber bundle models are one of the most fundamental modeling approaches for the investigation of the fracture of heterogeneous materials being able to capture a broad spectrum of damage mechanisms, loading conditions, and types of load sharing. In the framework of the fiber bundle model we introduce a kinetic Monte Carlo algorithm to investigate the thermally induced creep rupture of materials occurring under a constant external load. We demonstrate that the method overcomes several limitations of previous techniques and provides an efficient numerical framework at any load and temperature values. We show for both equal and localized load sharing that the computational time does not depend on the temperature; it is solely determined by the external load and the system size. In the limit of low load where the lifetime of the system diverges, the computational time saturates to a constant value. The method takes into account the secondary failures induced by subsequent load redistributions after breaking events, with the additional advantage that breaking avalanches always start from a single broken fiber.

We consider the low but nonzero-temperature regimes of the Glauber dynamics in a chain of Ising spins with first- and second-neighbor interactions $J_1,J_2$. For $0<-J_2/\left|J_1\right|<1$ it is known that at $T=0$ the dynamics is both metastable and noncoarsening, while being always ergodic and coarsening in the limit of $T\to 0^{+}$. Based on finite-size scaling analyses of relaxation times, here we argue that in that latter situation the asymptotic kinetics of small or weakly frustrated $-J_2/\left|J_1\right|$ ratios is characterized by an almost ballistic dynamic exponent $z\simeq 1.03\left(2\right)$ and arbitrarily slow velocities of growth. By contrast, for noncompeting interactions the coarsening length scales are estimated to be almost diffusive.

The dynamical responses of a complex system to external perturbations are of both fundamental interest and practical significance. Here, by the model of networked chaotic oscillators, we investigate how the synchronization behavior of a complex network is influenced by an externally added periodic driving. Interestingly, it is found that by a slight change of the properties of the external driving, e.g., the frequency or phase lag between its intrinsic oscillation and external driving, the network synchronizability could be significantly modified. We demonstrate this phenomenon by different network models and, based on the method of master stability function, give an analysis on the underlying mechanisms. Our studies highlight the importance of external perturbations on the collective behaviors of complex networks, and also provide an alternate approach for controlling network synchronization.

We investigate numerically Turing patterns in the Lengyel-Epstein model in three dimensions. In a bulk homogeneous system under periodic boundary conditions, we obtain not only lamellar, cylindrical, and spherical structures but also several interconnected periodic structures including the Schwartz P-surface structure. In order to examine Turing patterns in the conditions accessible experimentally, we consider inhomogeneous systems where a parameter in the reaction-diffusion equations depends on the space coordinate with either Dirichlet or Neumann boundary conditions. In this situation, we find that a perforated-lamellar structure and an $Fddd$ structure, both of which have a uniaxial symmetry, appear depending on the boundary conditions.

We present the nonautonomous discrete bright soliton solutions and their interactions in the discrete Ablowitz-Ladik (DAL) equation with variable coefficients, which possesses complicated wave propagation in time and differs from the usual bright soliton waves. The differential-difference similarity transformation allows us to relate the discrete bright soliton solutions of the inhomogeneous DAL equation to the solutions of the homogeneous DAL equation. Propagation and interaction behaviors of the nonautonomous discrete solitons are analyzed through the one- and two-soliton solutions. We study the discrete snaking behaviors, parabolic behaviors, and interaction behaviors of the discrete solitons. In addition, the interaction management with free functions and dynamic behaviors of these solutions is investigated analytically, which have certain applications in electrical and optical systems.

This paper presents a nonparametric modeling approach for forecasting stochastic dynamical systems on low-dimensional manifolds. The key idea is to represent the discrete shift maps on a smooth basis which can be obtained by the diffusion maps algorithm. In the limit of large data, this approach converges to a Galerkin projection of the semigroup solution to the underlying dynamics on a basis adapted to the invariant measure. This approach allows one to quantify uncertainties (in fact, evolve the probability distribution) for nontrivial dynamical systems with equation-free modeling. We verify our approach on various examples, ranging from an inhomogeneous anisotropic stochastic differential equation on a torus, the chaotic Lorenz three-dimensional model, and the Niño-3.4 data set which is used as a proxy of the El Niño Southern Oscillation.

Free-surface turbulence, albeit constrained onto a two-dimensional space, exhibits features that barely resemble predictions of simplified two-dimensional modeling. We demonstrate that, in a three-dimensional open channel flow, surface turbulence is characterized by upscale energy transfer, which controls the long-term evolution of the larger scales. We are able to associate downscale and upscale energy transfer at the surface with the two-dimensional divergence of velocity. We finally demonstrate that surface compressibility confirms the strongly three-dimensional nature of surface turbulence.

A numerical and experimental investigation of the acoustic streaming flow in the near field of a circular plane ultrasonic transducer in water is performed. The experimental domain is a parallelepipedic cavity delimited by absorbing walls to avoid acoustic reflection, with a top free surface. The flow velocities are measured by particle image velocimetry, leading to well-resolved velocity profiles. The theoretical model is based on a linear acoustic propagation model, which correctly reproduces the acoustic field mapped experimentally using a hydrophone, and an acoustic force term introduced in the Navier-Stokes equations under the plane-wave assumption. Despite the complexity of the acoustic field in the near field, in particular in the vicinity of the acoustic source, a good agreement between the experimental measurements and the numerical results for the velocity field is obtained, validating our numerical approach and justifying the planar wave assumption in conditions where it is a priori far from obvious. The flow structure is found to be correlated with the acoustic field shape. Indeed, the longitudinal profiles of the velocity present a wavering linked to the variations in acoustic intensity along the beam axis and transverse profiles exhibit a complex shape strongly influenced by the transverse variations of the acoustic intensity in the beam. Finally, the velocity in the jet is found to increase as the square root of the acoustic force times the distance from the origin of the jet over a major part of the cavity, after a strong short initial increase, where the velocity scales with the square of the distance from the upstream wall.

We study how the stability of spherical crystalline shells under external pressure is influenced by the defect structure. In particular, we compare stability for shells with a minimal set of topologically required defects to shells with extended defect arrays (grain boundary “scars” with nonvanishing net disclination charge). We perform both Monte Carlo and conjugate gradient simulations to compare how shells with and without scars deform quasistatically under external hydrostatic pressure. We find that the critical pressure at which shells collapse is lowered for scarred configurations that break icosahedral symmetry and raised for scars that preserve icosahedral symmetry. The particular shapes which arise from breaking of an initial icosahedrally symmetric shell depend on the Föppl–von Kármán number.

Scroll waves of reentrant activity and their interactions pose a serious threat to cardiac health. In experiments with the Belousov-Zhabotinsky reaction we demonstrate the interaction of scroll waves. We show that depending on their mutual orientation, two scroll rings can push each other away and rupture on touching the system boundary, or they can reconnect to form a single, large ring. Reconnection only occurs when the filaments lie within one core length of each other. The reconnected filament has extended lifetimes, which could have serious implications in systems where they occur. The experimental results are explained on the basis of a simple numerical model.

We study rigidity percolation transitions in two-dimensional central-force isostatic lattices, including the square and the kagome lattices, as next-nearest-neighbor bonds (“braces”) are randomly added to the system. In particular, we focus on the differences between regular lattices, which are perfectly periodic, and generic lattices with the same topology of bonds but whose sites are at random positions in space. We find that the regular square and kagome lattices exhibit a rigidity percolation transition when the number of braces is $\sim LlnL$, where $L$ is the linear size of the lattice. This transition exhibits features of both first-order and second-order transitions: The whole lattice becomes rigid at the transition, and a diverging length scale also exists. In contrast, we find that the rigidity percolation transition in the generic lattices occur when the number of braces is very close to the number obtained from Maxwell's law for floppy modes, which is $\sim L$. The transition in generic lattices is a very sharp first-order-like transition, at which the addition of one brace connects all small rigid regions in the bulk of the lattice, leaving only floppy modes on the edge. We characterize these transitions using numerical simulations and develop analytic theories capturing each transition. Our results relate to other interesting problems, including jamming and bootstrap percolation.

We present an efficient method for simulating a stationary Gaussian noise with an arbitrary covariance function, and then we study numerically the impact of time-correlated noise on the time evolution of a ($1+1$)-dimensional generalized Langevin equation by comparing also to analytical results. Finally, we apply our method to the generalized Langevin equation with an external harmonic and double-well potential.

We study the thermodynamic properties of the three-dimensional Blume-Capel model on the simple cubic lattice by means of computer simulations. In particular, we implement a parallelized variant of the multicanonical approach and perform simulations by keeping a constant temperature and crossing the phase boundary along the crystal-field axis. We obtain numerical data for several temperatures in both the first- and second-order regime of the model. Finite-size scaling analyses provide us with transition points and the dimensional scaling behavior in the numerically demanding first-order regime, as well as a clear verification of the expected Ising universality in the respective second-order regime. Finally, we discuss the scaling behavior in the vicinity of the tricritical point.

The paper considers heat conduction in a model chain of composite particles with hard core and elastic external shell. Such model mimics three main features of realistic interatomic potentials—hard repulsive core, quasilinear behavior in a ground state, and possibility of dissociation. It has become clear recently that this latter feature has crucial effect on convergence of the heat conduction coefficient in thermodynamic limit. We demonstrate that in one-dimensional chain of elastic particles with hard core the heat conduction coefficient also converges, as one could expect. Then we explore effect of dimensionality on the heat transport in this model. For this sake, longitudinal and transversal motions of the particles are allowed in a long narrow channel. With varying width of the channel, we observe sharp transition from “one-dimensional” to “two-dimensional” behavior. Namely, the heat conduction coefficient drops by about order of magnitude for relatively small widening of the channel. This transition is not unique for the considered system. Similar phenomenon of transition to quasi-1D behavior with growth of aspect ratio of the channel is observed also in a gas of densely packed hard (billiard) particles, both for two- and three-dimensional cases. It is the case despite the fact that the character of transition in these two systems is not similar, due to different convergence properties of the heat conductivity. In the billiard model, the divergence pattern of the heat conduction coefficient smoothly changes from logarithmic to power-like law with increase of the length.

We study a molecular engine constituted by a gas of $N\sim {10}^2$ molecules enclosed between a massive piston and a thermostat. The force acting on the piston and the temperature of the thermostat are cyclically changed with a finite period $\tau$. In the adiabatic limit $\tau \to \infty$, even for finite size $N$, the average work and heat reproduce the thermodynamic values, recovering the Carnot result for the efficiency. The system exhibits a stall time ${\tau }^{*}$ where the net work is zero: for $\tau <{\tau }^{*}$ it consumes work instead of producing it, acting as a refrigerator or as a heat sink. At $\tau >{\tau }^{*}$ the efficiency at maximum power is close to the Curzorn-Ahlborn limit. The fluctuations of work and heat display approximatively a Gaussian behavior. Based upon kinetic theory, we develop a three-variables Langevin model in which the piston's position and velocity are linearly coupled together with the internal energy of the gas. The model reproduces many of the system's features, such as the inversion of the work's sign, the efficiency at maximum power, and the approximate shape of the fluctuations. A further simplification in the model allows us to compute analytically the average work, explaining its nontrivial dependence on $\tau$.

By means of extensive contact dynamics simulations, we analyze the combined effects of polydispersity both in particle size and in particle shape, defined as the degree of shape irregularity, on the shear strength and microstructure of sheared granular materials composed of pentagonal particles. We find that the shear strength is independent of the size span, but unexpectedly, it declines with increasing shape polydispersity. At the same time, the solid fraction is an increasing function of both the size span and the shape polydispersity. Hence, the densest and loosest packings have the same shear strength. At the scale of the particles and their contacts, we analyze the connectivity of particles, force transmission, and friction mobilization as well as their anisotropies. We show that stronger forces are carried by larger particles and propped by an increasing number of small particles. The independence of shear strength with regard to size span is shown to be a consequence of contact network self-organization, with the falloff of contact anisotropy compensated by increasing force anisotropy.

We consider a model for two delay-coupled optoelectronic oscillators under positive delayed feedback as prototypical to study the conditions for synchronization of asymmetric square-wave oscillations, for which the duty cycle is not half of the period. We show that the scenario arising for positive feedback is much richer than with negative feedback. First, it allows for the coexistence of multiple in- and out-of-phase asymmetric periodic square waves for the same parameter values. Second, it is tunable: The period of all the square-wave periodic pulses can be tuned with the ratio of the delays, and the duty cycle of the asymmetric square waves can be changed with the offset phase while the total period remains constant. Finally, in addition to the multiple in- and out-of-phase periodic square waves, low-frequency periodic asymmetric solutions oscillating in phase may coexist for the same values of the parameters. Our analytical results are in agreement with numerical simulations and bifurcation diagrams obtained by using continuation techniques.

Experimental evidence of the existence of an intermediate region between a capacitively coupled plasma and the collisional space-charge region at its borders is presented. This proof is generated by monitoring—in an airplane carrying out parabolic flights—the trajectory of plasma-confined microparticles. Based on only primary data and without the need for a sophisticated model, our analysis concludes a sharply marked transition from the sheath region into another region with a significantly lower—yet nonzero—space-charge density, i.e., a region which is often called the presheath.

Aiming at the study of critical phenomena in the presence of boundaries with a nontrivial shape we discuss how lattices with an adaptive lattice spacing can be implemented. Since the parameters of the Hamiltonian transform nontrivially under changes of the length scale, adapting the lattice spacing is much more difficult than in the case of the numerical solution of partial differential equations, where this method is common practice.

Here we shall focus on the universality class of the three-dimensional Ising model. Our starting point is the improved Blume-Capel model on the simple cubic lattice. In our approach, the system is composed of sectors with lattice spacing $a,2a,4a,....$ We work out how parts of the lattice with lattice spacing $a$ and $2a$, respectively, can be coupled in a consistent way. Here, we restrict ourself to the case where the boundary between the sectors is perpendicular to one of the lattice axes. Based on the theory of defect planes one expects that it is sufficient to tune the coupling between these two regions. To this end we perform a finite-size scaling study. However, first numerical results show that slowly decaying corrections remain. It turns out that these corrections can be removed by adjusting the strength of the couplings within the boundary layers. As a benchmark, we simulate films with strongly symmetry breaking boundary conditions. We determine the magnetization profile and the thermodynamic Casimir force. For our largest thickness $L_0=64.5$, we find that results obtained for the homogeneous system are nicely reproduced.

We perform molecular dynamics (MD) simulations of the crystallization process in binary Lennard-Jones systems during heating and cooling to investigate atomic-scale crystallization kinetics in glass-forming materials. For the cooling protocol, we prepared equilibrated liquids above the liquidus temperature $T_l$ and cooled each sample to zero temperature at rate $R_c$. For the heating protocol, we first cooled equilibrated liquids to zero temperature at rate $R_p$ and then heated the samples to temperature $T>T_l$ at rate $R_h$. We measured the critical heating and cooling rates $R_h^{*}$ and $R_c^{*}$, below which the systems begin to form a substantial fraction of crystalline clusters during the heating and cooling protocols. We show that $R_h^{*}>R_c^{*}$ and that the asymmetry ratio $R_h^{*}/R_c^{*}$ includes an intrinsic contribution that increases with the glass-forming ability (GFA) of the system and a preparation-rate dependent contribution that increases strongly as $R_p\to R_c^{*}$ from above. We also show that the predictions from classical nucleation theory (CNT) can qualitatively describe the dependence of the asymmetry ratio on the GFA and preparation rate $R_p$ from the MD simulations and results for the asymmetry ratio measured in Zr- and Au-based bulk metallic glasses (BMG). This work emphasizes the need for and benefits of an improved understanding of crystallization processes in BMGs and other glass-forming systems.