We have proposed a new method for solving the problem of ship waves excited on the surface of a nonviscous liquid by a submerged object that moves at a variable speed. As a first application of this method, we have obtained a new solution to the classic problem of ship waves generated by a submerged ball that moves rectilinearly with constant velocity parallel to the equilibrium surface of the liquid. For this example, we have derived asymptotic expressions describing the vertical displacement of the liquid surface in the limit of small and large values of the Froude number. The exact solution is presented in the form of two terms, each of which is reduced to one-dimensional integrals. One term describes the “Bernoulli hump” and another term the “Kelvin wedge.” As a second example, we considered vertical oscillation of the submerged ball. In this case, the solution leads to the calculation of one-dimensional integral and describes surface waves propagating from the epicenter above the ball.

An analytical expression is derived for the secondary Bjerknes force experienced by two cylindrical bubbles in a microfluidic channel with planar elastic walls. The derived expression takes into account that the bubbles generate two types of scattered acoustic waves: bulk waves that propagate in the fluid gap with the speed of sound and Lamb-type surface waves that propagate at the fluid-wall interfaces with a much lower speed than that of the bulk waves. It is shown that the surface waves cause the bubbles to form a bound pair in which the equilibrium interbubble distance is determined by the wavelength of the surface waves, which is much smaller than the acoustic wavelength. Comparison of theoretical and experimental results demonstrates good agreement.

We present a method for applying a class of velocity-dependent forces within a multicomponent lattice Boltzmann equation simulation that is designed to recover continuum regime incompressible hydrodynamics. This method is applied to the problem, in two dimensions, of constraining to uniformity the tangential velocity of a vesicle membrane implemented within a recent multicomponent lattice Boltzmann simulation method, which avoids the use of Lagrangian boundary tracers. The constraint of uniform tangential velocity is carried by an additional contribution to an immersed boundary force, which we derive here from physical arguments. The result of this enhanced immersed boundary force is to apply a physically appropriate boundary condition at the interface between separated lattice fluids, defined as that region over which the phase-field varies most rapidly. Data from this enhanced vesicle boundary method are in agreement with other data obtained using related methods [e.g., T. Krüger, S. Frijters, F. Günther, B. Kaoui, and J. Harting, Eur. Phys. J. 222, 177 (2013)] and underscore the importance of a correct vesicle membrane condition.

We present a systematic study on the effects of tangential-type boundary condition discontinuities on the accuracy of the lattice Boltzmann equation (LBE) method for Dirichlet and Neumann problems in heat and mass transfer modeling. The second-order accurate boundary condition treatments for continuous Dirichlet and Neumann problems are directly implemented for the corresponding discontinuous boundary conditions. Results from three numerical tests, including both straight and curved boundaries, are presented to show the accuracy and order of convergence of the LBE computations. Detailed error assessments are conducted for the interior temperature or concentration (denoted as a scalar ϕ) and the interior derivatives of ϕ for both types of boundary conditions, for the boundary flux in the Dirichlet problem and for the boundary ϕ values in the Neumann problem. When the discontinuity point on the straight boundary is placed at the center of the unit lattice in the Dirichlet problem, it yields only first-order accuracy for the interior distribution of ϕ, first-order accuracy for the boundary flux, and zeroth-order accuracy for the interior derivatives compared with the second-order accuracy of all quantities of interest for continuous boundary conditions. On the lattice scale, the LBE solution for the interior derivatives near the singularity is largely independent of the resolution and correspondingly the local distribution of the absolute errors is almost invariant with the changing resolution. For Neumann problems, when the discontinuity is placed at the lattice center, second-order accuracy is preserved for the interior distribution of ϕ; and a “superlinear” convergence order of 1.5 for the boundary ϕ values and first-order accuracy for the interior derivatives are obtained. For straight boundaries with the discontinuity point arbitrarily placed within the lattice and curved boundaries, the boundary flux becomes zeroth-order accurate for Dirichlet problems; and all three quantities, including the interior and boundary ϕ values and the interior derivatives, are only first-order accurate for Neumann problems.

A universal characteristic of mesoscale turbulence in active suspensions is the emergence of a typical vortex length scale, distinctly different from the scale invariance of turbulent high-Reynolds number flows. Collective length-scale selection has been observed in bacterial fluids, endothelial tissue, and active colloids, yet the physical origins of this phenomenon remain elusive. Here, we systematically derive an effective fourth-order field theory from a generic microscopic model that allows us to predict the typical vortex size in microswimmer suspensions. Building on a self-consistent closure condition, the derivation shows that the vortex length scale is determined by the competition between local alignment forces, rotational diffusion, and intermediate-range hydrodynamic interactions. Vortex structures found in simulations of the theory agree with recent measurements in Bacillus subtilis suspensions. Moreover, our approach yields an effective viscosity enhancement (reduction), as reported experimentally for puller (pusher) microorganisms.

