Source author record

Shugo Yasuda

Shugo Yasuda 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

9works
7topics
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

9 published item(s)

preprint2026arXiv

Multiscale analysis of a kinetic equation for mechanotaxis

We present a new kinetic equation for cell migration driven by mechanical interactions with the substrate, an effect not previously captured in kinetic models, and essential for explaining observed collective behaviors such as those in bacterial colonies. The model introduces an acceleration term that accounts for the dynamics of motile cells undergoing mechanotaxis, where extracellular signals modulate the forces arising from cell-substrate interactions. From this formulation, we derive a family of macroscopic limit equations and analyze their principal properties. In particular, we examine linear stability and pattern formation ability through theoretical analysis, supported by numerical simulations.

preprint2022arXiv

Numerical study of the volcano e ect in chemotactic aggregation based on a kinetic transport equation with non-instantaneous tumbling

Aggregation of chemotactic bacteria under a unimodal distribution of chemical cues was investigated by Monte Carlo (MC) simulation based on a kinetic transport equation, which considers an internal adaptation dynamics as well as a finite tumbling duration. It was found that there exist two different regimes of the adaptation time, between which the effect of the adaptation time on the aggregation behavior is reversed; that is, when the adaptation time is as small as the running duration, the aggregation becomes increasingly steeper as the adaptation time increases, while, when the adaptation time is as large as the diffusion time of the population density, the aggregation becomes more diffusive as the adaptation time increases. Moreover, the aggregation profile becomes bimodal (volcano) at the large adaptation-time regime when the tumbling duration is sufficiently large while it is always unimodal at the small adaptation-time regime. A remarkable result of this study is the identification of the parameter regime and scaling for the volcano effect. That is, by comparing the results of MC simulations to the continuum-limit models obtained at each of the small and large adaptation-time scalings, it is clarified that the volcano effect arises due to the coupling of diffusion, adaptation, and finite tumbling duration, which occurs at the large adaptation-time scaling.

preprint2020arXiv

Numerical scheme for kinetic transport equation with internal state *

We investigate the numerical discretization of a two-stream kinetic system with an internal state, such system has been introduced to model the motion of cells by chemotaxis. This internal state models the intracellular methylation level. It adds a variable in the mathematical model, which makes it more challenging to simulate numerically. Moreover, it has been shown that the macroscopic or mesoscopic quantities computed from this system converge to the Keller-Segel system at diffusive scaling or to the velocity-jump kinetic system for chemotaxis at hyperbolic scaling. Then we pay attention to propose numerical schemes uniformly accurate with respect to the scaling parameter. We show that these schemes converge to some limiting schemes which are consistent with the limiting macroscopic or kinetic system. This study is illustrated with some numerical simulations and comparisons with Monte Carlo simulations.

preprint2015arXiv

Synchronized molecular-dynamics simulation for the thermal lubrication of a polymeric liquid between parallel plates

The Synchronized Molecular-Dynamics simulation which was recently proposed by authors [Phys. Rev. X {\bf 4}, 041011 (2014)] is applied to the analysis of polymer lubrication between parallel plates. The rheological properties, conformational change of polymer chains, and temperature rise due to the viscous heating are investigated with changing the values of thermal conductivity of the polymeric liquid. It is found that at a small applied shear stress on the plate, the temperature of polymeric liquid only slightly increases in inverse proportion to the thermal conductivity and the apparent viscosity of polymeric liquid is not much affected by changing the thermal conductivity. However, at a large shear stress, the transitional behaviors of the polymeric liquid occur due to the interplay of the shear deformation and viscous heating by changing the thermal conductivity. This transition is characterized by the Nahme-Griffith number $Na$ which is defined as the ratio of the viscous heating to the thermal conduction at a characteristic temperature. When the Nahme-Griffith number exceeds the unity, the temperature of polymeric liquid increases rapidly and the apparent viscosity also exponentially decreases as the thermal conductivity decreases. The conformation of polymer chains is stretched and aligned by the shear flow for $Na<1$, but the coherent structure becomes disturbed by the thermal motion of molecules for $Na>1$.

preprint2011arXiv

Dynamic rheology of a supercooled polymer melt in non-uniform oscillating flows in rapidly oscillating plates

The dynamic rheology of a polymer melt composed of short chains with ten beads between rapidly oscillating plates is investigated for various oscillation frequencies by using the hybrid simulation of the molecular dynamics and computational fluid dynamics. In the rapidly oscillating plates, non-uniform boundary layer flows are generated over the plate due to inertia of the fluid, and the local rheological properties of the melt are spatially varied according to the local flow fields. Far from the plate, the local strain decreases as the oscillation frequency increases such that the dynamic rheology deviates from the linear moduli at the low oscillation frequencies rather than high oscillation frequencies. On the contrary, near the plate, the local strain rate increases with the oscillation frequency such that the shear thinning is enhanced at high oscillation frequencies. In close vicinity to the plate, the dynamic viscosity is mostly independent of the oscillation frequency, and the shear thinning behavior becomes similar to that observed in steady shear flows. Non-linear rheological properties are also investigated by the spectrum analysis and the Lissajous-Bowditch curve. It is found that the fractional amplitude of the higher harmonics to the linear harmonics is suppressed within the boundary layer due to the non-slip boundary on the oscillating plate. We also find that the melt exhibits inter-cycle shear thinning between different positions but exhibits intra-cycle shear thickening at a fixed position in the vicinity of the plate.

