Researcher profile

J. Ravi Prakash

J. Ravi Prakash contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

30 published item(s)

preprint2022arXiv

Universality of dilute solutions of ring polymers in the thermal crossover region between $θ$ and athermal solvents

Due to their unique topology of having no chain ends, dilute solutions of ring polymers exhibit behaviour distinct from their linear chain counterparts. The universality of their static and dynamic properties, as a function of solvent quality $z$ in the thermal crossover regime between $θ$ and athermal solvents, is studied here using Brownian dynamics simulations. The universal ratio $U_{\text{RD}}$ of the radius of gyration $R_g$ to the hydrodynamic radius $R_H$ is determined, and a comparative study of the swelling ratio $α_g$ of the radius of gyration, the swelling ratio $α_H$ of the hydrodynamic radius, and the swelling ratio $α_X$ of the mean polymer stretch $X$ along the $x$-axis, for linear and ring polymers, is carried out. The ratio $U_{\text{RD}}$ for dilute ring polymer solutions is found to converge asymptotically to a constant value as $z \to \infty$, which is a major difference from the behaviour of solutions of linear chains, where no such asymptotic limit exists. Additionally, the ratio of the mean stretch along the $x$-axis to the hydrodynamic radius, $(X/R_H)$, is found to be independent of $z$ for polymeric rings, unlike in the case for linear polymers. These results indicate a fundamental difference in the scaling of static and dynamic properties of rings and linear chains in the thermal crossover regime.

preprint2021arXiv

How important are fluctuations in the treatment of internal friction in polymers?

The Rouse model with internal friction (RIF), a widely used theoretical framework to interpret the effects of internal friction on conformational transitions in biomolecules, is shown to be an approximate treatment that is based on preaveraging internal friction. By comparison with Brownian dynamics simulations of an exact coarse-grained model that incorporates fluctuations in internal friction, the accuracy of the preaveraged model predictions is examined both at and away from equilibrium. While the two models predict intrachain autocorrelations that approach each other for long enough chain segments, they differ in their predictions for shorter segments. Furthermore, the two models differ qualitatively in their predictions for the chain extension and viscosity in shear flow, which is taken to represent a prototypical out-of-equilibrium condition.

preprint2020arXiv

Wet and dry internal friction can be measured with the Jarzynski equality

The existence of two types of internal friction wet and dry is revisited, and a simple protocol is proposed for distinguishing between the two types and extracting the appropriate internal friction coefficient. The scheme requires repeatedly stretching a polymer molecule, and measuring the average work dissipated in the process by applying the Jarzynski equality. The internal friction coefficient is then estimated from the average dissipated work in the extrapolated limit of zero solvent viscosity. The validity of the protocol is established through analytical calculations on a one-dimensional free-draining Hookean spring-dashpot model for a polymer, and Brownian dynamics simulations of: (a) a single-mode nonlinear spring-dashpot model for a polymer, and (b) a finitely extensible bead-spring chain with cohesive intra-chain interactions, both of which incorporate fluctuating hydrodynamic interactions. Well-established single-molecule manipulation techniques, such as optical tweezer-based pulling, can be used to implement the suggested protocol experimentally.

preprint2019arXiv

Uniaxial extensional viscosity of semidilute DNA solutions

The extensional rheology of polymeric liquids has been extensively examined through experiments and theoretical predictions. However, a systematic study of the extensional rheology of polymer solutions in the semidilute regime, in terms of examining the effects of concentration and molecular weight, has not been carried out so far. Prior studies of the shear rheology of semidilute polymer solutions have demonstrated that their behaviour is distinctively different from that observed in the dilute and concentrated regimes. This difference in behaviour is anticipated to be even more pronounced in extensional flows. In this work, the extensional rheology of linear, double-stranded DNA molecules, spanning an order of magnitude of molecular weights (25 to 289 kilobasepairs) and concentrations (0.03 to 0.3 mg/ml), has been investigated. DNA solutions are now used routinely as model polymeric systems due to their near-perfect monodispersity. Measurements have been carried out with a filament stretching rheometer since it is the most reliable method for obtaining an estimate of the elongational stress growth of a polymer solution. Transient and steady-state uniaxial extensional viscosities of DNA dissolved in a solvent under excess salt conditions, with a high concentration of sucrose in order to achieve a sufficiently high solvent viscosity, have been determined in the semidilute regime at room temperature. The dependence of the steady state uniaxial extensional viscosity on molecular weight, concentration and extension rate is measured with a view to determining if data collapse can be observed with an appropriate choice of variables. Steady state shear viscosity measurements suggest that sucrose-DNA interactions might play a role in determining the observed rheological behaviour of semidilute DNA solutions with sucrose as a component in the solvent.

