BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ryoichi Yamamoto

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

ResearcherUnclaimed source record
cond-mat.softphysics.flu-dyncond-mat.stat-mechcond-mat.dis-nncond-mat.mtrl-sciastro-ph.COBiological PhysicsCell BehaviorMachine Learningmath-phmath.MPmath.OCphysics.comp-phphysics.data-an

Catalog footprint

What is connected

31works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

31 published item(s)

preprint2022arXiv

Dynamics of microswimmers near a soft penetrable interface

Few simulations exist for microswimmers near deformable interfaces. Here, we present numerical simulations of the hydrodynamic flows associated with a single microswimmer embedded in a binary fluid mixture. The two fluids demix, separated by a penetrable and deformable interface that we assume to be initially prepared in its planar ground-state. We find that the microswimmer can either penetrate the interface, move parallel to it or bounce back off it. We analyze how the trajectory depends on the swimmer type (pusher/puller) and the angle of incidence with respect to the interface. Our simulations are performed in a system with periodic boundary conditions, corresponding to an infinite array of fluid interfaces. A puller reaches a steady state in which it either swims parallel to the interface or selects a perpendicular orientation, repeatedly penetrating through the interface. In contrast, a pusher follows a bouncing trajectory between two interfaces. We discuss several examples in biology in which swimmers penetrate soft interfaces. Our work can be seen as a highly simplified model of such processes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nash Neural Networks : Inferring Utilities from Optimal Behaviour

We propose Nash Neural Networks ($N^3$) as a new type of Physics Informed Neural Network that is able to infer the underlying utility from observations of how rational individuals behave in a differential game with a Nash equilibrium. We assume that the dynamics for both the population and the individual are known, but not the payoff function, which specifies the cost per unit time of being in any particular state. We construct our network in such a way that the Euler-Lagrange equations of the corresponding optimal control problem are satisfied and the optimal control is self-consistently determined. In this way, we are able to learn the unknown payoff function in an unsupervised manner. We have applied the $N^3$ to study the optimal behaviour during epidemics, in which individuals can choose to socially distance depending on the state of the pandemic and the cost of being infected. Training our network against synthetic data for a simple SIR model, we showed that it is possible to accurately reproduce the hidden payoff function, in such a way that the game dynamics are respected. Our approach will have far-reaching applications, as it allows one to infer utilities from behavioural data, and can thus be applied to study a wide array of problems in science, engineering, economics and government planning.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Relation between dynamic heterogeneities observed in scattering experiments and four-body correlations

Dynamic heterogeneity is expected to be a key concept for understanding the origin of slow dynamics near the glass transition. In previous studies, quantitative evaluations of dynamic heterogeneity have been attempted using two different routes, i.e., the speckle patterns in scattering experiments or the four-body correlation functions of microscopic configuration data obtained from molecular dynamics simulations or real-space observations using confocal microscopy. However, the physical relationship between these dynamic heterogeneities obtained using different methods has not been clarified. This study proposes a connection between dynamic heterogeneities characterized based on speckle patterns and those obtained from four-body correlations. The validity of the relationship is also clarified.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Active gas features in three HSC-SSP CAMIRA clusters revealed by high angular resolution analysis of MUSTANG-2 SZE and XXL X-ray observations

We present results from simultaneous modeling of high angular resolution GBT/MUSTANG-2 90 GHz Sunyaev-Zel'dovich effect (SZE) measurements and XMM-XXL X-ray images of three rich galaxy clusters selected from the HSC-SSP Survey. The combination of high angular resolution SZE and X-ray imaging enables a spatially resolved multi-component analysis, which is crucial to understand complex distributions of cluster gas properties. The targeted clusters have similar optical richnesses and redshifts, but exhibit different dynamical states in their member galaxy distributions: a single-peaked cluster, a double-peaked cluster, and a cluster belonging to a supercluster. A large-scale residual pattern in both regular Compton-parameter $y$ and X-ray surface brightness distributions is found in the single-peaked cluster, indicating a sloshing mode. The double-peaked cluster shows an X-ray remnant cool core between two SZE peaks associated with galaxy concentrations. The temperatures of the two peaks reach $\sim20-30$ keV in contrast to the cool core component of $\sim2$ keV, indicating a violent merger. The main SZE signal for the supercluster is elongated along a direction perpendicular to the major axis of the X-ray core, suggesting a minor merger before core passage. The $S_X$ and $y$ distributions are thus perturbed at some level, regardless of the optical properties. We find that the integrated Compton $y$ parameter and the temperature for the major merger are boosted from those expected by the weak-lensing mass and those for the other two clusters show no significant deviations, which is consistent with predictions of numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamics of a chiral swimmer sedimenting on a flat plate

