Source author record

Hisao Hayakawa

Hisao Hayakawa appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

85works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

85 published item(s)

preprint2022arXiv

Geometrical Quantum Chemical Engine

We propose a geometrical engine undergoing an adiabatic (Thouless) pumping process for a small system connected to external isothermal reservoirs with the control of electrochemical potentials of the reservoirs and one parameter in the system Hamiltonian. Thanks to the geometrical nature of this process, the entropy production is characterized by the geometric metric tensor which is connected to the Fisher information and the Hessian of the density matrix in a nonequilibrium steady state. The existence of an inequality between the thermodynamic length and entropy production is established. We also establish that the work done on this system is characterized by a vector potential and is equivalent to the thermodynamic flux. To characterize the engine, we the introduce effective efficiency as the relation between the work and entropy production. Through the theoretical analysis of the quantum master equation for the Anderson model of a quantum dot within the wide-band approximation, we illustrate the explicit values of the work, thermodynamic length, and effective efficiency of the engine as functions of the phase difference of the externally controlled electrochemical potentials.

preprint2022arXiv

Softening and residual loss modulus of jammed grains under oscillatory shear in an absorbing state

From a theoretical study of the mechanical response of jammed materials comprising frictionless and overdamped particles under oscillatory shear, we find that the material becomes soft, and the loss modulus remains finite even in an absorbing state where any irreversible plastic deformation does not exist. The trajectories of the particles in this region exhibit hysteresis loops. We succeed in clarifying the origin of the softening of the material and the residual loss modulus with the aid of Fourier analysis. We also clarify the roles of the yielding point in the softening to distinguish the plastic deformation from reversible deformation in the absorbing state.

preprint2021arXiv

Geometrical Formulation of Adiabatic Pumping as a Heat Engine

We investigate a heat engine under an adiabatic (Thouless) pumping process. In this process, the extracted work and lower bound on dissipated availability are characterized by a vector potential and a Riemannian metric tensor, respectively. We derive a trade-off relation between the power and effective efficiency. We also explicitly calculate the trade-off relation as well as the power and effective efficiency for a spin-boson model coupled to two reservoirs.

preprint2021arXiv

Mpemba effect in inertial suspensions

The Mpemba effect (a counterintuitive thermal relaxation process where an initially hotter system may cool down to the steady state sooner than an initially colder system) is studied in terms of a model of inertial suspensions under shear. The relaxation to a common steady state of a suspension initially prepared in a quasi-equilibrium state is compared with that of a suspension initially prepared in a nonequilibrium sheared state. Two classes of Mpemba effect are identified, the normal and the anomalous one. The former is generic, in the sense that the kinetic temperature starting from a cold nonequilibrium sheared state is overtaken by the one starting from a hot quasi-equilibrium state, due to the absence of initial viscous heating in the latter, resulting in a faster initial cooling. The anomalous Mpemba effect is opposite to the normal one since, despite the initial slower cooling of the nonequilibrium sheared state, it can eventually overtake an initially colder quasi-equilibrium state. The theoretical results based on kinetic theory agree with those obtained from event-driven simulations for inelastic hard spheres. It is also confirmed the existence of the inverse Mpemba effect, which is a peculiar heating process, in these suspensions. More particularly, we find the existence of a mixed process in which both heating and cooling can be observed during relaxation.

preprint2020arXiv

Enskog kinetic theory of rheology for a moderately dense inertial suspension

The Enskog kinetic theory for moderately dense inertial suspensions under simple shear flow is considered as a model to analyze the rheological properties of the system. The influence of the background fluid on suspended particles is modeled via a viscous drag force plus a Langevin-like term defined in terms of the background temperature. In a previous paper [Hayakawa et al., Phys. Rev. E 96, 042903 (2017)], Grad's moment method with the aid of a linear shear-rate expansion was employed to obtain a theory which gave good agreement with the results of event-driven Langevin simulations of hard spheres for low densities and/or small shear rates. Nevertheless, the previous approach had a limitation of applicability to the high shear-rate and high density regime. Thus, in the present paper, we extend the previous work and develop Grad's theory including higher order terms in the shear rate. This improves significantly the theoretical predictions, a quantitative agreement between theory and simulation being found in the high-density region (volume fractions smaller than or equal to $0.4$).

preprint2020arXiv

Fluctuation Relations For Adiabatic Pumping

We derive an extended fluctuation relation for an open system coupled with two reservoirs under adiabatic one-cycle modulation. We confirm that the geometric phase caused by the Berry-Sintisyn-Nemenman curvature in the parameter space generates non-Gaussian fluctuations. This non-Gaussianity is enhanced for the instantaneous fluctuation relation when the bias between the two reservoirs disappears.

preprint2020arXiv

Impact-induced hardening in dense frictional suspensions

We numerically study the impact-induced hardening in dense suspensions. We employ the lattice Boltzmann method and perform simulations of dense suspensions under impacts, which incorporate the contact between suspended particles with the free surface of the suspension. Our simulation for a free-falling impactor on a dense suspension reproduces experimental results, where rebound takes place for frictional particles at high-speed impact and high volume fraction shortly after the impact before subsequently sinking. We found that the shear stress of the suspension is not affected by the impact, which clearly distinguishes the impact-induced hardening from the discontinuous shear thickening. Instead, we found the existence of a localized region with distinctively high value of normal stress corresponding to the dynamically jammed region. Our simulation indicates that the frictional interaction between suspended particles is important for the impact-induced hardening to maintain the dynamically jammed region. Furthermore, persistent homology analysis successfully elucidates the topological structure of force chains.

preprint2020arXiv

Kinetic theory of shear thickening for a moderately dense gas-solid suspension: from discontinuous thickening to continuous thickening

The Enskog kinetic theory for moderately dense gas-solid suspensions under simple shear flow is considered as a model to analyze the rheological properties of the system. The influence of the environmental fluid on solid particles is modeled via a viscous drag force plus a stochastic Langevin-like term. The Enskog equation is solved by means of two independent but complementary routes: (i) Grad's moment method and (ii) event-driven Langevin simulation of hard spheres. Both approaches clearly show that the flow curve (stress-strain rate relation) depends significantly on the volume fraction of the solid particles. In particular, as the density increases, there is a transition from the discontinuous shear thickening (observed in dilute gases) to the continuous shear thickening for denser systems. The comparison between theory and simulations indicate that while the theoretical predictions for the kinetic temperature agree well with simulations for densities $φ\lesssim 0.5$, the agreement for the other rheological quantities (the viscosity, the stress ratio and the normal stress differences) is limited to more moderate densities ($φ\lesssim 0.3$) if the inelasticity during collisions between particles is not large. [This paper has been published in Phys. Rev. E {\bf 96}, 42903 (2017) but we have realized that there are some typos and mistakes after its publication. So we add the Erratum which will be published in PRE in the top of this paper.]