preprint2019arXiv

Universal dynamics of dilute and semidilute solutions of flexible linear polymers

Experimental observations and computer simulations have recently revealed universal aspects of polymer solution behaviour that were previously unknown. This progress has been made possible due to developments in experimental methodologies and simulation techniques, and the formulation of scaling arguments that have identified the proper scaling variables with which to achieve data collapse, both for equilibrium and nonequilibrium properties. In this review, these advances in our understanding of universal polymer solution behaviour are discussed, along with a brief overview of prior experimental observations and theoretical predictions in order to provide a context for the discussion. Properties for which the existence of universality is still unclear, and those for which it has not yet been verified theoretically are highlighted, and theoretical predictions of universal behaviour that require experimental validation are pointed out.

preprint2019arXiv

Universality of the collapse transition of sticky polymers

The universality of the swelling of the radius of gyration of a homopolymer relative to its value in the $θ$ state, independent of polymer-solvent chemistry, in the crossover regime between $θ$ and athermal solvent conditions, is well known. Here we study, by Brownian dynamics, a polymer model where a subset of monomers is labelled as "stickers". The mutual interaction of the stickers is more attractive than those of the other ("backbone") monomers, and has the additional important characteristic of "functionality" $φ$, i.e., the maximum number of stickers that can locally bind to a given sticker. A saturated bond formed in this manner remains bound until it breaks due to thermal fluctuations, a requirement which can be viewed as an additional Boolean degree of freedom that describes the bonding. This, in turn, makes the question of the order of the collapse transition a non-trivial one. Nevertheless, for the parameters that we have studied (in particular, $φ=1$), we find a standard second-order $θ$ collapse, using a renormalised solvent quality parameter that takes into account the increased average attraction due to the presence of stickers. We examine the swelling of the radius of gyration of such a sticky polymer relative to its value in the altered $θ$ state, using a novel potential to model the various excluded volume interactions that occur between the monomers on the chain. We find that the swelling of such sticky polymers is identical to the universal swelling of homopolymers in the thermal crossover regime. Additionally, for our model, the Kuhn segment length under $θ$ conditions is found to be the same for chains with and without stickers.

preprint2018arXiv

Rheological consequences of wet and dry friction in a dumbbell model with hydrodynamic interactions and internal viscosity

The effect of fluctuating internal viscosity and hydrodynamic interactions on a range of rheological properties of dilute polymer solutions is examined using a finitely extensible dumbbell model for a polymer. Brownian dynamics simulations are used to compute both transient and steady state viscometric functions in shear flow. The results enable a careful differentiation of the influence, on rheological properties, of solvent-mediated friction from that of a dissipative mechanism that is independent of solvent viscosity. In particular, hydrodynamic interactions have a significant influence on the magnitude of the stress jump at the inception of shear flow, and on the transient viscometric functions, but a negligible effect on the steady state viscometric functions at high shear rates. Zero-shear rate viscometric functions of free-draining dumbbells remain essentially independent of the internal viscosity parameter, as predicted by the Gaussian approximation, but the inclusion of hydrodynamic interactions induces a dependence on both the hydrodynamic interaction and the internal viscosity parameter. Large values of the internal viscosity parameter lead to linear viscoelastic predictions that mimic the behavior of rigid dumbbell solutions. On the other hand, steady-shear viscometric functions at high shear rates differ in general from those for rigid dumbbells, depending crucially on the finite extensibility of the dumbbell spring.

