In a sheared and confined granular flow, the mean force and the force fluctuations on a rigid wall are studied by means of numerical simulations based on the discrete element method. An original periodic immersed-wall system is designed to investigate a wide range of confinement pressure and shearing velocity imposed at the top of the flow, considering different obstacle heights. The mean pressure on the wall relative to the confinement pressure is found to be a monotonic function of the boundary macroscopic inertial number which encapsulates the confinement pressure, the shearing velocity, and the thickness of the sheared layer above the wall. The one-to-one relation is slightly affected by the length of the granular system. The force fluctuations on the wall are quantified through the analysis of both the distributions of grain-wall contact forces and the autocorrelation of force time series. The distributions narrow as the boundary macroscopic inertial number decreases, moving from asymmetric log-normal shape to nearly Gaussian-type shape. That evolution of the grain-wall force distributions is accompanied at the lowest inertial numbers by the occurrence of a system memory in terms of the force transmitted to the wall, provided that the system length is not too large. Moreover, the distributions of grain-wall contact forces are unchanged when the inertial number is increased above a critical value. All those results allow to clearly identify the transitions from quasistatic to dense inertial, and from dense inertial to collisional, granular flow regimes.

We introduce a variant of the Hybrid Monte Carlo (HMC) algorithm to address large-deviation statistics in stochastic hydrodynamics. Based on the path-integral approach to stochastic (partial) differential equations, our HMC algorithm samples space-time histories of the dynamical degrees of freedom under the influence of random noise. First, we validate and benchmark the HMC algorithm by reproducing multiscale properties of the one-dimensional Burgers equation driven by Gaussian and white-in-time noise. Second, we show how to implement an importance sampling protocol to significantly enhance, by orders of magnitudes, the probability to sample extreme and rare events, making it possible to estimate moments of field variables of extremely high order (up to 30 and more). By employing reweighting techniques, we map the biased configurations back to the original probability measure in order to probe their statistical importance. Finally, we show that by biasing the system towards very intense negative gradients, the HMC algorithm is able to explore the statistical fluctuations around instanton configurations. Our results will also be interesting and relevant in lattice gauge theory since they provide unique insights into reweighting techniques.

A monolithic solver based on a diffuse-interface immersed-boundary (IB) approach for conjugate-heat-transfer (CHT) problems is presented. The IB strategy assumes that the solid which is “immersed” into the computational grid is occupied by a “virtual” fluid to facilitate construction of “unified” governing equations that are solved everywhere in the domain. A unified momentum equation is devised using the solid volume fraction that reduces to the Navier-Stokes equation outside of the solid and to the no-slip boundary condition inside of it. The “unified” energy equation is constructed in an analogous fashion reducing to a convective-diffusive equation in the fluid domain and a fully diffusive equation in the solid domain with different thermal conductivities (or diffusivities) for both domains. The resulting equations are solved in both domains simultaneously using a hybrid staggered and nonstaggered finite-volume (FV) framework for incompressible flows. The second-order accurate IB-FV solver is employed to carry out investigations for CHT problems in natural and forced convective regimes. Numerical studies for different fluid-to-solid conductivity ratios show that the monolithic IB-CHT solver is a fast, simple, and accurate framework for simulations of CHT problems for Boussinesq flows.

We consider magnetic friction between two square lattices of the ferromagnetic Ising model of finite thickness. We analyze the dependence on the boundary conditions and the sample thickness. Monte Carlo results indicate that the setup enables us to control the frictional force by magnetic fields on the boundaries. In addition, we confirm that the temperature derivative of the frictional force as well as that of the boundary energy has singularity at a velocity-dependent critical temperature.

The short-time relaxation critical behavior of the $XY$ model on a simple-cubic lattice is investigated within the scope of deterministic Hamiltonian dynamics. The Hamiltonian includes a first-neighbor interaction between planar vectors and a rotational kinetic term from which the motion equations are derived. The dynamical evolution from a fully ordered initial state is followed by employing a symplectic algorithm based on a high-order Trotter-Suzuki decomposition of the time-evolution operator. A finite-time scaling analysis is performed to provide accurate estimates of the critical energy density, the order-parameter relaxation exponent, and the dynamical critical exponent. The estimated critical exponents are consistent with prior theoretical and experimental values reported for the superfluid ${}^4\mathrm{He}$, extreme type-II superconducting, and Bose-Einstein condensation transitions.

