Conformational fluctuations of a tethered polymer in uniform flow

The decay of correlations in the conformational fluctuations of a tethered polymer subjected to a uniform flow is analyzed in terms of relaxation times and associated normal modes. These quantities are calculated numerically from Brownian dynamics simulations of several bead spring polymer models. In this way, the influence of different effects like a finite extensibility of the springs and excluded-volume as well as hydrodynamic interactions between the beads on the decay of fluctuations is identified. Moreover, by comparison of the simulation results to analytically tractable blob models with corresponding assumptions, the capability of the tensile-blob picture to predict relaxation times and modes is assessed. For excluded-volume and hydrodynamic interactions a crossover to Rouse-like behavior occurs when the flow velocity and hence the polymer deformation exceeds a certain value. For finitely extensible springs, in contrast, the relaxation times decrease monotonically with increasing polymer deformation. This latter behavior differs from assumptions often used in rheological modeling by dumbbells and is not captured by the blob model.Received: 4 April 2003, Published online: 12 August 2003PACS: 83.80.Rs Polymer solutions - 83.10.Mj Molecular dynamics, Brownian dynamics - 36.20.Ey Conformation (statistics and dynamics) - 47.50. + d Non-Newtonian fluid flows     

Crooks Fluctuation Theorem for Single Polymer Dynamics in Time-Dependent Flows: Understanding Viscoelastic Hysteresis

Nonequilibrium work relations have fundamentally advanced our understanding of molecular processes. In recent years, fluctuation theorems have been extensively applied to understand transitions between equilibrium steady-states, commonly described by simple control parameters such as molecular extension of a protein or polymer chain stretched by an external force in a quiescent fluid. Despite recent progress, far less is understood regarding the application of fluctuation theorems to processes involving nonequilibrium steady-states such as those described by polymer stretching dynamics in nonequilibrium fluid flows. In this work, we apply the Crooks fluctuation theorem to understand the nonequilibrium thermodynamics of dilute polymer solutions in flow. We directly determine the nonequilibrium free energy for single polymer molecules in flow using a combination of single molecule experiments and Brownian dynamics simulations. We further develop a time-dependent extensional flow protocol that allows for probing viscoelastic hysteresis over a wide range of flow strengths. Using this framework, we define quantities that uniquely characterize the coil-stretch transition for polymer chains in flow. Overall, generalized fluctuation theorems provide a powerful framework to understand polymer dynamics under far-from-equilibrium conditions.     

A new mixed subgrid-scale model for large eddy simulation of turbulent drag-reducing flows of viscoelastic fluids

A mixed subgrid-scale(SGS) model based on coherent structures and temporal approximate deconvolution(MCT) is proposed for turbulent drag-reducing flows of viscoelastic fluids. The main idea of the MCT SGS model is to perform spatial filtering for the momentum equation and temporal filtering for the conformation tensor transport equation of turbulent flow of viscoelastic fluid, respectively. The MCT model is suitable for large eddy simulation(LES) of turbulent dragreducing flows of viscoelastic fluids in engineering applications since the model parameters can be easily obtained. The LES of forced homogeneous isotropic turbulence(FHIT) with polymer additives and turbulent channel flow with surfactant additives based on MCT SGS model shows excellent agreements with direct numerical simulation(DNS) results. Compared with the LES results using the temporal approximate deconvolution model(TADM) for FHIT with polymer additives, this mixed SGS model MCT behaves better, regarding the enhancement of calculating parameters such as the Reynolds number.For scientific and engineering research, turbulent flows at high Reynolds numbers are expected, so the MCT model can be a more suitable model for the LES of turbulent drag-reducing flows of viscoelastic fluid with polymer or surfactant additives.     

Influence of the nematic order on the rheology and conformation of stretched comb-like liquid crystalline polymers