preprint2018arXiv

Shear thinning in dilute and semidilute solutions of polystyrene and DNA

The viscosity of dilute and semidilute unentangled DNA solutions, in steady simple shear flow, has been measured across a range of temperatures and concentrations. For polystyrene solutions, measurements of viscosity have been carried out in the semidilute unentangled regime, while results of prior experimental measurements in the dilute regime have been used for the purpose of data analysis, and for comparison with the behaviour of DNA solutions. Interpretation of the shear rate dependence of viscosity in terms of suitably defined non-dimensional variables, is shown to lead to master plots, independent of temperature and concentration, in each of the two concentration regimes. In the case of semidilute unentangled solutions, defining the Weissenberg number in terms of a concentration dependent large scale relaxation time is found not to lead to data collapse across different concentrations. On the other hand, the use of an alternative relaxation time, with the concentration dependence of a single correlation blob, suggests the existence of universal shear thinning behaviour at large shear rates.

preprint2016arXiv

Direct observation of DNA dynamics in semi-dilute solutions in extensional flow

The dynamic behavior of semi-dilute polymer solutions is governed by an interplay between solvent quality, concentration, molecular weight, and flow type. Semi-dilute solutions are characterized by large fluctuations in concentration, wherein polymer coils interpenetrate but may not be topologically entangled at equilibrium. In non-equilibrium flows, it is generally thought that polymer chains can self-entangle in semi-dilute solutions, thereby leading to entanglements in solutions that are nominally unentangled at equilibrium. Despite recent progress, we still lack a complete molecular-level understanding of these dynamics. In this work, we use single molecule techniques to study the dynamics of semi-dilute solutions of DNA in planar extensional flow, including polymer relaxation from high stretch, transient stretching dynamics in step-strain experiments, and steady-state stretching in flow. Our results are consistent with a power-law scaling of the longest polymer relaxation time in semi-dilute solutions, revealing an effective excluded volume exponent $ν$ = 0.56. We further studied the non-equilibrium stretching dynamics in extensional flow, with results showing a decrease in transient polymer stretch at moderate Weissenberg number and a milder coil-to-stretch transition in semi-dilute solutions compared to dilute solutions. Interestingly, a unique set of molecular conformations during the transient stretching process for single polymers in semi-dilute solutions is observed, which suggests transient stretching pathways for polymer chains in semi-dilute solutions are qualitatively different than dilute solutions due to intermolecular interactions. Taken together, this work provides a molecular framework for understanding the non-equilibrium stretching dynamics of semi-dilute solutions in strong flows.

preprint2016arXiv

Effect of stretching-induced changes in hydrodynamic screening on coil-stretch hysteresis of unentangled polymer solutions

Extensional rheometry and Brownian Dynamics simulations of flexible polymer solutions confirm predictions based on blob concepts that coil--stretch hysteresis in extensional flows increases with concentration, reaching a maximum at the critical overlap concentration $c^\ast$ before progressively vanishing in the semidilute regime. These observations demonstrate that chain stretching strengthens intermolecular hydrodynamic screening in dilute solutions, but weakens it in semidilute solutions. Flow can thus strongly modify the concentration dependence of viscoelastic properties of polymer solutions.

preprint2016arXiv

Equilibrium binding energies from fluctuation theorems and force spectroscopy simulations

Brownian dynamics simulations are used to study the detachment of a particle from a substrate. Although the model is simple and generic, we attempt to map its energy, length and time scales onto a specific experimental system, namely a bead that is weakly bound to a cell and then removed by an optical tweezer. The external driving force arises from the combined optical tweezer and substrate potentials, and thermal fluctuations are taken into account by a Brownian force. The Jarzynski equality and Crooks' fluctuation theorem are applied to obtain the equilibrium free energy difference between the final and initial states. To this end, we sample non--equilibrium work trajectories for various tweezer pulling rates. We argue that this methodology should also be feasible experimentally for the envisioned system. Furthermore, we outline how the measurement of a whole free energy profile would allow the experimentalist to retrieve the unknown substrate potential by means of a suitable deconvolution. The influence of the pulling rate on the accuracy of the results is investigated, and umbrella sampling is used to obtain the equilibrium probability of particle escape for a variety of trap potentials.