preprint2020arXiv

Nonadiabatic Control of Geometric Pumping

We study nonadiabatic effects of geometric pumping. With arbitrary choices of periodic control parameters, we go beyond the adiabatic approximation to obtain the exact pumping current. We find that a geometrical interpretation for the nontrivial part of the current is possible even in the nonadiabatic regime. The exact result allows us to find a smooth connection between the adiabatic Berry phase theory at low frequencies and the Floquet theory at high frequencies. We also study how to control the geometric current. Using the method of shortcuts to adiabaticity with the aid of an assisting field, we illustrate that it enhances the current.

preprint2020arXiv

Scaling law of the drag force in dense granular media

Making use of the system of pulling a spherical intruder in static three-dimensional granular media, we numerically study the scaling law for the drag force $F_{\rm drag}$ acting on the moving intruder under the influence of the gravity. Suppose if the intruder of diameter $D$ immersed in a granular medium consisting of grains of average diameter $d$ is located at a depth $h$ and moves with a speed $V$, we find that $F_{\rm drag}$ can be scaled as $(D+d)^{ϕ_μ} h^{α_μ}$ with two exponents $ϕ_μ$ and $α_μ$, which depend on the friction coefficient $μ$ and satisfy an approximate sum rule $ϕ_μ+α_μ\approx 3$. This scaling law is valid for the arbitrary Froude number (defined by $\mathrm{Fr}={2 V}\sqrt{{2D}/{g}}\big/(D+d)$), if $h$ is sufficiently deep. We also identify the existence of three regimes (quasistatic, linear, and quadratic) at least for frictional grains in the velocity dependence of drag force. The crossovers take place at $\mathrm{Fr}\approx 1$ between the quasistatic to the linear regimes and at $\mathrm{Fr}\approx 5$ between the linear to the quadratic regimes. We also observe that Froude numbers at which these crossovers between the regimes happen are independent of the depth $h$ and the diameter of the intruder $D$. We also report the numerical results on the average coordination number of the intruder and average contact angle as functions of intruder velocity.

preprint2020arXiv

Scaling laws for frictional granular materials confined by constant pressure under oscillatory shear

Herein, we numerically study the rheology of a two-dimensional frictional granular system confined by constant pressure under oscillatory shear. Several scaling laws for the storage and loss moduli against the scaled strain amplitude have been found. The scaling laws in plastic regime for large strain amplitude can be understood by the angular distributions of the contact force. The scaling exponents are estimated by considering the physical mechanism.

preprint2020arXiv

Shear jamming, discontinuous shear thickening, and fragile states in dry granular materials under oscillatory shear

We numerically study the linear response of two-dimensional frictional granular materials under oscillatory shear. The storage modulus $G'$ and the loss modulus $G''$ in the zero strain rate limit depend on the initial strain amplitude of the oscillatory shear before measurement. The shear jammed state (satisfying $G'>0$) can be observed at an amplitude greater than a critical initial strain amplitude. The fragile state is defined by the emergence of liquid-like and solid-like states depending on the form of the initial shear. In this state, the observed $G'$ after the reduction of the strain amplitude depends on the phase of the external shear strain. The loss modulus $G''$ exhibits a discontinuous jump corresponding to discontinuous shear thickening in the fragile state.

preprint2019arXiv

Simulation of dense non-Brownian suspensions with the lattice Boltzmann method: Shear jammed and fragile states

Dense non-Brownian suspensions including both the hydrodynamic interactions and the frictional contacts between particles are numerically studied under simple and oscillatory shears in terms of the lattice Boltzmann method. We successfully reproduce the discontinuous shear thickening (DST) under a simple shear for bulk three-dimensional systems. For our simulation of an oscillatory shear in a quasi-two-dimensional system, we measure the mechanical response when we reduce the strain amplitude after the initial oscillations with a larger strain amplitude. Here, we find the existence of the shear-jammed state under this protocol in which the storage modulus $G^{\prime}$ is only finite for high initial strain amplitude $γ_0^{I}$. We also find the existence of the fragile state in which both fluid-like and solid-like responses can be detected for an identical area fraction and an initial strain amplitude $γ_0^{I}$ depending on the initial phase $Θ$ (or the asymmetricity of the applied strain) of the oscillatory shear. We also observe the DST-like behavior under the oscillatory shear in the fragile state. Moreover, we find that the stress anisotropy becomes large in the fragile state. Finally, we confirm that the stress formula based on the angular distribution of the contact force recovers the contact contributions to the stress tensors for both simple and oscillatory shears with large strains.

preprint2016arXiv

Active cage model of glassy dynamics

We build up a phenomenological picture in terms of the effective dynamics of a tracer confined in a cage experiencing random hops to capture somec haracteristics of glassy systems. This minimal description exhibits scale invariance properties for the small-displacement distribution that echo experimental observations. We predict the existence of exponential tails as a cross-over between two Gaussian regimes. Moreover, we demonstrate that the onset of glassy behavior is controlled only by two dimensionless numbers: the number of hops occurring during the relaxation of the particle within a local cage, and the ratio of the hopping length to the cage size.

preprint2016arXiv

Efficiency at maximum power output for an engine with a passive piston

Efficiency at maximum power (MP) output for an engine with a passive piston without mechanical controls between two reservoirs is theoretically studied. We enclose a hard core gas partitioned by a massive piston in a temperature-controlled container and analyze the efficiency at MP under a heating and cooling protocol without controlling the pressure acting on the piston from outside. We find the following three results: (i) The efficiency at MP for a dilute gas is close to the Chambadal-Novikov-Curzon-Ahlborn (CNCA) efficiency if we can ignore the side wall friction and the loss of energy between a gas particle and the piston, while (ii) the efficiency for a moderately dense gas becomes smaller than the CNCA efficiency even when the temperature difference of reservoirs is small. (iii) Introducing the Onsager matrix for an engine with a passive piston, we verify that the tight coupling condition for the matrix of the dilute gas is satisfied, while that of the moderately dense gas is not satisfied because of the inevitable heat leak. We confirm the validity of these results using the molecular dynamics simulation and introducing an effective mean-field-like model which we call stochastic mean field model.

