Diffusion of linear polymer melts in shear and extensional flows

Micro-structure and tensile properties of microfluidic spinning konjac glucomannan and sodium alginate composite bio-fibers regulated by shear and elongational flow: experiment and multi-scale simulation

Microfluidic spinning has been widely used to produce bio-fibers with excellent tensile performances by regulating the conformation of biological macromolecules. However, the effect of channel shapes on fiber tensile performances is unclear. In this study, bio-fibers were prepared using konjac glucomannan and sodium alginate by five channels. The micro-morphology and tensile performance of fibers were characterized and measured. Then, the dynamical behaviours of macromolecule clusters in flow fields were simulated by multi-scale numerical methods. The results show that the elongational flow with increasing extension rates produced fibers with a tensile strength of 32.34 MPa and a tensile strain of 18.72 %, which were 1.37 and 1.55 times that for a shear flow, respectively. The difference in tensile performances was attributed to the micro-morphology regulated by flow fields. The continuously increasing extension rate of flow was more effective than the shear rate or the maximum extension rate for the stretching of macromolecule clusters. We conclude that the channel shapes significantly influence flow fields, dynamical behaviours of molecule clusters, the morphology of fibers, and tensile performances. This study provides a novel numerical method and understanding of microfluidic spinning, which will promote the optimization and applications of bio-fibers.

Fingerprinting Brownian Motions of Polymers under Flow

We present a quantitative approach to the self-dynamics of polymers under steady flow by employing a set of complementary reference frames and extending the spherical harmonic expansion technique to dynamic density correlations. Application of this method to nonequilibrium molecular dynamics simulations of polymer melts reveals a number of universal features. For both unentangled and entangled melts, the center-of-mass motions in the flow frame are described by superdiffusive, anisotropic Gaussian distributions, whereas the isotropic component of monomer self-dynamics in the center-of-mass frame is strongly suppressed. Spatial correlation analysis shows that the heterogeneity of monomer self-dynamics increases significantly under flow.

An anisotropic immerse precipitation process for the preparation of polymer membranes

We study the immerse precipitation process in a ternary polymer/solvent/non-solvent system by numerically solving the two-dimensional Cahn-Hilliard phase field equation. In particular, we introduce anisotropic mobility, namely the mobility of a polymer varies over different spatial directions, and focus on the porosity morphology of the obtained polymer membrane. Simulations reveal that as the anisotropy increases in the polymer mobility, the polymer pattern changes from nearly isotropic and random voids to strips parallel to the direction with smaller mobility. The influence of anisotropy quickly saturates. The anisotropic mobility model is also applied to a ternary system mimicking the preparation of a hollow fiber membrane, and shows strong effects on the membrane porosity pattern.

Prediction and validation of diffusion coefficients in a model drug delivery system using microsecond atomistic molecular dynamics simulation and vapour sorption analysis

Diffusion of small to medium sized molecules in polymeric medical device materials underlies a broad range of public health concerns related to unintended leaching from or uptake into implantable medical devices. However, obtaining accurate diffusion coefficients for such systems at physiological temperature represents a formidable challenge, both experimentally and computationally. While molecular dynamics simulation has been used to accurately predict the diffusion coefficients, D, of a handful of gases in various polymers, this success has not been extended to molecules larger than gases, e.g., condensable vapours, liquids, and drugs. We present atomistic molecular dynamics simulation predictions of diffusion in a model drug eluting system that represent a dramatic improvement in accuracy compared to previous simulation predictions for comparable systems. We find that, for simulations of insufficient duration, sub-diffusive dynamics can lead to dramatic over-prediction of D. We present useful metrics for monitoring the extent of sub-diffusive dynamics and explore how these metrics correlate to error in D. We also identify a relationship between diffusion and fast dynamics in our system, which may serve as a means to more rapidly predict diffusion in slowly diffusing systems. Our work provides important precedent and essential insights for utilizing atomistic molecular dynamics simulations to predict diffusion coefficients of small to medium sized molecules in condensed soft matter systems.

Modeling steady-state dynamics of macromolecules in exponential-stretching flow using multiscale molecular-dynamics-multiparticle-collision simulations