preprint2016arXiv

Parameter-free prediction of DNA dynamics in planar extensional flow of semidilute solutions

The dynamics of individual DNA molecules in semidilute solutions undergoing planar extensional flow is simulated using a multi-particle Brownian dynamics algorithm, which incorporates hydrodynamic and excluded volume interactions in the context of a coarse-grained bead-spring chain model for DNA. The successive fine-graining protocol [1, 2], in which simulation data acquired for bead-spring chains with increasing values of the number of beads $N_b$, is extrapolated to the number of Kuhn steps $N_\text{K}$ in DNA (while keeping key physical parameters invariant), is used to obtain parameter-free predictions for a range of Weissenberg numbers and Hencky strain units. A systematic comparison of simulation predictions is carried out with the experimental observations of [3], who have recently used single molecule techniques to investigate the dynamics of dilute and semidilute solutions of $λ$-phage DNA in planar extensional flow. In particular, they examine the response of individual chains to step-strain deformation followed by cessation of flow, thereby capturing both chain stretch and relaxation in a single experiment. The successive fine-graining technique is shown to lead to quantitatively accurate predictions of the experimental observations in the stretching and relaxation phases. Additionally, the transient chain stretch following a step strain deformation is shown to be much smaller in semidilute solutions than in dilute solutions, in agreement with experimental observations.

preprint2015arXiv

Significance of thermal fluctuations and hydrodynamic interactions in receptor-ligand mediated adhesive dynamics of a spherical particle in wall bound shear flow

The dynamics of adhesion of a spherical micro-particle to a ligand-coated wall, in shear flow, is studied using a Langevin equation that accounts for thermal fluctuations, hydrodynamic interactions and adhesive interactions. Contrary to the conventional assumption that thermal fluctuations play a negligible role at high P$\acute{e}$clet numbers, we find that for particles with low surface densities of receptors, rotational diffusion caused by fluctuations about the flow and gradient directions aids in bond formation, leading to significantly greater adhesion on average, compared to simulations where thermal fluctuations are completely ignored. The role of wall hydrodynamic interactions on the steady state motion of a particle, when the particle is close to the wall, has also been explored. At high P$\acute{e}$clet numbers, the shear induced force that arises due to the stresslet part of the Stokes dipole, plays a dominant role, reducing the particle velocity significantly, and affecting the states of motion of the particle. The coupling between the translational and rotational degrees of freedom of the particle, brought about by the presence of hydrodynamic interactions, is found to have no influence on the binding dynamics. On the other hand, the drag coefficient, which depends on the distance of the particle from the wall, plays a crucial role at low rates of bond formation. A significant difference in the effect of both the shear force and the position dependent drag force, on the states of motion of the particle, is observed when the P$\acute{e}$let number is small.

preprint2015arXiv

Size, shape and diffusivity of a single Debye-Hückel polyelectrolyte chain in solution

Brownian dynamics simulations of a coarse-grained bead-spring chain model, with Debye-Hückel electrostatic interactions between the beads, are used to determine the root-mean-square end-to-end vector, the radius of gyration, and various shape functions (defined in terms of eigenvalues of the radius of gyration tensor) of a weakly-charged polyelectrolyte chain in solution, in the limit of low polymer concentration. The long-time diffusivity is calculated from the mean square displacement of the centre of mass of the chain, with hydrodynamic interactions taken into account through the incorporation of the Rotne-Prager-Yamakawa tensor. Simulation results are interpreted in the light of the OSFKK blob scaling theory (R. Everaers, A. Milchev, and V. Yamakov, Eur. Phys. J. E 8, 3 (2002)) which predicts that all solution properties are determined by just two scaling variables--the number of electrostatic blobs $X$, and the reduced Debye screening length, $Y$. We identify three broad regimes, the ideal chain regime at small values of $Y$, the blob-pole regime at large values of $Y$, and the crossover regime at intermediate values of $Y$, within which the mean size, shape, and diffusivity exhibit characteristic behaviours. In particular, when simulation results are recast in terms of blob scaling variables, universal behaviour independent of the choice of bead-spring chain parameters, and the number of blobs $X$, is observed in the ideal chain regime and in much of the crossover regime, while the existence of logarithmic corrections to scaling in the blob-pole regime leads to non-universal behaviour.