preprint2016arXiv

Fulde-Ferrell-Larkin-Ovchinnikov states in a superconducting ring with magnetic fields: Phase diagram and the first-order phase transitions

We find the angular Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states (or the twisted kink crystals) in which a phase and an amplitude of a pair potential modulate simultaneously in a quasi-onedimensional superconducting ring with a static Zeeman magnetic field applied on the ring and static Aharonov-Bohm magnetic flux penetrating the ring. The superconducting ring with magnetic flux produces a persistent current, whereas the Zeeman split of Fermi energy results in the spatial modulation of the pair potential. We show that these two magnetic fields stabilize the FFLO phase in a large parameter region of the magnetic fields. We further draw the phase diagram with the two kinds of first-order phase transitions; one corresponds to phase slips separating the Aharonov-Bohm magnetic flux, and the other separates the number of peaks of the pair amplitude for the Zeeman magnetic field.

preprint2016arXiv

Granular rotor as a probe for a non-equilibrium bath

This study numerically and analytically investigates the dynamics of a rotor under viscous or dry friction as a non-equilibrium probe of a granular gas. In order to demonstrate the role of the rotor as a probe for a non-equilibrium bath, the molecular dynamics (MD) simulation of the rotor is performed under viscous or dry friction surrounded by a steady granular gas under gravity. A one- to-one map between the velocity distribution function (VDF) of the granular gas and the angular distribution function for the rotor is theoretically derived. The MD simulation demonstrates that the one-to-one map accurately infers the local VDF of the granular gas from the angular VDF of the rotor, and vice versa.

preprint2016arXiv

Kinetic theory for dilute cohesive granular gases with a square well potential

We develop the kinetic theory of dilute cohesive granular gases in which the attractive part is described by a square well potential. We derive the hydrodynamic equations from the kinetic theory with the microscopic expressions for the dissipation rate and the transport coefficients. We check the validity of our theory by performing the direct simulation Monte Carlo.

preprint2015arXiv

Divergence of Viscosity in Jammed Granular Materials: A Theoretical Approach

A theory for jammed granular materials is developed with the aid of a nonequilibrium steady-state distribution function. The approximate nonequilibrium steady-state distribution function is explicitly given in the weak dissipation regime by means of the relaxation time. The theory quantitatively agrees with the results of the molecular dynamics simulation on the critical behavior of the viscosity below the jamming point without introducing any fitting parameter.

preprint2015arXiv

Drag Law of Two Dimensional Granular Fluids

The drag force law acting on a moving circular disk in a two-dimensional granular medium is analyzed based on the discrete element method (DEM). It is remarkable that the drag force on the moving disk in moderate dense and pure two-dimensional granular medium can be well reproduced by a perfect fluid with separation from the surface of the tracer. A yield force, being independent of the moving speed of the disk, appears if a dry friction between the granular disks and the bottom plate exists. The perfect fluidity is violated in this case. The yield force and the drag force diverge at the jamming point.

preprint2015arXiv

Minimal Model of Stochastic Athermal Systems: Origin of Non-Gaussian Noise

For a wide class of stochastic athermal systems, we derive Langevin-like equations driven by non-Gaussian noise, starting from master equations and developing a new asymptotic expansion. We found an explicit condition whereby the non-Gaussian properties of the athermal noise become dominant for tracer particles associated with both thermal and athermal environments. Furthermore, we derive an inverse formula to infer microscopic properties of the athermal bath from the statistics of the tracer particle. We apply our formulation to a granular motor under viscous friction, and analytically obtain the angular velocity distribution function. Our theory demonstrates that the non-Gaussian Langevin equation is the minimal model of athermal systems.

preprint2015arXiv

Size distribution of particles in Saturn's rings from aggregation and fragmentation

Saturn's rings consist of a huge number of water ice particles, with a tiny addition of rocky material. They form a flat disk, as the result of an interplay of angular momentum conservation and the steady loss of energy in dissipative inter-particle collisions. For particles in the size range from a few centimeters to a few meters, a power-law distribution of radii, $\sim r^{-q}$ with $q \approx 3$, has been inferred; for larger sizes, the distribution has a steep cutoff. It has been suggested that this size distribution may arise from a balance between aggregation and fragmentation of ring particles, yet neither the power-law dependence nor the upper size cutoff have been established on theoretical grounds. Here we propose a model for the particle size distribution that quantitatively explains the observations. In accordance with data, our model predicts the exponent $q$ to be constrained to the interval $2.75 \le q \le 3.5$. Also an exponential cutoff for larger particle sizes establishes naturally with the cutoff-radius being set by the relative frequency of aggregating and disruptive collisions. This cutoff is much smaller than the typical scale of micro-structures seen in Saturn's rings.

preprint2014arXiv

Asymptotic derivation of Langevin-like equation with non-Gaussian noise and its analytical solution

We asymptotically derive a non-linear Langevin-like equation with non-Gaussian white noise for a wide class of stochastic systems associated with multiple stochastic environments, by developing the expansion method in our previous paper [K. Kanazawa et al., arXiv: 1407.5267 (2014)]. We further obtain a full-order asymptotic formula of the steady distribution function in terms of a large friction coefficient for a non-Gaussian Langevin equation with an arbitrary non-linear frictional force. The first-order truncation of our formula leads to the independent-kick model and the higher-order correction terms directly correspond to the multiple-kicks effect during relaxation. We introduce a diagrammatic representation to illustrate the physical meaning of the high-order correction terms. As a demonstration, we apply our formula to a granular motor under Coulombic friction and get good agreement with our numerical simulations.

preprint2014arXiv

Avalanche contribution to shear modulus of granular materials

Shear modulus of frictionless granular materials near the jamming transition under oscillatory shear is numerically investigated. It is found that the shear modulus $G$ satisfies a scaling law to interpolate between $G \sim (ϕ- ϕ_J)^{1/2}$ and $G \sim γ_0^{-1/2}(ϕ-ϕ_J)$ for a linear spring model of the elastic interaction between contacting grains, where $ϕ$, $ϕ_J$, and $γ_0$ are, respectively, the volume fraction of grains, the fraction at the jamming point, and the amplitude of the oscillatory shear. The linear relation between the shear modulus and $ϕ- ϕ_J$ can be understood by slip avalanches.

