• Bulletin of the Korean Chemical Society
• /
• v.21 no.9
• /
• pp.875-881
• /
• 2000

We study the rheological properties of concentrated polymer solution and the melt under simple shear and elon-gational flow using Brownian dynamicssimulation. In order to describe the anisotropic molecular motion, we modifiedthe Giesekus' mobility tensor by incorporating the finitely extensible non-linear elastic (FENE) spring force into dumbbell model. To elucidate the nature of this model, our simulation results are compared with the data of FENE-P ("P"standsfor the Perterin) dumbbell model and experiments. While in steady state both original FENE and FENE-P models exhibit a similar viscosity response,the growthof viscosity becomes dissimilar as the anisotropy decreases and the flowrate increases. The steady state viscosity obtained from the simulation well describes the experiments including the shear-thinning behavior in shear flow and viscosity-thinning behavior in elongational flow. But the growth of viscosity oforiginal FENE dumbbell model cannot describe the experimental results in both flow fields.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.8
• /
• pp.799-807
• /
• 2010

Direct numerical simulations (DNS) of turbulent channel flow for which the drag is reduced by using polymer additives have been performed by a pseudo-spectral method. The Reynolds number based on the friction velocity and half-channel height is 395, and the polymeric stresses due to the polymer additives are evaluated using the FENE-P (finitely extensible nonlinear elastic-Peterlin) model. The numerical results show that the drag reduction rate is significantly affected by the parameters used in the FENE-P model, such as the maximum extensibility and relaxation time of the polymer molecules. The turbulence data for both low- and high-drag reduction regimes are analyzed. In addition, the effects of FENE-P model parameters on the flow characteristics have been investigated for the same drag reduction rate due to the polymer additives. Finally, the present DNS results have been used to verify the correlation between rheological parameters and the extent of drag reduction, which was suggested by Li et al. (2006).

• Korea-Australia Rheology Journal
• /
• v.14 no.4
• /
• pp.161-174
• /
• 2002

A new technique for numerical calculation of viscoelastic flow based on the combination of Neural Net-works (NN) and Brownian Dynamics simulation or Stochastic Simulation Technique (SST) is presented in this paper. This method uses a "universal approximator" based on neural network methodology in combination with the kinetic theory of polymeric liquid in which the stress is computed from the molecular configuration rather than from closed form constitutive equations. Thus the new method obviates not only the need for a rheological constitutive equation to describe the fluid (as in the original Calculation Of Non-Newtonian Flows: Finite Elements St Stochastic Simulation Techniques (CONNFFESSIT) idea) but also any kind of finite element-type discretisation of the domain and its boundary for numerical solution of the governing PDE's. As an illustration of the method, the time development of the planar Couette flow is studied for two molecular kinetic models with finite extensibility, namely the Finitely Extensible Nonlinear Elastic (FENE) and FENE-Peterlin (FENE-P) models.P) models.

• Korea-Australia Rheology Journal
• /
• v.17 no.3
• /
• pp.131-140
• /
• 2005

In this work we show the results of our most recent Direct Numerical Simulations (DNS) of turbulent viscoelastic channel flow using spectral spatial approximations and a stabilizing artificial diffusion in the viscoelastic constitutive model. The Finite-Elasticity Non-Linear Elastic Dumbbell model with the Peterlin approximation (FENE-P) is used to represent the effect of polymer molecules in solution, The corresponding rheological parameters are chosen so that to get closer to the conditions corresponding to maximum drag reduction: A high extensibility parameter (60) and a moderate solvent viscosity ratio (0.8) are used with two different friction Weissenberg numbers (50 and 100). We then first find that the corresponding achieved drag reduction, in the range of friction Reynolds numbers used in this work (180-590), is insensitive to the Reynolds number (in accordance to previous work). The obtained drag reduction is at the level of $49\%\;and\;63\%$, for the friction Weissenberg numbers 50 and 100, respectively. The largest value is substantially higher than any of our previous simulations, performed at more moderate levels of viscoelasticity (i.e. higher viscosity ratio and smaller extensibility parameter values). Therefore, the maximum extensional viscosity exhibited by the modeled system and the friction Weissenberg number can still be considered as the dominant factors determining the levels of drag reduction. These can reach high values, even for of dilute polymer solution (the system modeled by the FENE-P model), provided the flow viscoelasticity is high, corresponding to a high polymer molecular weight (which translates to a high extensibility parameter) and a high friction Weissenberg number. Based on that and the changes observed in the turbulent structure and in the most prevalent statistics, as presented in this work, we can still rationalize for an increasing extensional resistance-based drag reduction mechanism as the most prevalent mechanism for drag reduction, the same one evidenced in our previous work: As the polymer elasticity increases, so does the resistance offered to extensional deformation. That, in turn, changes the structure of the most energy-containing turbulent eddies (they become wider, more well correlated, and weaker in intensity) so that they become less efficient in transferring momentum, thus leading to drag reduction. Such a continuum, rheology-based, mechanism has first been proposed in the early 70s independently by Metzner and Lamley and is to be contrasted against any molecularly based explanations.