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

Aleksandar Donev

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

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softphysics.flu-dynmath.NAphysics.comp-phcond-mat.mes-hallNumerical Analysismath-phmath.MGmath.MPmath.PR

Catalog footprint

What is connected

15works
11topics
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

15 published item(s)

preprint2022arXiv

The hydrodynamics of a twisting, bending, inextensible fiber in Stokes flow

In swimming microorganisms and the cell cytoskeleton, inextensible fibers resist bending and twisting, and interact with the surrounding fluid to cause or resist large-scale fluid motion. In this paper, we develop a novel numerical method for the simulation of cylindrical fibers by extending our previous work on inextensible bending fibers [Maxian et al., Phys. Rev. Fluids 6 (1), 014102] to fibers with twist elasticity. In our "Euler" model, twist is a scalar function that measures the deviation of the fiber cross section relative to a twist-free frame, the fiber exerts only torque parallel to the centerline on the fluid, and the perpendicular components of the rotational fluid velocity are discarded in favor of the translational velocity. In the first part of this paper, we justify this model by comparing it to another commonly-used "Kirchhoff" formulation where the fiber exerts both perpendicular and parallel torque on the fluid, and the perpendicular angular fluid velocity is required to be consistent with the translational fluid velocity. We then develop a spectral numerical method for the hydrodynamics of the Euler model. We define hydrodynamic mobility operators using integrals of the Rotne-Prager-Yamakawa tensor, and evaluate these integrals through a novel slender-body quadrature, which requires on the order of 10 points along the fiber to obtain several digits of accuracy. We demonstrate that this choice of mobility removes the unphysical negative eigenvalues in the translation-translation mobility associated with asymptotic slender body theories, and ensures strong convergence of the fiber velocity and weak convergence of the fiber constraint forces. We pair the spatial discretization with a semi-implicit temporal integrator to confirm the negligible contribution of twist elasticity to the relaxation dynamics of a bent fiber and study the instability of a twirling fiber.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An integral-based spectral method for inextensible slender fibers in Stokes flow

Every animal cell is filled with a cytoskeleton, a dynamic gel made of inextensible fibers, such as microtubules, actin fibers, and intermediate filaments, all suspended in a viscous fluid. Numerical simulation of this gel is challenging because the fiber aspect ratios can be as large as $10^4$. We describe a new method for rapidly computing the dynamics of inextensible slender filaments in periodically-sheared Stokes flow. The dynamics of the filaments are governed by a nonlocal slender body theory which we partially reformulate in terms of the Rotne-Prager-Yamakawa hydrodynamic tensor. To enforce inextensibility, we parameterize the space of inextensible fiber motions and strictly confine the dynamics to the manifold of inextensible configurations. To do this, we introduce a set of Lagrange multipliers for the tensile force densities on the filaments and impose the constraint of no virtual work in an $L^2$ weak sense. We augment this approach with a spectral discretization of the local and nonlocal slender body theory operators which is linear in the number of unknowns and gives improved spatial accuracy over approaches based on solving a line tension equation. For dynamics, we develop a second-order semi-implicit temporal integrator which requires at most a few evaluations of nonlocal hydrodynamics and a few block diagonal linear solves per time step. After demonstrating the improved accuracy and robustness of our approach through numerical examples, we apply our formulation to a permanently cross-linked actin mesh in a background oscillatory shear flow. We observe a characteristic frequency at which the network transitions from quasi-static, primarily elastic behavior to dynamic, primarily viscous behavior. We find that nonlocal hydrodynamics increases the viscous modulus by as much as 25%, even for semi-dilute fiber suspensions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Driven dynamics in dense suspensions of microrollers

We perform detailed computational and experimental measurements of the driven dynamics of a dense, uniform suspension of sedimented microrollers driven by a magnetic field rotating around an axis parallel to the floor. We develop a lubrication-corrected Brownian Dynamics method for dense suspensions of driven colloids sedimented above a bottom wall. The numerical method adds lubrication friction between nearby pairs of particles, as well as particles and the bottom wall, to a minimally-resolved model of the far-field hydrodynamic interactions. Our experiments combine fluorescent labeling with particle tracking to trace the trajectories of individual particles in a dense suspension, and to measure their propulsion velocities. Previous computational studies [B. Sprinkle et al., J. Chem. Phys., 147, 244103, 2017] predicted that at sufficiently high densities a uniform suspension of microrollers separates into two layers, a slow monolayer right above the wall, and a fast layer on top of the bottom layer. Here we verify this prediction, showing good quantitative agreement between the bimodal distribution of particle velocities predicted by the lubrication-corrected Brownian Dynamics and those measured in the experiments. The computational method accurately predicts the rate at which particles are observed to switch between the slow and fast layers in the experiments. We also use our numerical method to demonstrate the important role that pairwise lubrication plays in motility-induced phase separation in dense monolayers of colloidal microrollers, as recently suggested for suspensions of Quincke rollers [D. Geyer et al., Physical Review X, 9(3), 031043, 2019].

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Fluctuation-enhanced electric conductivity in electrolyte solutions

