Next Article in Journal
Preparation of p-Phenylenediamine Modified Graphene Foam/Polyaniline@Epoxy Composite with Superior Thermal and EMI Shielding Performance
Next Article in Special Issue
Modelling of Elongational Flow of HDPE Melts by Hierarchical Multi-Mode Molecular Stress Function Model
Previous Article in Journal
Silicone Rubber Composites Reinforced by Carbon Nanofillers and Their Hybrids for Various Applications: A Review
Previous Article in Special Issue
Simulation Approach for Hydrophobicity Replication via Injection Molding
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Influence of Viscoelastic Modeling in Fluid Flow Simulations of Gum Acrylonitrile Butadiene Rubber

by
Sebastian Stieger
1,*,
Evan Mitsoulis
2,
Matthias Walluch
3,
Catharina Ebner
4,
Roman Christopher Kerschbaumer
5,
Matthias Haselmann
5,
Mehdi Mostafaiyan
6,
Markus Kämpfe
6,
Ines Kühnert
6,
Sven Wießner
6,7 and
Walter Friesenbichler
1
1
Institute of Injection Moulding of Polymers, Montanuniversitaet Leoben, Otto Gloeckel-Straße 2, 8700 Leoben, Austria
2
School of Mining Engineering and Metallurgy, National Technical University of Athens, 157 80 Zografou, Greece
3
Anton Paar GmbH, Anton-Paar Straße 20, 8054 Graz, Austria
4
Institute of Chemistry of Polymeric Materials, Montanuniversitaet Leoben, Otto Gloeckel-Straße 2, 8700 Leoben, Austria
5
Polymer Competence Center Leoben GmbH, Roseggerstrasse 12, 8700 Leoben, Austria
6
Leibniz-Institut für Polymerforschung Dresden e.V., Hohe Straße 6, 01069 Dresden, Germany
7
Institute of Materials Science, Technische Universität Dresden, 01062 Dresden, Germany
*
Author to whom correspondence should be addressed.
Polymers 2021, 13(14), 2323; https://doi.org/10.3390/polym13142323
Submission received: 10 June 2021 / Revised: 8 July 2021 / Accepted: 9 July 2021 / Published: 15 July 2021

Abstract

:
Computational fluid dynamics (CFD) simulation is an important tool as it enables engineers to study different design options without a time-consuming experimental workload. However, the prediction accuracy of any CFD simulation depends upon the set boundary conditions and upon the applied rheological constitutive equation. In the present study the viscoelastic nature of an unfilled gum acrylonitrile butadiene rubber (NBR) is considered by applying the integral and time-dependent Kaye–Bernstein–Kearsley–Zapas (K-BKZ) rheological model. First, exhaustive testing is carried out in the linear viscoelastic (LVE) and non-LVE deformation range including small amplitude oscillatory shear (SAOS) as well as high pressure capillary rheometer (HPCR) tests. Next, three abrupt capillary dies and one tapered orifice die are modeled in Ansys POLYFLOW. The pressure prediction accuracy of the K-BKZ/Wagner model was found to be excellent and insensitive to the applied normal force in SAOS testing as well as to the relation of first and second normal stress differences, provided that damping parameters are fitted to steady-state rheological data. Moreover, the crucial importance of viscoelastic modeling is proven for rubber materials, as two generalized Newtonian fluid (GNF) flow models severely underestimate measured pressure data, especially in contraction flow-dominated geometries.

Graphical Abstract

1. Introduction

In order to solve macroscale flow problems like mold filling or extrusion, use is made of continuum mechanics by omitting microscopic discontinuities of the investigated fluid. Discretizing the area of interest by an appropriate mesh, conservation equations of mass, momentum, and energy, are applied, leading to a set of partial differential equations, which are typically solved by means of computational fluid dynamics (CFD) simulation. The prediction accuracy depends upon the set boundary conditions and upon the applied rheological constitutive equation. In most CFD simulations, the viscoelastic nature of polymers is not taken into account. Instead, generalized Newtonian fluid (GNF) flow models are used to describe the relation between the stress and rate of deformation tensors. However, these models fail to reflect important rheological properties (e.g., normal stress differences or transient data). Moreover, they intrinsically assume a ratio of three (Trouton ratio) between steady-state shear and uniaxial elongational viscosities. As a result, they fail to predict extrudate (die) swells or inlet vortices and massively underestimate pressure drops in contraction flow areas [1,2]. Hence, developing and improving viscoelastic constitutive equations is an eminent subject in polymer rheology research [3].
Numerous differential [4,5,6,7,8] and integral [9,10,11,12] models have been proposed to describe the viscoelastic behavior of polymeric fluids. Most of them are based on either a network [7,9,12] or reptation (tube) [6,13] theory. Although some differential equations such as the (extended) Pom-Pom [14] or exponential Phan Thien-Tanner (ePTT) [15] model are in a good agreement with rheological data, a conflict of interest arises as follows: An accurate description of a real polymer melt requires a multi-mode formulation. Each added number of mode increases the degrees of freedom (dof) and therefore the complexity from a numerical point of view [2].
The second group of models overcomes this problem by applying an integral formulation [2]. Inspired by the theory of rubber elasticity [11] Kaye [9], Bernstein, Kearsley and Zapas [12] proposed a temporary network model (K-BKZ), which has been continuously improved over the past 50 years [2,16]. In the current mathematical formulation implemented in Ansys POLYFLOW (Equations (1)–(3) and (6)) it is able to reflect steady-state shear η and steady-state uniaxial elongational η e viscosities, first N 1 and second N 2 normal stress differences, start-up shear η + and start-up uniaxial elongational η e + viscosities as well as storage G , loss G and relaxation G t moduli of unfilled polymer melts in the linear viscoelastic (LVE) and non-LVE region well [2]:
τ = = 1 1   θ t M t t H I C 1 ,   I C   C = t 1 t   +   θ   ·   C = t t d t
with the stress tensor τ = , the material constant θ , the present time t, the past time t , the memory function M t t , the damping function H I C 1 ,   I C , the Finger strain tensor C = t 1 t , and the Cauchy–Green tensor C = t t . The material constant θ is given by Equation (2):
N 2 N 1 = θ 1 θ  
In Equation (1) the stress results from an integral over the past time, where effects of strain are separated from those of time. The latter is described through a memory function:
M t t   = i = 1 N g i λ i · exp t t λ i  
with i representing the ith mode of the total number of modes N, the relaxation moduli g i and relaxation times λ i . The required relaxation spectrum can be easily obtained from small amplitude oscillatory shear (SAOS) tests and more precisely fitting G and G data ( ω is the angular frequency) applying Equations (4) and (5):
G ω   = i = 1 N g i   ω 2   λ i 2 1 + ω 2   λ i 2 ,
G ω   = i = 1 N g i   ω   λ i 1 + ω 2   λ i 2 .  
The dependency of rheological properties on the imposed strain is given by Wagner’s [17,18] damping function:
H I C 1 ,   I C   = exp ( α β I C 1 + 1 β I C 3 )  
where α and β are fitting parameters to be determined from shear and elongational data and I C 1 as well as I C are the first invariants of the Finger strain tensor and Cauchy–Green tensor, respectively.
Recently, different research groups had great success applying K-BKZ equations to reproduce benchmark polymer-flow phenomena such as inlet vortices [19] or extrudate swelling [20,21,22]. Moreover, the K-BKZ model was found to predict measured pressure drops with an accuracy that has not been reported for any other constitutive rheological model [21,23,24]. However, in all aforementioned cases the investigated fluid was an unfilled polyolefin melt.
In our recent studies [25,26] the applicability of the K-BKZ/Wagner model was tested to highly carbon black filled rubber compounds. This class of materials differs distinctively from polyolefin melts, as rubber compounds exhibit a material behavior dominated by elasticity. Moreover, they are processed at significant lower temperature levels (80 °C to 120 °C) and contain reinforcing filler particles. When comparing recorded pressure drops in various capillary dies to K-BKZ/Wagner predictions considerable deviations were observed [25,26]. One possible explanation is that the K-BKZ/Wagner model is not applicable to highly filled polymeric melts. The second possible explanation is that the K-BKZ/Wagner model is not applicable to rubber materials. In order to answer this open research question the present study tests for the first time the ability of the K-BKZ/Wagner model to correctly predict pressure drops of an unfilled gum rubber in CFD simulation.

