BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

A. Donev

A. Donev contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NAcond-mat.mes-hallcond-mat.stat-mechphysics.flu-dyncond-mat.softcond-mat.mtrl-scimath-phmath.MPphysics.chem-phphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
10topics
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

20 published item(s)

preprint2015arXiv

An Immersed Boundary Method for Rigid Bodies

We develop an immersed boundary (IB) method for modeling flows around fixed or moving rigid bodies that is suitable for a broad range of Reynolds numbers, including steady Stokes flow. The spatio-temporal discretization of the fluid equations is based on a standard staggered-grid approach. Fluid-body interaction is handled using Peskin's IB method; however, unlike existing IB approaches to such problems, we do not rely on penalty or fractional-step formulations. Instead, we use an unsplit scheme that ensures the no-slip constraint is enforced exactly in terms of the Lagrangian velocity field evaluated at the IB markers. Fractional-step approaches, by contrast, can impose such constraints only approximately. Imposing these constraints exactly requires the solution of a large linear system that includes the fluid velocity and pressure as well as Lagrange multiplier forces that impose the motion of the body. To solve this system efficiently, we develop a preconditioner for the constrained IB formulation that is based on an analytical approximation to the Schur complement. This approach is enabled by the near translational and rotational invariance of Peskin's IB method. We demonstrate that only a few cycles of a geometric multigrid method for the fluid equations are required in each application of the preconditioner, and we demonstrate robust convergence of the overall Krylov solver despite the approximations made in the preconditioner. We apply the method to a number of test problems at zero and finite Reynolds numbers, and we demonstrate first-order convergence of the method to several analytical solutions and benchmark computations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Brownian Dynamics without Green's Functions

We develop a Fluctuating Immersed Boundary (FIB) method for performing Brownian dynamics simulations of confined particle suspensions. Unlike traditional methods which employ analytical Green's functions for Stokes flow in the confined geometry, the FIB method uses a fluctuating finite-volume Stokes solver to generate the action of the response functions "on the fly". Importantly, we demonstrate that both the deterministic terms necessary to capture the hydrodynamic interactions among the suspended particles, as well as the stochastic terms necessary to generate the hydrodynamically-correlated Brownian motion, can be generated by solving the steady Stokes equations numerically only once per time step. This is accomplished by including a stochastic contribution to the stress tensor in the fluid equations consistent with fluctuating hydrodynamics. We develop novel temporal integrators that account for the multiplicative nature of the noise in the equations of Brownian dynamics and the strong dependence of the mobility on the configuration for confined systems. Notably, we propose a random finite difference approach to approximating the stochastic drift proportional to the divergence of the configuration-dependent mobility matrix. Through comparisons with analytical and existing computational results, we numerically demonstrate the ability of the FIB method to accurately capture both the static (equilibrium) and dynamic properties of interacting particles in flow.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fluctuating hydrodynamics of multispecies reactive mixtures

We formulate and study computationally the fluctuating compressible Navier-Stokes equations for reactive multi-species fluid mixtures. We contrast two different expressions for the covariance of the stochastic chemical production rate in the Langevin formulation of stochastic chemistry, and compare both of them to predictions of the chemical Master Equation for homogeneous well-mixed systems close to and far from thermodynamic equilibrium. We develop a numerical scheme for inhomogeneous reactive flows, based on our previous methods for non-reactive mixtures [K. Balakrishnan, A. L. Garcia, A. Donev and J. B. Bell, Phys. Rev. E 89:013017, 2014]. We study the suppression of non-equilibrium long-ranged correlations of concentration fluctuations by chemical reactions, as well as the enhancement of pattern formation by spontaneous fluctuations. Good agreement with available theory demonstrates that the formulation is robust and a useful tool in the study of fluctuations in reactive multi-species fluids. At the same time, several problems with Langevin formulations of stochastic chemistry are identified, suggesting that future work should examine combining Langevin and Master Equation descriptions of hydrodynamic and chemical fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inertial Coupling Method for particles in an incompressible fluctuating fluid