Three-dimensional simulations with fully resolved hydrodynamics are performed to study the dynamics of a single squirmer under gravity, in order to clarify its motion in the vicinity of a flat plate. Different dynamics emerge for different gravity strengths. In a moderate gravity regime, neutral squirmers and pullers eventually stop moving and reorient in a direction perpendicular to the plate; pushers, instead, exhibit continuous motion in a tilted direction. In the strong gravity regime, all types of squirmers sediment and reorient perpendicularly to the plate. In this study, the chirality is introduced to model realistic micro-swimmers, and its crucial effects on the swimmer dynamics are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Collective motion of cells crawling on a substrate: roles of cell shape and contact inhibition

Contact inhibition plays a crucial role in the motility of cells, the process of wound healing, and the formation of tumors. By mimicking the mechanical motion of calls crawling on a substrate using a pseudopod, we constructed a minimal model for migrating cells which gives rise to contact inhibition of locomotion (CIL) naturally. The model cell consists of two disks, one in the front (a pseudopod) and the other one in the back (cell body), connected by a finitely extensible spring. Despite the simplicity of the model, the cells' collective behavior is highly nontrivial, depending on the shape of cells and whether CIL is enabled or not. Cells with a small front circle (i.e. a narrow pseudopod) form immobile colonies. In contrast, cells with a large front circle (i.e. such as a lamellipodium) exhibit coherent migration without any explicit alignment mechanism being present in the model. This suggests that crawling cells often exhibit broad fronts because it helps them avoid clustering. Upon increasing the density, the cells develop density waves which propagate against the direction of cell migration and finally arrest at higher densities.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Purely hydrodynamic origin for swarming of swimming particles

Three-dimensional simulations with fully resolved hydrodynamics are performed to study the collective motion of model swimmers in confinement. We show that certain swimming mechanisms can lead to traveling wave-like collective motion even without any direct alignment mechanism. It is also shown that by varying the swimming mechanism, this collective motion can be suppressed, contrary to the perception that hydrodynamic effects are completely screened at high volume fraction. From an analysis of bulk systems, it is shown that this traveling wave-like motion, which can be characterized as a pseudo-acoustic mode, is mainly due to the intrinsic swimming property of the particles.

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Direct numerical simulations of rigid body dispersions. I. Mobility/Friction tensors of assemblies of spheres

An improved formulation of the Smoothed Profile method is introduced to perform direct numerical simulations of arbitrary rigid body dispersions in a Newtonian host solvent. Previous implementations of the method were restricted to spherical particles, severely limiting the types of systems that could be studied. The validity of the method is carefully examined by computing the friction/mobility tensors for a wide variety of geometries and comparing them to reference values obtained from accurate solutions to the Stokes-Equation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hydrodynamic Interactions of Self-Propelled Swimmers

The hydrodynamic interactions of a suspension of self-propelled particles are studied using a direct numerical simulation method which simultaneously solves for the host fluid and the swimming particles. A modified version of the "Smoothed Profile" method (SPM) is developed to simulate microswimmers as squirmers, which are spherical particles with a specified surface-tangential slip velocity between the particles and the fluid. This simplified swimming model allows one to represent different types of propulsion (pullers and pushers) and is thus ideal to study the hydrodynamic interactions among swimmers. We use the SPM to study the diffusive behavior which arises due to the swimming motion of the particles, and show that there are two basic mechanisms responsible for this phenomena: the hydrodynamic interactions caused by the squirming motion of the particles, and the particle-particle collisions. This dual nature gives rise to two distinct time- and length- scales, and thus to two diffusion coefficients, which we obtain by a suitable analysis of the swimming motion. We show that the collisions between swimmers can be interpreted in terms of binary collisions, in which the effective collision radius is reduced due to the collision dynamics of swimming particles in viscous fluids. At short time-scales, the dynamics of the swimmer is analogous to that of an inert tracer particle in a swimming suspension, in which the diffusive motion is caused by fluid-particle collisions. Our results, along with the simulation method we have introduced, will allow us to gain a better understanding of the complex hydrodynamic interactions of self-propelled swimmers.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Propagation of hydrodynamic interactions between particles in a compressible fluid