We have studied the rheology and the conformation of stretched comb-like liquid-crystalline polymers. Both the influence of the comb-like structure and the specific effect of the nematic interaction on the dynamics are investigated. For this purpose, two isomers of a comb-like polymetacrylate polymer, of well-defined molecular weights, were synthesized: one displays a nematic phase over a wide range of temperature, the other one has only an isotropic phase. Even with high degrees of polymerization N, between 40 and 1000, the polymer chains studied were not entangled. The stress-strain curves during the stretching and relaxation processes show differences between the isotropic and nematic comb-like polymers. They suggest that, in the nematic phase, the chain dynamics is more cooperative than for a usual linear polymer. Small-angle neutron scattering has been used in order to determine the evolution of the chain conformation after stretching, as a function of the duration of relaxation t _r. The conformation can be described with two parameters only: , the global deformation of the polymer chain, and p, the number of statistical units of locally relaxed sub-chains. For the comb-like polymer, the chain deformation is pseudo-affine: is always smaller than (the deformation ratio of the whole sample). In the isotropic phase, has a constant value, while pincreases as tr. This latter behavior is not that expected for non-entangled chains, in which p varies as t _r ^1/2 (Rouse model). In the nematic phase, decreases as a stretched exponential function of t _r, while p remains constant. The dynamics of the comb-like polymers is discussed in terms of living clusters from which junctions are produced by interactions between side chains. The nematic interaction increases the lifetime of these junctions and, strikingly, the relaxation is the same at all scales of the whole polymer chain. Received 5 May 1999 and Received in final form 18 October 1999     

Dynamics of a polymer chain confined in a membrane

We present a Brownian dynamics theory with full hydrodynamics (Stokesian dynamics) for a Gaussian polymer chain embedded in a liquid membrane which is surrounded by bulk solvent and walls. The mobility tensors are derived in Fourier space for the two geometries, namely, a free membrane embedded in a bulk fluid, and a membrane sandwiched by the two walls. Within the preaveraging approximation, a new expression for the diffusion coefficient of the polymer is obtained for the free-membrane geometry. We also carry out a Rouse normal mode analysis to obtain the relaxation time and the dynamical structure factor. For large polymer size, both quantities show Zimm-like behavior in the free-membrane case, whereas they are Rouse-like for the sandwiched membrane geometry. We use the scaling argument to discuss the effect of excluded-volume interactions on the polymer relaxation time.     

Molecular dynamics in thin grafted and spin-coated polymer layers

Broadband Dielectric Spectroscopy is employed to study the molecular dynamics in thin layers of grafted and spin-coated polydimethylsiloxane (PDMS) and of the type-A polymer cis-1,4-polyisoprene (PI). Strong resemblance with low-molecular-weight systems is found but additionally the conformation of the chain and its extension becomes important. For PI a confinement-induced relaxation is found. It is observable if the layer thickness is comparable to the end-to-end vector of the chain and it depends strongly on the separation between interfaces. In contrast, the dynamic (and hence the calorimetric) glass transition is not affected by the geometrical confinement.Received: 1 January 2003, Published online: 8 October 2003PACS: 64.70.Pf Glass transitions - 83.10.Mj Molecular dynamics, Brownian dynamics - 68.60.-p Physical properties of thin films, nonelectronic     

Study of hysteresis behavior in reactive sputtering of cylindrical magnetron plasma

In order to make sufficient use of reactive cylindrical magnetron plasma for depositing compound thin films, it is necessary to characterize the hysteresis behavior of the discharge. Cylindrical magnetron plasmas with different targets namely titanium and aluminium are studied in an argon/oxygen and an argon/nitrogen gas environment respectively. The aluminium and titanium emission lines are observed at different flows of reactive gases. The emission intensity is found to decrease with the increase of the reactive gas flow rate. The hysteresis behavior of reactive cylindrical magnetron plasma is studied by determining the variation of discharge voltage with increasing and then reducing the flow rate of reactive gas,while keeping the discharge current constant at 100 m A. Distinct hysteresis is found to be formed for the aluminium target and reactive gas oxygen. For aluminium/nitrogen, titanium/oxygen and titanium/nitrogen, there is also an indication of the formation of hysteresis; however, the characteristics of variation from metallic to reactive mode are different in different cases. The hysteresis behaviors are different for aluminium and titanium targets with the oxygen and nitrogen reactive gases, signifying the difference in reactivity between them. The effects of the argon flow rate and magnetic field on the hysteresis are studied and explained.     

Solitary coherent structures in viscoelastic shear flow: computation and mechanism

Starting from stationary bifurcations in Couette-Dean flow, we compute nontrivial stationary solutions in inertialess viscoelastic circular Couette flow. These solutions are strongly localized vortex pairs, exist at arbitrarily large wavelengths, and show hysteresis in the Weissenberg number, similar to experimentally observed "diwhirl" patterns. Based on the computed velocity and stress fields, we elucidate a heuristic, fully nonlinear mechanism for these flows. We propose that these localized, fully nonlinear structures comprise fundamental building blocks for complex spatiotemporal dynamics in the flow of elastic liquids.     

On the role of pressure in elasto-inertial turbulence