Micron-sized colloidal spheres that are dispersed in an isotropic-nematic biphasic host suspension of charged rods (fd-virus particles) are shown to spontaneously form dimers, which exhibit a synchronized oscillatory motion. Dimer formation is not observed in the monophase of isotropic and nematic suspensions. The synchronized oscillations of dimers are connected to the inhomogeneous state of the host suspension of charged rods (fd viruses) where nematic domains are in coexistence with isotropic regions. The synchronization of oscillations occurs in bulk states, in the absence of an external field. With a low field strength of an applied electric field, the synchronization is rather reduced, but it recovers again when the field is turned off. In this Rapid Communication, we report this observation as an example of the strange attractor, occurring in the mixture of PS (polystyrene) dimers in an isotropic-nematic coexistence biphasic fd-virus network. Furthermore, we highlight that the synchronization of PS-dimer oscillations is the result of a global bifurcation diagram, driven by a delicate balance between the short-attractive “twisted” interaction of PS dimers and long-ranged electrostatic repulsive interactions of charged fd rods. The interest is then in the local enhancement of “twist-nematic” elasticity in reorientation of the dimer oscillations. An analysis of image-time correlations is provided with the data movies and Fourier transforms of averaged orientations for the synchronized oscillations of dimers in the biphasic $I-N$ coexistence concentration of charged fd-virus suspensions.

We investigate the influence of time-varying environmental noise, i.e., temporal disorder, on the nonequilibrium phase transition of the contact process. Combining a real-time renormalization group, scaling theory, and large scale Monte-Carlo simulations in one and two dimensions, we show that the temporal disorder gives rise to an exotic critical point. At criticality, the effective noise amplitude diverges with increasing time scale, and the probability distribution of the density becomes infinitely broad, even on a logarithmic scale. Moreover, the average density and survival probability decay only logarithmically with time. This infinite-noise critical behavior can be understood as the temporal counterpart of infinite-randomness critical behavior in spatially disordered systems, but with exchanged roles of space and time. We also analyze the generality of our results, and we discuss potential experiments.

The dynamics of domain walls in optical bistable systems with pump and loss is considered. It is shown that an oscillating component of the pump affects the average drift velocity of the domain walls. The cases of harmonic and biharmonic pumps are considered. It is demonstrated that in the case of biharmonic pulse the velocity of the domain wall can be controlled by the mutual phase of the harmonics. The analogy between this phenomenon and the ratchet effect is drawn. Synchronization of the moving domain walls by the oscillating pump in discrete systems is studied and discussed.

Surface and interfacial creases induced by biological growth are common types of instability in soft biological tissues. This study focuses on the criteria for the onset of surface and interfacial creases as well as their morphological evolution in a growing bilayer soft tube within a confined environment. Critical growth ratios for triggering surface and interfacial creases are investigated both analytically and numerically. Analytical interpretations provide preliminary insights into critical stretches and growth ratios for the onset of instability and formation of both surface and interfacial creases. However, the analytical approach cannot predict the evolution pattern of the model after instability; therefore nonlinear finite element simulations are carried out to replicate the poststability morphological patterns of the structure. Analytical and computational simulation results demonstrate that the initial geometry, growth ratio, and shear modulus ratio of the layers are the most influential factors to control surface and interfacial crease formation in this soft tubular bilayer. The competition between the stretch ratios in the free and interfacial surfaces is one of the key driving factors to determine the location of the first crease initiation. These findings may provide some fundamental understanding in the growth modeling of tubular biological tissues such as esophagi and airways as well as offering useful clues into normal and pathological functions of these tissues.

We study the spreading of a bacterial colony undergoing turbulentlike collective motion. We present two minimalistic models to investigate the interplay between population growth and coherent structures arising from turbulence. Using direct numerical simulation of the proposed models we find that turbulence has two prominent effects on the spatial growth of the colony: (a) the front speed is enhanced, and (b) the front gets crumpled. Both these effects, which we highlight by using statistical tools, are markedly different in our two models. We also show that the crumpled front structure and the passive scalar fronts in random flows are related in certain regimes.