preprint2015arXiv

Viscoelastic fluid flow in a 2D channel bounded above by a deformable finite thickness elastic wall

The steady flow of three viscoelastic fluids (Oldroyd-B, FENE-P, and Owens model for blood) in a two-dimensional channel, partly bound by a deformable, finite thickness neo-Hookean solid, is computed. The limiting Weissenberg number beyond which computations fail to converge is found to increase with increasing dimensionless solid elasticity parameter Γ, following the trend Owens > FENE- P > Oldroyd-B. The highly shear thinning nature of Owens model leads to the elastic solid always collapsing into the channel, for the wide range of values of Γ considered here. In the case of the FENE-P and Oldroyd-B models, however, the fluid-solid interface can be either within the channel, or bulge outwards, depending on the value of Γ. This behaviour differs considerably from predictions of earlier models that treat the deformable solid as a zero-thickness membrane, in which case the membrane always lies within the channel. The capacity of the solid wall to support both pressure and shear stress, in contrast to the zero-thickness membrane that only responds to pressure, is responsible for the observed difference. Compar- ison of the stress and velocity fields in the channel for the three viscoelastic fluids, with the predictions for a Newtonian fluid, reveals that shear thinning rather than elasticity is the key source of the observed differences in behaviour.

preprint2014arXiv

Brownian dynamics simulations of planar mixed flows of polymer solutions at finite concentrations

Periodic boundary conditions for planar mixed flows are implemented in the context of a multi-chain Brownian dynamics simulation algorithm. The effect of shear rate $\dotγ$, and extension rate $\dotε$, on the size of polymer chains, $\left<R_e^2\right>$, and on the polymer contribution to viscosity, $η$, is examined for solutions of FENE dumbbells at finite concentrations, with excluded volume interactions between the beads taken into account. The influence of the mixedness parameter, $χ$, and flow strength, $\dotΓ$, on $\left<R_e^2\right>$ and $η$, is also examined, where $χ\rightarrow 0$ corresponds to pure shear flow, and $χ\rightarrow 1$ corresponds to pure extensional flow. It is shown that there exists a critical value, $χ_\text{c}$, such that the flow is shear dominated for $χ< χ_\text{c}$, and extension dominated for $χ> χ_\text{c}$.

preprint2014arXiv

The viscosity radius in dilute polymer solutions: Universal behaviour from DNA rheology and Brownian dynamics simulations

The swelling of the viscosity radius, $α_η$, and the universal viscosity ratio, $U_{ηR}$, have been determined experimentally for linear DNA molecules in dilute solutions with excess salt, and numerically by Brownian dynamics simulations, as a function of the solvent quality. In the latter instance, asymptotic parameter free predictions have been obtained by extrapolating simulation data for finite chains to the long chain limit. Experiments and simulations show a universal crossover for $α_η$ and $U_{ηR}$ from $θ$ to good solvents in line with earlier observations on synthetic polymer-solvent systems. The significant difference between the swelling of the dynamic viscosity radius from the observed swelling of the static radius of gyration, is shown to arise from the presence of hydrodynamic interactions in the non-draining limit. Simulated values of $α_η$ and $U_{ηR}$ are in good agreement with experimental measurements in synthetic polymer solutions reported previously, and with the measurements in linear DNA solutions reported here.

preprint2013arXiv

Universal solvent quality crossover of the zero shear rate viscosity of semidilute DNA solutions