We develop an inertial coupling method for modeling the dynamics of point-like 'blob' particles immersed in an incompressible fluid, generalizing previous work for compressible fluids. The coupling consistently includes excess (positive or negative) inertia of the particles relative to the displaced fluid, and accounts for thermal fluctuations in the fluid momentum equation. The coupling between the fluid and the blob is based on a no-slip constraint equating the particle velocity with the local average of the fluid velocity, and conserves momentum and energy. We demonstrate that the formulation obeys a fluctuation-dissipation balance, owing to the non-dissipative nature of the no-slip coupling. We develop a spatio-temporal discretization that preserves, as best as possible, these properties of the continuum formulation. In the spatial discretization, the local averaging and spreading operations are accomplished using compact kernels commonly used in immersed boundary methods. We find that the special properties of these kernels make the discrete blob a particle with surprisingly physically-consistent volume, mass, and hydrodynamic properties. We develop a second-order semi-implicit temporal integrator that maintains discrete fluctuation-dissipation balance, and is not limited in stability by viscosity. Furthermore, the temporal scheme requires only constant-coefficient Poisson and Helmholtz linear solvers, enabling a very efficient and simple FFT-based implementation on GPUs. We numerically investigate the performance of the method on several standard test problems...

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Low Mach Number Fluctuating Hydrodynamics of Binary Liquid Mixtures

Continuing on our previous work [ArXiv:1212.2644], we develop semi-implicit numerical methods for solving low Mach number fluctuating hydrodynamic equations appropriate for modeling diffusive mixing in isothermal mixtures of fluids with different densities and transport coefficients. We treat viscous dissipation implicitly using a recently-developed variable-coefficient Stokes solver [ArXiv:1308.4605]. This allows us to increase the time step size significantly compared to the earlier explicit temporal integrator. For viscous-dominated flows, such as flows at small scales, we develop a scheme for integrating the overdamped limit of the low Mach equations, in which inertia vanishes and the fluid motion can be described by a steady Stokes equation. We also describe how to incorporate advanced higher-order Godunov advection schemes in the numerical method, allowing for the treatment of fluids with high Schmidt number including the vanishing mass diffusion coefficient limit. We incorporate thermal fluctuations in the description in both the inertial and overdamped regimes. We apply our algorithms to model the development of giant concentration fluctuations during the diffusive mixing of water and glycerol, and compare numerical results with experimental measurements. We find good agreement between the two, and observe propagative (non-diffusive) modes at small wavenumbers (large spatial scales), not reported in published experimental measurements of concentration fluctuations in fluid mixtures. Our work forms the foundation for developing low Mach number fluctuating hydrodynamics methods for miscible multi-species mixtures of chemically reacting fluids.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Low Mach Number Fluctuating Hydrodynamics of Multispecies Liquid Mixtures

We develop a low Mach number formulation of the hydrodynamic equations describing transport of mass and momentum in a multispecies mixture of incompressible miscible liquids at specified temperature and pressure that generalizes our prior work on ideal mixtures of ideal gases and binary liquid mixtures. In this formulation we combine and extend a number of existing descriptions of multispecies transport available in the literature. The formulation applies to non-ideal mixtures of arbitrary number of species, without the need to single out a 'solvent' species, and includes contributions to the diffusive mass flux due to gradients of composition, temperature and pressure. Momentum transport and advective mass transport are handled using a low Mach number approach that eliminates fast sound waves (pressure fluctuations) from the full compressible system of equations and leads to a quasi-incompressible formulation. Thermal fluctuations are included in our fluctuating hydrodynamics description following the principles of nonequilibrium thermodynamics. We extend the semi-implicit staggered-grid finite-volume numerical method developed in our prior work on binary liquid mixtures, and use it to study the development of giant nonequilibrium concentration fluctuations in a ternary mixture subjected to a steady concentration gradient. We also numerically study the development of diffusion-driven gravitational instabilities in a ternary mixture, and compare our numerical results to recent experimental measurements in a Hele-Shaw cell. We find that giant nonequilibrium fluctuations can trigger the instability but are eventually dominated by the deterministic growth of the unstable mode, in both quasi two-dimensional (Hele-Shaw), and fully three-dimensional geometries used in typical shadowgraph experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiscale temporal integrators for fluctuating hydrodynamics