We study large deviations of the time-averaged size of stochastic populations described by a continuous-time Markov jump process. When the expected population size $N$ in the steady state is large, the large deviation function (LDF) of the time-averaged population size can be evaluated by using a Wentzel-Kramers-Brillouin (WKB) method, applied directly to the master equation for the Markov process. For a class of models that we identify, the direct WKB method predicts a giant disparity between the probabilities of observing an unusually small and an unusually large values of the time-averaged population size. The disparity results from a qualitative change in the “optimal” trajectory of the underlying classical mechanics problem. The direct WKB method also predicts, in the limit of $N\to \infty$, a singularity of the LDF, which can be interpreted as a second-order dynamical phase transition. The transition is smoothed at finite $N$, but the giant disparity remains. The smoothing effect is captured by the van-Kampen system size expansion of the exact master equation near the attracting fixed point of the underlying deterministic model. We describe the giant disparity at finite $N$ by developing a different variant of WKB method, which is applied in conjunction with the Donsker-Varadhan large-deviation formalism and involves subleading-order calculations in $1/N$.

This article describes high-field nuclear magnetic resonance experiments using pulsed gradient spin echo methods to investigate water mobility inside homogeneously swollen polyacrylamide hydrogel beads. The NMR data permitted determining the solvent self-diffusion coefficient dependence on either the swelling time or the polymer volume fraction. The sensitivity to the structure of the hydrogel network could be tuned to a certain extent via controlling the diffusion probing time. Relaxation measurements have helped substantiate a theoretical description for these self-diffusion dependencies based upon a fast-exchange two-site model. A tortuosity-porosity dependence is extracted from the model and then compared to other known tortuosity regimes. The diffusion dependence on the degree of swelling obtained here corroborates the existence of a gel swelling-front observed during inhomogeneous solvent imbibition experiments. Finally, we stress that the NMR experimental results reported here are of great value on performing theoretical modeling of gel swelling via solvent imbibition.

A fluid particle changes its dynamics from diffusive to oscillatory as the system density increases up to the melting density. Hence the notion of the Frenkel line was introduced to demarcate the fluid region into rigid and nonrigid liquid subregions based on the collective particle dynamics. In this work, we apply a topological framework to locate the Frenkel lines of the soft-sphere and the hard-sphere models relying on the system configurations. The topological characteristics of the ideal gas and the maximally random jammed state are first analyzed, then the classification scheme designed in our earlier work is applied to soft-sphere and hard-sphere fluids. The dependence of the classification result on the bulk density is understood based on the theory of fluid polyamorphism. The percolation behavior of solid-like clusters is described based on the fraction of solid-like molecules in an integrated manner. The crossover densities are obtained by examining the percolation of solid-like clusters. The resultant crossover densities of soft-sphere fluids converge to that of hard-sphere fluid. Hence the topological method successfully highlights the generality of the Frenkel line.

An efficient route to the displacement field around a rigid spherical inclusion in an infinitely extended homogeneous elastic medium is presented in a slightly alternative way when compared to some common textbook methods. Moreover, two Faxén relations of next-higher order beyond the stresslet are calculated explicitly for compressible media. They quantify higher-order moments involving the force distribution on a rigid spherical particle in a deformed elastic medium. As a consequence, additional contributions to the distortions of the deformed elastic medium are identified that are absent to lower order. Furthermore, the displaceability and rotateability matrices for an ensemble of rigid spheres are calculated up to (including) sixth order in inverse particle separation distance. These matrices describe the interactions mediated between the rigid embedded particles by the elastic environment. In this way, additional coupling effects are identified that are absent to lower order, particularly when rotations and torques are involved. All methods and results can formally be transferred to the corresponding case of incompressible hydrodynamic low-Reynolds-number Stokes flow by considering the limit of an incompressible environment. The roles of compressibility of the embedding medium and of the here additionally derived higher-order contributions are highlighted by some selected example configurations.

