The self-consistent probabilistic approach has proven itself powerful in studying the percolation behavior of interdependent or multiplex networks without tracking the percolation process through each cascading step. In order to understand how directed dependency links impact criticality, we employ this approach to study the percolation properties of networks with both undirected connectivity links and directed dependency links. We find that when a random network with a given degree distribution undergoes a second-order phase transition, the critical point and the unstable regime surrounding the second-order phase transition regime are determined by the proportion of nodes that do not depend on any other nodes. Moreover, we also find that the triple point and the boundary between first- and second-order transitions are determined by the proportion of nodes that depend on no more than one node. This implies that it is maybe general for multiplex network systems, some important properties of phase transitions can be determined only by a few parameters. We illustrate our findings using Erdős-Rényi networks.

We study the relative translation of two arbitrarily shaped objects, caused by their hydrodynamic interaction as they are forced through a viscous fluid in the limit of zero Reynolds number. It is well known that in the case of two rigid spheres in an unbounded fluid, the hydrodynamic interaction does not produce relative translation. More generally, such an effective pair-interaction vanishes in configurations with spatial inversion symmetry; for example, an enantiomorphic pair in mirror image positions has no relative translation. We show that the breaking of inversion symmetry by boundaries of the system accounts for the interactions between two spheres in confined geometries, as observed in experiments. The same general principle also provides new predictions for interactions in other object configurations near obstacles. We examine the time-dependent relative translation of two self-aligning objects, extending the numerical analysis of our preceding publication [Goldfriend, Diamant, and Witten, Phys. Fluids 27, 123303 (2015)]. The interplay between the orientational interaction and the translational one, in most cases, leads over time to repulsion between the two objects. The repulsion is qualitatively different for self-aligning objects compared to the more symmetric case of uniform prolate spheroids. The separation between the two objects increases with time $t$ as $t^{1/3}$ in the former case, and more strongly, as $t$, in the latter.

The shearless mixing layer is generated from the interaction of two homogeneous isotropic turbulence (HIT) fields with different integral scales ${\ell }_1$ and ${\ell }_2$ and different turbulent kinetic energies $E_1$ and $E_2$. In this study, the sensitivity of temporal evolutions of two-dimensional, incompressible shearless mixing layers to the parametric variations of ${\ell }_1/{\ell }_2$ and $E_1/E_2$ is investigated. The sensitivity methodology is based on the nonintrusive approach; using direct numerical simulation and generalized polynomial chaos expansion. The analysis is carried out at ${\mathrm{Re}}_{{\ell }_1}=90$ for the high-energy HIT region and different integral length scale ratios $1/4\le {\ell }_1/{\ell }_2\le 4$ and turbulent kinetic energy ratios $1\le E_1/E_2\le 30$. It is found that the most influential parameter on the variability of the mixing layer evolution is the turbulent kinetic energy while variations of the integral length scale show a negligible influence on the flow field variability. A significant level of anisotropy and intermittency is observed in both large and small scales. In particular, it is found that large scales have higher levels of intermittency and sensitivity to the variations of ${\ell }_1/{\ell }_2$ and $E_1/E_2$ compared to the small scales. Reconstructed response surfaces of the flow field intermittency and the turbulent penetration depth show monotonic dependence on ${\ell }_1/{\ell }_2$ and $E_1/E_2$. The mixing layer growth rate and the mixing efficiency both show sensitive dependence on the initial condition parameters. However, the probability density function of these quantities shows relatively small solution variations in response to the variations of the initial condition parameters.

We use microscopic density functional theory to study filling transitions in systems with long-ranged wall-fluid and short-ranged fluid-fluid forces occurring in a right-angle wedge. By changing the strength of the wall-fluid interaction we can induce both wetting and filling transitions over a wide range of temperatures and study the order of these transitions. At low temperatures we find that both wetting and filling transitions are first order in keeping with predictions of simple local effective Hamiltonian models. However close to the bulk critical point the filling transition is observed to be continuous even though the wetting transition remains first order and the wetting binding potential still exhibits a small activation barrier. The critical singularities for adsorption for the continuous filling transitions depend on whether retarded or nonretarded wall-fluid forces are present and are in excellent agreement with predictions of effective Hamiltonian theory even though the change in the order of the transition was not anticipated.

We study coarse-grained entropy production in an asymmetric random walk system on a periodic one-dimensional lattice. In coarse-grained systems, the original dynamics are unavoidably destroyed, but the coarse-grained entropy production is not hidden below the critical time-scale separation. The hidden entropy production is rapidly increasing near the critical time-scale separation.

Optimizing information transmission across a network is an essential task for controlling and manipulating generic information-processing systems. Here, we show how topological amplification effects in scale-free networks of signaling devices are optimally enhanced when the impulse transmitted by periodic external signals (time integral over two consecutive zeros) is maximum. This is demonstrated theoretically by means of a star-like network of overdamped bistable systems subjected to generic zero-mean periodic signals and confirmed numerically by simulations of scale-free networks of such systems. Our results show that the enhancer effect of increasing values of the signal's impulse is due to a correlative increase of the energy transmitted by the periodic signals, while it is found to be resonant-like with respect to the topology-induced amplification mechanism.