Following on our previous work [S. Delong and B. E. Griffith and E. Vanden-Eijnden and A. Donev, Phys. Rev. E, 87(3):033302, 2013], we develop temporal integrators for solving Langevin stochastic differential equations that arise in fluctuating hydrodynamics. Our simple predictor-corrector schemes add fluctuations to standard second-order deterministic solvers in a way that maintains second-order weak accuracy for linearized fluctuating hydrodynamics. We construct a general class of schemes and recommend two specific schemes: an explicit midpoint method, and an implicit trapezoidal method. We also construct predictor-corrector methods for integrating the overdamped limit of systems of equations with a fast and slow variable in the limit of infinite separation of the fast and slow timescales. We propose using random finite differences to approximate some of the stochastic drift terms that arise because of the kinetic multiplicative noise in the limiting dynamics. We illustrate our integrators on two applications involving the development of giant nonequilibrium concentration fluctuations in diffusively-mixing fluids. We first study the development of giant fluctuations in recent experiments performed in microgravity using an overdamped integrator. We then include the effects of gravity, and find that we also need to include the effects of fluid inertia, which affects the dynamics of the concentration fluctuations greatly at small wavenumbers.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A reversible mesoscopic model of diffusion in liquids: from giant fluctuations to Fick's law

We study diffusive mixing in the presence of thermal fluctuations under the assumption of large Schmidt number. In this regime we obtain a limiting equation that contains a diffusive thermal drift term with diffusion coefficient obeying a Stokes-Einstein relation, in addition to the expected advection by a random velocity. The overdamped limit correctly reproduces both the enhanced diffusion in the ensemble-averaged mean and the long-range correlated giant fluctuations in individual realizations of the mixing process, and is amenable to efficient numerical solution. Through a combination of Eulerian and Lagrangian numerical methods we demonstrate that diffusion in liquids is not most fundamentally described by Fick's irreversible law; rather, diffusion is better modeled as reversible random advection by thermal velocity fluctuations. We find that the diffusion coefficient is effectively renormalized to a value that depends on the scale of observation. Our work reveals somewhat unexpected connections between flows at small scales, dominated by thermal fluctuations, and flows at large scales, dominated by turbulent fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamic Density Functional Theory with hydrodynamic interactions and fluctuations

We derive a closed equation for the empirical concentration of colloidal particles in the presence of both hydrodynamic and direct interactions. The ensemble average of our functional Langevin equation reproduces known deterministic Dynamic Density Functional Theory (DDFT) [Rex and Löwen, Phys. Rev. Lett., 101(14):148302, 2008], and, at the same time, it also describes the microscopic fluctuations around the mean behavior. We suggest separating the ideal (non-interacting) contribution from additional corrections due to pairwise interactions. We find that, for an incompressible fluid and in the absence of direct interactions, the mean concentration follows Fick's law just as for uncorrelated walkers. At the same time, the nature of the stochastic terms in fluctuating DDFT is shown to be distinctly different for hydrodynamically-correlated and uncorrelated walkers. This leads to striking differences in the behavior of the fluctuations around Fick's law, even in the absence of pairwise interactions. We connect our own prior work [A. Donev, T. G. Fai, E. Vanden-Eijnden, J. Stat. Mech., P04004, 2014] on fluctuating hydrodynamics of diffusion in liquids to the DDFT literature, and demonstrate that the fluid cannot easily be eliminated from consideration if one wants to describe the collective diffusion in colloidal suspensions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Efficient Variable-Coefficient Finite-Volume Stokes Solvers