The scaling behaviour of the zero shear rate viscosity of semidilute unentangled DNA solutions, in the double crossover regime driven by temperature and concentration, is mapped out by systematic experiments. The viscosity is shown to have a power law dependence on the scaled concentration $c/c^*$, with an effective exponent that depends on the solvent quality parameter $z$. The determination of the form of this universal crossover scaling function requires the estimation of the $θ$ temperature of dilute DNA solutions in the presence of excess salt, and the determination of the solvent quality parameter at any given molecular weight and temperature. The $θ$ temperature is determined to be $T_θ\approx 15^\circ$ C using static light scattering, and the solvent quality parameter has been determined by dynamic light scattering.

preprint2012arXiv

Dynamic crossover scaling in polymer solutions

The crossover region in the phase diagram of polymer solutions, in the regime above the overlap concentration, is explored by Brownian Dynamics simulations, to map out the universal crossover scaling functions for the gyration radius and the single-chain diffusion constant. Scaling considerations, our simulation results, and recently reported data on the polymer contribution to the viscosity obtained from rheological measurements on DNA systems, support the assumption that there are simple relations between these functions, such that they can be inferred from one another.

preprint2012arXiv

Fluid-Structure Interaction in Deformable Microchannels

A microfluidic device is constructed from PDMS with a single channel having a short section that is a thin flexible membrane, in order to investigate the complex fluid-structure interaction that arises between a flowing fluid and a deformable wall. Experimental measurements of membrane deformation and pressure drop are compared with predictions of two-dimensional and three-dimensional computational models which numerically solve the equations governing the elasticity of the membrane coupled with the equations of motion for the fluid. It is shown that the two-dimensional model, which assumes a finite thickness elastic beam that is infinitely wide, approximates reasonably well the three-dimensional model, and is in excellent agreement with experimental observations of the profile of the membrane, when the width of the membrane is beyond a critical thickness, determined to be roughly twice the length of the membrane.

preprint2012arXiv

Optimisation of a Brownian dynamics algorithm for semidilute polymer solutions

Simulating the static and dynamic properties of semidilute polymer solutions with Brownian dynamics (BD) requires the computation of a large system of polymer chains coupled to one another through excluded-volume and hydrodynamic interactions. In the presence of periodic boundary conditions, long-ranged hydrodynamic interactions are frequently summed with the Ewald summation technique. By performing detailed simulations that shed light on the influence of several tuning parameters involved both in the Ewald summation method, and in the efficient treatment of Brownian forces, we develop a BD algorithm in which the computational cost scales as O(N^{1.8}), where N is the number of monomers in the simulation box. We show that Beenakker&#39;s original implementation of the Ewald sum, which is only valid for systems without bead overlap, can be modified so that θ-solutions can be simulated by switching off excluded-volume interactions. A comparison of the predictions of the radius of gyration, the end-to-end vector, and the self-diffusion coefficient by BD, at a range of concentrations, with the hybrid Lattice Boltzmann/Molecular Dynamics (LB/MD) method shows excellent agreement between the two methods. In contrast to the situation for dilute solutions, the LB/MD method is shown to be significantly more computationally efficient than the current implementation of BD for simulating semidilute solutions. We argue however that further optimisations should be possible.

preprint2010arXiv

Collapse dynamics of copolymers in a poor solvent: Influence of hydrodynamic interactions and chain sequence

We investigate the dynamics of the collapse of a single copolymer chain, when the solvent quality is suddenly quenched from good to poor. We employ Brownian dynamics simulations of a bead-spring chain model and incorporate fluctuating hydrodynamic interactions via the Rotne-Prager-Yamakawa tensor. Various copolymer architectures are studied within the framework of a two-letter HP model, where monomers of type H (hydrophobic) attract each other, while all interactions involving P (polar or hydrophilic) monomers are purely repulsive. The hydrodynamic interactions are found to assist the collapse. Furthermore, the chain sequence has a strong influence on the kinetics and on the compactness and energy of the final state. The dynamics is typically characterised by initial rapid cluster formation, followed by coalescence and final rearrangement to form the compact globule. The coalescence stage takes most of the collapse time, and its duration is particularly sensitive to the details of the architecture. Long blocks of type P are identified as the main bottlenecks to find the globular state rapidly.