The influence of an applied shear on the planar crystal-melt interface is modeled by a nonlinear stochastic partial differential equation of the interface fluctuations. A feature of this theory is the asymmetric destruction of interface fluctuations due to advection of the crystal protrusions on the liquid side of the interface only. We show that this model is able to qualitatively reproduce the nonequilibrium coexistence line found in simulations. The impact of shear on spherical clusters is also addressed.

Microbeam resonant x-ray scattering experiments recently revealed the sequential emergence of electric-field-induced subphases (stable states) with exceptionally large unit cells consisting of 12 and 15 smectic layers. We explain the emergence of the field-induced subphases by the quasimolecular model based on the Emelyanenko-Osipov long-range interlayer interactions (LRILIs) together with our primitive way of understanding the frustration in clinicity using the $q_E$ number defined as $q_E=|\left[R\right]-\left[L\right]|/(\left[R\right]+\left[L\right])$; here $\left[R\right]$ and $\left[L\right]$ refer to the numbers of smectic layers with directors tilted to the right and to the left, respectively, in the unit cell of a field-induced subphase. We show that the model actually stabilizes the field-induced subphases with characteristic composite unit cells consisting of several blocks, each of which is originally a ferrielectric three-layer unit cell stabilized by the LRILIs, but some of which would be modified to become ferroelectric by an applied electric field. In a similar line of thought, we also try to understand the puzzling electric-field-induced birefringence data in terms of the LRILIs.

The laminar-turbulent (LT) transition of dilute polymer solutions is of great interest not only for the complex transition dynamics itself, but also for its potential link to the maximum drag reduction (MDR) phenomenon. We present an in-depth investigation of the edge state (ES), an asymptotic solution on the LT boundary, in viscoelastic channel flow. For given $\text{Re}$ and simulation domain size, mean flow statistics of the ES do not vary with the introduction of polymers, proving that there is a region of turbulent states not susceptible to polymer drag reduction effects. The dynamics of the ES features low-frequency fluctuations and in the longer domains we studied it is nearly periodic with regular bursts of turbulent activities separated by extended quiescent periods. Its flow field is dominated by elongated vortices and streaks, with very weak extensional and rotational flow motions. Polymer stretching is almost exclusively contributed by the mean shear and polymer-turbulence interaction is minimal. Flow structures and the kinematics of the ES match hibernating turbulence, an MDR-like phase intermittently occurring in turbulent dynamics. Its observation now seems to result from recurrent visits to certain parts of the ES. The ES offers explanations for the existence and universality of MDR, the quantitative magnitude of which, however, still remains unsolved.

In Taylor's theory, the longitudinal dispersion in turbulent pipe flows approaches, on long time scales, a diffusive behavior with a constant diffusivity $K_L$, which depends empirically on the Reynolds number Re. We show that the dependence on Re can be determined from the turbulent energy spectrum. By using the intimate connection between the friction factor and the longitudinal dispersion in wall-bounded turbulence, we predict different asymptotic scaling laws of $K_L$(Re) depending on the different turbulent cascades in two-dimensional turbulence. We also explore numerically the $K_L$(Re) dependence in turbulent channel flows with smooth and rough walls using a lattice Boltzmann method.

Scalar transport by turbulent flows is best described in terms of Lagrangian parcel motions. Here we measure the Eulerian distance travel along Lagrangian trajectories in a simple point vortex flow to determine the probabilistic impulse response function for scalar transport in the absence of molecular diffusion. As expected, the mean squared Eulerian displacement scales ballistically at very short times and diffusively for very long times, with the displacement distribution at any given time approximating that of a random walk. However, significant deviations in the displacement distributions from Rayleigh are found. The probability of long distance transport is reduced over inertial range time scales due to spatial and temporal intermittency. This can be modeled as a series of trapping events with durations uniformly distributed below the Eulerian integral time scale. The probability of long distance transport is, on the other hand, enhanced beyond that of the random walk for both times shorter than the Lagrangian integral time and times longer than the Eulerian integral time. The very short-time enhancement reflects the underlying Lagrangian velocity distribution, while that at very long times results from the spatial and temporal variation of the flow at the largest scales. The probabilistic impulse response function, and with it the expectation value of the scalar concentration at any point in space and time, can be modeled using only the evolution of the lowest spatial wave number modes (the mean and the lowest harmonic) and an eddy based constrained random walk that captures the essential velocity phase relations associated with advection by vortex motions. Preliminary examination of Lagrangian tracers in three-dimensional homogeneous isotropic turbulence suggests that transport in that setting can be similarly modeled.