2. Materials and Methods

2.1. Material

The investigated material is a butadiene acrylonitrile copolymer labeled PERBUNAN® 3965 F. This industrial grade was obtained from ARLANXEO (ARLANXEO Deutschland GmbH, Dormagen, Germany) containing no filler particles as well as no curing system (Figure A1). The Mooney viscosity ML (1 + 4) [27] is according to product specifications [28] without pre-treatment 65 ± 7 Mooney units (MUs), the acrylonitrile content is 39 ± 1.0 wt %, and specific gravity 0.99. It was delivered in bales (Figure A1a) and will be referred to as “gum NBR” in the following chapters.

2.2. Rheological Testing

2.2.1. SAOS

LVE rheological properties were detected performing small amplitude oscillatory shear tests applying the Modular Compact Rheometer (MCR) 501 (Anton Paar GmbH, Graz, Austria) in parallel-plate configuration with serrated surfaces [25]. For all experiments with the MCR 501 device, samples of 25 mm diameter and approximately 1 mm thickness were prepared. As a second measurement device the rotorless Rubber Process Analyzer (RPA) D-MDR 3000 (MonTech Werkstoffprüfmaschinen GmbH, Buchen, Germany) was tested, employing a closed, bi-conical, and grooved measurement chamber. In order to specify the LVE deformation range of gum NBR, first, runs with constant angular frequency and step-wise amplitude increase of the sinusoidal strain (“amplitude sweeps”) were done. One exemplary curve is displayed in Figure A2, where a strain amplitude of γ 0 = 3.5% was selected for further frequency sweep measurements, which is well within the linear viscoelastic deformation range of gum NBR. This procedure was performed for each individual measurement setting. Second, the influence of sample preparation (roller milled vs. compression molded specimens), normal force (Fn), pre-shearing as well as the influence of the measurement device (MCR 501 vs. RPA) was analyzed by keeping the selected strain amplitude constant and varying the angular frequency level (“frequency sweep”). Third, storage G and loss G moduli were recorded with the MCR 501 device at four different temperature levels (60, 80, 100, and 120 °C) and a master curve was derived by time-temperature shifting (TTS) to the reference temperature of 100 °C. The mastered LVE moduli were subsequently used to fit the memory function in Equation (3).

2.2.2. High Pressure Capillary Rheometry

The pressure-driven flow of gum NBR was studied employing a high pressure capillary rheometer (HPCR) of the type Rheograph 50 (GÖTTFERT Werkstoff-Prüfmaschinen GmbH, Buchen, Germany). In order to analyze possible slippage at the wall, first, tests with two abrupt (entrance angle ϕ = 180°) capillary dies of same length-to-diameter (L/D) ratio (L/D = 10/1, 20/2) were carried out. The entrance pressure loss was measured in both cases with tapered (ϕ = 90°) orifice dies (L/D = 0.2/1, 0.2/2). Second, the pressure linearity in dependence of the capillary length (“Bagley plot”) was analyzed by measuring pressure drops in two additional abrupt capillary dies (L/D = 5/1, 20/1). Third, the steady-state shear viscosity η was calculated for the 1 mm diameter die set applying Bagley and Weißenberg–Rabinowitsch corrections. This material property (η) was used to fit the damping parameter α in Equation (6) as well as GNF models. Finally, Binding’s model [29] was selected to estimate the steady-state uniaxial elongational viscosity η e from entrance pressure losses and subsequently used to fit the final damping parameter β in Equation (6). All tests were performed at a reference temperature of 100 °C.

2.3. Constitutive and Numerical Modeling

In order to test K-BKZ/Wagner’s applicability to correctly predict pressure drops of gum NBR, the commercial software package Ansys POLYFLOW (ANSYS Inc., Canonsburg, CA, USA) was selected to simulate the pressure-driven flow of the high pressure capillary rheometer (HPCR) experiment.
A two-dimensional axisymmetric finite element model was built for one tapered orifice die (L/D = 0.2/1) and three abrupt capillary dies (L/D = 5/1, 10/1, 20/1). Figure 1a illustrates a schematic drawing of the numerical setup including boundary conditions (BCs), geometrical dimensions, as well as the position of the pressure evaluation (p), which corresponds to the position of the pressure transducer in the HPCR. The mesh design was based on [20,24], who used a denser grid moving toward the stress singularity at the entrance of the die.
Mesh independency was checked in one of our previous studies [25] by comparing predicted pressure drops using the meshes displayed in Figure 1b to denser ones. The applied meshes in the present study (Figure 1b) consisted of 936 elements (1099 nodes) for the tapered orifice die (L/D = 0.2/1), 450 elements (501 nodes) for the short abrupt capillary die (L/D = 5/1), and 600 elements (666 nodes) for the two remaining abrupt capillary dies (L/D = 10/1, 20/1).
In all calculations steady-state and isothermal flow conditions were assumed leading for an incompressible fluid to the following governing equations:
  ·   v ¯ = 0 ,
