A new general class of exact solutions is presented for the time evolution of a bubble of arbitrary initial shape in a Hele-Shaw cell when surface tension effects are neglected. These solutions are obtained by conformal mapping the viscous flow domain to an annulus in an auxiliary complex plane. It is then demonstrated that the only stable fixed point (attractor) of the nonsingular bubble dynamics corresponds precisely to the selected pattern. This thus shows that, contrary to the established theory, bubble selection in a Hele-Shaw cell does not require surface tension. The solutions reported here significantly extend previous results for a simply connected geometry (finger) to a doubly connected one (bubble). We conjecture that the same selection rule without surface tension holds for Hele-Shaw flows of arbitrary connectivity.

Semiflexible polymers can generate a range of filamentous networks significantly different in structure from those seen in conventional polymer solutions. Our coarse-grained simulations with an implicit cross-linker potential show that networks of branching bundles, knotted morphologies, and structural chirality can be generated by a generalized approach independent of specific cross-linkers. Network structure depends primarily on filament flexibility and separation, with significant connectivity increase after percolation. Results should guide the design of engineered semiflexible polymers.

In this article, we present a framework for investigating the order-disorder transition in lipid droplets using the standard Ising model. While a single lipid droplet is itself a complex system whose constituent cholesteryl esters each possesses many degrees of freedom, we present justification for using this effective approach to isolate the underlying physics. It is argued that the behavior of the esters confined within lipid droplets is significantly different from that of a bulk system of similar esters, which is adequately described by continuum mean-field theory in the thermodynamic limit. When the droplet's shell is modeled as an elastic membrane, a simple picture emerges for a transition between two ordered phases within the core which is tuned by the strength of interactions between the esters. Triglyceride concentration is proposed as a variable which strongly influences the strength of interactions between cholesteryl esters within droplets. The possible relevance of this mechanism to the well known atherogenic nature of small low-density lipoprotein particles is discussed in detail.

Propulsion at microscopic scales is often achieved through propagating traveling waves along hairlike organelles called flagella. Taylor's two-dimensional swimming sheet model is frequently used to provide insight into problems of flagellar propulsion. We derive numerically the large-amplitude wave form of the two-dimensional swimming sheet that yields optimum hydrodynamic efficiency: the ratio of the squared swimming speed to the rate-of-working of the sheet against the fluid. Using the boundary element method, we show that the optimal wave form is a front-back symmetric regularized cusp that is $25\%$ more efficient than the optimal sine wave. This optimal two-dimensional shape is smooth, qualitatively different from the kinked form of Lighthill's optimal three-dimensional flagellum, not predicted by small-amplitude theory, and different from the smooth circular-arc-like shape of active elastic filaments.

Random walks (RW) of particles adsorbed in the internal walls of porous deposits produced by ballistic-type growth models are studied. The particles start at the external surface of the deposits and enter their pores in order to simulate an external flux of a species towards a porous solid. For short times, the walker concentration decays as a stretched exponential of the depth $z$, but a crossover to long-time normal diffusion is observed in most samples. The anomalous concentration profile remains at long times in very porous solids if the walker steps are restricted to nearest neighbors and is accompanied with subdiffusion features. These findings are correlated with a decay of the explored area with $z$. The study of RW of tracer particles left at the internal part of the solid rules out an interpretation by diffusion equations with position-dependent coefficients. A model of RW in a tube of decreasing cross section explains those results by showing long crossovers from an effective subdiffusion regime to an asymptotic normal diffusion. The crossover position and density are analytically calculated for a tube with area decreasing exponentially with $z$ and show good agreement with numerical data. The anomalous decay of the concentration profile is interpreted as a templating effect of the tube shape on the total number of diffusing particles at each depth, while the volumetric concentration in the actually explored porous region may not have significant decay. These results may explain the anomalous diffusion of metal atoms in porous deposits observed in recent works. They also confirm the difficulty in interpreting experimental or computational data on anomalous transport reported in recent works, particularly if only the concentration profiles are measured.

Bubble monolayers are model systems for experiments and simulations of two-dimensional packing problems of deformable objects. We explore the relation between the distributions of the number of bubble sides (topology) and the bubble areas (geometry) in the low liquid fraction limit. We use a statistical model [M. Durand, Europhys. Lett. 90, 60002 (2010)] which takes into account Plateau laws. We predict the correlation between geometrical disorder (bubble size dispersity) and topological disorder (width of bubble side number distribution) over an extended range of bubble size dispersities. Extensive data sets arising from shuffled foam experiments, surface evolver simulations, and cellular Potts model simulations all collapse surprisingly well and coincide with the model predictions, even at extremely high size dispersity. At moderate size dispersity, we recover our earlier approximate predictions [M. Durand, J. Kafer, C. Quilliet, S. Cox, S. A. Talebi, and F. Graner, Phys. Rev. Lett. 107, 168304 (2011)]. At extremely low dispersity, when approaching the perfectly regular honeycomb pattern, we study how both geometrical and topological disorders vanish. We identify a crystallization mechanism and explore it quantitatively in the case of bidisperse foams. Due to the deformability of the bubbles, foams can crystallize over a larger range of size dispersities than hard disks. The model predicts that the crystallization transition occurs when the ratio of largest to smallest bubble radii is $1.4$.