We investigate several robust preconditioners for solving the saddle-point linear systems that arise from spatial discretization of unsteady and steady variable-coefficient Stokes equations on a uniform staggered grid. Building on the success of using the classical projection method as a preconditioner for the coupled velocity-pressure system [B. E. Griffith, J. Comp. Phys., 228 (2009), pp. 75657595], as well as established techniques for steady and unsteady Stokes flow in the finite-element literature, we construct preconditioners that employ independent generalized Helmholtz and Poisson solvers for the velocity and pressure subproblems. We demonstrate that only a single cycle of a standard geometric multigrid algorithm serves as an effective inexact solver for each of these subproblems. Contrary to traditional wisdom, we find that the Stokes problem can be solved nearly as efficiently as the independent pressure and velocity subproblems, making the overall cost of solving the Stokes system comparable to the cost of classical projection or fractional step methods for incompressible flow, even for steady flow and in the presence of large density and viscosity contrasts. Two of the five preconditioners considered here are found to be robust to GMRES restarts and to increasing problem size, making them suitable for large-scale problems. Our work opens many possibilities for constructing novel unsplit temporal integrators for finite-volume spatial discretizations of the equations of low Mach and incompressible flow dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Low Mach Number Fluctuating Hydrodynamics of Diffusively Mixing Fluids

We formulate low Mach number fluctuating hydrodynamic equations appropriate for modeling diffusive mixing in isothermal mixtures of fluids with different density and transport coefficients. These equations eliminate the fluctuations in pressure associated with the propagation of sound waves by replacing the equation of state with a local thermodynamic constraint. We demonstrate that the low Mach number model preserves the spatio-temporal spectrum of the slower diffusive fluctuations. We develop a strictly conservative finite-volume spatial discretization of the low Mach number fluctuating equations in both two and three dimensions and construct several explicit Runge-Kutta temporal integrators that strictly maintain the equation of state constraint. The resulting spatio-temporal discretization is second-order accurate deterministically and maintains fluctuation-dissipation balance in the linearized stochastic equations. We apply our algorithms to model the development of giant concentration fluctuations in the presence of concentration gradients, and investigate the validity of common simplifications such as neglecting the spatial non-homogeneity of density and transport properties. We perform simulations of diffusive mixing of two fluids of different densities in two dimensions and compare the results of low Mach number continuum simulations to hard-disk molecular dynamics simulations. Excellent agreement is observed between the particle and continuum simulations of giant fluctuations during time-dependent diffusive mixing.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Reversible Diffusion by Thermal Fluctuations

A model for diffusion in liquids that couples the dynamics of tracer particles to a fluctuating Stokes equation for the fluid is investigated in the limit of large Schmidt number. In this limit, the concentration of tracers is shown to satisfy a closed-form stochastic advection-diffusion equation that is used to investigate the collective diffusion of hydrodynamically-correlated tracers through a combination of Eulerian and Lagrangian numerical methods. This analysis indicates that transport in liquids is quite distinct from the traditional Fickian picture of diffusion. While the ensemble-averaged concentration follows Fick's law with a diffusion coefficient that obeys the Stokes-Einstein relation, each instance of the diffusive mixing process exhibits long-ranged giant fluctuations around its average behavior. We construct a class of mesoscopic models for diffusion in liquids at different observation scales in which the renormalized diffusion coefficient depends on this scale. This indicates that the Fickian diffusion coefficient in liquids is not a material constant, but rather, changes with the scale at which experimental measurements are performed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Minimally-Resolved Immersed Boundary Model for Reaction-Diffusion Problems

We develop an immersed-boundary approach to modeling reaction-diffusion processes in dispersions of reactive spherical particles, from the diffusion-limited to the reaction-limited setting. We represent each reactive particle with a minimally-resolved 'blob' using many fewer degrees of freedom per particle than standard discretization approaches. More complicated or more highly resolved particle shapes can be built out of a collection of reactive blobs. We demonstrate numerically that the blob model can provide an accurate representation at low to moderate packing densities of the reactive particles, at a cost not much larger than solving a Poisson equation in the same domain. Unlike multipole expansion methods, our method does not require analytically-computed Green's functions, but rather, computes regularized discrete Green's functions on the fly by using a standard grid-based discretization of the Poisson equation. This allows for great flexibility in implementing different boundary conditions, coupling to fluid flow or thermal transport, and the inclusion of other effects such as temporal evolution and even nonlinearities. We develop multigrid-based preconditioners for solving the linear systems that arise when using implicit temporal discretizations or studying steady states. In the diffusion-limited case the resulting linear system is a saddle-point problem, the efficient solution of which remains a challenge for suspensions of many particles. We validate our method by comparing to published results on reaction-diffusion in ordered and disordered suspensions of reactive spheres.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Stokes-Einstein Relation at Moderate Schmidt Number