preprint2014arXiv

Energetics of active fluctuations in living cells

The nonequilibrium activity taking place in a living cell can be monitored with a tracer embedded in the medium. While microrheology experiments based on optical manipulation of such probes have become increasingly standard, we put forward a number of experiments with alternative protocols that, we claim, will provide new insight into the energetics of active fluctuations. These are based on either performing thermodynamic--like cycles in control-parameter space, or on determining response to external perturbations of the confining trap beyond simple translation. We illustrate our proposals on an active itinerant Brownian oscillator modeling the dynamics of a probe embedded in a living medium.

preprint2014arXiv

Energy pumping in electrical circuits under avalanche noise

We theoretically study energy pumping processes in an electrical circuit with avalanche diodes, where non-Gaussian athermal noise plays a crucial role. We show that a positive amount of energy (work) can be extracted by an external manipulation of the circuit in a cyclic way, even when the system is spatially symmetric. We discuss the properties of the energy pumping process for both quasi-static and fnite-time cases, and analytically obtain formulas for the amounts of the work and the power. Our results demonstrate the significance of the non-Gaussianity in energetics of electrical circuits.

preprint2014arXiv

Geometric pumping induced by shear flow in dilute liquid crystalline polymer solutions

We investigate nonlinear rheology of dilute liquid crystalline polymer solutions under time dependent two-directional shear flow. We analyze the Smoluchowski equation, which describes the dynamics of the orientation of a liquid crystalline polymer, by employing technique of the full counting statistics. In the adiabatic limit, we derive the expression for time integrated currents generated by a Berry-like curvature. Using this expression, it is shown that the expectation values of the time-integrated angular velocity of a liquid crystalline polymer and the time-integrated stress tensor are generally not zero even if the time average of the shear rate is zero. The validity of the theoretical calculations is confirmed by direct numerical simulations of the Smoluchowski equation. Nonadiabatic effects are also investigated by simulations and it is found that the time-integrated stress tensor depends on the speed of the modulation of the shear rate if we adopt the isotropic distribution as an initial state.

preprint2014arXiv

Non-adiabatic effect in quantum pumping for a spin-boson system

We clarify the role of non-adiabatic effects in a quantum pumping for a spin-boson system. When we sinusoidally control the temperatures of two reservoirs with π/2 phase difference, we find that the pumping current strongly depends on the initial condition, and thus, the current deviates from that predicted by the adiabatic treatment. We also analytically obtain the contribution of non-adiabatic effects in the pumping current proportional to Ω^3 where Ωis the angular frequency of the temperature control. The validity of the analytic expression is verified by our numerical calculation. Moreover, we extend the steady heat fluctuation theorem to the case for slowly modulated temperatures and large transferred energies.

preprint2014arXiv

Non-divergent representation of non-Hermitian operator near the exceptional point with application to a quantum Lorentz gas

We propose a non-singular representation for a non-Hermitian operator even if the parameter space contains exceptional points (EPs), at which the operator cannot be diagonalized and the usual spectral representation ceases to exist. Our representation has a generalized Jordan block form and is written in terms of extended pseudo-eigenstates. Our method is free from a divergence in the spectral representation at EPs, at which multiple eigenvalues and eigenvectors coalesce and the eigenvectors cannot be normalized. Our representation improves the accuracy of numerical calculations of physical quantities near EPs. We also find that our method is applicable to various problems related to EPs in the parameter space of non-Hermitian operators. We demonstrate the usefulness of our representation by investigating Boltzmann's collision operator in a one-dimensional quantum Lorentz gas in the weak coupling approximation.

preprint2014arXiv

Rheology of Dense Sheared Granular Liquids

The rheology of dense sheared granular liquids is investigated based on the mode-coupling theory (MCT). This extended MCT includes correlations for the density-current mode as well as the density-density correlation mode, and a self-consistent coupling equation for the energy balance condition. The extended MCT exhibits disappearance of the two-step relaxation of the density-density correlation function, and also successfully reproduces the density dependence of the shear viscosity for volume fractions between 0.50 and 0.60, if we shift the density. However, it predicts unphysical tendency for the granular temperature. The cause of this drawback and the possibilities of its amendment are discussed.

preprint2014arXiv

Roles of Dry Friction in Fluctuating Motion of Adiabatic Piston

The motion of an adiabatic piston under dry friction is investigated to clarify the roles of dry friction in non-equilibrium steady states. We clarify that dry friction can reverse the direction of the piston motion and causes a discontinuity or a cusp-like singularity for velocity distribution functions of the piston. We also show that the heat fluctuation relation is modified under dry friction.

preprint2014arXiv

Simulation of cohesive fine powders under a plane shear

Three dimensional molecular dynamics simulations of cohesive dissipative powders under a plane shear are performed. We find the various phases depending on the dimensionless shear rate and the dissipation rate as well as the density. We also find that the shape of clusters depends on the initial condition of velocities of particles when the dissipation is large. Our simple stochastic model reproduces the non-Gaussian velocity distribution function appearing in the coexistence phase of a gas and a plate.

preprint2013arXiv

Analytical expression of geometrical pumping for a quantum dot based on quantum master equation

We analytically investigate a non-equilibrium quantum pumping for a single quantum dot connected to external leads on the basis of the quantum master equation (QME). We show that the Coulomb interaction associated with the spin effect in the dot induces the Berry-like phase in the parameter space and this phase results in the excess charge transfer for the cyclic modulation of parameters in leads. We obtain an analytical expression of the curvature of the phase and that for the pumped currents.

preprint2013arXiv

Breakdown of a renormalized perturbation expansion around mode-coupling theory of the glass transition

We analyze a renormalized perturbation expansion around the mode-coupling theory of the glass transition. We focus on the long-time limit of the irreducible memory function. We discuss a renormalized diagrammatic expansion for this function and re-sum two infinite classes of diagrams. We show that the resulting contributions to the irreducible memory function diverge at the mode-coupling transition. A further re-summation of ladder diagrams constructed by iterating these divergent contributions gives a finite result which cancels the mode-coupling theory's expression for the irreducible memory function.

preprint2013arXiv

Constitutive Relations of Jammed Frictionless Granular Materials under Oscillatory Shear