preprint2009arXiv

Implicit and explicit solvent models for the simulation of a single polymer chain in solution: Lattice Boltzmann vs Brownian dynamics

We present a comparative study of two computer simulation methods to obtain static and dynamic properties of dilute polymer solutions. The first approach is a recently established hybrid algorithm based upon dissipative coupling between Molecular Dynamics and lattice Boltzmann (LB), while the second is standard Brownian Dynamics (BD) with fluctuating hydrodynamic interactions. Applying these methods to the same physical system (a single polymer chain in a good solvent in thermal equilibrium) allows us to draw a detailed and quantitative comparison in terms of both accuracy and efficiency. It is found that the static conformations of the LB model are distorted when the box length L is too small compared to the chain size. Furthermore, some dynamic properties of the LB model are subject to an $L^{-1}$ finite size effect, while the BD model directly reproduces the asymptotic $L \to \infty$ behavior. Apart from these finite size effects, it is also found that in order to obtain the correct dynamic properties for the LB simulations, it is crucial to properly thermalize all the kinetic modes. Only in this case, the results are in excellent agreement with each other, as expected. Moreover, Brownian Dynamics is found to be much more efficient than lattice Boltzmann as long as the degree of polymerization is not excessively large.

preprint2007arXiv

Coil-Stretch Transition and the Break Down of Continuum Models

The breakdown of finite element (FEM) computations for the flow of an Oldroyd-B fluid around a cylinder confined between parallel plates, at Weissenberg numbers Wi = O(1), is shown to arise due to a coil-stretch transition experienced by polymer molecules traveling along the centerline in the wake of the cylinder. With increasing Wi, the coil-stretch transition leads to an unbounded growth in the stress maximum in the cylinder wake. Finite element computations for a FENE-P fluid reveal that, although polymer molecules undergo a coil-stretch transition in the cylinder wake, the mean extension of the molecules saturates to a value close to the fully extended length, leading to bounded stresses with increasing Wi. The existence of a coil-stretch transition has been deduced by examining the behavior of ultra-dilute Oldroyd-B and FENE-P fluids. In this case, the solution along the centerline in the cylinder wake can be obtained exactly since the velocity field is uncoupled from the stress and conformation tensor fields. Estimation of the number of finite elements required to achieve convergence reveals the in-feasibility of obtaining solutions for the Oldroyd-B model for Wi > 1.

preprint2006arXiv

Gaussian approximation for finitely extensible bead-spring chains with hydrodynamic interaction

The Gaussian Approximation, proposed originally by Ottinger [J. Chem. Phys., 90 (1) : 463-473, 1989] to account for the influence of fluctuations in hydrodynamic interactions in Rouse chains, is adapted here to derive a new mean-field approximation for the FENE spring force. This &#34;FENE-PG&#34; force law approximately accounts for spring-force fluctuations, which are neglected in the widely used FENE-P approximation. The Gaussian Approximation for hydrodynamic interactions is combined with the FENE-P and FENE-PG spring force approximations to obtain approximate models for finitely-extensible bead-spring chains with hydrodynamic interactions. The closed set of ODE&#39;s governing the evolution of the second-moments of the configurational probability distribution in the approximate models are used to generate predictions of rheological properties in steady and unsteady shear and uniaxial extensional flows, which are found to be in good agreement with the exact results obtained with Brownian dynamics simulations. In particular, predictions of coil-stretch hysteresis are in quantitative agreement with simulations&#39; results. Additional simplifying diagonalization-of-normal-modes assumptions are found to lead to considerable savings in computation time, without significant loss in accuracy.

preprint2005arXiv

The origin of dynamic scaling in dilute polymer solutions

The hydrodynamic radius of a polymer chain, obtained using Brownian dynamics simulations of the continuum Edwards model, is found to obey a crossover in the excluded volume parameter z, which is significantly different from that observed for the radius of gyration. It is shown that this difference arises from contributions due to dynamic correlations to the diffusivity, which are ignored in the commonly used definition of hydrodynamic radius based on the Kirkwood expression. The swelling of the hydrodynamic radius from the theta-state, obtained from simulations, shows remarkable agreement with experimental measurements.

