The lack of a first-principles derivation has made the hydrodynamic boundary condition a classical issue for the past century. The fact that the fluid can have interfacial structures adds additional complications and ambiguities to the problem. Here we report the use of molecular dynamics to identify from equilibrium thermal fluctuations the hydrodynamic modes in a fluid confined by solid walls, thereby extending the application of the fluctuation-dissipation theorem to yield not only the accurate location of the hydrodynamic boundary at the molecular scale, but also the relevant parameter value(s) for the description of the macroscopic boundary condition. We present molecular dynamics results on two examples to illustrate the application of this approach—one on the hydrophilic case and one on the hydrophobic case. It is shown that the use of the orthogonality condition of the modes can uniquely locate the hydrodynamic boundary to be inside the fluid in both cases, separated from the molecular solid-liquid interface by a small distance $\mathrm{\Delta }$ that is a few molecules in size. The eigenvalue equation of the hydrodynamic modes directly yields the slip length, which is about equal to $\mathrm{\Delta }$ in the hydrophilic case but is larger than $\mathrm{\Delta }$ in the hydrophobic case. From the decay time we also obtain the bulk viscosity which is in good agreement with the value obtained from dynamic simulations. To complete the picture, we derive the Green-Kubo relation for a finite fluid system and show that the boundary fluctuations decouple from the bulk only in the infinite-fluid-channel limit; and in that limit we recover the interfacial fluctuation-dissipation theorem first presented by Bocquet and Barrat. The coupling between the bulk and the boundary fluctuations provides both the justification and the reason for the effectiveness of the present approach, which promises broad utility for probing the hydrodynamic boundary conditions relevant to structured or elastic interfaces, as well as two-phase immiscible flows.

Mixing at low Reynolds numbers, especially in the framework of confined flows occurring in Hele-Shaw cells, porous media, and microfluidic devices, has attracted considerable attention lately. Under such circumstances, enhanced mixing is limited due to the lack of turbulence, and absence of sizable inertial effects. Recent studies, performed in rectangular Hele-Shaw cells, have demonstrated that the combined action of viscous fluid fingering and alternating injection can dramatically improve mixing efficiency. In this work, we revisit this important fluid mechanical problem, and analyze it in the context of radial Hele-Shaw flows. The development of radial fingering instabilities under alternating injection conditions is investigated by intensive numerical simulations. We focus on the impact of the relevant physical parameters of the problem (Péclet number $\mathrm{Pe}$, viscosity contrast $A$, and injection time interval $\mathrm{\Delta }t)$ on fluid mixing performance.

The microscopic quantum theory of plasma nonlinear interaction with the coherent shortwave electromagnetic radiation of arbitrary intensity is developed. The Liouville–von Neumann equation for the density matrix is solved analytically considering a wave field exactly and a scattering potential of plasma ions as a perturbation. With the help of this solution we calculate the nonlinear inverse-bremsstrahlung absorption rate for a grand canonical ensemble of electrons. The latter is studied in Maxwellian, as well as in degenerate quantum plasma for x-ray lasers at superhigh intensities and it is shown that one can achieve the efficient absorption coefficient in these cases.

We investigate site percolation on a weighted planar stochastic lattice (WPSL), which is a multifractal and whose dual is a scale-free network. Percolation is typically characterized by a threshold value $p_c$ at which a transition occurs and by a set of critical exponents $\beta , \gamma , \nu$ which describe the critical behavior of the percolation probability $P\left(p\right)$, mean cluster size $S\left(p\right)$, and the correlation length $\xi$. Besides, the exponent $\tau$ characterizes the cluster size distribution function $n_s\left(p_c\right)$ and the fractal dimension $d_f$ characterizes the spanning cluster. We numerically obtain the value of $p_c$ and of all the exponents. These results suggest that the percolation on WPSL belong to a separate universality class than on all other planar lattices.

Connectivity and conductivity of two-dimensional fracture networks (FNs), as an important type of continuous random networks, are examined systematically through Monte Carlo simulations under a variety of conditions, including different power law distributions of the fracture lengths and domain sizes. The simulation results are analyzed using analogies of the percolation theory for discrete random networks. With a characteristic length scale and conductivity scale introduced, we show that the connectivity and conductivity of FNs can be well described by universal scaling solutions. These solutions shed light on previous observations of scale-dependent FN behavior and provide a powerful method for quantifying effective bulk properties of continuous random networks.

We consider linear arrays of cells of volume $V_{\mathrm{c}}$ populated by monodisperse rods of size $\sigma V_{\mathrm{c}},\sigma =1,2,...$, subject to hardcore exclusion interaction. Each rod experiences a position-dependent external potential. In one application we also examine effects of contact forces between rods. We employ two distinct methods of exact analysis with complementary strengths and different limits of spatial resolution to calculate profiles of pressure and density on mesoscopic and microscopic length scales at thermal equilibrium. One method uses density functionals and the other statistically interacting vacancy particles. The applications worked out include gravity, power-law traps, and hard walls. We identify oscillations in the profiles on a microscopic length scale and show how they are systematically averaged out on a well-defined mesoscopic length scale to establish full consistency between the two approaches. The continuum limit, realized as $V_{\mathrm{c}}\to 0,\sigma \to \infty$ at nonzero and finite $\sigma V_{\mathrm{c}}$, connects our highest-resolution results with known exact results for monodisperse rods in a continuum. We also compare the pressure profiles obtained from density functionals with the average microscopic pressure profiles derived from the pair distribution function.

