Development of an integral viscoelastic flow solver in OpenFOAM®

The usual high cost of commercial codes, and some technical limitations, clearly limits the employment of numerical modelling tools in both industry and academia. Consequently, the number of companies that use numerical code is limited and there a lot of effort put on the development and maintenance of in-house academic based codes. Having in mind the potential of using numerical modelling tools as a design aid, of both products and processes, different research teams have been contributing to the development of open source codes/libraries. In this framework, any individual can take advantage of the available code capabilities and/or implement additional features based on his specific needs. These type of codes are usually developed by large communities, which provide improvements and new features in their specific fields of research, thus increasing significantly the code development process. Among others, OpenFOAM® multi-physics computational library, developed by a very large and dynamic community, nowadays comprises several features usually only available in their commercial counterparts; e.g. dynamic meshes, large diversity of complex physical models, parallelization, multiphase models, to name just a few. This computational library is developed in C++ and makes use of most of all language capabilities to facilitate the implementation of new functionalities. Concerning the field of computational rheology, OpenFOAM® solvers were recently developed to deal with the most relevant differential viscoelastic rheological models, and stabilization techniques are currently being verified. This work describes the implementation of a new solver in OpenFOAM® library, able to cope with integral viscoelastic models based on the deformation field method. The implemented solver is verified through the comparison of the predicted results with analytical solutions, results published in the literature and by using the Method of Manufactured Solutions.

A new integral viscoelastic flow solver in OpenFOAM®

The usual high cost of commercial codes, and some technical limitations, clearly limits the employment of numerical modelling tools in both industry and academia. Consequently, the number of companies that use numerical code is limited and there a lot of effort put on the development and maintenance of in-house academic based codes . Having in mind the potential of using numerical modelling tools as a design aid, of both products and processes, different research teams have been contributing to the development of open source codes/libraries. In this framework, any individual can take advantage of the available code capabilities and/or implement additional features based on his specific needs. These type of codes are usually developed by large communities, which provide improvements and new features in their specific fields of research, thus increasing significantly the code development process. Among others, OpenFOAM® multi-physics computational library, developed by a very large and dynamic community, nowadays comprises several features usually only available in their commercial counterparts; e.g. dynamic meshes, large diversity of complex physical models, parallelization, multiphase models, to name just a few. This computational library is developed in C++ and makes use of most of all language capabilities to facilitate the implementation of new functionalities. Concerning the field of computational rheology, OpenFOAM® solvers were recently developed to deal with the most relevant differential viscoelastic rheological models, and stabilization techniques are currently being verified. This work describes the implementation of a new solver in OpenFOAM® library, able to cope with integral viscoelastic models based on the deformation field method. The implemented solver is verified through the comparison of the predicted results with analytical solutions, results published in the literature and by using the Method of Manufactured Solutions

Analytical and numerical study of the electro-osmotic annular flow of viscoelastic fluids

In this work we present semi-analytical solutions for the electro-osmotic annular flow of viscoelastic fluids modeled by the Linear and Exponential PTT models. The viscoelastic fluid flows in the axial direction between two concentric cylinders under the combined influences of electrokinetic and pressure forcings. The analysis invokes the Debye-Hückel approximation and includes the limit case of pure electro-osmotic flow. The solution is valid for both no slip and slip velocity at the walls and the chosen slip boundary condition is the linear Navier slip velocity model. The combined effects of fluid rheology, electro-osmotic and pressure gradient forcings on the fluid velocity distribution are also discussed.

Annular flow of viscoelastic fluids: Analytical and numerical solutions

This work provides analytical and numerical solutions for the linear, quadratic and exponential Phan–Thien–Tanner (PTT) viscoelastic models, for axial and helical annular fully-developed flows under no slip and slip boundary conditions, the latter given by the linear and nonlinear Navier slip laws. The rheology of the three PTT model functions is discussed together with the influence of the slip velocity upon the flow velocity and stress fields. For the linear PTT model, full analytical solutions for the inverse problem (unknown velocity) are devised for the linear Navier slip law and two different slip exponents. For the linear PTT model with other values of the slip exponent and for the quadratic PTT model, the polynomial equation for the radial location (β) of the null shear stress must be solved numerically. For both models, the solution of the direct problem is given by an iterative procedure involving three nonlinear equations, one for β, other for the pressure gradient and another for the torque per unit length. For the exponential PTT model we devise a numerical procedure that can easily compute the numerical solution of the pure axial flow problem