Hydrodynamic interactions are transmitted by viscous diffusion and sound propagation: the temporal evolution of hydrodynamic interactions by both mechanisms is studied by direct numerical simulation in this paper. The hydrodynamic interactions for a system of two particles in a fluid are estimated by the velocity correlation of the particles. In an incompressible fluid, hydrodynamic interactions propagate instantaneously at the infinite speed of sound, followed by the temporal evolution of viscous diffusion. On the other hand, in a compressible fluid, sound propagates at a finite speed, which affects the temporal evolution of the hydrodynamic interactions by the order of magnitude relation between the time scales of viscous diffusion and sound propagation. The hydrodynamic interactions are characterized by introducing the ratio of these time scales as an interactive compressibility factor.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Velocity relaxation of a particle in a confined compressible fluid

The velocity relaxation of an impulsively forced spherical particle in a fluid confined by two parallel plane walls is studied using a direct numerical simulation approach. During the relaxation process, the momentum of the particle is transmitted in the ambient fluid by viscous diffusion and sound wave propagation, and the fluid flow accompanied by each mechanism has a different character and affects the particle motion differently. Because of the bounding walls, viscous diffusion is hampered, and the accompanying shear flow is gradually diminished. However, the sound wave is repeatedly reflected and spreads diffusely. As a result, the particle motion is governed by the sound wave and backtracks differently in a bulk fluid. The time when the backtracking of the particle occurs changes non-monotonically with respect to the compressibility factor and is minimized at the characteristic compressibility factor. This factor depends on the wall spacing, and the dependence is different at small and large wall spacing regions based on the different mechanisms causing the backtracking.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

2D Lattice Liquid Models

A family of novel models of liquid on a 2D lattice (2D lattice liquid models) have been proposed as primitive models of soft-material membrane. As a first step, we have formulated them as single-component, single-layered, classical particle systems on a two-dimensional surface with no explicit viscosity. Among the family of the models, we have shown and constructed two stochastic models, a vicious walk model and a flow model, on an isotropic regular lattice and on the rectangular honeycomb lattice of various sizes. In both cases, the dynamics is governed by the nature of the frustration of the particle movements. By simulations, we have found the approximate functional form of the frustration probability, and peculiar anomalous diffusions in their time-averaged mean square displacements in the flow model. The relations to other existing statistical models and possible extensions of the models are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Acoustic wave propagation through a supercooled liquid: A normal mode analysis

The mechanism of acoustic wave propagation in supercooled liquids is not yet fully understood since the vibrational dynamics of supercooled liquids are strongly affected by their amorphous inherent structures. In this paper, the acoustic wave propagation in a supercooled model liquid is studied by using normal mode analysis. Due to the highly disordered inherent structure, a single acoustic wave is decomposed into many normal modes in broad frequency range. This causes the rapid decay of the acoustic wave and results in anomalous wavenumber dependency of the dispersion relation and the rate of attenuation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Direct numerical simulation of dispersed particles in a compressible fluid

We present a direct numerical simulation method for investigating the dynamics of dispersed particles in a compressible solvent fluid. The validity of the simulation is examined by calculating the velocity relaxation of an impulsively forced spherical particle with a known analytical solution. The simulation also gives information about the fluid motion, which provides some insight into the particle motion. Fluctuations are also introduced by random stress, and the validity of this case is examined by comparing the calculation results with the fluctuation-dissipation theorem.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mechanical Responses and Stress Fluctuations of a Supercooled Liquid in a Sheared Non-Equilibrium State

A steady shear flow can drive supercooled liquids into a non-equilibrium state. Using molecular dynamics simulations under steady shear flow superimposed with oscillatory shear strain for a probe, non-equilibrium mechanical responses are studied for a model supercooled liquid composed of binary soft spheres. We found that even in the strongly sheared situation, the supercooled liquid exhibits surprisingly isotropic responses to oscillating shear strains applied in three different components of the strain tensor. Based on this isotropic feature, we successfully constructed a simple two-mode Maxwell model that can capture the key features of the storage and loss moduli, even for highly non-equilibrium state. Furthermore, we examined the correlation functions of the shear stress fluctuations, which also exhibit isotropic relaxation behaviors in the sheared non-equilibrium situation. In contrast to the isotropic features, the supercooled liquid additionally demonstrates anisotropies in both its responses and its correlations to the shear stress fluctuations. Using the constitutive equation (a two-mode Maxwell model), we demonstrated that the anisotropic responses are caused by the coupling between the oscillating strain and the driving shear flow. We measured the magnitude of this violation in terms of the effective temperature. It was demonstrated that the effective temperature is notably different between different components, which indicates that a simple scalar mapping, such as the concept of an effective temperature, oversimplifies the true nature of supercooled liquids under shear flow. An understanding of the mechanism of isotropies and anisotropies in the responses and fluctuations will lead to a better appreciation of these violations of the FDT, as well as certain consequent modifications to the concept of an effective temperature.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamical heterogeneity in a highly supercooled liquid: Consistent calculations of correlation length, intensity, and lifetime

