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.

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. (...)

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.

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

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

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.

Effect of binder on performance of intumescent coatings

This study investigates the role of the polymeric binder on the properties and performance of an intumescent coating. Waterborne resins of different types (vinylic, acrylic, and styrene-acrylic) were incorporated in an intumescent paint formulation, and characterized extensively in terms of thermal degradation behavior, intumescence thickness, and thermal insulation. Thermal microscopy images of charred foam development provided further information on the particular performance of each type of coating upon heating. The best foam expansion and heat protection results were obtained with the vinyl binders. Rheological measurements showed a complex evolution of the viscoelastic characteristics of the materials with temperature. As an example, the vinyl binders unexpectedly hardened significantly after thermal degradation. The values of storage moduli obtained at the onset of foam blowing (melamine decomposition) were used to explain different intumescence expansion behaviors.

Mechanical analysis of asphalt mixtures produced with waste plastic modified binders

This work compares the viscoelastic properties of an asphalt binder (70/100 pen) modified with different waste plastics and the mechanical properties of the resultant asphalt mixtures. Two different plastic wastes were used, namely recycled HDPE and EVA. Three different polymer modified binders were produced with these plastic wastes: i) 5% HDPE modified binder (P5); ii) 5% EVA modified binder (E5) and; iii) a modified binder with 4% of EVA and 2% HDPE (E4P2). Asphalt mixtures were produced with these modified binders, and their mechanical properties were analysed and compared with a conventional mixture produced with a 30/50 pen bitumen. It was possible to conclude that these recycled polymers are able to improve the mechanical performance of the asphalt mixtures used in road paving.

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

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.

Multiphasic, multistructured and hierarchical strategies for cartilage regeneration

Cartilage tissue is a complex nonlinear, viscoelastic, anisotropic, and multiphasic material with a very low coefficient of friction, which allows to withstand millions of cycles of joint loading over decades of wear. Upon damage, cartilage tissue has a low self-reparative capacity due to the lack of neural connections, vascularization, and a latent pool of stem/chondroprogenitor cells. Therefore, the healing of articular cartilage defects remains a significant clinical challenge, affecting millions of people worldwide. A plethora of biomaterials have been proposed to fabricate devices for cartilage regeneration, assuming a wide range of forms and structures, such as sponges, hydrogels, capsules, fibers, and microparticles. In common, the fabricated devices were designed taking in consideration that to fully achieve the regeneration of functional cartilage it is mandatory a well-orchestrated interplay of biomechanical properties, unique hierarchical structures, extracellular matrix (ECM), and bioactive factors. In fact, the main challenge in cartilage tissue engineering is to design an engineered device able to mimic the highly organized zonal architecture of articular cartilage, specifically its spatiomechanical properties and ECM composition, while inducing chondrogenesis, either by the proliferation of chondrocytes or by stimulating the chondrogenic differentiation  of stem/chondro-progenitor cells. In this chapter we present the recent advances in the development of innovative and complex biomaterials that fulfill the required structural key elements for cartilage regeneration. In particular, multiphasic, multiscale, multilayered, and hierarchical strategies composed by single or multiple biomaterials combined in a welldefined structure will be addressed. Those strategies include biomimetic scaffolds mimicking the structure of articular cartilage or engineered scaffolds as models of research to fully understand the biological mechanisms that influence the regeneration of cartilage tissue.