We numerically investigate the rheological properties of jammed frictionless granular materials under an oscillatory shear in terms of our simulation of the distinct element method. It is demonstrated that the constitutive relation between the shear stress and the shear strain of grains strongly depends on the amplitude and the frequency of the oscillatory shear. For a small amplitude region, the granular material behaves as a visco-elastic material characterized by the Kelvin-Voigt model, while it behaves as a yield stress liquid for a large amplitude region. We find that the rheology of grains under the large amplitude oscillatory shear can be described by a phenomenological constitutive model.

preprint2013arXiv

Effect of Elastic Vibrations on Normal Head on Collisions of Isothermal Spheres

We numerically investigate head on collisions of isothermal visco-elastic spheres. We find that the restitution coefficient oscillates against the impact speed if the solid viscosity inside the sphere is small enough. We confirm that the oscillation arises from the resonance between the duration of contact and the eigen-frequencies of the sphere. This oscillation disappears if there exists the strong solid viscosity in spheres. We also find that a sinusoidal behavior of the restitution coefficient against the initial phase in the eigenmodes for collisions between a thermally activated sphere and a flat wall. As a result, the restitution coefficient can exceed unity if the impact speed of the colliding sphere is nearly equal to or slower than the thermal speed. We have confirmed the existence of the fluctuation theorem for impact processes through our simulation.

preprint2013arXiv

Geometrical Excess Entropy Production in Nonequilibrium Quantum Systems

For open systems described by the quantum Markovian master equation, we study a possible extension of the Clausius equality to quasistatic operations between nonequilibrium steady states (NESSs). We investigate the excess heat divided by temperature (i.e., excess entropy production) which is transferred into the system during the operations. We derive a geometrical expression for the excess entropy production, which is analogous to the Berry phase in unitary evolution. Our result implies that in general one cannot define a scalar potential whose difference coincides with the excess entropy production in a thermodynamic process, and that a vector potential plays a crucial role in the thermodynamics for NESSs. In the weakly nonequilibrium regime, we show that the geometrical expression reduces to the extended Clausius equality derived by Saito and Tasaki (J. Stat. Phys. {\bf 145}, 1275 (2011)). As an example, we investigate a spinless electron system in quantum dots. We find that one can define a scalar potential when the parameters of only one of the reservoirs are modified in a non-interacting system, but this is no longer the case for an interacting system.

preprint2013arXiv

Geometrical pumping in spin coupled double quantum dots

We analytically investigate the non-equilibrium pumping for a double quantum dots system on the basis of the quantum master equation (QME), where the double quantum dots are connected to two external leads by the spin coupling. Each of leads has two tunable parameters, temperature and chemical potential. Using QME formalism, we obtain an analytical expression of the Berry-like phase for eigenstate of the QME in the parameter space. We show that the Berry-like curvature is non zero in the whole region in the parameter space. We also show that the Berry-like curvature vanishes in the case of single quantum dot and the case of isolated two dots.

preprint2013arXiv

Heat conduction induced by non-Gaussian athermal fluctuations

We study the properties of heat conduction induced by non-Gaussian noises from athermal environments. We find that new terms should be added to the conventional Fourier law and the fluctuation theorem for the heat current, where its average and fluctuation are determined not only by the noise intensities but also by the non-Gaussian nature of the noises. Our results explicitly show the absence of the zeroth law of thermodynamics in athermal systems.

preprint2013arXiv

Jet-induced jammed states of granular jet impacts

The impacts of granular jets for both frictional and frictionless grains in two dimensions are numerically investigated. A dense flow with a dead zone emerges during the impact. From our two-dimensional simulation, we evaluate the equations of state and the con- stitutive equations of the flow. The asymptotic divergences of pressure and shear stress similar to the situation near the jamming transition appear for the frictionless case, while their exponents are smaller than those of the sheared granular systems, and are close to the extrapolation from the kinetic theoretical regime. In a similar manner to the jam- ming for frictional grains, the critical density decreases as the friction constant of grains increases. For bi-disperse systems, the effective friction constant defined as the ratio of shear stress to normal stress, monotonically increases from near zero, as the strain rate increases. On the other hand, the effective friction constant has two metastable branches for mono-disperse systems because of the coexistence of a crystallized state and a liquid state.

preprint2013arXiv

Nonequilibrium identities and response theory for dissipative particles

We derive some nonequilibrium identities such as the integral fluctuation theorem and the Jarzynski equality starting from a nonequilibrium state for dissipative classical systems. Thanks to the existence of the integral fluctuation theorem we can naturally introduce an entropy-like quantity for dissipative classical systems in far from equilibrium states. We also derive the generalized Green-Kubo formula as a nonlinear response theory for a steady dynamics around a nonequilibrium state. We numerically verify the validity of the derived formulas for sheared frictionless granular particles.

preprint2013arXiv

Nonequilibrium mode-coupling theory for uniformly sheared underdamped systems

Mode-coupling theory for uniformly sheared underdamped systems with an isothermal condition is presented. As a result of the isothermal condition, it is shown that the shear stress exhibits significant relaxation at the alpha-relaxation regime due to the cooling effect which is accompanied by the growth of the current fluctuation. This indicates that nonequilibrium underdamped MCT is not equivalent to the corresponding overdamped MCT, even at long time-scales after the early stage of the beta relaxation. It is verified that the effect of correlations to the density-current modes is, however, negligible in sheared thermostated underdamped systems.

preprint2013arXiv

Numerical analysis of impact processes of granular jets

The rheology of a three-dimensional granular jet during an impact is investigated numerically. The cone-like scattering pattern and the sheet-like pattern observed in an experiment [X. Cheng, et al. Phys. Rev. Lett. 99, 188001 (2007)] can be reproduced through our calculation. We discuss the constitutive equation for granular jet impact in terms of our simulation. From the analysis of an effective friction constant, which is the ratio between the shear stress and the pressure the assumption of the zero yield stress would be natural in our setup and the shear visocity is not small in contrast to the suggestion by the experiment.

preprint2013arXiv

The Legacy of Hideki Yukawa, Sin-itiro Tomonaga, and Shoichi Sakata: Some Aspects from their Archives

Hideki Yukawa, Sin-itiro Tomonaga and Shoichi Sakata pioneered nuclear and particle physics and left enduring legacies. Their friendly collaboration and severe competition laid the foundation to bring up the active postwar generation of nuclear and particle physicists in Japan. In this presentation we illustrate milestones of nuclear and particle physics in Japan from 1930's to mid-1940's which have been clarified in Yukawa Hall Archival Library, Tomonaga Memorial Room and Sakata Memorial Archival Library.