In this letter we analyze the effects of an externally applied electric field on thermal fluctuations for a fluid containing charged species. We show in particular that the fluctuating Poisson-Nernst-Planck equations for charged multispecies diffusion coupled with the fluctuating fluid momentum equation, result in enhanced charge transport. Although this transport is advective in nature, it can macroscopically be represented as electrodiffusion with renormalized electric conductivity. We calculate the renormalized electric conductivity by deriving and integrating the structure factor coefficients of the fluctuating quantities and show that the renormalized electric conductivity and diffusion coefficients are consistent although they originate from different noise terms. In addition, the fluctuating hydrodynamics approach recovers the electrophoretic and relaxation corrections obtained by Debye-Huckel-Onsager theory, and provides a quantitative theory that predicts a non-zero cross-diffusion Maxwell-Stefan coefficient that agrees well with experimental measurements. Finally, we show that strong applied electric fields result in anisotropically enhanced velocity fluctuations and reduced fluctuations of salt concentrations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Low Mach Number Fluctuating Hydrodynamics for Electrolytes

We formulate and study computationally the low Mach number fluctuating hydrodynamic equations for electrolyte solutions. We are interested in studying transport in mixtures of charged species at the mesoscale, down to scales below the Debye length, where thermal fluctuations have a significant impact on the dynamics. Continuing our previous work on fluctuating hydrodynamics of multicomponent mixtures of incompressible isothermal miscible liquids (A. Donev, et al., Physics of Fluids, 27, 3, 2015), we now include the effect of charged species using a quasielectrostatic approximation. Localized charges create an electric field, which in turn provides additional forcing in the mass and momentum equations. Our low Mach number formulation eliminates sound waves from the fully compressible formulation and leads to a more computationally efficient quasi-incompressible formulation. We demonstrate our ability to model saltwater (NaCl) solutions in both equilibrium and nonequilibrium settings. We show that our algorithm is second-order in the deterministic setting, and for length scales much greater than the Debye length gives results consistent with an electroneutral/ambipolar approximation. In the stochastic setting, our model captures the predicted dynamics of equilibrium and nonequilibrium fluctuations. We also identify and model an instability that appears when diffusive mixing occurs in the presence of an applied electric field.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Brownian Dynamics of Confined Rigid Bodies

We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust, space efficient, and easy to accumulate. We construct a system of overdamped Langevin equations in the quaternion representation that accounts for hydrodynamic effects, preserves the unit-norm constraint on the quaternion, and is time reversible with respect to the Gibbs-Boltzmann distribution at equilibrium. We introduce two schemes for temporal integration of the overdamped Langevin equations of motion, one based on the Fixman midpoint method and the other based on a random finite difference approach, both of which ensure the correct stochastic drift term is captured in a computationally efficient way. We study several examples of rigid colloidal particles diffusing near a no-slip boundary, and demonstrate the importance of the choice of tracking point on the measured translational mean square displacement (MSD). We examine the average short-time as well as the long-time quasi-two-dimensional diffusion coefficient of a rigid particle sedimented near a bottom wall due to gravity. For several particle shapes we find a choice of tracking point that makes the MSD essentially linear with time, allowing us to estimate the long-time diffusion coefficient efficiently using a Monte Carlo method. However, in general such a special choice of tracking point does not exist, and numerical techniques for simulating long trajectories, such as the ones we introduce here, are necessary to study diffusion on long timescales.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coupling a nano-particle with isothermal fluctuating hydrodynamics: Coarse-graining from microscopic to mesoscopic dynamics