We have investigated dynamical heterogeneity in a highly supercooled liquid using molecular-dynamics simulations in three dimensions. Dynamical heterogeneity can be characterized by three quantities: correlation length $ξ_4$, intensity $χ_4$, and lifetime $τ_{\text{hetero}}$. We evaluated all three quantities consistently from a single order parameter. In a previous study (H. Mizuno and R. Yamamoto, Phys. Rev. E {\bf 82}, 030501(R) (2010)), we examined the lifetime $τ_{\text{hetero}}(t)$ in two time intervals $t=τ_α$ and $τ_{\text{ngp}}$, where $τ_α$ is the $α$-relaxation time and $τ_{\text{ngp}}$ is the time at which the non-Gaussian parameter of the Van Hove self-correlation function is maximized. In the present study, in addition to the lifetime $τ_{\text{hetero}}(t)$, we evaluated the correlation length $ξ_4(t)$ and the intensity $χ_4(t)$ from the same order parameter used for the lifetime $τ_{\text{hetero}}(t)$. We found that as the temperature decreases, the lifetime $τ_{\text{hetero}}(t)$ grows dramatically, whereas the correlation length $ξ_4(t)$ and the intensity $χ_4(t)$ increase slowly compared to $τ_{\text{hetero}}(t)$ or plateaus. Furthermore, we investigated the lifetime $τ_{\text{hetero}}(t)$ in more detail. We examined the time-interval dependence of the lifetime $τ_{\text{hetero}}(t)$ and found that as the time interval $t$ increases, $τ_{\text{hetero}}(t)$ monotonically becomes longer and plateaus at the relaxation time of the two-point density correlation function. At the large time intervals for which $τ_{\text{hetero}}(t)$ plateaus, the heterogeneous dynamics migrate in space with a diffusion mechanism, such as the particle density.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Reentrant transition in the shear viscosity of dilute rigid rod dispersions

The intrinsic viscosity of a dilute dispersion of rigid rods is studied using a recently developed direct numerical simulation (DNS) method for particle dispersions. A reentrant transition from shear-thinning to the 2nd Newtonian regime is successfully reproduced in the present DNS results around a Peclet number ${\rm Pe}=150$, which is in good agreement with our theoretical prediction of ${\rm Pe}=143$, at which the dynamical crossover from Brownian to non-Brownian behavior takes place in the rotational motion of the rotating rod. The viscosity undershoot is observed in our simulations before reaching the 2nd Newtonian regime. The physical mechanisms behind these behaviors are analyzed in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Implementation of Lees-Edwards periodic boundary conditions for direct numerical simulations of particle dispersions under shear flow

A general methodology is presented to perform direct numerical simulations of particle dispersions in a shear flow with Lees-Edwards periodic boundary conditions. The Navier-Stokes equation is solved in oblique coordinates to resolve the incompatibility of the fluid motions with the sheared geometry, and the force coupling between colloidal particles and the host fluid is imposed by using a smoothed profile method. The validity of the method is carefully examined by comparing the present numerical results with experimental viscosity data for particle dispersions in a wide range of volume fractions and shear rates including nonlinear shear-thinning regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Lifetime of dynamical heterogeneity in a highly supercooled liquid

We numerically examine dynamical heterogeneity in a highly supercooled three-dimensional liquid via molecular-dynamics simulations. To define the local dynamics, we consider two time intervals, $τ_α$ and $τ_{\text{ngp}}$. $τ_α$ is the $α$ relaxation time, and $τ_{\text{ngp}}$ is the time at which non-Gaussian parameter of the van Hove self-correlation function is maximized. We determine the lifetimes of the heterogeneous dynamics in these two different time intervals, $τ_{\text{hetero}}(τ_α)$ and $τ_{\text{hetero}}(τ_{\text{ngp}})$, by calculating the time correlation function of the particle dynamics, i.e., the four-point correlation function. We find that the difference between $τ_{\text{hetero}}(τ_α)$ and $τ_{\text{hetero}}(τ_{\text{ngp}})$ increases with decreasing temperature. At low temperatures, $τ_{\text{hetero}}(τ_α)$ is considerably larger than $τ_α$, while $τ_{\text{hetero}}(τ_{\text{ngp}})$ remains comparable to $τ_α$. Thus, the lifetime of the heterogeneous dynamics depends strongly on the time interval.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Tumbling motion of a single chain in shear flow: a crossover from Brownian to non-Brownian behavior

