In this research, the surface waves of a horizontal fluid layer open to air under gravity field and vertical temperature gradient effects are studied. The governing equations of this model are reformulated and converted to a nonlinear evolution equation, the perturbed Korteweg–de Vries (pKdV) equation. We investigate the latter equation, which includes dispersion, diffusion, and instability effects, in order to examine the evolution of long surface waves in a convective fluid. Dispersion relation of the pKdV equation and its properties are discussed. The Painlevé analysis is applied not only to check the integrability of the pKdV equation but also to establish the Bäcklund transformation form. In addition, traveling wave solutions and a general form of the multiple-soliton solutions of the pKdV equation are obtained via Bäcklund transformation, the simplest equation method using Bernoulli, Riccati, and Burgers' equations as simplest equations, and the factorization method.

Biochemical oscillations are prevalent in living organisms. Systems with a small number of constituents cannot sustain coherent oscillations for an indefinite time because of fluctuations in the period of oscillation. We show that the number of coherent oscillations that quantifies the precision of the oscillator is universally bounded by the thermodynamic force that drives the system out of equilibrium and by the topology of the underlying biochemical network of states. Our results are valid for arbitrary Markov processes, which are commonly used to model biochemical reactions. We apply our results to a model for a single KaiC protein and to an activator-inhibitor model that consists of several molecules. From a mathematical perspective, based on strong numerical evidence, we conjecture a universal constraint relating the imaginary and real parts of the first nontrivial eigenvalue of a stochastic matrix.

The shear modulus of jammed frictional granular materials with harmonic repulsive interaction under an oscillatory shear is numerically investigated. It is confirmed that the storage modulus, the real part of the shear modulus, for frictional grains with sufficiently small strain amplitude ${\gamma }_0$ discontinuously emerges at the jamming transition point. The storage modulus for small ${\gamma }_0$ differs from that of frictionless grains even in the zero friction limit, whereas they are almost identical with each other for sufficiently large ${\gamma }_0$, where the transition becomes continuous. The stress-strain curve exhibits a hysteresis loop even for a small strain, which connects a linear region for sufficiently small strain to another linear region for larger strain. We propose a scaling law to interpolate between the states of small and large ${\gamma }_0$.

Multimegabar laser-driven shock waves are unique tools for studying matter under extreme conditions. Accurate characterization of shocked matter is for instance necessary for measurements of equation of state data or opacities. This paper reports experiments performed at the LULI facility on the diagnosis of shock waves, using x-ray-absorption radiography. Radiographs are analyzed using standard Abel inversion. In addition, synthetic radiographs, which also take into account the finite size of the x-ray source, are generated using density maps produced by hydrodynamic simulations. Reported data refer to both plane cylindrical targets and hemispherical targets. Evolution and deformation of the shock front could be followed using hydrodynamic simulations.

We present a simple method to estimate the central charge of the conformal field theory corresponding to a critical point of a two-dimensional lattice model from Monte Carlo simulations. The main idea is to use the Wang-Landau flat-histogram algorithm, which allows us to obtain the free energy of a lattice model on a torus as a function of torus radii. The central charge is calculated with good precision from a free-energy scaling at the critical point. We apply the method to the Ising, tricritical Ising (Blume-Capel), Potts, and site-diluted Ising models, and we also discuss an estimation of the conformal weights.

Flow and heat transfer in periodic structures are of great interest for many applications. In this paper, we carefully examine the periodic features of fully developed periodic incompressible thermal flows, and incorporate them in the lattice Boltzmann method (LBM) for flow and heat transfer simulations. Two numerical approaches, the distribution modification (DM) approach and the source term (ST) approach, are proposed; and they can both be used for periodic thermal flows with constant wall temperature (CWT) and surface heat flux boundary conditions. However, the DM approach might be more efficient, especially for CWT systems since the ST approach requires calculations of the streamwise temperature gradient at all lattice nodes. Several example simulations are conducted, including flows through flat and wavy channels and flows through a square array with circular cylinders. Results are compared to analytical solutions, previous studies, and our own LBM calculations using different simulation techniques (i.e., the one-module simulation vs. the two-module simulation, and the DM approach vs. the ST approach) with good agreement. These simple, however, representative simulations demonstrate the accuracy and usefulness of our proposed LBM methods for future thermal periodic flow simulations.

One of the most basic characterizations of the relationship between two random variables, $X$ and $Y$, is the value of their mutual information. Unfortunately, calculating it analytically and estimating it empirically are often stymied by the extremely large dimension of the variables. One might hope to replace such a high-dimensional variable by a smaller one that preserves its relationship with the other. It is well known that either $X$ (or $Y)$ can be replaced by its minimal sufficient statistic about $Y$ (or $X)$ while preserving the mutual information. While intuitively reasonable, it is not obvious or straightforward that both variables can be replaced simultaneously. We demonstrate that this is in fact possible: the information $X$'s minimal sufficient statistic preserves about $Y$ is exactly the information that $Y$'s minimal sufficient statistic preserves about $X$. We call this procedure information trimming. As an important corollary, we consider the case where one variable is a stochastic process' past and the other its future. In this case, the mutual information is the channel transmission rate between the channel's effective states. That is, the past-future mutual information (the excess entropy) is the amount of information about the future that can be predicted using the past. Translating our result about minimal sufficient statistics, this is equivalent to the mutual information between the forward- and reverse-time causal states of computational mechanics. We close by discussing multivariate extensions to this use of minimal sufficient statistics.

