Researcher profile

Satoshi Takada

Satoshi Takada contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

8 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.

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.