The dynamics of elasto-inertial turbulence is investigated numerically from the perspective of the coupling between polymer dynamics and flow structures. In particular, direct numerical simulations of channel flow with Reynolds numbers ranging from 1000 to 6000 are used to study the formation and dynamics of elastic instabilities and their effects on the flow. Based on the splitting of the pressure into inertial and polymeric contributions, it is shown that the polymeric pressure is a non-negligible component of the total pressure fluctuations, although the rapid inertial part dominates. Unlike Newtonian flows, the slow inertial part is almost negligible in elasto-inertial turbulence. Statistics on the different terms of the Reynolds stress transport equation also illustrate the energy transfers between polymers and turbulence and the redistributive role of pressure. Finally, the trains of cylindrical structures around sheets of high polymer extension that are characteristics of elasto-inertial turbulence are shown to be correlated with the polymeric pressure fluctuations.     

A higher-order macroscopic model for pedestrian flows

This paper develops a higher-order macroscopic model of pedestrian crowd dynamics derived from fluid dynamics that consists of two-dimensional Euler equations with relaxation. The desired directional motion of pedestrians is determined by an Eikonal-type equation, which describes a problem that minimizes the instantaneous total walking cost from origin to destination. A linear stability analysis of the model demonstrates its ability to describe traffic instability in crowd flows. The algorithm to solve the macroscopic model is composed of a splitting technique introduced to treat the relaxation terms, a second-order positivity-preserving central-upwind scheme for hyperbolic conservation laws, and a fast-sweeping method for the Eikonal-type equation on unstructured meshes. To test the applicability of the model, we study a challenging pedestrian crowd flow problem of the presence of an obstruction in a two-dimensional continuous walking facility. The numerical results indicate the rationality of the model and the effectiveness of the computational algorithm in predicting the flux or density distribution and the macroscopic behavior of the pedestrian crowd flow. The simulation results are compared with those obtained by the two-dimensional Lighthill-Whitham-Richards pedestrian flow model with various model parameters, which further shows that the macroscopic model is able to correctly describe complex phenomena such as “stop-and-go waves” observed in empirical pedestrian flows.     

Dynamics of hairpin vortices and polymer-induced turbulent drag reduction

It has been known for over six decades that the dissolution of minute amounts of high molecular weight polymers in wall-bounded turbulent flows results in a dramatic reduction in turbulent skin friction by up to 70%. First principles simulations of turbulent flow of model polymer solutions can predict the drag reduction (DR) phenomenon. However, the essential dynamical interactions between the coherent structures present in turbulent flows and polymer conformation field that lead to DR are poorly understood. We examine this connection via dynamical simulations that track the evolution of hairpin vortices, i.e., counter-rotating pairs of quasistreamwise vortices whose nonlinear autogeneration and growth, decay and breakup are centrally important to turbulence stress production. The results show that the autogeneration of new vortices is suppressed by the polymer stresses, thereby decreasing the turbulent drag.     

纳米尺度流动速度计算的新方法

采用流体力学中流量与流速的计算和分子动力学相结合的方法,模拟液态氩在纳米通道内的三维Poiseuille流动和驱动方腔流动,计算流体流速.结果表明:平板形纳米通道内,该方法求得的流速与传统分子动力学方法所求流速基本吻合,可以用该方法计算不同壁面情况下的流速;对于纳米方腔通道内的流体,在不同模型下两种方法计算出的流速分布大致相同,但是其边界速度差别很大,在边界的速度计算方面新方法的精确度更高,收敛速度比传统方法快.     

A Lattice Boltzmann Kinetic Model for Microflow and Heat Transfer

In this paper, we propose a lattice Boltzmann BGK model for simulation of micro flows with heat transfer based on kinetic theory and the thermal lattice Boltzmann method (He et al., J. Comp. Phys. 146:282, 1998). The relaxation times are redefined in terms of the Knudsen number and a diffuse scattering boundary condition (DSBC) is adopted to consider the velocity slip and temperature jump at wall boundaries. To check validity and potential of the present model in modelling the micro flows, two two-dimensional micro flows including thermal Couette flow and thermal developing channel flow are simulated and numerical results obtained compare well with previous studies of the direct simulation Monte Carlo (DSMC), molecular dynamics (MD) approaches and the Maxwell theoretical analysis     

Different pathways in mechanical unfolding/folding cycle of a single semiflexible polymer