We extend the exchange fluctuation theorem for energy exchange between thermal quantum systems beyond the assumption of molecular chaos, and describe the nonequilibrium exchange dynamics of correlated quantum states. The relation quantifies how the tendency for systems to equilibrate is modified in high-correlation environments. In addition, a more abstract approach leads us to a “correlation fluctuation theorem”. Our results elucidate the role of measurement disturbance for such scenarios. We show a simple application by finding a semiclassical maximum work theorem in the presence of correlations. We also present a toy example of qubit-qudit heat exchange, and find that non-classical behaviour such as deterministic energy transfer and anomalous heat flow are reflected in our exchange fluctuation theorem.

Equilibrium behavior of a single chain of dipolar spheres is investigated by the method of molecular dynamics in a wide range of the dipolar coupling constant $\lambda$. Two cases are considered: rodlike and flexible chains. In the first case, particle centers are immovably fixed on one axis, but their magnetic moments retain absolute orientational freedom. It has been found that at $\lambda \gtrsim 1.5$ particle moments are chiefly aligned parallel to the chain axis, but the total moment of the chain continuously changes its sign with some mean frequency, which exponentially decreases with the growth of $\lambda$. Such behavior of the rodlike chain is analogous to the Néel relaxation of a superparamagnetic particle with a finite energy of magnetic anisotropy. In the flexible chain particles are able to move in the three-dimensional space, but the distance between centers of the first-nearest neighbors never exceeds a given limiting value $r_{max}$. If $r_{max}\simeq d$ $(d$ is the particle diameter) then the most probable shape of the chain of five or more particles at $\lambda \gtrsim 6$ is that of a ring. The behavior of chains with $r_{max}\ge 2d$ is qualitatively different: At $\lambda \simeq 4$ long chains collapse into dense quasispherical globules and at $\lambda \gtrsim 8$ these globules take toroidal configuration with a spontaneous azimuthal ordering of magnetic dipoles. With the increase of $r_{max}$ to larger values ($r_{max}>10d)$ globules expand and break down to form separate rings.

We provide a comprehensive physical description of the vaporization, self-assembly, agglomeration, and buckling kinetics of sessile nanofluid droplets pinned on a hydrophobic substrate. We have deciphered five distinct regimes of the droplet life cycle. Regimes I–III consists of evaporation-induced preferential agglomeration that leads to the formation of a unique dome-shaped inhomogeneous shell with a stratified varying-density liquid core. Regime IV involves capillary-pressure-initiated shell buckling and stress-induced shell rupture. Regime V marks rupture-induced cavity inception and growth. We demonstrate through scaling arguments that the growth of the cavity (which controls the final morphology or structure) can be described by a universal function.

Liquid crystalline elastomers combine the ordering properties of liquid crystals with elasticity of crosslinked polymer networks. In monodomain (permanently aligned) elastomers, altering the orientational (nematic) order causes changes in the equilibrium sample length, which is the basis of the famous effect of large-amplitude reversible mechanical actuation. The stimulus for this effect could be a change in temperature, or illumination by light in photosensitized elastomers, but equally the nematic order changes by mixing with a solvent. This work theoretically investigates a competition between the spontaneous contraction on swelling of a monodomain nematic elastomer and the externally imposed stretching. We find that this competition leads to bistability in the system and allows a two-phase separation between a nematic state with lower swelling and an isotropic state with higher solvent concentration. We calculated the conditions in which the instability occurs as well as the mechanical and geometric parameters of equilibrium states. Being able to predict how this instability arises will provide opportunities for exploiting nematic elastomer filaments.

Integrin receptors mediate interaction between the cellular actin-cytoskeleton and extracellular matrix. Based on their activation properties, we propose a reaction-diffusion model where the kinetics of the two-state receptors is modulated by their lipidic environment. This environment serves as an activator variable, while a second variable plays the role of a scaffold protein and controls the self-sustained activation of the receptors. Due to receptor diffusion which couples dynamically the activator and the inhibitor, our model connects major classes of reaction diffusion systems for excitable media. Spot and rosette solutions, characterized by receptor clustering into localized static or dynamic structures, are organized into a phase diagram. It is shown that diffusion and kinetics of receptors determines the dynamics and the stability of these structures. We discuss this model as a precursor model for cell signaling in the context of podosomes forming actoadhesive metastructures, and we study how generic signaling defects influence their organization.

The impact of connectivity and individual dynamics on the basin stability of the burst synchronization regime in small-world networks consisting of chaotic slow-fast oscillators is studied. It is shown that there are rewiring probabilities corresponding to the largest basin stabilities, which uncovers a reason for finding small-world topologies in real neuronal networks. The impact of coupling density and strength as well as the nodal parameters of relaxation or excitability are studied. Dynamic mechanisms are uncovered that most strongly influence basin stability of the burst synchronization regime.