Below the roughening transition, crystal surfaces have macroscopic plateaus, facets, whose evolution is driven by the microscale dynamics of steps. A long-standing puzzle was how to reconcile discrete effects in facet motion with fully continuum approaches. We propose a resolution of this issue via connecting, through a  jump condition, the continuum-scale surface chemical potential away from the facet, characterized by variations of the continuum surface free energy, with a chemical potential originating from the decay of atomic steps on top of the facet. The proposed condition accounts for step flow inside a discrete boundary layer near the facet. To validate this approach, we implement in a radial geometry a hybrid discrete-continuum scheme in which the continuum theory is coupled with only a few, minimally three, steps in diffusion-limited kinetics with conical initial data.

We have investigated the behavior of bouncing drops on a liquid surface by using particle image velocimetry analysis. A drop on an oscillating liquid surface is observed to not coalesce with the liquid and to travel along the surface if the oscillation is strong enough. A streaming vortex pair, induced by the alternatively distorted liquid surface, shows up below a bouncing drop. The time-averaged flow fields of the vortices are measured. In our quasi-one-dimensional setup, there are three stable distances for the drops, which can be characterized by the Faraday wavelength. The interactions of the vortex-mediated bouncing drops are deduced from the streamlines in the liquid bulk. We further show that a three-dimensional vortex ring is induced by a bouncing drop in a square cell.

Electrically discharged plasma channels can guide laser pulses, extending the laser-plasma interaction length to many Rayleigh ranges. In applications such as the laser-plasma accelerator, the laser group velocity in the channel plays a critical role. The laser travel time (and thus the averaged group velocity) was measured through two-pulse frequency-domain interferometry and was found to depend on the on-axis plasma density and laser spot size. The data is in agreement with theory.

For a shell model of the fully developed turbulence and the incompressible Navier-Stokes equations in the Fourier space, when a Gaussian white noise is artificially added to the equation of each mode, an expression of the mean linear response function in terms of the velocity correlation functions is derived by applying the method developed for nonequilibrium Langevin systems [Harada and Sasa, Phys. Rev. Lett. 95, 130602 (2005)]. We verify numerically for the shell-model case that the derived expression of the response function, as the noise tends to zero, converges to the response function of the noiseless shell model.

We study the rock-paper-scissors game in structured populations, where the invasion rates determine individual payoffs that govern the process of strategy change. The traditional version of the game is recovered if the payoffs for each potential invasion stem from a single pairwise interaction. However, the transformation of invasion rates to payoffs also allows the usage of larger interaction ranges. In addition to the traditional pairwise interaction, we therefore consider simultaneous interactions with all nearest neighbors, as well as with all nearest and next-nearest neighbors, thus effectively going from single pair to group interactions in games of cyclic dominance. We show that differences in the interaction range affect not only the stationary fractions of strategies but also their relations of dominance. The transition from pairwise to group interactions can thus decelerate and even revert the direction of the invasion between the competing strategies. Like in evolutionary social dilemmas, in games of cyclic dominance, too, the indirect multipoint interactions that are due to group interactions hence play a pivotal role. Our results indicate that, in addition to the invasion rates, the interaction range is at least as important for the maintenance of biodiversity among cyclically competing strategies.

Molecular dynamics simulation is used to examine the thermal slip of liquids at rough solid surfaces as characterized by fractal Cantor structures. The temperature profiles, potential energy distributions, thermal slip, and interfacial thermal resistance are investigated and evaluated for a variety of surface topographies. In addition, the effects of liquid-solid interaction, surface stiffness, and boundary condition on thermal slip length are presented. Our results indicate that the presence of roughness expands the low potential energy regions in adjacent liquids, enhances the energy transfer at liquid-solid interface, and decreases the thermal slip. Interestingly, the thermal slip length and thermal resistance for liquids in contact with solid surfaces depends not only on the statistical roughness height, but also on the fractal dimension (i.e., topographical spectrum).