The Stokes-Einstein relation for the self-diffusion coefficient of a spherical particle suspended in an incompressible fluid is an asymptotic result in the limit of large Schmidt number, that is, when momentum diffuses much faster than the particle. When the Schmidt number is moderate, which happens in most particle methods for hydrodynamics, deviations from the Stokes-Einstein prediction are expected. We study these corrections computationally using a recently-developed minimally-resolved method for coupling particles to an incompressible fluctuating fluid in both two and three dimensions. We find that for moderate Schmidt numbers the diffusion coefficient is reduced relative to the Stokes-Einstein prediction by an amount inversely proportional to the Schmidt number in both two and three dimensions. We find, however, that the Einstein formula is obeyed at all Schmidt numbers, consistent with linear response theory. The numerical data is in good agreement with an approximate self-consistent theory, which can be used to estimate finite-Schmidt number corrections in a variety of methods. Our results indicate that the corrections to the Stokes-Einstein formula come primarily from the fact that the particle itself diffuses together with the momentum. Our study separates effects coming from corrections to no-slip hydrodynamics from those of finite separation of time scales, allowing for a better understanding of widely observed deviations from the Stokes-Einstein prediction in particle methods such as molecular dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Staggered Schemes for Fluctuating Hydrodynamics

We develop numerical schemes for solving the isothermal compressible and incompressible equations of fluctuating hydrodynamics on a grid with staggered momenta. We develop a second-order accurate spatial discretization of the diffusive, advective and stochastic fluxes that satisfies a discrete fluctuation-dissipation balance, and construct temporal discretizations that are at least second-order accurate in time deterministically and in a weak sense. Specifically, the methods reproduce the correct equilibrium covariances of the fluctuating fields to third (compressible) and second (incompressible) order in the time step, as we verify numerically. We apply our techniques to model recent experimental measurements of giant fluctuations in diffusively mixing fluids in a micro-gravity environment [A. Vailati et. al., Nature Communications 2:290, 2011]. Numerical results for the static spectrum of non-equilibrium concentration fluctuations are in excellent agreement between the compressible and incompressible simulations, and in good agreement with experimental results for all measured wavenumbers.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Diffusive Transport Enhanced by Thermal Velocity Fluctuations

We study the contribution of advection by thermal velocity fluctuations to the effective diffusion coefficient in a mixture of two indistinguishable fluids. The enhancement of the diffusive transport depends on the system size L and grows as \ln(L/L_0) in quasi two-dimensional systems, while in three dimensions it scales as L_0^{-1}-L^{-1}, where L_0 is a reference length. The predictions of a simple fluctuating hydrodynamics theory are compared to results from particle simulations and a finite-volume solver and excellent agreement is observed. Our results conclusively demonstrate that the nonlinear advective terms need to be retained in the equations of fluctuating hydrodynamics when modeling transport in small-scale finite systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Enhancement of Diffusive Transport by Nonequilibrium Thermal Fluctuations

We study the contribution of advection by thermal velocity fluctuations to the effective diffusion coefficient in a mixture of two identical fluids. The steady-state diffusive flux in a finite system subject to a concentration gradient is enhanced because of long-range correlations between concentration fluctuations and fluctuations of the velocity parallel to the concentration gradient. The enhancement of the diffusive transport depends on the system size L and grows as \ln(L/L_{0}) in quasi-two dimensional systems, while in three dimensions it grows as L_{0}^{-1}-L^{-1}, where L_{0} is a reference length. The predictions of a simple fluctuating hydrodynamics theory, closely related to second-order mode-mode coupling analysis, are compared to results from particle simulations and a finite-volume solver and excellent agreement is observed. We elucidate the direct connection to the long-time tail of the velocity autocorrelation function in finite systems, as well as finite-size corrections employed in molecular dynamics calculations. Our results conclusively demonstrate that the nonlinear advective terms need to be retained in the equations of fluctuating hydrodynamics when modeling transport in small-scale finite systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A First-Passage Kinetic Monte Carlo Algorithm for Complex Diffusion-Reaction Systems