The Richtmyer-Meshkov instability of small perturbed single-mode interfaces between an elastic-plastic solid and an inviscid liquid is investigated by theoretical analysis and numerical simulation in this work. A modified model including the Atwood number effect is proposed to describe the long-term behaviors of small perturbations at the solid-liquid interface. In contrast to an effective theoretical model at the solid-vacuum interface, this model is appropriate at different Atwood numbers. Owing to the effect of elastic-plastic characteristics and the density ratio, the evolution of the spike amplitude exhibits nonlinear mechanical behavior. As the absolute value of the Atwood number decreases, the maximum spike amplitude also decreases. To validate this model, an Eulerian finite-difference multicomponent code is adopted to study the time evolution of the spike amplitude at different Atwood numbers. The model coefficients are obtained by analyzing the relevant characteristic statistics collected from the numerical results. Under different initial conditions such as Atwood number and shock strength, the applicability of this modified model is verified by comparing the numerical results with the model profile. The consistency in results signifies that the modified model is not only suitable for specific shock intensity and Atwood number, but also adaptable within a certain range.

In this paper, we investigate the first-passage-time distribution (FPTD) within a time-dependent parabolic potential in the presence of roughness with two methods: the Kramers theory and a nonsingular integral equation. By spatially averaging, the rough potential is equivalent to the combination of an effective smooth potential and an effective diffusion coefficient. Based on the Kramers theory, we first obtain Kramers approximations (KAs) of FPTD for both smooth and rough potentials. As expected, KA is valid only for high barriers and small external forces, and generally applicable for high barriers in rough potentials. To overcome the shortcoming of KA, a probability asymptotic approximation (PAA) based on an integral equation is proposed, which uses the transient probability density function (PDF) of the natural boundary conditions instead of the absorbing boundary conditions. We find that PAA fits very well even for large external forces. This enables us to analytically solve the FPTD for large external forces and low barriers as a strong extension to KA. In addition, we show that in the presence of a rough potential, the PAA of FPTD is in good agreement with numerical simulations for low barrier potentials. The PAA makes it possible to investigate the first-passage problem with ultrafast varying potentials and short exiting time. Thus, KA and PAA are complementary in determining the FPTD both for various barriers and external forces. Finally, the mean first-passage time (MFPT) is studied, which illustrates that the PAA of MFPT is effective almost in the whole range of external forces, while the KA of MFPT is valid only for small external forces.

Surface growth is a crucial component of many natural and artificial processes, from cell proliferation to additive manufacturing. In elastic systems surface growth is usually accompanied by the development of geometrical incompatibility, leading to residual stresses and triggering various instabilities. In a recent paper [G. Zurlo and L. Truskinovsky, Phys. Rev. Lett. 119, 048001 (2017)] we presented a linearized elasticity theory of incompatible surface growth, which provides a quantitative link between deposition protocols and postgrowth states of stress. Here we extend this analysis to account for both physical and geometrical nonlinearities of an elastic solid. This development reveals the shortcomings of the linearized theory, in particular its inability to describe kinematically confined surface growth and to account for growth-induced elastic instabilities.

Expressions of pressure in warm and hot dense matter using the average-atom model are presented. They are based on the stress-tensor approach. Nonrelativistic and relativistic cases are considered. The obtained formulas are simple and can be easily implemented in an average-atom model code. Comparisons with experimental data and quantum molecular dynamics and path integral Monte Carlo simulations are shown. The present formalism agrees well with experimental results for a large variety of elements in the warm dense matter regime and with ab initio simulations in the warm and hot dense matter regime for aluminum.

The difficulty of obtaining appreciable quantities of biologically important molecules in thermodynamic equilibrium has long been identified as an obstacle to life's emergence, and determining the specific nonequilibrium conditions that might have given rise to life is challenging. To address these issues, we investigate how the concentrations of life's building blocks change as a function of the distance from equilibrium on average, in two example settings: (i) the synthesis of heavy amino acids and (ii) their polymerization into peptides. We find that relative concentrations of the heaviest amino acids can be boosted by four orders of magnitude, and concentrations of the longest peptide chains can be increased by hundreds of orders of magnitude. The average nonequilibrium distribution does not depend on the details of how the system was driven from equilibrium, indicating that environments might not have to be fine-tuned to support life.