We experimentally investigated the noise-induced dynamics of a prototypical thermoacoustic system undergoing a subcritical Hopf bifurcation to limit cycle oscillations. The study was performed prior to the bistable regime. Analysis of the characteristics of pressure oscillations in the combustor and fluctuations in the heat release rate from the flame—the two physical entities involved in thermoacoustic coupling—at increasing levels of noise indicated precursors to the Hopf bifurcation. These precursors were further identified to be a result of coherence resonance.

Ultrasonic melt processing significantly improves the properties of metallic materials. However, this promising technology has not been successfully transferred to the industry because of difficulties in treating large volumes of melt. To circumvent these difficulties, a fundamental understanding of the efficiency of ultrasonic treatment of liquid metals is required. In this endeavor, the dynamics of two interacting hydrogen bubbles in liquid aluminum are studied to determine the effect of a strong acoustic field on their behavior. It is shown that coalescence readily occurs at low frequencies in the range of 16 to 20 kHz; forcing frequencies at these values are likely to promote degassing. Emitted acoustic pressures from relatively isolated bubbles that resonate with the driving frequency are in the megapascal range and these cavitation shock waves are presumed to promote grain refinement by disrupting the growth of the solidification front.

We present fluid velocity measurements in a modified Taylor-Couette device operated in the quasi-Keplerian regime, where it is observed that nearly ideal flows exhibit self-similarity under scaling of the Reynolds number. In contrast, nonideal flows show progressive departure from ideal Couette as the Reynolds number is increased. We present a model that describes the observed departures from ideal Couette rotation as a function of the fluxes of angular momentum across the boundaries, capturing the dependence on Reynolds number and boundary conditions.

The dynamics in the thin boundary layers of temperature and velocity is the key to a deeper understanding of turbulent transport of heat and momentum in thermal convection. The velocity gradient at the hot and cold plates of a Rayleigh-Bénard convection cell forms the two-dimensional skin friction field and is related to the formation of thermal plumes in the respective boundary layers. Our analysis is based on a direct numerical simulation of Rayleigh-Bénard convection in a closed cylindrical cell of aspect ratio $\mathrm{\Gamma }=1$ and focused on the critical points of the skin friction field. We identify triplets of critical points, which are composed of two unstable nodes and a saddle between them, as the characteristic building block of the skin friction field. Isolated triplets as well as networks of triplets are detected. The majority of the ridges of linelike thermal plumes coincide with the unstable manifolds of the saddles. From a dynamical Lagrangian perspective, thermal plumes are formed together with an attractive hyperbolic Lagrangian coherent structure of the skin friction field. We also discuss the differences from the skin friction field in turbulent channel flows from the perspective of the Poincaré-Hopf index theorem for two-dimensional vector fields.

We investigated ion acceleration by an electrostatic shock in an exploded target irradiated by an ultrashort, circularly polarized laser pulse by means of one- and three-dimensional particle-in-cell simulations. We discovered that the laser field penetrating via relativistic transparency (RT) rapidly heated the upstream electron plasma to enable the formation of a high-speed electrostatic shock. Owing to the RT-based rapid heating and the fast compression of the initial density spike by a circularly polarized pulse, a new regime of the shock ion acceleration driven by an ultrashort (20–40 fs), moderately intense (1–1.4 PW) laser pulse is envisaged. This regime enables more efficient shock ion acceleration under a limited total pulse energy than a linearly polarized pulse with crystal laser systems of $\lambda \sim 1\mu \mathrm{m}$.

An alternative application of duality relations of stochastic processes is demonstrated. Although conventional usages of the duality relations need analytical solutions for the dual processes, here I employ numerical solutions of the dual processes and investigate the usefulness. As a demonstration, estimation problems of hidden variables in stochastic differential equations are discussed. Employing algebraic probability theory, a little complicated birth-death process is derived from the stochastic differential equations, and an estimation method based on the ensemble Kalman filter is proposed. As a result, the possibility for making faster computational algorithms based on the duality concepts is shown.

An interconnected network features a structural transition between two regimes [F. Radicchi and A. Arenas, Nat. Phys. 9, 717 (2013)]: one where the network components are structurally distinguishable and one where the interconnected network functions as a whole. Our exact solution for the coupling threshold uncovers network topologies with unexpected behaviors. Specifically, we show conditions that superdiffusion, introduced by Gómez et al. [Phys. Rev. Lett. 110, 028701 (2013)], can occur despite the network components functioning distinctly. Moreover, we find that components of certain interconnected network topologies are indistinguishable despite very weak coupling between them.

Networks of fast nonlinear elements may display slow fluctuations if interactions are strong. We find a transition in the long-term variability of a sparse recurrent network of perfect integrate-and-fire neurons at which the Fano factor switches from zero to infinity and the correlation time is minimized. This corresponds to a bifurcation in a linear map arising from the self-consistency of temporal input and output statistics. More realistic neural dynamics with a leak current and refractory period lead to smoothed transitions and modified critical couplings that can be theoretically predicted.