We develop an asynchronous event-driven First-Passage Kinetic Monte Carlo (FPKMC) algorithm for continuous time and space systems involving multiple diffusing and reacting species of spherical particles in two and three dimensions. The FPKMC algorithm presented here is based on the method introduced in [Phys. Rev. Lett., 97:230602, 2006] and is implemented in a robust and flexible framework. Unlike standard KMC algorithms such as the n-fold algorithm, FPKMC is most efficient at low densities where it replaces the many small hops needed for reactants to find each other with large first-passage hops sampled from exact time-dependent Green's functions, without sacrificing accuracy. We describe in detail the key components of the algorithm, including the event-loop and the sampling of first-passage probability distributions, and demonstrate the accuracy of the new method. We apply the FPKMC algorithm to the challenging problem of simulation of long-term irradiation of metals, relevant to the performance and aging of nuclear materials in current and future nuclear power plants. The problem of radiation damage spans many decades of time-scales, from picosecond spikes caused by primary cascades, to years of slow damage annealing and microstructure evolution. Our implementation of the FPKMC algorithm has been able to simulate the irradiation of a metal sample for durations that are orders of magnitude longer than any previous simulations using the standard Object KMC or more recent asynchronous algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Thermodynamically-Consistent Non-Ideal Stochastic Hard-Sphere Fluid

A grid-free variant of the Direct Simulation Monte Carlo (DSMC) method is proposed, named the Isotropic DSMC (I-DSMC) method, that is suitable for simulating dense fluid flows at molecular scales. The I-DSMC algorithm eliminates all grid artifacts from the traditional DSMC algorithm; it is Galilean invariant and microscopically isotropic. The stochastic collision rules in I-DSMC are modified to yield a non-ideal structure factor that gives consistent compressibility, as first proposed in [Phys. Rev. Lett. 101:075902 (2008)]. The resulting Stochastic Hard Sphere Dynamics (SHSD) fluid is empirically shown to be thermodynamically identical to a deterministic Hamiltonian system of penetrable spheres interacting with a linear core pair potential, well-described by the hypernetted chain (HNC) approximation. We apply a stochastic Enskog kinetic theory for the SHSD fluid to obtain estimates for the transport coefficients that are in excellent agreement with particle simulations over a wide range of densities and collision rates. The fluctuating hydrodynamic behavior of the SHSD fluid is verified by comparing its dynamic structure factor against theory based on the Landau-Lifshitz Navier-Stokes equations. We also study the Brownian motion of a nano-particle suspended in an SHSD fluid and find a long-time power-law tail in its velocity autocorrelation function consistent with hydrodynamic theory and molecular dynamics calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

First-Passage Kinetic Monte Carlo method

We present a new efficient method for Monte Carlo simulations of diffusion-reaction processes. First introduced by us in [Phys. Rev. Lett., 97:230602, 2006], the new algorithm skips the traditional small diffusion hops and propagates the diffusing particles over long distances through a sequence of super-hops, one particle at a time. By partitioning the simulation space into non-overlapping protecting domains each containing only one or two particles, the algorithm factorizes the N-body problem of collisions among multiple Brownian particles into a set of much simpler single-body and two-body problems. Efficient propagation of particles inside their protective domains is enabled through the use of time-dependent Green's functions (propagators) obtained as solutions for the first-passage statistics of random walks. The resulting Monte Carlo algorithm is event-driven and asynchronous; each Brownian particle propagates inside its own protective domain and on its own time clock. The algorithm reproduces the statistics of the underlying Monte-Carlo model exactly. Extensive numerical examples demonstrate that for an important class of diffusion-reaction models the new algorithm is efficient at low particle densities, where other existing algorithms slow down severely.

0 citations0 reviewsNo code linkedOpen access