p +   ·   τ = = 0
where v ¯ is the velocity vector, p is the pressure, and τ = is the stress tensor. Effects of inertia and gravity were neglected due to the high viscosity of gum NBR. In case of viscoelastic modeling the stress tensor τ = was calculated with Equations (1)–(3) and (6), which represent the K-BKZ/Wagner model. Viscous stresses were computed as:
τ = = 2   η γ ˙   D =
where D = is the rate of deformation tensor, with two different GNF flow models describing the shear dependence of the viscosity η γ ˙ :
η γ ˙ = K λ   ·   γ ˙ n 1 .
The power-law model in Equation (10) is arguably the simplest and most common GNF model in rubber rheology with K as the consistency index, λ as the natural time, and n as the power-law index. It assumes a power-law behavior of the steady-state shear viscosity even in the LVE region. However, most polymer melts exhibit a constant viscosity at low shear rates (“Newtonian plateau”). Consequently, a second GNF model was applied able to reflect the aforementioned material behavior properly:
η γ ˙   =   η 0 / 1 + λ   · γ ˙ m .
The Cross model in Equation (11) consists of the three model parameters η 0 , λ , and m , which represent the zero-shear-rate viscosity, the natural time, and the Cross-law index (m ≙ 1 −   n ), respectively.
Finally, all boundary conditions applied in this study are listed in Table 1 with corresponding surfaces illustrated in Figure 1a.

3. Results and Discussion

The overall goal of the study is to answer the open research question whether the K-BKZ/Wagner model is able to correctly predict pressure drops of an unfilled gum rubber in CFD simulation. In order to address this objective, a specific research approach was designed, which is depicted in Figure 2 in a flow chart schematic.
First, rheological tests are performed in the LVE (SAOS) and non-LVE (HPCR) deformation range of gum NBR. The material properties obtained from these measurements are storage (G’) and loss (G”) moduli, both in dependence of the angular frequency (ω) as well as steady-state shear (η) and uniaxial elongational ( η e ) viscosities, which are in dependence of the shear ( γ ˙ ) and strain ( ε ˙ ) rate, respectively. These material properties are used to determine the model parameters of the K-BKZ/Wagner, Cross (GNF), and power-law (GNF) models. Additionally, the influence of normal force applied on the sample in SAOS testing as well as the influence of the model parameter θ (Equation (2)) is analyzed by performing in total three K-BKZ fits. Finally, CFD simulations are performed aiming to compare predicted with recorded pressure data of the high pressure capillary rheometer experiment.

3.1. Rheological Testing and Constitutive Modeling

In rubber rheology a testing instrument called Rubber Process Analyzer is widely employed to determine viscoelastic properties before, during, and after curing. This rotational shear rheometer consists of a closed and sealed die system, where the tested rubber is compressed by a pre-defined clamping pressure leading to normal forces on the sample during the measurement. This setup aims to prevent any slippage and differs from conventional dynamic rheometers like the MCR 501 device, where test specifications recommend keeping normal forces at a minimum. Comparing the complex viscosity ( η * ) of gum NBR, normal forces applied on the sample did not affect the shape of the curve but shifted measured data to higher values (Figure 3a). Applying a normal force of 10 N using the MCR 501 device, η * corresponded to RPA data. On the other hand, neither sample preparation (roller milled vs. compression molded) nor pre-shearing the sample for 240 s with a strain amplitude of 42% and constant angular frequency of 31.42 rad/s (recommended by Fasching [30] for highly filled rubber compounds) had any effect on η * (Figure 3). The detected insensitivity of gum NBR validated the listed set of pre-shearing parameters, as they did not damage macromolecular chains and may serve for highly filled rubber compounds as a fixed precondition to break down filler networks and minimize their impact on measured LVE properties [25]. The mean measurement deviations of the four settings displayed in Figure 3 (RPA; RPA, presheared; MCR, Fn < 1 N; MCR 501, Fn = 10 N) are 0.3, 0.3, 1.8, and 3.4%, respectively. In the present study, both MCR setups (Fn < 1 N, Fn = 10 N) were further used to determine G’ and G” aiming to fit the memory function in Equation (3) and finally analyze the influence of normal force in SAOS testing to viscoelastic modeling and pressure prediction in CFD simulation.
Next, time-temperature superposition was applied, to extend the frequency range of LVE modeling. From a thermo-rheological point of view, rubber materials are complex fluids, where vertical shifting may be necessary in order to construct a master curve for G’ and G”. Consequently, a guideline presented in [25] was followed, where horizontal shift factors ( a t ) were obtained by mastering first the loss factor ( tan δ ) , which is intrinsically invariant to vertical shifting ( b t ):
tan δ = b t · G a t   ·   ω b t · G a t   ·   ω .
Figure 4a displays the mastered dissipation factor at the reference temperature of 100 °C. Moreover, natural logarithms of a t were plotted against the inverse temperature, resulting in the so-called “Arrhenius Plot”.
The calculated activation energy (Ea) of 68 kJ/mol corresponded well with values obtained for NBR and HNBR compounds [25,26] and was insensitive to Fn (Figure A3a). Second, vertical shifting was applied to G’ and G’’ measured at 60, 80, and 120 °C, minimizing any deviations to the reference temperature segment (Figure 4b). Over the whole measured angular frequency range, G’ exceeded G” (no crossover point), proving a material response that was dominated by elasticity. Finally, LVE moduli were used to determine relaxation moduli g i of the memory function in Equation (3). A fit was performed, which minimized the sum of squared differences according to Equations (4) and (5), keeping pre-defined relaxation times λ i and modes (N = 8) constant. To give the fit a physical meaning, a constraint was added, which ensured that relaxation moduli would not decrease with increasing corresponding relaxation time ( g i g i + 1 ). The obtained material parameters (Table A1) well described the LVE material behavior of gum NBR in a frequency range of almost five decades (Figure 4b). The mean measurement deviations of the three recorded properties G’, G”, and tan (δ) were 2.6, 3.1, and 0.6%, respectively.
Next, the exact same procedure was carried out for the second rheological characterization setup (Fn = 10 N). Table A2 lists corresponding material parameters and Figure A3b compares model predictions with measured data. The mean measurement deviations of G’, G”, and tan (δ) with an applied normal force of 10 N were 3.9, 3.9, and 0.8%, respectively. In the present study “K-BKZ fit 1” represented LVE data measured according to conventional test specifications (Fn < 1 N) and “K-BKZ fit 2” represented LVE data, where a normal force of 10 N was applied on the sample during the measurement. The latter was similar to conditions in a Rubber Process Analyzer for gum NBR.
A key assumption in rheological testing is no slip at the wall, which indicates adhesion of the investigated fluid to the capillary die or to the plates, in cases of SAOS testing. However, as reviewed by Hatzikiriakos [31], this classic boundary condition is not always valid for polymer melts, especially at higher wall shear stresses. A popular method for detection was introduced by Mooney [32], who proved that for a wall slipping fluid the observed flow curves become dependent on the diameter of the capillary die. Comparing flow curves recorded with two dies of same L/D ratio but different geometry, gum NBR exhibits no clear dependency on the die diameter (Figure 5).
If in fact slip was present in HPCR experiments, one would expect the wall shear stress in the L/D = 20/2 die to be larger than in the L/D = 10/1 die, especially at higher apparent shear rate levels, and the calculated steady-state shear viscosity to be substantially lower than the complex viscosity. Comparing the aforementioned material properties (Figure A4), the steady-state shear viscosity of gum NBR even exceeded the complex viscosity (Fn < 1 N). For the second SAOS testing setup, where a normal force of 10 N was applied on the sample during the measurement, the Cox–Merz rule [33] ( η * ω = η( γ ˙ )   ω   = γ ˙ ) was obeyed. Moreover, at the highest apparent shear rate level the wall shear stress in the L/D = 20/2 capillary die was lower than in the L/D = 10/1 die. These observations indicated no slippage at the wall for gum NBR.
Furthermore, strong extrudate swelling was observed even at low apparent shear rate levels (Figure A5). This non-linear flow phenomenon arises from stored energy of uncoiled polymer molecules due to shear and elongational stresses in the capillary die. These uncoiled molecules strive to reach the former state of higher entropy resulting in additional forces that press the fluid in the capillary against the wall and may prevent slippage. Thus, a no slip BC was applied in CFD simulations.
Next, the detected pressure drops of the 1 mm die set were plotted in dependence of the capillary length (Bagley plot [34]), with mean measurement deviations of 1.5, 3.1, 1.1, and 1.4% for the L/D = 0.2, 5, 10, and 20 capillary dies, respectively. The apparent shear rate was calculated applying Equation (13):
γ ˙ a = 32 Q π D 3
where Q is the volumetric flow rate and D the diameter of the capillary die.
All apparent shear rate levels exhibited high linearity of the pressure drop, despite some observed pressure fluctuations (Figure A6), with coefficients of determination (R2) exceeding 98% (Figure 6).
A highly linear Bagley plot indicated that pressure and temperature dependencies of the viscosity effectively leveled off each other [25,26]. Since Ansys POLYFLOW offers no rheological model, which accounts for the pressure dependency of the viscosity, the best option was to consider neither of them. Consequently, isothermal flow conditions were assumed in CFD simulations, and η was modeled as a simple shear rate dependent function.
Next to the two damping parameters α and β , the material constant θ also needed to be determined to complete the K-BKZ modeling. For most polymeric melts small negative values are reported for θ [35]. Since no experimental data of the second normal stress difference were available for gum NBR, this study tested the sensitivity of θ with respect to pressure prediction in CFD simulation. Thus, a third fit (K-BKZ fit 3) was performed, which maintained the relaxation spectrum of K-BKZ fit 1 and set θ to –0.15 (Table A3). K-BKZ fit 1 (Table A1) and K-BKZ fit 2 (Table A2) both assumed θ to be zero.
In a simple shear flow, Equation (6) depends only on one parameter, so the steady-state shear viscosity was used next to determine α . The final model parameter β was obtained from uniaxial elongational viscosity data calculated from entrance pressure losses according to Binding [29]. In general, extensional rheometers developed by Sentmanat [36], Meissner [37,38], or Münstedt [39] are recommended to characterize the elongational behavior of polymeric melts. However, in the present study we preferred material properties obtained with the HPCR to fit damping parameters α and β , since recorded pressure drops of the exact same measurement device are aimed to be predicted in CFD simulations. All fits were performed in MATLAB minimizing the sum of squared differences to measured data.
Despite applying different relaxation spectra (fit 1 vs. fit 2) and different values for the material constant θ (fit 1 vs. fit 3), model predictions of all fits described the steady-state shear viscosity of gum NBR well (Figure 7).
A higher damping parameter α compensated for the vertical shifted LVE data of K-BKZ fit 2 (Fn = 10 N). However, all fits overestimated the uniaxial elongational viscosity at low strain rates and underestimated η e at higher ones. To improve K-BKZ modeling of elongational properties, Luo and Tanner [40] proposed an approach that assigns a different value of β to each mode ( β i ). As multiple betas are not implemented in Ansys POLYFLOW, a best-fitted single value of β was employed despite the observed deviations.
Finally, coefficients (Table A4 and Table A5) of the two selected GNF models were determined applying the fitting tool implemented in Ansys POLYFLOW. Figure A7 proves the ability of both models to describe the steady-state shear viscosity of gum NBR well. However, the Cross model assumes a Newtonian plateau, resulting in differences between the two models in the lower shear rate region.