In many complex diffusion processes the drift of random walkers is not caused by an external force, as in the case of Brownian motion, but by local variations of fitness perceived by the random walkers. In this paper, a simple but general framework is presented that describes such a type of random motion and may be of relevance in different problems, such as opinion dynamics, cultural spreading, and animal movement. To this aim, we study the problem of a random walker in $d$ dimensions moving in the presence of a local heterogeneous attracting factor expressed in terms of an assigned position-dependent “attractiveness function.” At variance with standard Brownian motion, the attractiveness function introduced here regulates both the advection and diffusion of the random walker, thus providing testable predictions for a specific form of fluctuation-relations. We discuss the relation between the drift-diffusion equation based on the attractiveness function and that describing standard Brownian motion, and we provide some explicit examples illustrating its relevance in different fields, such as animal movement, chemotactic diffusion, and social dynamics.

We develop a two-species exclusion process with a distinct pair of entry and exit sites for each species of rigid rods. The relatively slower forward stepping of the rods in an extended bottleneck region, located in between the two entry sites, controls the extent of interference of the codirectional flow of the two species of rods. The relative positions of the sites of entry of the two species of rods with respect to the location of the bottleneck are motivated by a biological phenomenon. However, the primary focus of the study here is to explore the effects of the interference of the flow of the two species of rods on their spatiotemporal organization and the regulations of this interference by the extended bottleneck. By a combination of mean-field theory and computer simulation, we calculate the flux of both species of rods and their density profiles as well as the composite phase diagrams of the system. If the bottleneck is sufficiently stringent, then some of the phases become practically unrealizable, although not ruled out on the basis of any fundamental physical principle. Moreover, the extent of suppression of flow of the downstream entrants by the flow of the upstream entrants can also be regulated by the strength of the bottleneck. We speculate on the possible implications of the results in the context of the biological phenomenon that motivated the formulation of the theoretical model.

We present an analytic method to find the areas of nonmixing regions in orientation-preserving spherical piecewise isometries (PWIs), and apply it to determine the mixing efficacy of a class of spherical PWIs derived from granular flow in a biaxial tumbler. We show that mixing efficacy has a complex distribution across the protocol space, with local minima in mixing efficacy, termed resonances, that can be determined analytically. These resonances are caused by the interaction of two mode-locking-like phenomena.

The paper derives a relation for the viscosity of undercooled liquids on the basis of the pragmatical model concept of Eshelby relaxations with a finite lifetime. From accurate shear relaxation data in the literature, one finds that slightly less than half of the internal stresses relax directly via single Eshelby relaxations; the larger part dissolves at the terminal lifetime, which is a combined effect of many Eshelby relaxations.

Picosecond time-resolved x-ray spectroscopy is used to measure the spectral line shift of the $1s2p\text{−}1s^2$ transition in He-like Al ions as a function of the instantaneous plasma conditions. The plasma temperature and density are inferred from the Al $\mathrm{H}{\mathrm{e}}_{\alpha }$ complex using a nonlocal-thermodynamic-equilibrium atomic physics model. The experimental spectra show a linearly increasing redshift for electron densities of $1-5\times {10}^{23}\mathrm{c}{\mathrm{m}}^{-3}$. The measured line shifts are broadly consistent with a generalized analytic line-shift model based on calculations of a self-consistent field ion-sphere model.

Chemical potential, as an important thermodynamic quantity, has been popularly used in thermodynamic modeling for complex systems, especially for those involving the phase transitions and chemical reactions. Here we present a chemical-potential-based multiphase lattice Boltzmann model, in which the nonideal force is directly evaluated by a chemical potential. The numerical computation is more efficient than the pressure-tensor-based model [Wen et al. Europhys. Lett. 112, 44002 (2015)] because the calculations of the pressure tensor and its divergence are avoided. We have derived several chemical potentials of the popular equations of state from the free-energy density function. The theoretical analyses and numerical results support that the present model satisfies thermodynamics and Galilean invariance. An effective chemical-potential boundary condition is also implemented to investigate the wettability of a solid surface, and the contact angle can be linearly tuned by the surface chemical potential.