We derive a coarse-grained description of the dynamics of a nanoparticle immersed in an isothermal simple fluid by performing a systematic coarse graining of the underlying microscopic dynamics. As coarse-grained or relevant variables we select the position of the nanoparticle and the \emph{total} mass and momentum density field of the fluid, which are locally conserved slow variables because they are defined to include the contribution of the nanoparticle. The theory of coarse graining based on the Zwanzing projection operator leads us to a system of stochastic \emph{ordinary} differential equations (SODEs) that are closed in the relevant variables. We demonstrate that our discrete coarse-grained equations are consistent with a Petrov-Galerkin finite-element discretization of a system of formal stochastic \emph{partial} differential equations (SPDEs) which resemble previously-used phenomenological models based on fluctuating hydrodynamics. Under suitable approximations we obtain \emph{closed} approximations of the coarse-grained dynamics in a manner which gives them a clear physical interpretation, and provides \emph{explicit} microscopic expressions for all of the coefficients appearing in the closure. Our work leads to a model for dilute nanocolloidal suspensions that can be simulated effectively using feasibly short molecular dynamics simulations as input to a FEM fluctuating hydrodynamic solver.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dynamic scaling for the growth of non-equilibrium fluctuations during thermophoretic diffusion in microgravity

Diffusion processes are widespread in biological and chemical systems, where they play a fundamental role in the exchange of substances at the cellular level and in determining the rate of chemical reactions. Recently, the classical picture that portrays diffusion as random uncorrelated motion of molecules has been revised, when it was shown that giant non-equilibrium fluctuations develop during diffusion processes. Under microgravity conditions and at steady-state, non-equilibrium fluctuations exhibit scale invariance and their size is only limited by the boundaries of the system. In this work, we investigate the onset of non-equilibrium concentration fluctuations induced by thermophoretic diffusion in microgravity, a regime not accessible to analytical calculations but of great relevance for the understanding of several natural and technological processes. A combination of state of the art simulations and experiments allows us to attain a fully quantitative description of the development of fluctuations during transient diffusion in microgravity. Both experiments and simulations show that during the onset the fluctuations exhibit scale invariance at large wave vectors. In a broader range of wave vectors simulations predict a spinodal-like growth of fluctuations, where the amplitude and length-scale of the dominant mode are determined by the thickness of the diffuse layer.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Finite element discretization of non-linear diffusion equations with thermal fluctuations

We present a finite element discretization of a non-linear diffusion equation used in the field of critical phenomena and, more recently, in the context of Dynamic Density Functional Theory. The discretized equation preserves the structure of the continuum equation. Specifically, it conserves the total number of particles and fulfills an H-theorem as the original partial differential equation. Guided by the Theory of Coarse-Graining, we discuss the inclusion of thermal fluctuations in the non-linear diffusion equation. This sheds light on the meaning of such a fluctuating hydrodynamics equation and to the limitations of the approximations involved. The methodology proposed for the introduction of thermal fluctuations in finite element methods is general and valid for both regular and irregular grids in arbitrary dimensions. We focus here on simulations of the Ginzburg-Landau free energy functional using both regular and irregular 1D grids. Convergence of the numerical results is obtained for the static and dynamic structure factors as the resolution of the grid is increased.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Slowing-down of non-equilibrium concentration fluctuations in confinement

Fluctuations in a fluid are strongly affected by the presence of a macroscopic gradient making them long-ranged and enhancing their amplitude. While small-scale fluctuations exhibit diffusive lifetimes, larger-scale fluctuations live shorter because of gravity, as theoretically and experimentally well-known. We explore here fluctuations of even larger size, comparable to the extent of the system in the direction of the gradient, and find experimental evidence of a dramatic slowing-down in their dynamics. We recover diffusive behaviour for these strongly-confined fluctuations, but with a diffusion coefficient that depends on the solutal Rayleigh number. Results from dynamic shadowgraph experiments are complemented by theoretical calculations and numerical simulations based on fluctuating hydrodynamics, and excellent agreement is found. The study of the dynamics of non-equilibrium fluctuations allows to probe and measure the competition of physical processes such as diffusion, buoyancy and confinement.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fluctuating hydrodynamics of multi-species, non-reactive mixtures