3.2. Numerical Modeling and Evaluation

The overall goal of this study was to test the ability of the K-BKZ/Wagner model to correctly predict pressure drops of an unfilled gum rubber in CFD simulation. Thus, four different capillary dies applied in HPCR tests with varying L/D ratios (0.2/1, 5/1, 10/1, 20/1) were modeled in Ansys POLYFLOW (Figure 1). At the inlet a fully developed normal velocity profile with a pre-defined volumetric flow rate Q was imposed (BC 1). Moreover, a normal and tangential velocity of zero was assigned at the capillary walls (BC 2). Next, use was made of the axisymmetric die geometry by applying a tangential force and normal velocity of zero at the axis of symmetry (BC 3). Since extruded strands of gum NBR clearly exhibited die swelling at all apparent shear rate levels (Figure A5), a “zero force” boundary condition (BC 4) was applied. This BC considered non-zero normal stresses at the outlet and, as a consequence, exit pressure dropped associated with die swelling.
First, the influence of normal forces applied on the sample in SAOS testing (K-BKZ fit 2) as well as the influence of the model parameter θ (K-BKZ fit 3) were analyzed. For this purpose, Figure 8 compares recorded HPCR data with CFD simulation results, which proved the excellent pressure prediction accuracy of all three K-BKZ fits for gum NBR.
Mean deviations of K-BKZ fit 1, fit 2, and fit 3 to recorded pressure data were 8.6, 10.1, and 8.3%, respectively, indicating an insensitivity to the applied normal force and to the model parameter θ .
The tapered orifice (L/D = 0.2/1) represented an extreme condition, which was dominated by the contraction flow from the larger reservoir (   = 15 mm) in the small die (   = 1 mm). Thus, its measured pressure level was strongly affected by both entrance and exit pressure drops and consequently by the viscoelasticity of the material. The K-BKZ/Wagner model was able to correctly predict these extreme conditions for gum NBR with mean deviations (compared to recorded pressure drops of the L/D = 0.2/1 die) of 7.5, 5.9, and 9.3% for K-BKZ fit 1, fit 2, and fit 3, respectively. In contrast, pressure drops in the L/D = 20/1 die were dominated by the capillary flow and consequently by shear deformation. However, the K-BKZ model predicted also the second condition well with mean deviations (compared to recorded pressure drops of the L/D = 20/1 die) of 8.4, 7.6, and 8.0% for K-BKZ fit 1, fit 2, and fit 3, respectively. These results proved that the K-BKZ/Wagner was in fact able to correctly predict pressure drops of an unfilled gum rubber in both contraction and shear-dominated geometries. Thus, the observed deviations in our recent studies for two highly carbon black filled rubber compounds [25,26] can now be attributed to the inability of the K-BKZ/Wagner model to reflect the decreasing molecular mobility with increasing filler content. As a consequence, the K-BKZ/Wagner model is not applicable to highly filled polymer systems. The research groups of Walter Friesenbichler, Evan Mitsoulis, Ines Kühnert, and Sven Wießner have agreed to address this issue, aiming to develop a holistic viscoelastic–plastic constitutive rheological model for particle-filled, multiphase polymer systems based on the original K-BKZ equations.
Finally, the influence of viscoelastic modeling was assessed by comparing CFD simulation results of K-BKZ fit 1 to pressure predictions of two generalized Newtonian fluid flow models (Figure 9).
Both GNF models (BCs listed in Table 1) underestimated measured pressure data distinctly. Deviations increased with rising apparent shear rate level and decreasing L/D ratio. Any GNF model considers viscous stresses only and intrinsically assumes a ratio of three between steady-state shear and uniaxial elongational viscosities. However, gum NBR exhibited a material response dominated by elasticity (Figure 4b), leading to an additional pressure drop at the exit and a ratio between η and η e , clearly exceeding three (Figure 7), leading to an enhanced pressure level at the entrance. With decreasing L/D ratio these two factors became more and more important and mean deviations increased accordingly.
Moreover, Figure 7 displays an increasing ratio between η and η e with rising apparent shear rate level ( η e / η Cross = 10.4 γ ˙ , ε ˙ = 26   s 1 ;   η e / η Cross = 19 γ ˙ , ε ˙ = 207   s 1 ). Thus, the largest deviations between measured data and GNF predictions were expected and found to be true for the orifice die (L/D = 0.2) at the highest apparent shear rate level. Considering all dies as well as apparent shear rate levels, mean deviations to recorded pressure data were 8.6, 49, and 39.3% for K-BKZ, Cross, and power-law models, respectively (Figure 9). These values highlighted the importance of viscoelastic modeling to correctly predict pressure drops of an unfilled gum rubber especially in contraction flow-dominated geometries.