Slip flows of Newtonian and viscoelastic fluids in a 4:1 contraction

This work presents a numerical study of the 4:1 planar contraction flow of a viscoelastic fluid described by the simplified Phan-Thien–Tanner model under the influence of slip boundary conditions at the channel walls. The linear Navier slip law was considered with the dimensionless slip coefficient varying in the range ½0; 4500. The simulations were carried out for a small constant Reynolds number of 0.04 and Deborah numbers (De) varying between 0 and 5. Convergence could not be achieved for higher values of the Deborah number, especially for large values of the slip coefficient, due to the large stress gradients near the singularity of the reentrant corner. Increasing the slip coefficient leads to the formation of two vortices, a corner and a lip vortex. The lip vortex grows with increasing slip until it absorbs the corner vortex, creating a single large vortex that continues to increase in size and intensity. In the range De = 3–5 no lip vortex was formed. The flow is characterized in detail for De ¼ 1 as function of the slip coefficient, while for the remaining De only the main features are shown for specific values of the slip coefficient.

Weld lines in extrusion: Understanding the role of the flow conditions

This work reports the implemen tation and verification of a new so lver in OpenFOAM® open source computational library, able to cope w ith integral viscoelastic models based on the integral upper-convected Maxwell model. The code is verified through the comparison of its predictions with anal ytical solutions and numerical results obtained with the differential upper-convected Maxwell model

Dispersion of fillers in viscoelastic polymer matrices using Eulerian-Lagrangian approaches

Understanding the behavior of c omplex composite materials using mixing procedures is fundamental in several industrial processes. For instance, polymer composites are usually manufactured using dispersion of fillers in polymer melt matrices. The success of the filler dispersion depends both on the complex flow patterns generated and on the polymer melt rheological behavior. Consequently, the availability of a numerical tool that allow to model both fluid and particle would be very useful to increase the process insight. Nowadays there ar e computational tools that allow modeling the behavior of filled systems, taking into account both the behavior of the fluid (Computational Rheology) and the particles (Discrete Element Method). One example is the DPMFoam solver of the OpenFOAM ® framework where the averaged volume fraction momentum and mass conservation equations are used to describe the fluid (continuous phase) rheology, and the Newton’s second law of motion is used to compute the particles (discrete phase) movement. In this work the refer red solver is extended to take into account the elasticity of the polymer melts for the continuous phase. The solver capabilities will be illustrated by studying the effect of the fluid rheology on the filler dispersion, taking into account different fluid types (generalized Newtonian or viscoelastic) and particles volume fraction and size. The results obtained are used to evaluate the relevance of considering the fluid complex rheology for the prediction of the composites morphology

Innovative method for the continuous monitoring of concrete viscoelastic properties since early ages: concept and pilot experiments

The experimental evaluation of viscoelastic properties of concrete is traditionally made upon creep tests that consist in the application of sustained loads either in compression or in tension. This kind of testing demands for specially devised rigs and requires careful monitoring of the evolution of strains, whereas assuring proper load constancy. The characterization of creep behaviour at early ages offers additional challenges due to the strong variations in viscoelastic behaviour of concrete during such stages, demanding for several testing ages to be assessed. The present research work aims to assist in reducing efforts for continuous assessment of viscoelastic properties of concrete at early ages, by application of a dynamic testing technique inspired in methodologies used in polymer science: Dynamic Mechanical Analyses. This paper briefly explains the principles of the proposed methodology and exhibits the first results obtained in a pilot application. The results are promising enough to encourage further developments.

Dispersion and re-agglomeration of graphite nanoplates in polypropylene melts under controlled flow conditions