We introduce a method to obtain steady-state uniaxial exponential-stretching flow of a fluid (akin to extensional flow) in the incompressible limit, which enables us to study the response of suspended macromolecules to the flow by computer simulations. The flow field in this flow is defined by v(x) = εx, where v(x) is the velocity of the fluid and ε is the stretch flow gradient. To eliminate the effect of confining boundaries, we produce the flow in a channel of uniform square cross section with periodic boundary conditions in directions perpendicular to the flow, but simultaneously maintain uniform density of fluid along the length of the tube. In experiments a perfect elongational flow is obtained only along the axis of symmetry in a four-roll geometry or a filament-stretching rheometer. We can reproduce flow conditions very similar to extensional flow near the axis of symmetry by exponential-stretching flow; we do this by adding the right amounts of fluid along the length of the flow in our simulations. The fluid particles added along the length of the tube are the same fluid particles which exit the channel due to the flow; thus mass conservation is maintained in our model by default. We also suggest a scheme for possible realization of exponential-stretching flow in experiments. To establish our method as a useful tool to study various soft matter systems in extensional flow, we embed (i) spherical colloids with excluded volume interactions (modeled by the Weeks-Chandler potential) as well as (ii) a bead-spring model of star polymers in the fluid to study their responses to the exponential-stretched flow and show that the responses of macromolecules in the two flows are very similar. We demonstrate that the variation of number density of the suspended colloids along the direction of flow is in tune with our expectations. We also conclude from our study of the deformation of star polymers with different numbers of arms f that the critical flow gradient ε(c) at which the star undergoes the coil-to-stretch transition is independent of f for f = 2,5,10, and 20.

Generalized extended Navier-Stokes theory: correlations in molecular fluids with intrinsic angular momentum

The extended Navier-Stokes theory accounts for the coupling between the translational and rotational molecular degrees of freedom. In this paper, we generalize this theory to non-zero frequencies and wavevectors, which enables a new study of spatio-temporal correlation phenomena present in molecular fluids. To discuss these phenomena in detail, molecular dynamics simulations of molecular chlorine are performed for three different state points. In general, the theory captures the behavior for small wavevector and frequencies as expected. For example, in the hydrodynamic regime and for molecular fluids with small moment of inertia like chlorine, the theory predicts that the longitudinal and transverse intrinsic angular velocity correlation functions are almost identical, which is also seen in the molecular dynamics simulations. However, the theory fails at large wavevector and frequencies. To account for the correlations at these scales, we derive a phenomenological expression for the frequency dependent rotational viscosity and wavevector and frequency dependent longitudinal spin viscosity. From this we observe a significant coupling enhancement between the molecular angular velocity and translational velocity for large frequencies in the gas phase; this is not observed for the supercritical fluid and liquid state points.

Transient-time correlation function applied to mixed shear and elongational flows

The transient-time correlation function (TTCF) method is used to calculate the nonlinear response of a homogeneous atomic fluid close to equilibrium. The TTCF response of the pressure tensor subjected to a time-independent planar mixed flow of shear and elongation is compared to directly averaged non-equilibrium molecular dynamics (NEMD) simulations. We discuss the consequence of noise in simulations with a small rate of deformation. The generalized viscosity for planar mixed flow is also calculated with TTCF. We find that for small rates of deformation, TTCF is far more efficient than direct averages of NEMD simulations. Therefore, TTCF can be applied to fluids with deformation rates which are much smaller than those commonly used in NEMD simulations. Ultimately, TTCF applied to molecular systems is amenable to direct comparison between NEMD simulations and experiments and so in principle can be used to study the rheology of polymer melts in industrial processes.

Anisotropic mobility model for polymers under shear and its linear response functions

We propose a simple dynamic model of polymers under shear with an anisotropic mobility tensor. We calculate the shear viscosity, the rheo-dielectric response function, and the parallel relaxation modulus under shear flow deduced from our model. We utilize recently developed linear response theories for nonequilibrium systems to calculate linear response functions. Our results are qualitatively consistent with experimental results. We show that our anisotropic mobility model can reproduce essential dynamical nature of polymers under shear qualitatively. We compare our model with other models or theories such as the convective constraint release model or nonequilibrium linear response theories.

Viscosity kernel of molecular fluids: butane and polymer melts

The wave-vector dependent shear viscosities for butane and freely jointed chains have been determined. The transverse momentum density and stress autocorrelation functions have been determined by equilibrium molecular dynamics in both atomic and molecular hydrodynamic formalisms. The density, temperature, and chain length dependencies of the reciprocal and real-space viscosity kernels are presented. We find that the density has a major effect on the shape of the kernel. The temperature range and chain lengths considered here have by contrast less impact on the overall normalized shape. Functional forms that fit the wave-vector-dependent kernel data over a large density and wave-vector range have also been tested. Finally, a structural normalization of the kernels in physical space is considered. Overall, the real-space viscosity kernel has a width of roughly 3-6 atomic diameters, which means that generalized hydrodynamics must be applied in predicting the flow properties of molecular fluids on length scales where the strain rate varies sufficiently in the order of these dimensions (e.g., nanofluidic flows).