4. Summary and Conclusions

The present study tested the applicability of the viscoelastic K-BKZ/Wagner model in CFD simulations to correctly predict pressure drops of an unfilled gum acrylonitrile butadiene rubber (NBR). First, exhaustive rheological testing was performed in the linear viscoelastic (LVE) and non-LVE deformation range, including small amplitude oscillatory shear (SAOS) as well as high pressure capillary rheometer (HPCR) tests. The obtained material properties were subsequently used to fit the K-BKZ/Wagner and two generalized Newtonian fluid (GNF) flow models (Cross, power-law). Moreover, the influence of normal forces applied on the sample in SAOS testing as well as the relation of first and second normal stress differences to viscoelastic modeling and CFD simulations were analyzed by performing in total three K-BKZ/Wagner fits. Second, capillary dies applied in HPCR tests with varying length-to-diameter (L/D) ratios were modeled in Ansys POLYFLOW. Comparing CFD simulation results to recorded pressure data proved (i) an excellent pressure prediction accuracy for the K-BKZ/Wagner model and consequently its applicability to unfilled rubber materials; (ii) an insensitivity of the K-BKZ/Wagner model to the applied normal force in SAOS testing as well as to the relation of first and second normal stress differences, provided that damping parameters are fitted to steady-state rheological data; and (iii) the crucial importance of viscoelastic modeling for unfilled gum rubber materials, as both GNF flow models severely underestimated measured pressure data, especially in contraction flow-dominated geometries.

Author Contributions

Conceptualization and methodology, S.S.; software, S.S. and M.H.; validation, S.S. and E.M.; formal analysis, S.S. and C.E.; investigation, S.S. and E.M.; resources, W.F. and M.W.; data curation, S.S., M.K., M.M., R.C.K. and M.W.; writing—original draft preparation, S.S.; writing—review and editing, S.S., E.M. and R.C.K.; visualization S.S., M.K., M.M. and C.E.; supervision, W.F., E.M., I.K. and S.W.; project administration, W.F., I.K. and S.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research work was supported by the Austrian Research Promotion Agency (FFG) as part of the “RubExject II” project (corresponding project number 855873) and the company partners SKF Sealing Solutions Austria GmbH (Judenburg, Austria), IB Steiner (Spielberg, Austria), and ELMET Elastomere Produktions- und Dienstleistungs-GmbH (Oftering, Austria).

Data Availability Statement

All data presented are only available after request from the corresponding author assuming a formal approval of our company partners.

Acknowledgments

The authors thank Ivica Duretek and Stephan Schuschnigg for their technical support and Erwin Mach Gummitechnik GmbH, Hirm, Austria, for providing the material.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix

Figure A1. (a) Image of gum NBR (PERBUNAN® 3965 F) delivered as bale; (b) scanning electron microscope (SEM) image of an exemplary compression molded sample ( = 8 mm, thickness = 1 mm) of gum NBR.
Figure A1. (a) Image of gum NBR (PERBUNAN® 3965 F) delivered as bale; (b) scanning electron microscope (SEM) image of an exemplary compression molded sample ( = 8 mm, thickness = 1 mm) of gum NBR.
Polymers 13 02323 g0a1
Figure A2. Exemplary amplitude sweep measurement performed with the MCR 501 device at 100 °C to detect the linear viscoelastic deformation range of gum NBR, with the selected shear strain rate of 3.5% for the frequency sweep measurements.
Figure A2. Exemplary amplitude sweep measurement performed with the MCR 501 device at 100 °C to detect the linear viscoelastic deformation range of gum NBR, with the selected shear strain rate of 3.5% for the frequency sweep measurements.
Polymers 13 02323 g0a2
Figure A3. (a) Time-temperature superposition of the phase angle δ and corresponding Arrhenius plot; (b) K-BKZ model predictions of storage G and loss G moduli compared to experimental data at the reference temperature Tref; K-BKZ fit 2 represents fitted LVE data measured at Fn = 10 N.
Figure A3. (a) Time-temperature superposition of the phase angle δ and corresponding Arrhenius plot; (b) K-BKZ model predictions of storage G and loss G moduli compared to experimental data at the reference temperature Tref; K-BKZ fit 2 represents fitted LVE data measured at Fn = 10 N.
Polymers 13 02323 g0a3
Figure A4. Complex viscosities η * compared to steady-state shear viscosity η.
Figure A4. Complex viscosities η * compared to steady-state shear viscosity η.
Polymers 13 02323 g0a4
Table A1. K-BKZ/Wagner model parameters (fit 1) for gum NBR at the reference temperature of 100 °C ( θ = 0); Fit 1 represents fitted LVE data measured at Fn < 1 N.
Table A1. K-BKZ/Wagner model parameters (fit 1) for gum NBR at the reference temperature of 100 °C ( θ = 0); Fit 1 represents fitted LVE data measured at Fn < 1 N.
Mode
λ i   s
g i   Pa
α   /
β   /
11.0000 × 10−51.6671 × 1050.12520.25
27.1969 × 10−51.6671 × 105
35.1795 × 10−31.6671 × 105
43.7276 × 10−31.6671 × 105
52.6827 × 10−21.6671 × 105
61.9307 × 10−17.4994 × 104
71.3895 × 10−02.9212 × 104
81.0000 × 1012.3313 × 104
Table A2. K-BKZ/Wagner model parameters (fit 2) for gum NBR at the reference temperature of 100 °C ( θ = 0); Fit 2 represents fitted LVE data measured at Fn = 10 N.
Table A2. K-BKZ/Wagner model parameters (fit 2) for gum NBR at the reference temperature of 100 °C ( θ = 0); Fit 2 represents fitted LVE data measured at Fn = 10 N.
Mode
λ i   s
g i   Pa
α   /
β   /
11.0000 × 10−52.3535 × 1050.18500.14
27.1969 × 10−52.3535 × 105
35.1795 × 10−32.3535 × 105
43.7276 × 10−32.3535 × 105
52.6827 × 10−22.3535 × 105
61.9307 × 10−11.0004 × 104
71.3895 × 1003.9518 × 104
81.0000 × 1013.0996 × 104
Table A3. K-BKZ/Wagner model parameters (fit 3) for gum NBR at the reference temperature of 100 °C ( θ = –0.15); Fit 3 represents fitted LVE data measured at Fn < 1 N.
Table A3. K-BKZ/Wagner model parameters (fit 3) for gum NBR at the reference temperature of 100 °C ( θ = –0.15); Fit 3 represents fitted LVE data measured at Fn < 1 N.
Mode
λ i   s
g i   Pa
α   /
β   /
11.0000 × 10−51.6671 × 1050.12600.20
27.1969 × 10−51.6671 × 105
35.1795 × 10−31.6671 × 105
43.7276 × 10−31.6671 × 105
52.6827 × 10−21.6671 × 105
61.9307 × 10−17.4994 × 104
71.3895 × 1002.9212 × 104
81.0000 × 1012.3313 × 104
Figure A5. Light microscope images of strands extruded through an abrupt capillary die with a length-to-diameter (L/D) ratio of 10/1 at an apparent shear rate level of (a) 27 s−1; (b) 30906 s−1.
Figure A5. Light microscope images of strands extruded through an abrupt capillary die with a length-to-diameter (L/D) ratio of 10/1 at an apparent shear rate level of (a) 27 s−1; (b) 30906 s−1.
Polymers 13 02323 g0a5
Figure A6. Recorded pressure in the L/D = 10/1 capillary die in dependence of the measurement time with apparent shear rate levels displayed ranging from 27 s−1 to 3906 s−1.
Figure A6. Recorded pressure in the L/D = 10/1 capillary die in dependence of the measurement time with apparent shear rate levels displayed ranging from 27 s−1 to 3906 s−1.
Polymers 13 02323 g0a6
Figure A7. Steady-state shear viscosity of gum NBR compared to power-law and Cross model predictions.
Figure A7. Steady-state shear viscosity of gum NBR compared to power-law and Cross model predictions.
Polymers 13 02323 g0a7
Table A4. Cross model parameters for gum NBR at the reference temperature of 100 °C.
Table A4. Cross model parameters for gum NBR at the reference temperature of 100 °C.
η 0 Pa   ·   s λ   s m   /
5.076 × 1041.31 × 10−19.095 × 10−1
Table A5. Power-law model parameters for gum NBR at the reference temperature of 100 °C.
Table A5. Power-law model parameters for gum NBR at the reference temperature of 100 °C.
K   Pa   ·   s λ   s n   /
8.467 × 1021.55 × 10−31.254 × 10−1