Kinetics of conformational change of a semiflexible polymer under mechanical external field were investigated with Langevin dynamics simulations. It is found that a semiflexible polymer exhibits large hysteresis in mechanical folding/unfolding cycle even with a slow operation, whereas in a flexible polymer, the hysteresis almost disappears at a sufficiently slow operation. This suggests that the essential features of the structural transition of a semiflexible polymer should be interpreted at least on a two-dimensional phase space. The appearance of such large hysteresis is discussed in relation to different pathways in the loading and unloading processes. By using a minimal two-variable model, the hysteresis loop is described in terms of different pathways on the transition between two stable states.     

Cyclic motion of a grafted polymer under shear flow

The long-time dynamics of a single end-tethered chain under shear flow are studied using molecular and Brownian dynamics simulations of a flexible polymer. As observed in previous experiments with tethered DNA [Phys. Rev. Lett. 84, 4769 (2000)], under a flow sheared at constant rate the chain performs a cyclic motion. But, contrary to what has been previously suggested, a well-defined characteristic period exists and it is clearly revealed in the cross spectra of the chain extension along flow and gradient directions. The main cycling time scales like the time needed to stretch the polymer by convection, being about 10 times the relaxation time of the chain in flow. This coherent recursive motion introduces long memory in the fluid and suggests resonance effects under periodic external forcing.     

Large-eddy simulations of a forced homogeneous isotropic turbulence with polymer additives

Large-eddy simulations （LES） based on the temporal approximate deconvolution model were performed for a forced homogeneous isotropic turbulence （FHIT） with polymer additives at moderate Taylor Reynolds number. Finitely extensible nonlinear elastic in the Peterlin approximation model was adopted as the constitutive equation for the filtered conformation tensor of the polymer molecules. The LES results were verified through comparisons with the direct numerical simulation results. Using the LES database of the FHIT in the Newtonian fluid and the polymer solution flows, the polymer effects on some important parameters such as strain, vorticity, drag reduction, and so forth were studied. By extracting the vortex structures and exploring the flatness factor through a high-order correlation function of velocity derivative and wavelet analysis, it can be found that the small-scale vortex structures and small-scale intermittency in the FHIT are all inhibited due to the existence of the polymers. The extended self-similarity scaling law in the polymer solution flow shows no apparent difference from that in the Newtonian fluid flow at the currently simulated ranges of Reynolds and Weissenberg numbers.     

Polymer dynamics in chaotic flows with a strong shear component

We consider the dynamics of a polymer molecule injected into a chaotic flow with a strong mean shear component. The polymer experiences aperiodic tumbling in such flows. We consider a simplified model of the chaotic velocity field given by the superposition of a steady shear flow and a large-scale isotropic short-correlated random component. In the framework of this model, we present a detailed study of the statistical properties of single-polymer dynamics. We obtain the stationary probability distribution function of the polymer orientation, find the distribution of time periods between consequent events of tumbling, and find the tails of the polymer size distribution. The text was submitted by the author in English.     

The dynamics of copolymers in homopolymer matrices

The ability of a polymer chain to relax when it is deformed, and the ease to which this relaxation occurs, dramatically influences the macroscopic properties of the polymeric material. However, the local segmental relaxation processes in multi-component polymer systems are not well understood. The impact of the dynamics of one component on the dynamics of the other is not simply proportional to the relative amounts of each component, as one might expect. This paper discusses recent results using neutron techniques and Monte Carlo simulation that monitor the dynamic properties of a copolymer in a homopolymer matrix. In particular, the results indicate that altering either copolymer sequence distribution or copolymer composition will dramatically impact the dynamics of the copolymer that is surrounded by homopolymers. These results provide important fundamental information on the coupling of the dynamics of two components in a multi-component polymer system. This data also offer insight into the local segmental relaxation processes in multi-component polymer systems, which are not well understood and yet influence the ultimate properties of these mixtures     

Modeling viscoelastic flow with discrete methods

The hydrodynamics of viscoelastic materials (for example, polymer melts and solutions) presents interesting and complex phenomena, for example, instabilities and turbulent flow at very low Reynolds numbers due to normal stress effects and the existence of a finite stress relaxation time. This present work is motivated by renewed interest in instabilities in polymer flow. The majority of currently used numerical methods discretize a constitutive equation on a grid with finite difference or similar methods. We present work in progress in which we simulate viscoelastic flow with dissipative particle dynamics. The advantage of this approach is that many of the numerical instabilities of conventional methods can be avoided, and that the model gives clear physical insight into the origins of many viscoelastic flow instabilities.