We provide a review of our recent work on the development of a multi-scale simulation methodology to calculate the rheology and flow of wormlike micelles. There is a great need for understanding the link between the detailed chemistry of surfactants, forming wormlike micelles, and their macroscopic rheological properties. We show how this link may be explored through particle simulations. First, we calculate the mechanical properties of small units of wormlike micelles from atomistic molecular dynamics simulations. These mechanical properties are subsequently used in a coarse-grained Brownian dynamics model, where the persistence length is the smallest length scale. We show that the non-Newtonian rheology of wormlike micellar systems can be accurately determined from our simple mesoscopic simulation model. Finally, we show that this mesoscopic model can be used to study the flow of wormlike micelles in contraction-expansion geometries.

We use molecular dynamics simulations to study phase separation of a 50:50 (by volume) fluid mixture in a confined and curved (Taylor-Couette) geometry, consisting of two concentric cylinders. The inner cylinder may be rotated to achieve a shear flow. In nonsheared systems we observe that, for all cases under consideration, the final equilibrium state has a stacked structure. Depending on the lowest free energy in the geometry the stack may be either flat, with its normal in the z direction, or curved, with its normal in the r or θ direction. In sheared systems we make several observations. First, when starting from a prearranged stacked structure, we find that sheared gradient and vorticity stacks retain their character for the durations of the simulation, even when another configuration is preferred (as found when starting from a randomly mixed configuration). This slow transition to another configuration is attributed to a large free energy barrier between the two states. In case of stacks with a normal in the gradient direction, we find interesting interfacial waves moving with a prescribed angular velocity in the flow direction. Because such a wave is not observed in simulations with a flat geometry at similar shear rates, the curvature of the wall is an essential ingredient of this phenomenon. Second, when starting from a randomly mixed configuration, stacks are also observed, with an orientation that depends on the applied shear rate. Such transitions to other orientations are similar to observations in microphase separated diblock copolymer melts. At higher shear rates complex patterns emerge, accompanied by deviations from a homogeneous flow profile. The transition from steady stacks to complex patterns takes place around a shear rate 1/ τdv, where τdv is the crossover time from diffusive to viscous dominated growth of phase-separated domains, as measured in equilibrium simulations.

Recently there has been a great deal of attention, from researchers both in academia and in industry, focused on the rheological properties of solutions of viscoelastic wormlike micelles formed by surfactants. It is particularly vital to understand the properties of these solutions with regard to their flow in porous media, given their application to the recovery of hydrocarbons from subterranean formations. In this study a realistic mesoscopic Brownian dynamics model has been utilized to investigate the flow of viscoelastic surfactant (VES) fluid through individual pores with sizes of around one micron. In particular the influence of micelle size, pore geometry and flow rate on the ability of worms to pass through the pores was studied. The ways in which these parameters influence the conformational properties of the worms and the spatial distribution of micelles inside the simulation cell was also investigated. Despite the observation that the density and length distributions became non-uniform at higher scission energy, the distribution of breaking and fusion events remained spatially uniform.

A coarse-grained simulation model was used to investigate the transient and steady-state nonlinear flow properties of unentangled and moderately entangled polyethylene melts. It was found that the model's predictions agree very well with various rheological experiments for shearing flow. This enabled to test the stress-optical rule and it was observed that it was valid, but only for low and moderate shear rates.

Molecular dynamics (MD) simulations of a polyethylene melt were performed to investigate the validity of the Rouse model predictions for a comprehensive set of correlation functions. It was found that the chains behave like Rouse theory predicts, but only on length scales larger than the segment length.

A soft-sphere discrete particle model was used to simulate mixing behavior of solid substrate particles in a slow rotating drum for solid-state fermentation. In this approach, forces acting on and subsequent motion of individual particles can be predicted. The (2D) simulations were qualitatively and quantitatively validated by mixing experiments using video and image analysis techniques. It was found that the simulations successfully predicted the mixing progress as a function of the degree of filling and size of the drum. It is shown that only relatively large, straight baffles perpendicular to the drum wall (67% of the drum radius) increase the mixing performance of the rotating drum. Considering the different aspects of mixing dealt with in this work, it is concluded that the soft sphere discrete particle model can serve as a valuable tool for investigating mixing of solid substrate particles. Finally, it is expected that this model may evolve into a potential tool for design and scale-up of mixed solid-state fermenters.