The kinetics of GnP dispersion in polypropylene melt was studied using a prototype small scale modular extensional mixer. Its modular nature enabled the sequential application of a mixing step, melt relaxation, and a second mixing step. The latter could reproduce the flow conditions on the first mixing step, or generate milder flow conditions. The effect of these sequences of flow constraints upon GnP dispersion along the mixer length was studied for composites with 2 and 10 wt.% GnP. The samples collected along the first mixing zone showed a gradual decrease of number and size of GnP agglomerates, at a rate that was independent of the flow conditions imposed to the melt, but dependent on composition. The relaxation zone induced GnP re-agglomeration, and the application of a second mixing step caused variable dispersion results that were largely dependent on the hydrodynamic stresses generated.

Implementation of integral viscoelastic constitutive models in OpenFOAM® computational library

This work reports the implementation and verification of a new so lver in OpenFOAM® open source computational library, able to cope with integral viscoelastic models based on the integral upper-convected Maxwell model. The code is verified through the comparison of its predictions with analytical solutions and numerical results obtained with the differential upper-convected Maxwell model

Modelling integral viscoelastic flows in OpenFOAM

[Excerpt] A large number of constitutive equations were developed for viscoelastic fluids, some empirical and other with strong physical foundations. The currently available macroscopic constitutive equations can be divided in two main types: differential and integral. Some of the constitutive equations, e.g. Maxwell are available both in differential and integral types. However, relevant in tegral models, like K - BKZ, just possesses the integral form. (...)

Understanding the dispersion of fillers in polymer matrices

Understanding the mixing process of complex composite materials is fundamental in several industrial processes. For instance, the dispersion of fillers in polymer melt matrices is commonly employed to manufacture polymer composites, using a twin-screw extruder. The effectiveness of the filler dispersion depends not only on the complex flow patterns generated, but also on the polymer melt rheological behavior. Therefore, the availability of a numerical tool able to predict mixing, taking into account both fluid and particles phases would be very useful to increase the process insight, and thus provide useful guidelines for its optimization. In this work, a new Eulerian-Lagrangian numerical solver is developed OpenFOAM® computational library, and used to better understand the mechanisms determining the dispersion of fillers in polymer matrices. Particular attention will be given to the effect of the rheological model used to represent the fluid behavior, on the level of dispersion obtained. For the Eulerian phase the averaged volume fraction governing equations (conservation of mass and linear momentum) are used to describe the fluid behavior. In the case of the Lagrangian phase, Newton’s second law of motion is used to compute the particles trajectories and velocity. To study the effect of fluid behavior on the filler dispersion, several systems are modeled considering different fluid types (generalized Newtonian or viscoelastic) and particles volume fraction and size. The results obtained are used to correlate the fluid and particle characteristics on the effectiveness of mixing and morphology obtained.

Improvign mixing processes with the help of numerical tools

In several industrial applications, highly complex behaviour materials are used together with intricate mixing processes, which difficult the achievement of the desired properties for the produced materials. This is the case of the well-known dispersion of nano-sized fillers in a melt polymer matrix, used to improve the nanocomposite mechanical and/or electrical properties. This mixing is usually performed in twin-screw extruders, that promote complex flow patterns, and, since an in loco analysis of the material evolution and mixing is difficult to perform, numerical tools can be very useful to predict the evolution and behaviour of the material. This work presents a numerical based study to improve the understanding of mixing processes. Initial numerical studies were performed with generalized Newtonian fluids, but, due to the null relaxation time that characterize this type of fluids, the assumption of viscoelastic behavior was required. Therefore, the polymer melt was rheologically characterized, and, a six mode Phan-Thien-Tanner and Giesekus models were used to fit the rheological data. These viscoelastic rheological models were used to model the process. The conclusions obtained in this work provide additional and useful data to correlate the type and intensity of the deformation history promoted to the polymer nanocomposite and the quality of the mixing obtained.

Measurement of colour flow with the jet pull angle in t(t)over-bar events using the ATLAS detector at root s=8 TeV

The distribution and orientation of energy inside jets is predicted to be an experimental handle on colour connections between the hard--scatter quarks and gluons initiating the jets. This Letter presents a measurement of the distribution of one such variable, the jet pull angle. The pull angle is measured for jets produced in tt¯ events with one W boson decaying leptonically and the other decaying to jets using 20.3 fb−1 of data recorded with the ATLAS detector at a centre--of--mass energy of s√=8 TeV at the LHC. The jet pull angle distribution is corrected for detector resolution and acceptance effects and is compared to various models.