In this paper we discuss the formulation of the fuctuating Navier-Stokes (FNS) equations for multi-species, non-reactive fluids. In particular, we establish a form suitable for numerical solution of the resulting stochastic partial differential equations. An accurate and efficient numerical scheme, based on our previous methods for single species and binary mixtures, is presented and tested at equilibrium as well as for a variety of non-equilibrium problems. These include the study of giant nonequilibrium concentration fluctuations in a ternary mixture in the presence of a diffusion barrier, the triggering of a Rayleigh-Taylor instability by diffusion in a four-species mixture, as well as reverse diffusion in a ternary mixture. Good agreement with theory and experiment demonstrates that the formulation is robust and can serve as a useful tool in the study of thermal fluctuations for multi-species fluids. The extension to include chemical reactions will be treated in a sequel paper.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Metropolis Integration Schemes for Self-Adjoint Diffusions

We present explicit methods for simulating diffusions whose generator is self-adjoint with respect to a known (but possibly not normalizable) density. These methods exploit this property and combine an optimized Runge-Kutta algorithm with a Metropolis-Hastings Monte-Carlo scheme. The resulting numerical integration scheme is shown to be weakly accurate at finite noise and to gain higher order accuracy in the small noise limit. It also permits to avoid computing explicitly certain terms in the equation, such as the divergence of the mobility tensor, which can be tedious to calculate. Finally, the scheme is shown to be ergodic with respect to the exact equilibrium probability distribution of the diffusion when it exists. These results are illustrated on several examples including a Brownian dynamics simulation of DNA in a solvent. In this example, the proposed scheme is able to accurately compute dynamics at time step sizes that are an order of magnitude (or more) larger than those permitted with commonly used explicit predictor-corrector schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modeling Multi-phase Flow using Fluctuating Hydrodynamics

Fluctuating hydrodynamics provides a model for fluids at mesoscopic scales where thermal fluctuations can have a significant impact on the behavior of the system. Here we investigate a model for fluctuating hydrodynamics of a single component, multiphase flow in the neighborhood of the critical point. The system is modeled using a compressible flow formulation with a van der Waals equation of state, incorporating a Korteweg stress term to treat interfacial tension. We present a numerical algorithm for modeling this system based on an extension of algorithms developed for fluctuating hydrodynamics for ideal fluids. The scheme is validated by comparison of measured structure factors and capillary wave spectra with equilibrium theory. We also present several non-equilibrium examples to illustrate the capability of the algorithm to model multi-phase fluid phenomena in a neighborhood of the critical point. These examples include a study of the impact of fluctuations on the spinodal decomposition following a rapid quench, as well as the piston effect in a cavity with supercooled walls. The conclusion in both cases is that thermal fluctuations affect the size and growth of the domains in off-critical quenches.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Superdense Crystal Packings of Ellipsoids

Particle packing problems have fascinated people since the dawn of civilization, and continue to intrigue mathematicians and scientists. Resurgent interest has been spurred by the recent proof of Kepler's conjecture: the face-centered cubic lattice provides the densest packing of equal spheres with a packing fraction $ϕ\approx0.7405$ \cite{Kepler_Hales}. Here we report on the densest known packings of congruent ellipsoids. The family of new packings are crystal (periodic) arrangements of nearly spherically-shaped ellipsoids, and always surpass the densest lattice packing. A remarkable maximum density of $ϕ\approx0.7707$ is achieved for both prolate and oblate ellipsoids with aspect ratios of $\sqrt{3}$ and $1/\sqrt{3}$, respectively, and each ellipsoid has 14 touching neighbors. Present results do not exclude the possibility that even denser crystal packings of ellipsoids could be found, and that a corresponding Kepler-like conjecture could be formulated for ellipsoids.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Generalized von Smoluchowski model of reaction rates, with reacting particles and a mobile trap

We study diffusion-limited coalescence, A+A<-->A, in one dimension, in the presence of a diffusing trap. The system may be regarded as a generalization of von Smoluchowski's model for reaction rates, in that: (a) it includes reactions between the particles surrounding the trap, and (b) the trap is mobile -- both considerations which render the model more physically relevant. As seen from the trap's frame of reference, the motion of the particles is highly correlated, because of the motion of the trap. An exact description of the long -time asymptotic limit is found using the IPDF method, and exploiting a "shielding" property of reversible coalescence that was discovered recently. In the case where the trap also acts as a source -- giving birth to particles -- the shielding property breaks down, but we find an "equivalence principle": Trapping and diffusion of the trap may be compensated by an appropriate rate of birth, such that the steady state of the system is identical with the equilibrium state in the absence of a trap.

0 citations0 reviewsNo code linkedOpen access