preprint2004arXiv

Universal consequences of the presence of excluded volume interactions in dilute polymer solutions undergoing shear flow

The role of solvent quality in determining the universal material properties of dilute polymer solutions undergoing steady simple shear flow is examined. A bead-spring chain representation of the polymer molecule is used, and the influence of solvent molecules on polymer conformations is modelled by a narrow Gaussian excluded volume potential that acts pair-wise between the beads of the chain. Brownian dynamics simulations data, acquired for chains of finite length, and extrapolated to the limit of infinite chain length, are shown to be model independent. This feature of the narrow Gaussian potential, which leads to results identical to a $δ$-function repulsive potential, enables the prediction of both universal crossover scaling functions and asymptotic behavior in the excluded volume limit. Universal viscometric functions, obtained by this procedure, are found to exhibit increased shear thinning with increasing solvent quality. In the excluded volume limit, they are found to obey power law scaling with the characteristic shear rate $β$, in close agreement with previously obtained renormalization group results. The presence of excluded volume interactions is also shown to lead to a weakening of the alignment of the polymer chain with the flow direction.

preprint2003arXiv

Equilibrium swelling and universal ratios in dilute polymer solutions: Exact Brownian dynamics simulations for a delta function excluded volume potential

A narrow Gaussian excluded volume potential, which tends to a delta-function repulsive potential in the limit of a width parameter d* going to zero, has been used to examine the universal consequences of excluded volume interactions on the equilibrium and linear viscoelastic properties of dilute polymer solutions. Brownian dynamics simulations data, acquired for chains of finite length, has been extrapolated to the limit of infinite chain length to obtain model independent predictions. The success of the method in predicting well known aspects of static solution properties suggests that it can be used as a systematic means by which the influence of solvent quality on both equilibrium and non-equilibrium properties can be studied.

preprint1999arXiv

Analysis of radial segregation of granular mixtures in a rotating drum

This paper considers the segregation of a granular mixture in a rotating drum. Extending a recent kinematic model for grain transport on sandpile surfaces to the case of rotating drums, an analysis is presented for radial segregation in the rolling regime, where a thin layer is avalanching down while the rest of the material follows rigid body rotation. We argue that segregation is driven not just by differences in the angle of repose of the species, as has been assumed in earlier investigations, but also by differences in the size and surface properties of the grains. The cases of grains differing only in size (slightly or widely) and only in surface properties are considered, and the predictions are in qualitative agreement with observations. The model yields results inconsistent with the assumptions for more general cases, and we speculate on how this may be corrected.

preprint1998arXiv

Viscometric Functions for a Dilute Solution of Polymers in a Good Solvent

A dilute polymer solution is modeled as a suspension of non-interacting Hookean dumbbells and the effect of excluded volume is taken into account by incorporating a narrow Gaussian repulsive potential between the beads of each dumbbell. The narrow Gaussian potential is a means of regularising a delta-function potential---it tends to the delta-function potential in the limit of the width parameter going to zero. Exact predictions of viscometric functions in simple shear flow are obtained with the help of a retarded motion expansion and by Brownian dynamics simulations. It is shown that for relatively small non-zero values of the width parameter, the presence of excluded volume causes a swelling of the dumbbell at equilibrium, and shear thinning in simple shear flow. On the other hand, a delta function excluded volume potential does not lead to either swelling or to shear thinning. Approximate viscometric functions, obtained by assuming that the bead-connector vector is described by a Gaussian non-equilibrium distribution function, are found to be accurate above a threshold value of the width parameter, for a given value of the strength of excluded volume interaction. A first order perturbation expansion reveals that the Gaussian approximation is exact to first order in the strength of excluded volume interaction. The predictions of an alternative quadratic excluded volume potential suggested earlier by Fixman (J. Chem. Phys., 1966, 45, 785; 793) are also compared with those of the narrow Gaussian potential.