We construct a fluid-dynamical model for the flow of a solution with a free surface at which surface tension acts. This model can describe both classical surfactants, which decrease the surface tension of the solution relative to that of the pure solvent, and antisurfactants (such as many salts when added to water, and small amounts of water when added to alcohol) which increase it. We demonstrate the utility of the model by considering the linear stability of an infinitely deep layer of initially quiescent fluid. In particular, we predict the occurrence of an instability driven by surface-tension gradients, which occurs for antisurfactant, but not for surfactant, solutions.

In this paper we study asymptotic behavior of a dynamical functional for an $\alpha$-stable autoregressive fractionally integrated moving average (ARFIMA) process. We find an analytical formula for this important statistics and show its usefulness as a diagnostic tool for ergodic properties. The obtained results point to the very fast convergence of the dynamical functional and show that even for short trajectories one may obtain reliable conclusions on the ergodic properties of the ARFIMA process. Moreover we use the obtained theoretical results to illustrate how the dynamical functional statistics can be used in the verification of the proper model for an analysis of some biophysical experimental data.

An irreversible Markov-chain Monte Carlo (MCMC) algorithm with skew detailed balance conditions originally proposed by Turitsyn et al. is extended to general discrete systems on the basis of the Metropolis-Hastings scheme. To evaluate the efficiency of our proposed method, the relaxation dynamics of the slowest mode and the asymptotic variance are studied analytically in a random walk on one dimension. It is found that the performance in irreversible MCMC methods violating the detailed balance condition is improved by appropriately choosing parameters in the algorithm.

What kind of one-dimensional modulated nematic structures (ODMNS) can form nonchiral and chiral bent-core and dimeric materials? Here, using the Landau–de Gennes theory of nematics, extended to account for molecular steric polarization, we study a possibility of formation of ODMNS, both in nonchiral and intrinsically chiral liquid crystalline materials. Besides nematic and cholesteric phases, we find four bulk ODMNS for nonchiral materials, two of which, to the best of our knowledge, have not been reported so far. These two structures are longitudinal $\left(N_{\text{LP}}\right)$ and transverse $\left(N_{\text{TP}}\right)$ periodic waves where the polarization field being periodic in one dimension stays parallel and perpendicular, respectively, to the wave vector. The other two phases are the twist-bend nematic phase $\left(N_{\text{TB}}\right)$ and the splay-bend nematic phase $\left(N_{\text{SB}}\right)$, but their fine structure appears more complex than that considered so far. The presence of molecular chirality converts nonchiral $N_{\text{TP}}$ and $N_{\text{SB}}$ into new $N_{\text{TB}}$ phases. Surprisingly, the nonchiral $N_{\text{LP}}$ phase can stay stable even in the presence of intrinsic chirality.

The transformative influence of compressibility on the Kelvin-Helmholtz instability (KHI) at the interface between two fluid streams of different velocities is explicated. When the velocity difference is small (subsonic), shear effects dominate the interface flow dynamics causing monotonic roll-up of vorticity and mixing between the two streams leading to the KHI. We find that at supersonic speed differentials, compressibility forces the dominance of dilatational (acoustic) rather than shear dynamics at the interface. Within this dilatational interface layer, traveling pressure waves cause the velocity perturbations to become oscillatory. We demonstrate that the oscillatory fluid motion reverses vortex roll-up and segregates the two streams leading to KHI suppression. Analysis and illustrations of the compressibility-induced suppression mechanism are presented.

We investigate the large-time scaling regimes arising from a variety of metastable structures in a chain of Ising spins with both first- and second-neighbor couplings while subject to Kawasaki dynamics. Depending on the ratio and sign of these former, different dynamic exponents are suggested by finite-size scaling analyses of relaxation times. At low but nonzero temperatures these are calculated via exact diagonalizations of the evolution operator in finite chains under several activation barriers. In the absence of metastability the dynamics is always diffusive.

Although the set of permutation symmetries of a complex network could be very large, few of them give rise to stable synchronous patterns. Here we present a general framework and develop techniques for controlling synchronization patterns in complex network of coupled chaotic oscillators. Specifically, according to the network permutation symmetry, we design a small-size and weighted network, namely the control network, and use it to control the large-size complex network by means of pinning coupling. We argue mathematically that for any of the network symmetries, there always exists a critical pinning strength beyond which the unstable synchronous pattern associated to this symmetry can be stabilized. The feasibility of the control method is verified by numerical simulations of both artificial and real-world networks and demonstrated experimentally in systems of coupled chaotic circuits. Our studies show the controllability of synchronous patterns in complex networks of coupled chaotic oscillators.

We consider the mechanism of formation of isolated localized wave structures in the diatomic Fermi-Pasta-Ulam (FPU) model. Using a singular multiscale asymptotic analysis in the limit of high mass mismatch between the alternating elements, we obtain the typical slow-fast time scale separation and formulate the Fredholm orthogonality condition approximating a sequence of mass ratios supporting the formation of solitary waves in the general type of diatomic FPU models. This condition is made explicit in the case of a diatomic Toda lattice. Results of numerical integration of the full diatomic Toda lattice equations confirm the formation of these genuinely localized wave structures at special values of the mass ratio that are close to the analytical predictions when the ratio is sufficiently small.