Entropic lattice Boltzmann methods have been developed to alleviate intrinsic stability issues of lattice Boltzmann models for under-resolved simulations. Its reliability in combination with moving objects was established for various laminar benchmark flows in two dimensions in our previous work [B. Dorschner, S. Chikatamarla, F. Bösch, and I. Karlin, J. Comput. Phys. 295, 340 (2015)] as well as for three-dimensional one-way coupled simulations of engine-type geometries in $\mathrm{B}$. Dorschner, F. Bösch, S. Chikatamarla, K. Boulouchos, and I. Karlin [J. Fluid Mech. 801, 623 (2016)] for flat moving walls. The present contribution aims to fully exploit the advantages of entropic lattice Boltzmann models in terms of stability and accuracy and extends the methodology to three-dimensional cases, including two-way coupling between fluid and structure and then turbulence and deforming geometries. To cover this wide range of applications, the classical benchmark of a sedimenting sphere is chosen first to validate the general two-way coupling algorithm. Increasing the complexity, we subsequently consider the simulation of a plunging SD7003 airfoil in the transitional regime at a Reynolds number of $\mathrm{Re}=40000$ and, finally, to access the model's performance for deforming geometries, we conduct a two-way coupled simulation of a self-propelled anguilliform swimmer. These simulations confirm the viability of the new fluid-structure interaction lattice Boltzmann algorithm to simulate flows of engineering relevance.

Modeling of multiphase compositional hydrodynamics at nanoscale is performed by means of density functional hydrodynamics (DFH). DFH is the method based on density functional theory and continuum mechanics. This method has been developed by the authors over 20 years and used for modeling in various multiphase hydrodynamic applications. In this paper, DFH was further extended to encompass phenomena inherent in liquids at nanoscale. The new DFH extension is based on the introduction of external potentials for chemical components. These potentials are localized in the vicinity of solid surfaces and take account of the van der Waals forces. A set of numerical examples, including disjoining pressure, film precursors, anomalous rheology, liquid in contact with heterogeneous surface, capillary condensation, and forward and reverse osmosis, is presented to demonstrate modeling capabilities.

Entropy measures are widely applied to quantify the complexity of dynamical systems in diverse fields. However, the practical application of entropy methods is challenging, due to the variety of entropy measures and estimators and the complexity of real-world time series, including nonstationarities and long-range correlations (LRC). We conduct a systematic study on the performance, bias, and limitations of three basic measures (entropy, conditional entropy, information storage) and three traditionally used estimators (linear, kernel, nearest neighbor). We investigate the dependence of entropy measures on estimator- and process-specific parameters, and we show the effects of three types of nonstationarities due to artifacts (trends, spikes, local variance change) in simulations of stochastic autoregressive processes. We also analyze the impact of LRC on the theoretical and estimated values of entropy measures. Finally, we apply entropy methods on heart rate variability data from subjects in different physiological states and clinical conditions. We find that entropy measures can only differentiate changes of specific types in cardiac dynamics and that appropriate preprocessing is vital for correct estimation and interpretation. Demonstrating the limitations of entropy methods and shedding light on how to mitigate bias and provide correct interpretations of results, this work can serve as a comprehensive reference for the application of entropy methods and the evaluation of existing studies.

When in an ensemble of oscillatory units the interaction occurs through a diffusion-like manner, the intrinsic oscillations can be quenched through two structurally different scenarios: amplitude death (AD) and oscillation death (OD). Unveiling the underlying principles of stable rhythmic activity against AD and OD is a challenging issue of substantial practical significance. Here, by developing a low-pass filter (LPF) to track the output signals of the local system in the coupling, we show that it can revoke both AD and OD, and even the AD to OD transition, thereby giving rise to oscillations in coupled nonlinear oscillators under diverse death scenarios. The effectiveness of the local LPF is proven to be valid in an arbitrary network of coupled oscillators with distributed propagation delays. The constructive role of the local LPF in revoking deaths provides a potential dynamic mechanism of sustaining a reliable rhythmicity in real-world systems.

A rattleback is a rigid, semielliptic toy which exhibits unintuitive behavior; when it is spun in one direction, it soon begins pitching and stops spinning, then it starts to spin in the opposite direction, but in the other direction, it seems to spin just steadily. This puzzling behavior results from the slight misalignment between the principal axes for the inertia and those for the curvature; the misalignment couples the spinning with the pitching and the rolling oscillations. It has been shown that under the no-slip condition and without dissipation the spin can reverse in both directions, and Garcia and Hubbard obtained the formula for the time required for the spin reversal $t_r$ [Proc. R. Soc. Lond. A 418, 165 (1988)]. In this work, we reformulate the rattleback dynamics in a physically transparent way and reduce it to a three-variable dynamics for spinning, pitching, and rolling. We obtain an expression of the Garcia-Hubbard formula for $t_r$ by a simple product of four factors: (1) the misalignment angle, (2) the difference in the inverses of inertia moment for the two oscillations, (3) that in the radii for the two principal curvatures, and (4) the squared frequency of the oscillation. We perform extensive numerical simulations to examine validity and limitation of the formula, and find that (1) the Garcia-Hubbard formula is good for both spinning directions in the small spin and small oscillation regime, but (2) in the fast spin regime especially for the steady direction, the rattleback may not reverse and shows a rich variety of dynamics including steady spinning, spin wobbling, and chaotic behavior reminiscent of chaos in a dissipative system.