preprint2011arXiv

Multiscale Simulations for Polymeric Flow

Multiscale simulation methods have been developed based on the local stress sampling strategy and applied to three flow problems with different difficulty levels: (a) general flow problems of simple fluids, (b) parallel (one-dimensional) flow problems of polymeric liquids, and (c) general (two- or three-dimensional) flow problems of polymeric liquids. In our multiscale methods, the local stress of each fluid element is calculated directly by performing microscopic or mesoscopic simulations according to the local flow quantities instead of using any constitutive relations. For simple fluids (a), such as the Lenard-Jones liquid, a multiscale method combining MD and CFD simulations is developed based on the local equilibrium assumption without memories of the flow history. (b), the multiscale method is extended to take into account the memory effects that arise in hydrodynamic stress due to the slow relaxation of polymer-chain conformations. The memory of polymer dynamics on each fluid element is thus resolved by performing MD simulations in which cells are fixed at the mesh nodes of the CFD simulations.For general (two- or three-dimensional) flow problems of polymeric liquids (c), it is necessary to trace the history of microscopic information such as polymer-chain conformation, which carries the memories of past flow history, along the streamline of each fluid element. A Lagrangian-based CFD is thus implemented to correctly advect the polymer-chain conformation consistently with the flow. On each fluid element, coarse-grained polymer simulations are carried out to consider the dynamics of entangled polymer chains that show extremely slow relaxation compared to microscopic time scales.

preprint2009arXiv

Multiscale modeling and simulation for polymer melt flows between parallel plates

The flow behaviors of polymer melt composed of short chains with ten beads between parallel plates are simulated by using a hybrid method of molecular dynamics and computational fluid dynamics. Three problems are solved: creep motion under a constant shear stress and its recovery motion after removing the stress, pressure-driven flows, and the flows in rapidly oscillating plates. In the creep/recovery problem, the delayed elastic deformation in the creep motion and evident elastic behavior in the recovery motion are demonstrated. The velocity profiles of the melt in pressure-driven flows are quite different from those of Newtonian fluid due to shear thinning. Velocity gradients of the melt become steeper near the plates and flatter at the middle between the plates as the pressure gradient increases and the temperature decreases. In the rapidly oscillating plates, the viscous boundary layer of the melt is much thinner than that of Newtonian fluid due to the shear thinning of the melt. Three different rheological regimes, i.e., the viscous fluid, visco-elastic liquid, and visco-elastic solid regimes, form over the oscillating plate according to the local Deborah numbers. The melt behaves as a viscous fluid in a region for $ωτ^R\lesssim 1$, and the crossover between the liquid-like and solid-like regime takes place around $ωτ^α\simeq 1$ (where $ω$ is the angular frequency of the plate and $τ^R$ and $τ^α$ are Rouse and $α$ relaxation time, respectively).

preprint2009arXiv

Rheology of a supercooled polymer melt near an oscillating plate: an application of multiscale modeling

The behavior of a supercooled polymer melt composed of short chains with ten beads near an oscillating plate are simulated by using a hybrid simulation of molecular dynamics (MD) and computational fluid dynamics (CFD). In the method, the macroscopic dynamics are solved by using CFD, but, instead of using any constitutive equations, a local stress is calculated by using a non-equilibrium MD simulation associated at each lattice node in the CFD calculation. It is seen that the local rheology of the melt varies considerably in a thin viscous diffusion layer that arises near an oscillating plate. It is also found that the local rheology of the melt is divided into the three different regimes, i.e., the viscous fluid, viscoelastic liquid, and viscoelastic solid regimes, according to the local Deborah number $De$, which is defined with the Rouse or $α$ relaxation time, $τ_R$ or $τ_α$, and the angular frequency of the plate $ω$ as $De^R$=$ωτ_R$ or $De^α$=$ωτ_α$. The melt behaves as a viscous fluid when $De^R\lesssim 1$, and the crossover between the liquid-like and solid-like regime takes place around $De^α\simeq 1$.

preprint2008arXiv

A Model for Hybrid Simulations of Molecular Dynamics and CFD

We propose a method for multi-scale hybrid simulations of molecular dynamics (MD) and computational fluid dynamics (CFD). In the method, usual lattice-mesh based simulations are applied for CFD level, but each lattice is associated with a small MD cell which generates a "local stress" according to a "local flow field" given from CFD instead of using any constitutive functions at CFD level. We carried out the hybrid simulations for some elemental flow problems of simple Lennard-Jones liquids and compared the results with those obtained by usual CFDs with a Newtonian constitutive relation in order to examine the validity of our hybrid simulation method. It is demonstrated that our hybrid simulations successfully reproduced the correct flow behavior obtained from usual CFDs as far as the mesh size $Δx$ and the time-step $Δt$ of CFD are not too large comparing to the system size $l_{\rm MD}$ and the sampling duration $t_{\rm MD}$ of MD simulations performed at each time step of CFDs. Otherwise, simulations are affected by large fluctuations due to poor statistical averages taken in the MD part. Properties of the fluctuations are analyzed in detail.