We discuss the application of fluorescence correlation spectroscopy (FCS) for characterization of anomalous diffusion of tracer particles in crowded environments. While the fact of anomaly may be detected by the standard fitting procedure, the value of the exponent $\alpha$ of anomalous diffusion may be not reproduced correctly for non-Gaussian anomalous diffusion processes. The important information is however contained in the asymptotic behavior of the fluorescence autocorrelation function at long and at short times. Thus, analysis of the short-time behavior gives reliable values of $\alpha$ and of lower moments of the distribution of particles' displacement, which allows us to confirm or reject its Gaussian nature. The method proposed was tested on the FCS data obtained in artificial crowded fluids and in living cells.

In irreversible aggregation processes droplets or polymers of microscopic size successively coalesce until a large cluster of macroscopic scale forms. This gelation transition is widely believed to be self-averaging, meaning that the order parameter (the relative size of the largest connected cluster) attains well-defined values upon ensemble averaging with no sample-to-sample fluctuations in the thermodynamic limit. Here, we report on anomalous gelation transition types. Depending on the growth rate of the largest clusters, the gelation transition can show very diverse patterns as a function of the control parameter, which includes multiple stochastic discontinuous transitions, genuine non-self-averaging and ultraslow convergence of the transition point. Our framework may be helpful in understanding and controlling gelation.

The emergence of clustering and coarsening in crowded ensembles of self-propelled agents is studied using a lattice model in one dimension. The persistent exclusion process, where particles move at directions that change randomly at a low tumble rate $\alpha$, is extended allowing sites to be occupied by more than one particle, with a maximum $n_{\text{max}}$ per site. Three phases are distinguished. For $n_{\text{max}}=1$ a gas of clusters form, with sizes distributed exponentially and no coarsening takes place. For $n_{\text{max}}\ge 3$ and small values of $\alpha$, coarsening takes place and few large clusters appear, with a large fraction of the total number of particles in them. In the same range of $n_{\text{max}}$ but for larger values of $\alpha$, a gas phase where a negligible fraction of particles takes part of clusters. Finally, $n_{\text{max}}=2$ corresponds to a crossover phase. The character of the transitions between phases is studied extending the model to allow $n_{\text{max}}$ to take real values and jumps to an occupied site are probabilistic. The transition from the gas of clusters to the coarsening phase is continuous and the mass of the large clusters grows continuously when varying the maximum occupancy, and the crossover found corresponds to values close to the transition. The second transition, from the coarsening to the gaseous phase, can be either continuous or discontinuous depending on the parameters, with a critical point separating both cases.

One possibility to adjust material properties to a specific need is to embed units of one substance into a matrix of another substance. Even materials that are readily tunable during operation can be generated in this way. In (visco)elastic substances, both the matrix material as well as the inclusions and/or their immediate environment can be dynamically deformed. If the typical dynamic response time of the inclusions and their surroundings approach the macroscopic response time, their deformation processes need to be included into a dynamic macroscopic characterization. Along these lines, we present a hydrodynamic description of (visco)elastic composite materials. For this purpose, additional strain variables reflect the state of the inclusions and their immediate environment. These additional strain variables in general are not set by a coarse-grained macroscopic displacement field. Apart from that, during our derivation, we also include the macroscopic variables of relative translations and relative rotations that were previously introduced in different contexts. As a central point, our approach reveals and classifies the importance of a macroscopic variable termed relative strains. We analyze two simplified minimal example geometries as an illustration.

We quantify the disorder-induced heating (DIH) of ultracold neutral plasmas (UCNPs) created from cold atoms in optical lattices with partial filling fractions, using a conservation of energy model involving the spatial correlations of the initial state and the equation of state in thermal equilibrium for a one-component plasma. We show, for experimentally achievable filling fractions, that the ionic Coulomb coupling parameter could be increased to a degree comparable to other proposed DIH-mitigation schemes. Molecular dynamics simulations were performed with compensation for finite-size and periodic boundary effects, which agree with calculations using the model. Reduction of DIH using optical lattices will allow for the study of strongly coupled plasma physics using low-density, low-temperature, laboratory-based plasmas, and lead to improved brightness in UCNP-based cold electron and ion beams, where DIH is otherwise a fundamental limitation to beam focal sizes and diffraction imaging capability.

We study the phenomenon of weak ergodicity breaking for a class of globally correlated random walk dynamics defined over a finite set of states. The persistence in a given state or the transition to another one depends on the whole previous temporal history of the system. A set of waiting time distributions, associated to each state, sets the random times between consecutive steps. Their mean value is finite for all states. The probability density of time-averaged observables is obtained for different memory mechanisms. This statistical object explicitly shows departures between time and ensemble averages. While the residence time in each state may have a divergent mean value, we demonstrate that this condition is in general not necessary for breaking ergodicity. Hence, we conclude that global memory effects are an alternative mechanism able to induce ergodicity breaking without involving power-law statistics. Analytical and numerical calculations support these results.