In this paper, we consider the complex modified Korteweg–de Vries (mKdV) equation as a model of few-cycle optical pulses. Using the Lax pair, we construct a generalized Darboux transformation and systematically generate the first-, second-, and third-order rogue wave solutions and analyze the nature of evolution of higher-order rogue waves in detail. Based on detailed numerical and analytical investigations, we classify the higher-order rogue waves with respect to their intrinsic structure, namely, fundamental pattern, triangular pattern, and ring pattern. We also present several new patterns of the rogue wave according to the standard and nonstandard decomposition. The results of this paper explain the generalization of higher-order rogue waves in terms of rational solutions. We apply the contour line method to obtain the analytical formulas of the length and width of the first-order rogue wave of the complex mKdV and the nonlinear Schrödinger equations. In nonlinear optics, the higher-order rogue wave solutions found here will be very useful to generate high-power few-cycle optical pulses which will be applicable in the area of ultrashort pulse technology.

Network of coupled oscillators has long been employed as an important approach to explore the complicated dynamics in spatially extended systems. Here we show how this approach can be used to the analysis of turbulence pinning control. Specifically, by use of a model of two-dimensional drift-wave plasma turbulence, we investigate how the performance of the turbulence control is influenced by the spatial distribution of the pinning strength. It is found that the dynamics of pinned turbulence can be well captured by a simple model of networked modes, based on which the dependence of the control performance on the pinning distribution can be analytically obtained. In particular, the model predicts that as the distribution of the pinning strength becomes more nonuniform, the performance of turbulence control will be gradually decreased. This theoretical prediction is in good agreement with the results of numerical simulations, including the sinusoidal and localized pinning distributions. Our studies provide a new viewpoint to the mechanism of mode couplings in drift-wave turbulence, as well as be constructive to the design of new schemes for controlling turbulence in realistic systems.

Estimating trajectories and parameters of dynamical systems from observations is a problem frequently encountered in various branches of science; geophysicists for example refer to this problem as data assimilation. Unlike as in estimation problems with exchangeable observations, in data assimilation the observations cannot easily be divided into separate sets for estimation and validation; this creates serious problems, since simply using the same observations for estimation and validation might result in overly optimistic performance assessments. To circumvent this problem, a result is presented which allows us to estimate this optimism, thus allowing for a more realistic performance assessment in data assimilation. The presented approach becomes particularly simple for data assimilation methods employing a linear error feedback (such as synchronization schemes, nudging, incremental 3DVAR and 4DVar, and various Kalman filter approaches). Numerical examples considering a high gain observer confirm the theory.

In microfluidic applications involving high-frequency acoustic waves over a solid boundary, the Stokes boundary-layer thickness $\delta$ is so small that some non-negligible slip may occur at the fluid-solid interface. This paper assesses the impact of this slip by revisiting the classical problem of steady acoustic streaming over a flat boundary, replacing the no-slip boundary condition with the Navier condition ${u|}_{y=0}=L_{\mathrm{s}}{\partial }_y{u|}_{y=0}$, where $u$ is the velocity tangent to the boundary $y=0$, and the parameter $L_{\mathrm{s}}$ is the slip length. A general expression is obtained for the streaming velocity across the boundary layer as a function of the dimensionless parameter $L_{\mathrm{s}}/\delta$. The limit outside the boundary layer provides an effective slip velocity satisfied by the interior mean flow. Particularizing to traveling and standing waves shows that the boundary slip respectively increases and decreases the streaming velocity.

A surface reaction boundary condition in multicomponent lattice Boltzmann simulations is developed. The method is applied to a test case with nonlinear reaction rates and nonlinear density profiles. The results are compared to the corresponding analytical solution, which shows that the error of the method scales with the square of the lattice spacing.

Exascale computing presents a challenge for the scientific community as new algorithms must be developed to take full advantage of the new computing paradigm. Atomistic simulation methods that offer full fidelity to the underlying potential, i.e., molecular dynamics (MD) and parallel replica dynamics, fail to use the whole machine speedup, leaving a region in time and sample size space that is unattainable with current algorithms. In this paper, we present an extension of the parallel replica dynamics algorithm [A. F. Voter, Phys. Rev. B 57, R13985 (1998)] by combining it with the synchronous sublattice approach of Shim and Amar [ and , Phys. Rev. B 71, 125432 (2005)], thereby exploiting event locality to improve the algorithm scalability. This algorithm is based on a domain decomposition in which events happen independently in different regions in the sample. We develop an analytical expression for the speedup given by this sublattice parallel replica dynamics algorithm and compare it with parallel MD and traditional parallel replica dynamics. We demonstrate how this algorithm, which introduces a slight additional approximation of event locality, enables the study of physical systems unreachable with traditional methodologies and promises to better utilize the resources of current high performance and future exascale computers.