The time that a diffusing particle spends in a certain region of space is known as the occupation time, or the residence time. Recently, the joint occupation-time statistics of an ensemble of noninteracting particles was addressed using the single-particle statistics. Here we employ the macroscopic fluctuation theory (MFT) to study the occupation-time statistics of many interacting particles. We find that interactions can significantly change the statistics and, in some models, even cause a singularity of the large-deviation function describing these statistics. This singularity can be interpreted as a dynamical phase transition. We also point out to a close relation between the MFT description of the occupation-time statistics of noninteracting particles and the level 2 large deviation formalism which describes the occupation-time statistics of a single particle.

A renormalized perturbation method is developed for quantum maps of periodically kicked rotor models to study the tunneling effect in the nearly integrable regime. Integrable Hamiltonians closely approximating the nonintegrable quantum map are systematically generated by the Baker-Hausdorff-Campbell (BHC) expansion for symmetrized quantum maps. The procedure results in an effective integrable renormalization, and the unrenormalized residual part is treated as the perturbation. If a sufficiently high-order BHC expansion is used as the base of perturbation theory, the lowest order perturbation well reproduces tunneling characteristics of the quasibound eigenstates, including the transition from the instanton tunneling to a noninstanton one. This approach enables a comprehensive understanding of the purely quantum mechanisms of tunneling in the nearly integrable regime. In particular, the staircase structure of tunneling probability dependence on quantum number can be clearly explained as the successive transition among multiquanta excitation processes. The transition matrix elements of the residual interaction have resonantly enhanced invariant components that are not removed by the renormalization. Eigenmodes coupled via these invariant components form noninstanton (NI) tunneling channels of two types contributing to the two regions of each step of the staircase structure: one type of channel is inside the separatrix, and the other goes across the separatrix. The amplitude of NI tunneling across the separatrix is insensitive to the Planck constant but shows an essentially singular dependence upon the nonintegrablity parameter. Its relation to the Melnikov integral, which characterizes the scale of classical chaos emerging close to the saddle on the potential top, is discussed.

We study through numerical simulation the Vicsek model for very low speeds and densities. We consider scalar noise in two and three dimensions and vector noise in three dimensions. We focus on the behavior of the critical noise with density and speed, trying to clarify seemingly contradictory earlier results. We find that, for scalar noise, the critical noise is a power law in both density and speed, but although we confirm the density exponent in two dimensions, we find a speed exponent different from earlier reports (we consider lower speeds than previous studies). On the other hand, for the vector noise case we find that the dependence of the critical noise cannot be separated as a product of power laws in speed and density. Finally, we study the dependence of the relaxation time with speed. At the critical point we find a power law, with the same exponent in two and three dimensions.

In this article, we formulate a set of the separate-phase governing equations at the representative-elementary-volume scale and develop its double-distribution function lattice Boltzmann (LB) algorithm to describe liquid-vapor two-phase flows with or without phase change in porous media. Different from those previous studies, the mathematical description in this article involves the Darcy force, viscous force, and pressure gradient, and the resulting LB simulations can well describe two-phase flows and mass transfer throughout porous media under the compounding effects of these forces. The LB algorithm was validated by simulating single-phase flows in porous media. Its results are in good agreement with those available analytical solutions. We also applied it to model water flows through a semi-infinite porous region bounded by a heated solid wall, where liquid-vapor phase change takes place. The numerical simulations recover the previous results in the limit of the zero Darcy number. Significantly, it reveals much richer two-phase flow and mass transfer characteristics in porous media adjacent to solid walls. The separate-phase model and its lattice Boltzmann algorithm in this article are effective means to gain more profound and clearer understandings of complex two-phase transport processes in a porous system.

Axons are linear structures of nerve cells that can range from a few tens of micrometers up to meters in length. In addition to external cues, the length of an axon is also regulated by unknown internal mechanisms. Molecular motors have been suggested to generate oscillations with an axon-length-dependent frequency that could be used to measure an axon's extension. Here, we present a mechanism for determining the axon length that couples the mechanical properties of an axon to the spectral decomposition of the oscillatory signal.

We theoretically investigate the nonlinear behavior of a buckled tip near the bifurcation point under external stress. We present a mechanical model for the buckled tip and derive the governing equation that describes the “buckling-to-flipping” nonlinear transition of the tip motion. Our minimal mechanistic model fully captures the velocity-dependent flipping phenomena, in which the flip position of the tip varies with the speed of the surface motion, as consistently observed in previous experiments. The present study could be applicable for sensitive detection of directional surface motion such as seismic waves.