References

  1. Mitsoulis, E.; Battisti, M.; Neunhäuserer, A.; Perko, L.; Friesenbichler, W. Flow behavior of PP-polymer nanocomposites in injection molding hyperbolical dies. Adv. Polym. Technol. 2018, 37, 3045–3055. [Google Scholar] [CrossRef]
  2. Mitsoulis, E. 50 Years of the K-BKZ Constitutive Relation for Polymers. ISRN Polym. Sci. 2013, 2013, 952379. [Google Scholar] [CrossRef]
  3. Larson, R.G. Constitutive Equations for Polymer Melts and Solutions; Butterworth-Heinemann: Boston, MA, USA, 1988; ISBN 9781483162867. [Google Scholar]
  4. Giesekus, H. A simple constitutive equation for polymer fluids based on the concept of deformation-dependent tensorial mobility. J. Non-Newton. Fluid Mech. 1982, 11, 69–109. [Google Scholar] [CrossRef]
  5. Oldroyd, J.G. On the formulation of rheological equations of state. Proc. R. Soc. Lond. A 1950, 200, 523–541. [Google Scholar] [CrossRef]
  6. McLeish, T.C.B.; Larson, R.G. Molecular constitutive equations for a class of branched polymers: The pom-pom polymer. J. Rheol. 1998, 42, 81–110. [Google Scholar] [CrossRef]
  7. Phan-Thien, N. A Nonlinear Network Viscoelastic Model. J. Rheol. 1978, 22, 259–283. [Google Scholar] [CrossRef]
  8. Leonov, A.I. On the rheology of filled polymers. J. Rheol. 1990, 34, 1039–1068. [Google Scholar] [CrossRef]
  9. Kaye, A. Non-Newtonian Flow in Incompressible Fluids; Aerospace Engineering Reports; College of Aeronautics Cranfield: Cranfield, UK, 1962. [Google Scholar]
  10. Doi, M.; Edwards, S.F. The Theory of Polymer Dynamics; Clarendon Press: Oxford, UK, 1986; ISBN 9780198519768. [Google Scholar]
  11. Lodge, A.S. A network theory of flow birefringence and stress in concentrated polymer solutions. Trans. Faraday Soc. 1956, 52, 120. [Google Scholar] [CrossRef]
  12. Bernstein, B.; Kearsley, E.A.; Zapas, L.J. A Study of Stress Relaxation with Finite Strain. Trans. Soc. Rheol. 1963, 7, 391–410. [Google Scholar] [CrossRef]
  13. de Gennes, P.G. Reptation of a Polymer Chain in the Presence of Fixed Obstacles. J. Chem. Phys. 1971, 55, 572–579. [Google Scholar] [CrossRef]
  14. Verbeeten, W.M.H.; Peters, G.W.M.; Baaijens, F.P.T. Differential constitutive equations for polymer melts: The extended Pom–Pom model. J. Rheol. 2001, 45, 823–843. [Google Scholar] [CrossRef] [Green Version]
  15. Tanner, R.I.; Nasseri, S. Simple constitutive models for linear and branched polymers. J. Non-Newton. Fluid Mech. 2003, 116, 1–17. [Google Scholar] [CrossRef]
  16. Tanner, R.I. From A to (BK)Z in Constitutive Relations. J. Rheol. 1988, 32, 673–702. [Google Scholar] [CrossRef]
  17. Wagner, M.H. Analysis of time-dependent non-linear stress-growth data for shear and elongational flow of a low-density branched polyethylene melt. Rheol. Acta 1976, 15, 136–142. [Google Scholar] [CrossRef]
  18. Wagner, M.H. Zur Netzwerktheorie von Polymer-Schmelzen. Rheol. Acta 1979, 18, 33–50. [Google Scholar] [CrossRef]
  19. Mitsoulis, E. Effect of Viscoelasticity in Fountain Flow of Polyethylene Melts. Int. Polym. Process. 2009, 24, 439–451. [Google Scholar] [CrossRef] [Green Version]
  20. Konaganti, V.K.; Ansari, M.; Mitsoulis, E.; Hatzikiriakos, S.G. The effect of damping function on extrudate swell. J. Non-Newton. Fluid Mech. 2016, 236, 73–82. [Google Scholar] [CrossRef]
  21. Ansari, M.; Mitsoulis, E.; Hatzikiriakos, S.G. Capillary Extrusion and Swell of a HDPE Melt Exhibiting Slip. Adv. Polym. Technol. 2013, 32, E369–E385. [Google Scholar] [CrossRef]
  22. Konaganti, V.K.; Derakhshandeh, M.; Ebrahimi, M.; Mitsoulis, E.; Hatzikiriakos, S.G. Non-isothermal extrudate swell. Phys. Fluids 2016, 28, 123101. [Google Scholar] [CrossRef]
  23. Ansari, M.; Zisis, T.; Hatzikiriakos, S.G.; Mitsoulis, E. Capillary flow of low-density polyethylene. Polym. Eng. Sci. 2012, 52, 649–662. [Google Scholar] [CrossRef]
  24. Ansari, M.; Hatzikiriakos, S.G.; Mitsoulis, E. Slip effects in HDPE flows. J. Non-Newton. Fluid Mech. 2011. [Google Scholar] [CrossRef]
  25. Stieger, S.; Kerschbaumer, R.C.; Mitsoulis, E.; Fasching, M.; Berger-Weber, G.R.; Friesenbichler, W.; Sunder, J. Contraction and capillary flow of a carbon black filled rubber compound. Polym. Eng. Sci. 2020, 60, 32–43. [Google Scholar] [CrossRef] [Green Version]
  26. Stieger, S. Development of the Exjection Process for Rubbers Taking into Account the Viscoelastic Material Behavior. Ph.D. Thesis, Montanuniversität Leoben, Leoben, Austria, 2020. [Google Scholar]
  27. ISO 289-1:2015. Rubber, Unvulcanized—Determinations Using a Shearing-Disc Viscometer—Part 1: Determination of Mooney Viscosity, 4th ed.; International Organization for Standardization: Geneva, Switzerland, 2015. [Google Scholar]
  28. Datasheet PERBUNAN® 3965 F; Arlanxeo Deutschland GmbH: Dormagen, Germany, 2019.
  29. Binding, D.M. An approximate analysis for contraction and converging flows. J. Non-Newton. Fluid Mech. 1988, 27, 173–189. [Google Scholar] [CrossRef]
  30. Fasching, M. Robust Processing in Rubber Injection Molding Using Advanced Simulation Methods and Material Data. Ph.D. Thesis, Montanuniversität Leoben, Leoben, Austria, 2015. [Google Scholar]
  31. Hatzikiriakos, S.G. Wall slip of molten polymers. Prog. Polym. Sci. 2012, 37, 624–643. [Google Scholar] [CrossRef]
  32. Mooney, M. Explicit Formulas for Slip and Fluidity. J. Rheol. 1931, 2, 210–222. [Google Scholar] [CrossRef]
  33. Cox, W.P.; Merz, E.H. Correlation of dynamic and steady flow viscosities. J. Polym. Sci. 1958, 28, 619–622. [Google Scholar] [CrossRef]
  34. Bagley, E.B. End Corrections in the Capillary Flow of Polyethylene. J. Appl. Phys. 1957, 28, 624–627. [Google Scholar] [CrossRef]
  35. Maklad, O.; Poole, R.J. A review of the second normal-stress difference; its importance in various flows, measurement techniques, results for various complex fluids and theoretical predictions. J. Non-Newton. Fluid Mech. 2021, 292, 104522. [Google Scholar] [CrossRef]
  36. Sentmanat, M.L. Miniature universal testing platform: From extensional melt rheology to solid-state deformation behavior. Rheol. Acta 2004, 43, 657–669. [Google Scholar] [CrossRef]
  37. Meißner, J. Rheometer zur Untersuchung der deformationsmechanischen Eigenschaften von Kunststoff-Schmelzen unter definierter Zugbeanspruchung. Rheol. Acta 1969, 8, 78–88. [Google Scholar] [CrossRef]
  38. Meissner, J.; Hostettler, J. A new elongational rheometer for polymer melts and other highly viscoelastic liquids. Rheol. Acta 1994, 33, 1–21. [Google Scholar] [CrossRef]
  39. Münstedt, H. New Universal Extensional Rheometer for Polymer Melts. Measurements on a Polystyrene Sample. J. Rheol. 1979, 23, 421–436. [Google Scholar] [CrossRef]
  40. Luo, X.-L.; Tanner, R.I. Finite element simulation of long and short circular die extrusion experiments using integral models. Int. J. Numer. Methods Eng. 1988, 25, 9–22. [Google Scholar] [CrossRef]