We present numerical results for the dynamics of a single chain in steady shear flow. The chain is represented by a bead-spring model, and the smoothed profile method is used to accurately account for the effects of thermal fluctuations and hydrodynamic interactions acting on beads due to host fluids. It is observed that the chain undergoes tumbling motions and that its dimensionless frequency F depends only on the Peclet number Pe with a power law. The exponent of Pe clearly changes from 2/3 to 1 around the critical Peclet number, indicating that the crossover reflects the competition of thermal fluctuation and shear flow. The presented numerical results agree well with our theoretical analysis based on Jeffery's work.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Direct numerical simulations for non-Newtonian rheology of concentrated particle dispersions

The non-Newtonian behavior of a monodisperse concentrated dispersion of spherical particles was investigated using a direct numerical simulation method, that takes into account hydrodynamic interactions and thermal fluctuations accurately. Simulations were performed under steady shear flow with periodic boundary conditions in the three directions. The apparent shear viscosity of the dispersions was calculated at volume fractions ranging from 0.31 to 0.56. Shear-thinning behavior was clearly observed at high volume fractions. The low- and high-limiting viscosities were then estimated from the apparent viscosity by fitting these data into a semi-empirical formula. Furthermore, the short-time motions were examined for Brownian particles fluctuating in concentrated dispersions, for which the fluid inertia plays an important role. The mean square displacement was monitored in the vorticity direction at several different Peclet numbers and volume fractions so that the particle diffusion coefficient is determined from the long-time behavior of the mean square displacement. Finally, the relationship between the non-Newtonian viscosity of the dispersions and the structural relaxation of the dispersed Brownian particles is examined.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
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$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Simulating (electro)hydrodynamic effects in colloidal dispersions: smoothed profile method

Previously, we have proposed a direct simulation scheme for colloidal dispersions in a Newtonian solvent [Phys.Rev.E 71,036707 (2005)]. An improved formulation called the ``Smoothed Profile (SP) method'' is presented here in which simultaneous time-marching is used for the host fluid and colloids. The SP method is a direct numerical simulation of particulate flows and provides a coupling scheme between the continuum fluid dynamics and rigid-body dynamics through utilization of a smoothed profile for the colloidal particles. Moreover, the improved formulation includes an extension to incorporate multi-component fluids, allowing systems such as charged colloids in electrolyte solutions to be studied. The dynamics of the colloidal dispersions are solved with the same computational cost as required for solving non-particulate flows. Numerical results which assess the hydrodynamic interactions of colloidal dispersions are presented to validate the SP method. The SP method is not restricted to particular constitutive models of the host fluids and can hence be applied to colloidal dispersions in complex fluids.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Velocity autocorrelation function of fluctuating particles in incompressible fluids. Toward direct numerical simulation of particle dispersions

Motions of fluctuating Brownian particles in an incompressible viscous fluid have been studied by coupled simulations of Brownian particles and host fluid. We calculated the velocity autocorrelation functions of Brownian particles and compared them with the theoretical results. Extensive discussions have been made on the time scales for which our numerical model is valid.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

A Smooth Interface Method for Simulating Liquid Crystal Colloid Dispersions

A new method is presented for mesoscopic simulations of particle dispersions in liquid crystal solvents. It allows efficient first-principle simulations of the dispersions involving many particles with many-body interactions mediated by the solvents. Demonstrations have been performed for the aggregation of colloid dispersions in two-dimensional nematic and smectic-C* solvents neglecting hydrodynamic effects, which will be taken into account in the near future.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Kinetic Heterogeneities in a Highly Supercooled Liquid

We study a highly supercooled two-dimensional fluid mixture via molecular dynamics simulation. We follow bond breakage events among particle pairs, which occur on the scale of the $α$ relaxation time $τ_α$. Large scale heterogeneities analogous to the critical fluctuations in Ising systems are found in the spatial distribution of bonds which are broken in a time interval with a width of order $0.05τ_α$. The structure factor of the broken bond density is well approximated by the Ornstein-Zernike form. The correlation length is of order $100 σ_1$ at the lowest temperature studied, $σ_1$ being the particle size. The weakly bonded regions thus identified evolve in time with strong spatial correlations.

0 citations0 reviewsNo code linkedOpen access