preprint2012arXiv

Critical scaling of jammed system after quench of temperature

Critical behavior of soft repulsive particles after quench of temperature near the jamming trasition is numerically investigated. It is found that the plateau of the mean square displacement of tracer particles and the pressure satisfy critical scaling laws. The critical density for the jamming transition depends on the protocol to prepare the system, while the values of the critical exponents which are consistent with the prediction of a phenomenology are independent of the protocol.

preprint2012arXiv

Geometrical Pumping in Quantum Transport: Quantum Master Equation Approach

For an open quantum system, we investigate the pumped current induced by a slow modulation of control parameters on the basis of the quantum master equation and full counting statistics. We find that the average and the cumulant generating function of the pumped quantity are characterized by the geometrical Berry-phase-like quantities in the parameter space, which is associated with the generator of the master equation. From our formulation, we can discuss the geometrical pumping under the control of the chemical potentials and temperatures of reservoirs. We demonstrate the formulation by spinless electrons in coupled quantum dots. We show that the geometrical pumping is prohibited for the case of non-interacting electrons if we modulate only temperatures and chemical potentials of reservoirs, while the geometrical pumping occurs in the presence of an interaction between electrons.

preprint2012arXiv

Instabilities and turbulence-like dynamics in an oppositely driven binary particle mixture

Using extensive particle-based simulations, we investigate out-of-equilibrium pattern dynamics in an oppositely driven binary particle system in two dimensions. A surprisingly rich dynamical behavior including lane formation, jamming, oscillation and turbulence-like dynamics is found. The ratio of two friction coefficients is a key parameter governing the stability of lane formation. When the friction coefficient transverse to the external force direction is sufficiently small compared to the longitudinal one, the lane structure becomes unstable to shear-induced disturbances, and the system eventually exhibits a dynamical transition into a novel turbulence-like phase characterized by random convective flows. We numerically construct an out-of-equilibrium phase diagram. Statistical analysis of complex spatio-temporal dynamics of the fully nonlinear turbulence-like phase suggests its apparent reminiscence to the swarming dynamics in certain active matter systems.

preprint2012arXiv

Nonequilibrium mode-coupling theory for uniformly sheared underdamped systems

Nonequilibrium mode-coupling theory (MCT) for uniformly sheared underdamped systems is developed, starting from the microscopic thermostatted SLLOD equation, and the corresponding Liouville equation. Special attention is paid to the translational invariance in the sheared frame, which requires an appropriate definition of the transient time-correlators. The derived MCT equation satisfies the alignment of the wavevectors, and is manifestly translationally invariant. Isothermal condition is implemented by the introduction of the current fluctuation in the dissipative coupling to the thermostat. This current fluctation grows in the $α$-relaxation regime, which generates a deviation of the yield stress in the glassy phase from the overdamped case. Response to a perturbation of the shear rate demonstrates an inertia effect which is not observed in the overdamped case. Our theory turns out to be a non-trivial extension of the theory by Fuchs and Cates (J. Rheol. 53(4), 2009) to underdamped systems. Since our starting point is identical to that by Chong and Kim (Phys. Rev. E 79, 2009), the contradictions between Fuchs-Cates and Chong-Kim are resolved.

preprint2012arXiv

Origin of rebounds with a restitution coefficient larger than unity in nanocluster collisions

We numerically investigate the mechanism of super rebounds for head-on collisions between nanoclusters in which the restitution coefficient is larger than unity. It is confirmed that the temperature and the entropy of the nanocluters decrease after the super rebounds by our molecular dynamics simulations. It is also found that the initial metastable structure plays a key role for the emergence of the super rebounds.

preprint2012arXiv

Pattern dynamics of cohesive granular particles under a plane shear

We perform three dimensional molecular dynamics simulations of cohesive granular particles under a plane shear. From the simulations, we found that the granular temperature of the system abruptly decreases to zero after reaching the critical temperature, where the characteristic time $t_{\rm cl}$ is approximately represented by $t_{\rm cl} \propto (ζ-ζ_{\rm cr})^{-β}$ with the dissipation rate $ζ$, the critical dissipation rate $ζ_{\rm cr}$ and the exponent $β\simeq 0.8$. We also found that there exist a variety types of clusters depending on the initial density and the dissipation rate.

preprint2012arXiv

Phase transition in peristaltic transport of frictionless granular particles

Flow of dissipative particles driven by peristaltic motion of a tube is numerically studied. A transition from slow unjammed flow to fast jammed flow is found through the observation of the mass flux if the minimum width of the peristaltic tube is smaller than a critical value. It is also found that average and fluctuation of the transition time, and the peak value of the second moment of the mass flux exhibit power-law divergences at the critical point and these variables satisfy scaling relationships near the critical point. Dependences of the critical width and exponents on peristaltic velocity and on density are also discussed.

preprint2012arXiv

Rheology of sheared granular particles near jamming transition

We investigate the rheology of sheared granular materials near the jamming transition point. We numerically determine the values of the critical fraction and the exponents for the jamming transition using a finite size scaling and the nonlinear minimization method known as the Levenberg-Marquardt algorithm. The exponents are close to our previous theoretical prediction, but there is a small discrepancy, if the critical point is independently determined.

preprint2012arXiv

Simulation of granular jet: Is granular flow really a "perfect fluid?"

We perform three-dimensional simulations of a granular jet impact for both frictional and frictionless grains. Small shear stress observed in the experiment[X. Cheng et al., Phys. Rev. Lett. 99, 188001 (2007) ] is reproduced through our simulation. However, the fluid state after the impact is far from a perfect fluid, and thus, similarity between granular jets and quark gluon plasma is superficial, because the observed viscosity is finite and its value is consistent with the prediction of the kinetic theory.

preprint2012arXiv

Simulation of pattern dynamics of cohesive granular particles under a plane shear

We have performed three-dimensional molecular dynamics simulation of cohesive granular particles under a plane shear. From the simulation, we found the existence of three distinct phases in steady states: (I) a uniform shear phase, (II) a coexistent phase of a shear band and a gas region and (III) a crystal phase. We also found that the critical line between (II) and (III) is approximately represented by $ζ\propto \exp(β\dotγL_y)$, where $ζ$, $β$, $\dotγ$, $L_y$ are the dissipation rate, an unimportant constant, the shear rate and the system size of the velocity gradient direction, respectively.

preprint2012arXiv

Steplike electric conduction in a classical two-dimensional electron system through a narrow constriction in a microchannel

Using molecular dynamics simulation, we investigate transport properties of a classical two-dimensional electron system confined in a microchannel with a narrow constriction. As a function of the confinement strength of the constriction, the calculated conductance in the simulations exhibits steplike increases as reported in a recent experiment [D. G. Rees et al., Phys. Rev. Lett. 106, 026803 (2011)]. It is confirmed that the number of the steps corresponds to the number of stream lines of electrons through the constriction. We verify that density fluctuation plays a major role in smoothing the steps in the conductance.

preprint2012arXiv

Stochastic Energetics for Non-Gaussian Processes

By introducing a new stochastic integral, we investigate the energetics of classical stochastic systems driven by non-Gaussian white noises. In particular, we introduce a decomposition of the total-energy difference into the work and the heat for each trajectory, and derive a formula to calculate the heat from experimental data on the dynamics. We apply our formulation and results to a Langevin system driven by a Poisson noise.

preprint2012arXiv

Time Dependent Ginzburg-Landau Equation for Sheared Granular Flow

The time dependent Ginzburg-Landau equation for a two-dimensional granular shear flow is numerically solved, where we study both the transient dynamics and the steady state of the order parameter. The structural changes of the numerical solutions are qualitatively similar to the shear bands observed in the discrete element method (DEM) simulation of the two-dimensional granular shear flow.

preprint2011arXiv

Critical scaling near jamming transition for frictional granular particles

The critical rheology of sheared frictional granular materials near jamming transition is numer- ically investigated. It is confirmed that there exist a true critical density which characterizes the onset of the yield stress, and two fictitious critical densities which characterize the scaling laws of rheological properties. We find the existence of a hysteresis loop between two of the critical densities for each friction coefficient. It is noteworthy that the critical scaling law for frictionless jamming transition seems to be still valid even for frictional jamming despite using fictitious critical density values.

preprint2011arXiv

Geometrical Expression of Excess Entropy Production

We derive a geometrical expression of the excess entropy production for quasi-static transitions between nonequilibrium steady states of Markovian jump processes, which can be exactly applied to nonlinear and nonequilibrium situations. The obtained expression is geometrical; the excess entropy production depends only on a trajectory in the parameter space, analogous to the Berry phase in quantum mechanics. Our results imply that vector potentials are needed to construct thermodynamics of nonequilibrium steady states.

preprint2010arXiv

Behavior of pressure and viscosity at high densities for two-dimensional hard and soft granular materials

The pressure and the viscosity in two-dimensional sheared granular assemblies are investigated numerically. The behavior of both pressure and viscosity is smoothly changing qualitatively when starting from a mono-disperse hard-disk system without dissipation and moving towards a system of (i) poly-disperse, (ii) soft particles with (iii) considerable dissipation. In the rigid, elastic limit of mono-disperse systems, the viscosity is approximately inverse proportional to the area fraction difference from $ϕ_η \simeq 0.7$, but the pressure is still finite at $ϕ_η$. In moderately soft, dissipative and poly-disperse systems, on the other hand, we confirm the recent theoretical prediction that both scaled pressure (divided by the kinetic temperature $T$) and scaled viscosity (divided by $\sqrt{T}$) diverge at the same density, i.e., the jamming transition point $ϕ_J > ϕ_η$, with the exponents -2 and -3, respectively. Furthermore, we observe that the critical region of the jamming transition becomes invisible as the restitution coefficient approaches unity, i.e. for vanishing dissipation. In order to understand the conflict between these two different predictions on the divergence of the pressure and the viscosity, the transition from soft to hard particles is studied in detail and the dimensionless control parameters are defined as ratios of various time-scales. We introduce a dimensionless number, i.e. the ratio of dissipation rate and shear rate, that can identify the crossover from the scaling of very hard, i.e. rigid disks to the scaling in the soft, jamming regime.

preprint2010arXiv

Representation of the nonequilibrium steady-state distribution function for sheared granular systems

We derive the representation of the nonequilibrium steady-state distribution function which is expressed in terms of the excess free energy production. This representation resembles the one derived recently by Komatsu and Nakagawa [Phys. Rev. Lett. 100 (2008), 030601] resting on the use of microscopic time-reversal symmetry, but our representation applies also to sheared granular systems in which such a symmetry is broken.

preprint2010arXiv

Unified description of long-time tails and long-range correlation functions for sheared granular liquids

Unified description on the long-time tail of velocity autocorrelation function and the long-range correlation for the equal-time spatial correlation functions is developed based on the generalized fluctuating hydrodynamics. The cross-over of the long-time tail from $t^{-3/2}$ to $t^{-5/2}$ is predicted independent of the density, and the equal-time spatial density correlation function and the equal-time spatial velocity correlation function respectively satisfy $r^{-11/3}$ and $r^{-5/3}$ for large $r$ limit.

preprint2009arXiv

Critical behaviors of sheared frictionless granular materials near jamming transition

Critical behaviors of sheared dense and frictionless granular materials in the vicinity of the jamming transition are numerically investigated. From the extensive molecular dynamics simulation, we verify the validity of the scaling theory near the jamming transition proposed by Otsuki and Hayakawa (Prog. Theor. Phys., 121, 647 (2009)). We also clarify the critical behaviors of the shear viscosity and the pair correlation function based on both a phenomenology and the simulation.

preprint2009arXiv

Generalized Green-Kubo formula for a dissipative quantum system

A generalized Green-Kubo formula is derived for a quantum dissipative system of driven Brownian particle, in which the coupling between the system and the environment is linear. The structure is essentially the same as that for the generalized Green-Kubo formula for driven granular particles. It is demonstrated that the correction to the conventional Green-Kubo formula is zero for a free Brownian particle.

preprint2009arXiv

Long-time tails in sheared fluids

The long-time behaviors of the velocity autocorrelation function $C(t)$ for sheared fluids ^M are investigated theoretically and numerically. It is found the existence of the cross-overs of^M $C(t)$ from $t^{-d/2}$ to $t^{-d}$ ^M in sheared fluids of elastic particles without any thermostat, and from $t^{-d/2}$ to $t^{-(d+2)/2}$ ^M in both sheared fluids of elastic particles with a thermostat and sheared granular fluids,^M where $d$ is the spatial dimension. ^M The validity of the predictions has been confirmed by^M our numerical simulations.

preprint2009arXiv

Molecular transport through a bottleneck driven by external force

The transport phenomena of Lennard-Jones molecules through a structural bottleneck driven by an external force are investigated by molecular dynamics simulations. We observe two distinct molecular flow regimes distinguished by a critical external force $F_{c}$ and find scaling behaviors between external forces and flow rates. Below the threshold $F_{c}$, molecules are essentially stuck in the bottleneck due to the attractive interaction between the molecules, while above $F_{c}$, molecules can smoothly move in the pipe. A critical flow rate $q_{c}$ corresponding to $F_{c}$ satisfies a simple relationship with angles and the value of $q_{c}$ can be estimated by a simple argument. We further clarify the role of the temperature dependence in the molecular flows through the bottleneck.

preprint2008arXiv

Super-elastic collisions in a thermally activated system

Impact phenomena of small clusters subject to thermal fluctuations are numerically investigated. From the molecular dynamics simulation for colliding two identical clusters, it is found that the restitution coefficient for head-on collisions can exceed unity when the colliding speed is smaller than the thermal velocity of one ``atom'' of the clusters. The averaged behavior can be understood by the quasi-static theory of impact phenomena. It is also confirmed that our result is governed by the fluctuation theorem.

preprint2008arXiv

Time correlation function of the shear stress in sheared particle systems

The long-time behaviors of the time correlation function $C_η(t)$ of the shear stress for three-dimensional sheared fluids are investigated theoretically. It is found that there are the cross-overs in $C_η(t)$ from $t^{-3/2}$ to $t^{-2}$ for sheared fluids of elastic particles without any thermostat, and from $t^{-3/2}$ to $t^{-5/2}$ for isothermal sheared fluids including granular fluids.

preprint2004arXiv

Kinetic Theory of a Dilute Gas System under Steady Heat Conduction

The velocity distribution function of the steady-state Boltzmann equation for hard-core molecules in the presence of a temperature gradient has been obtained explicitly to second order in density and the temperature gradient. Some thermodynamical quantities are calculated from the velocity distribution function for hard-core molecules and compared with those for Maxwell molecules and the steady-state Bhatnagar-Gross-Krook(BGK) equation. We have found qualitative differences between hard-core molecules and Maxwell molecules in the thermodynamical quantities, and also confirmed that the steady-state BGK equation belongs to the same universality class as Maxwell molecules.

preprint2004arXiv

Kinetic Theory of dilute gases under nonequilibrium conditions

The significance of the recent finding of the velocity distribution function of the steady-state Boltzmann equation under a steady heat current obtained by Kim and Haykawa (J. Phys. Soc. Jpn. {\bf 72}, 1904 (2003)) is discussed. Through the stability analysis, it seems that the steady solution is stable. One of possible applications to the nonequiliburium Knudsen effect in which one cell at equilibrium is connected to another cell under the steady heat conduction is discussed. This solution apparently shows that steady-state thermodynamics proposed by Sasa and Tasaki cannot be used in a naive setup. The preliminary result of our simulation based on molecular-dynamics for nonequilibrium Knudsen effect is also presented to verify the theoretical argument.

preprint2002arXiv

Hydrodynamics for inelastic Maxwell model

Hydrodynamic equations for an inelastic Maxwell model are derived from the inelastic Boltzmann equation based on a systematic Chapman-Enskog perturbative scheme. Transport coefficients appear in Navier-Stokes order have been determined as a function of the restitution coefficient $e$, which cannot be defined for small $e$ as a result of the high energy tail of the velocity distribution function obeying a power law. The dispersion relations for the linearized equation around a homogeneous cooling state have been obtained.

preprint2001arXiv

Simulation and Theory of the Impact of Two-Dimensional Elastic Disks

The impact of a two-dimensional elastic disk with a wall is numerically studied. It is clarified that the coefficient of restitution (COR) decreases with the impact velocity. The result is not consistent with the recent quasi-static theory of inelastic collisions even for very slow impact. This suggests that the elastic model cannot be used in the quasi-static limit. A new quasi-static theory of impacts is proposed, in which the effect of thermal diffusion is dominant. The abrupt decrease of COR has been found due to the plastic deformation of the disk, which is assisted by the initial internal motion. (This paper is to be published in Chemical Engineering Science.)

preprint2001arXiv

Statistical properties of inelastic Lorentz gas

The inelastic Lorentz gas in cooling states is studied. It is found that the inelastic Lorentz gas is localized and that the mean square displacement of the inelastic Lorentz gas obeys a power of a logarithmic function of time. It is also found that the scaled position distribution of the inelastic Lorentz gas has an exponential tail, while the distribution is close to the Gaussian near the peak. Using a random walk model, we derive an analytical expression of the mean square displacement as a function of time and the restitution coefficient, which well agrees with the data of our simulation. The exponential tail of the scaled position distribution function is also obtained by the method of steepest descent.

preprint1995arXiv

Dynamical Simulation of Fluidized Beds --- Hydrodynamically Interacting Granular Particles

A numerical simulation of a gas-fluidized bed is performed without introduction of any empirical parameters. Realistic bubbles and slugs are observed in our simulation. It is found that the convective motion of particles is important for the bubbling phase and there is no convection in the slugging phase. From the simulation results, non-Gaussian distributions are found in the particle velocities and the relation between the deviation from Gaussian and the local density of particles is suggested. It is also shown that the power spectra of particle velocities obey power laws. A brief explanation on the relationship between the simulation results and the Kolmogorov scaling argument is discussed.

preprint1995arXiv

Hydrodynamic Description of Granular Convection

We present a hydrodynamic model that captures the essence of granular dynamics in a vibrating bed. We carry out the linear stability analysis and uncover the instability mechanism that leads to the appearance of the convective rolls via a supercritical bifurcation of a bouncing solution. We also explicitly determine the onset of convection as a function of control parameters and confirm our picture by numerical simulations of the continuum equations.

preprint1995arXiv

Statistical Theory of Sedimentation of Disordered Suspensions

An analytical treatment for the sedimentation rate of disordered suspensions is presented in the context of a resistance problem. From the calculation it is confirmed that the lubrication effect is important in contrast to the suggestion by Brady and Durlofsky (Phys.Fluids Vol.31, 717 (1988)). The calculated sedimentation rate agrees well with the experimental results in all range of the volume fraction.