Figure 1. (a) Computational fluid dynamic (CFD) simulation setup including boundary conditions (BCs) 1 to 4, dimensions, and position of the pressure transducer (p); (b) Meshes for the tapered orifice (L/D = 0.2/1; ϕ = 90°) and three abrupt capillary dies (L/D = 5/1, 10/1, 20/1; ϕ = 180°) utilized in numerical calculations.
Figure 1. (a) Computational fluid dynamic (CFD) simulation setup including boundary conditions (BCs) 1 to 4, dimensions, and position of the pressure transducer (p); (b) Meshes for the tapered orifice (L/D = 0.2/1; ϕ = 90°) and three abrupt capillary dies (L/D = 5/1, 10/1, 20/1; ϕ = 180°) utilized in numerical calculations.
Polymers 13 02323 g001
Figure 2. Flow chart displaying the research design of the present study.
Figure 2. Flow chart displaying the research design of the present study.
Polymers 13 02323 g002
Figure 3. (a) Effect of applied normal force (Fn) on the sample during the measurement and influence of the measurement devices (RPA vs. MCR 501) on the complex viscosity of gum NBR; (b) No influence of sample preparation detected on the complex viscosity of gum NBR.
Figure 3. (a) Effect of applied normal force (Fn) on the sample during the measurement and influence of the measurement devices (RPA vs. MCR 501) on the complex viscosity of gum NBR; (b) No influence of sample preparation detected on the complex viscosity of gum NBR.
Polymers 13 02323 g003
Figure 4. (a) Time-temperature superposition of the dissipation factor ( tan δ ) and corresponding Arrhenius plot; (b) K-BKZ model predictions of storage G’ and loss G” moduli compared to experimental data at the reference temperature Tref; K-BKZ fit 1 represents LVE data measured at a normal force (Fn) < 1 N.
Figure 4. (a) Time-temperature superposition of the dissipation factor ( tan δ ) and corresponding Arrhenius plot; (b) K-BKZ model predictions of storage G’ and loss G” moduli compared to experimental data at the reference temperature Tref; K-BKZ fit 1 represents LVE data measured at a normal force (Fn) < 1 N.
Polymers 13 02323 g004
Figure 5. Entrance pressure corrected flow curves of two capillary dies with same length-to-diameter (L/D) ratios but different capillary diameters.
Figure 5. Entrance pressure corrected flow curves of two capillary dies with same length-to-diameter (L/D) ratios but different capillary diameters.
Polymers 13 02323 g005
Figure 6. Measured pressure drops for various apparent shear rates ( γ ˙ a ) in dependence of the capillary die length (Bagley plot) including linear approximations (Lin. approx.) with coefficients of determination (R²).
Figure 6. Measured pressure drops for various apparent shear rates ( γ ˙ a ) in dependence of the capillary die length (Bagley plot) including linear approximations (Lin. approx.) with coefficients of determination (R²).
Polymers 13 02323 g006
Figure 7. Steady-state shear (η) and uniaxial elongational ( η e ) viscosities of gum NBR compared to K-BKZ model predictions; K-BKZ fit 1 represents fitted LVE data measured at Fn < 1 N ( θ = 0), K-BKZ fit 2 represents fitted LVE data measured at Fn = 10 N ( θ = 0), and K-BKZ fit 3 represents fitted LVE data measured at Fn < 1 N ( θ = −0.15).
Figure 7. Steady-state shear (η) and uniaxial elongational ( η e ) viscosities of gum NBR compared to K-BKZ model predictions; K-BKZ fit 1 represents fitted LVE data measured at Fn < 1 N ( θ = 0), K-BKZ fit 2 represents fitted LVE data measured at Fn = 10 N ( θ = 0), and K-BKZ fit 3 represents fitted LVE data measured at Fn < 1 N ( θ = −0.15).
Polymers 13 02323 g007
Figure 8. Measured pressure data compared to CFD simulation results; K-BKZ fit 1 represents fitted LVE data measured at Fn < 1 N ( θ = 0), K-BKZ fit 2 represents fitted LVE data measured at Fn = 10 N ( θ = 0), and K-BKZ fit 3 represents fitted LVE data measured at Fn < 1 N ( θ = −0.15).
Figure 8. Measured pressure data compared to CFD simulation results; K-BKZ fit 1 represents fitted LVE data measured at Fn < 1 N ( θ = 0), K-BKZ fit 2 represents fitted LVE data measured at Fn = 10 N ( θ = 0), and K-BKZ fit 3 represents fitted LVE data measured at Fn < 1 N ( θ = −0.15).
Polymers 13 02323 g008
Figure 9. Measured pressure data compared to CFD simulation results of the (a) tapered orifice die (L/D = 0.2/1); (b) abrupt capillary die with a length-to-diameter (L/D) ratio of 5/1; (c) abrupt capillary die with L/D = 10/1; (d) abrupt capillary die with L/D = 20/1.
Figure 9. Measured pressure data compared to CFD simulation results of the (a) tapered orifice die (L/D = 0.2/1); (b) abrupt capillary die with a length-to-diameter (L/D) ratio of 5/1; (c) abrupt capillary die with L/D = 10/1; (d) abrupt capillary die with L/D = 20/1.
Polymers 13 02323 g009
Table 1. Boundary conditions (BCs) applied in CFD simulations.
Table 1. Boundary conditions (BCs) applied in CFD simulations.
BCDescription
BC 1Fully developed velocity profile (inlet)
BC 2Normal and tangential velocities are zero (no slip at the wall)
BC 3Tangential force and normal velocity are zero (axis of symmetry)
BC 4Normal force and tangential velocity are zero (viscous outlet)
BC 4Normal and tangential forces are zero “zero force BC” (viscoelastic outlet)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Stieger, S.; Mitsoulis, E.; Walluch, M.; Ebner, C.; Kerschbaumer, R.C.; Haselmann, M.; Mostafaiyan, M.; Kämpfe, M.; Kühnert, I.; Wießner, S.; et al. On the Influence of Viscoelastic Modeling in Fluid Flow Simulations of Gum Acrylonitrile Butadiene Rubber. Polymers 2021, 13, 2323. https://doi.org/10.3390/polym13142323

AMA Style

Stieger S, Mitsoulis E, Walluch M, Ebner C, Kerschbaumer RC, Haselmann M, Mostafaiyan M, Kämpfe M, Kühnert I, Wießner S, et al. On the Influence of Viscoelastic Modeling in Fluid Flow Simulations of Gum Acrylonitrile Butadiene Rubber. Polymers. 2021; 13(14):2323. https://doi.org/10.3390/polym13142323

Chicago/Turabian Style

Stieger, Sebastian, Evan Mitsoulis, Matthias Walluch, Catharina Ebner, Roman Christopher Kerschbaumer, Matthias Haselmann, Mehdi Mostafaiyan, Markus Kämpfe, Ines Kühnert, Sven Wießner, and et al. 2021. "On the Influence of Viscoelastic Modeling in Fluid Flow Simulations of Gum Acrylonitrile Butadiene Rubber" Polymers 13, no. 14: 2323. https://doi.org/10.3390/polym13142323

APA Style

Stieger, S., Mitsoulis, E., Walluch, M., Ebner, C., Kerschbaumer, R. C., Haselmann, M., Mostafaiyan, M., Kämpfe, M., Kühnert, I., Wießner, S., & Friesenbichler, W. (2021). On the Influence of Viscoelastic Modeling in Fluid Flow Simulations of Gum Acrylonitrile Butadiene Rubber. Polymers, 13(14), 2323. https://doi.org/10.3390/polym13142323

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop