Source author record

Satoshi Takada

Satoshi Takada 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

13works
5topics
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

13 published item(s)

preprint2023arXiv

Revisiting stress propagation in a two-dimensional elastic circular disk under diametric loading

In this paper, we present a comprehensive investigation of stress propagation in a two-dimensional elastic circular disk. To accurately describe the displacements and stress fields within the disk, we employ a scalar and vector potential approach, representing them as sums of Bessel functions. The determination of the coefficients for these expansions is accomplished in the Laplace space, where we compare the boundary conditions. By converting the inverse Laplace transforms into complex integrals using residue calculus, we successfully derive explicit expressions for the displacements and stress fields. Notably, these expressions encompass primary, secondary, and surface waves, providing a thorough characterization of the stress propagation phenomena within the disk. Our findings contribute to the understanding of mechanical behavior in disk-shaped components and can be valuable in the design and optimization of such structures across various engineering disciplines.

preprint2022arXiv

Drag of an elliptic intruder in a two-dimensional granular environment

The drag of an elliptic intruder in a two-dimensional granular environment is numerically studied. The movement parallel to the major axis of the intruder is found to be unstable. The drag law is given by the sum of the yield force and the dynamic term, the latter of which is approximately reproduced by a simple collision model. The flow field around the intruder for sufficiently larger drag force is well fitted by the streamlines obtained from the perfect fluid. The stress fields around the intruder are also investigated when the movement of the intruder is balanced with interactions with the surrounding particles. The Airy stress function is found to well reproduce the stress fields once the stress on the surface of the intruder is given.

preprint2022arXiv

Drag of Two Cylindrical Intruders in a Two-Dimensional Granular Environment

The drag of two cylindrical intruders in a two-dimensional granular environment is numerically studied by the discrete element method. We find the yield force, below which the intruders cannot move because of interactions with the surrounding particles. Above the yield force, on the other hand, the intruders can move at a constant speed. We investigate the relationship between the drag force and the steady speed of the intruders, where the speed becomes higher as the distance between the intruders decreases. We confirm that the origin of the yield is the Coulombic friction between the particles and the bottom plate by changing the value of the friction coefficient. We also find that the yield force is almost proportional to the friction coefficient, which means that the number of particles determining the yield force is almost constant. On the other hand, the two-dimensional elasticity is applicable to determine the stress fields around the intruders. We confirm that fields asymmetric with respect to the drag direction are reproduced by using the information of the stresses on the surfaces of the intruders by introducing bipolar coordinates.

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

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

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.

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

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

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.

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.

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

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.