Next Article in Journal
Effect of Photoinitiator on Precursory Stability and Curing Depth of Thiol-Ene Clickable Gelatin
Next Article in Special Issue
Simulation Approach for Hydrophobicity Replication via Injection Molding
Previous Article in Journal
Study of Polyvinyl Alcohol-SiO2 Nanoparticles Polymeric Membrane in Wastewater Treatment Containing Zinc Ions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Experimental Investigation of Viscoelastic Flow in a Contraction Channel

Key Laboratory of Rubber-Plastics, Ministry of Education/Shandong Provincial Key Laboratory of Rubber-Plastics, Qingdao University of Science and Technology, Qingdao 266042, China
*
Author to whom correspondence should be addressed.
Polymers 2021, 13(11), 1876; https://doi.org/10.3390/polym13111876
Submission received: 10 May 2021 / Revised: 31 May 2021 / Accepted: 2 June 2021 / Published: 4 June 2021

Abstract

:
In order to assess the predictive capability of the S–MDCPP model, which may describe the viscoelastic behavior of the low-density polyethylene melts, a planar contraction flow benchmark problem is calculated in this investigation. A pressure-stabilized iterative fractional step algorithm based on the finite increment calculus (FIC) method is adopted to overcome oscillations of the pressure field due to the incompressibility of fluids. The discrete elastic viscous stress splitting (DEVSS) technique in combination with the streamline upwind Petrov-Galerkin (SUPG) method are employed to calculate the viscoelastic flow. The equal low-order finite elements interpolation approximations for velocity-pressure-stress variables can be applied to calculate the viscoelastic contraction flows for LDPE melts. The predicted velocities agree well with the experimental results of particle imagine velocity (PIV) method, and the pattern of principal stress difference calculated by the S-MDCPP model has good agreement with the results measured by the flow induced birefringence (FIB) device. Numerical and experimental results show that the S-MDCPP model is capable of accurately capturing the rheological behaviors of branched polymers in complex flow.

Graphical Abstract

1. Introduction

During the past two decades with the developments of numerical algorithm and constitutive model, computational rheology has achieved great progress as an important tool. Many investigators have utilized computational rheology to obtain insight into the complex rheological behaviors of polymers, so as to deeply understand the flow status [1] and processing mechanisms, such as extrusion swell [2,3], flow instabilities [4] and secondary flow.
In order to explore the complex behavior of viscoelastic flow, the predictive capability of the constitutive model is of significance for numerical simulation. As for polymer melts, many rheologists put forward several constitutive models, such as the K-BKZ model, PTT model, and Gieseku model. Notably, for better understanding the rheological behaviors of branched polymer melts, McLeish and Larson [5] first proposed the Pom-Pom model that originated from tube theory [6]. The main contribution is to describe the rheological responses of macromolecules in terms of two independent relaxation timescales for orientation and stretch, respectively. After that, several researchers [7,8,9,10,11] developed the Pom-Pom model to remove some defects, so as to reproduce the real rheological behaviors of branched polymer melts. With the help of these constitutive models, we can acquire the orientation and stretch behaviors and reveal the viscoelastic flow mechanism. To validate the numerical results and assess the performance of the constitutive model, some experiments [12,13] are performed. Verbeeten et al. [14,15] calculated viscoelastic flows using the XPP model and obtained good agreement with the experimental results. Wang et al. [16] evaluated the performance of the double-convected Pom-Pom mode (DCPP) model using the particle imagine velocity (PIV) method. They observed that the LDPE melt appears as large corner vortex in the contraction entrance, which has good agreement with the numerical results predicted by Polyflow software. However, they gained less satisfactory for HDPE melt using the same method.
Numerous studies [1,13,15,17] were carried out both experimentally and numerically for contraction flow, which is a benchmark problem in computational rheology, because there exists a complex combination flow of both shear flow near the channel wall and extension flow along the centerline in the contraction channel. In practice, similar contraction flows also widely appear in industrial polymer processing, such as extrusion flow, injection molding, fiber spinning process, etc. By means of this complex flow field, the rheological responses of the constitutive model may be computed in a numerical way and verified by stress and velocity fields using the birefringence experiment and PIV measurement respectively. Aguayo et al. [18,19] calculated a 4:1 planar rounded-corner and sharp contraction flows governed by the SXPP model and indicated the influence of the number of dangling arms at the end of the pom-pom molecule on the flow behaviors with the aid of the hybrid finite-element/volume scheme. Clemeur et al. [13] exhibited the good performance of the DCPP model with the simulation and experiment results of contraction flow.
In order to further reveal the rheological behavior of the S-MDCPP model [10], we performed experiments using the 4:1 contraction flow channel for LDPE melt and carried out a numerical study for the extrusion flow. Moreover, to verify the numerical results, the PIV system and flow induced birefringence (FIB) device respectively are utilized to measure the velocity and the principal stress difference of the flow field in this work.
The objective of this investigation is to evaluate the predictive capability of the S-MDCPP model. The 4:1 planar abrupt contraction flow is studied numerically and experimentally, respectively. This paper is organized as follows: Section 2 introduces the material features and experimental set-up. Furthermore, the computational method and model description are also briefly depicted in Section 2. The computational and experimental results together with discussion are described in Section 3. The conclusions of this investigation are drawn in Section 4.

2. Materials and Methods

2.1. Materials

Low-density polyethylene (DSM Stamylan LD 2008 XC43) was used in the present study. The experimental temperature of the LDPE was 170 °C. See Table 1 for the main properties of the LDPE melt. The shear and uniaxial extension viscosities behaviors of the LDPE melts with various shear rates and tensile rates were depicted in our previous study [20]. The viscoelastic rheological behavior of the LDPE melts is described by the S-MDCPP model [10]. Hence, the model parameters of the S-MDCPP model are solely listed in Table 2 [20].

2.2. Experimental Set-Up

(1) Flow Induced Birefringence (FIB) System
The optics measurement system consists of two line polarizers, two quarter wave polymer retarders and a high-speed CCD camera. The schematic illustration of the FIB system is depicted in Figure 1. A laser device made by Beijing Laserwave Optoelectronics Technology Co., Ltd. is utilized to produce green light of wavelength λ = 532 nm. The flow channel device shown in Figure 2 has a width of 8 mm and height of 45 mm in the upstream channel and width of 2 mm and height of 45 mm in the downstream channel, and the thickness of the whole flow channel is 100 mm.
(2) Particle imagine velocity (PIV) system
The PIV measurement system made by Beijing Cubetiandi Science and Technology Development Co. Ltd. (Beijing, China) is used to obtain the velocities of the contraction flow field. The PIV system composes of a high-speed CCD camera, a laser with green light of wavelength λ = 532 nm and imaging analysis software.

2.3. Computational Method and Model Description

In this investigation, the pressure stabilized iterative fractional step algorithm adopted in our previous studies [10,20,21,22] will be applied to calculate the LDPE melts flow behavior in the contraction channel. It is worth noting that the finite increment calculus (FIC) method [23,24] is employed to reformulate the stabilized form of mass conservation equation to achieve the stable solution of the pressure field. Accordingly, the equal low-order triangle elements may be applied to discretize the flow field. Meanwhile, the discrete elastic viscous stress splitting (DEVSS) method of Guénette and Fortin [25] in combination with the streamline upwind Petrov-Galerkin (SUPG) approach [26] is applied to calculate the viscoelastic flow. The computational method is described more fully in our previous studies [10,20,21].
In the LDPE melt extrusion process, when the melt leaves the extruder die, it is squeezed into a home-made flow channel with two glass windows, as shown in Figure 2. The flow device consists of three plates, which create a 4:1 contraction channel. Because the ratio of thickness (depth) to width is above 10, the three-dimensional channel may be simplified to a planar model. So, a planar 4:1 contraction geometry is constructed in this study for calculating the viscoelastic flow. Because the contraction flow as a benchmark problem [1,13,15,17] is usually used to evaluate the capability of constitutive models and numerical algorithms. Accordingly, the planar 4:1 contraction viscoelastic flow is also calculated in this study. The finite element mesh model is displayed in Figure 3. The velocities at entrance and exit are prescribed by a fully developed velocity profile, respectively. The average velocity of the entrance is 1.275 mm/s. The vertical velocity at the centerline of the channel is set as zero. The no-slip boundary conditions are considered on the channel wall. The pressure of the symmetric point at the exit is specified zero. All components of the viscoelastic extra stress are initially applied by the computational results achieved from the UCM fluid which possesses the same Maxwell parameters as those used for the S–MDCPP model.

3. Results and Discussion

Figure 4 shows the sochromatic fringe patterns of the principal stress difference (PSD), that is,   PSD = τ x x τ y y 2 + 4 τ x y 2 . It can be seen that our experimental result agrees well with the PSD patterns measured by Verbeeten [15,27]. Because the channel contraction ratio of the experimental device used by Verbeeten is 3.29:1, but the contraction ratio of this study is 4:1, the size of the PSD pattern has a slight difference in the contraction region. However, the number and location of the PSD field pattern are consistent with the results reported by Verbeeten. Moreover, the computational results using our calculation scheme have good quantitative agreement with the measured PSD pattern by means of our FIB device, so the stability and accuracy of the computational method adopted in this study are verified. By comparison between numerical and experimental results, we observe the S-MDCPP model has good predictive capability and may reproduce the complex rheological behavior of branched polymer melts. The vortex behavior and butterfly pattern of PSD are well-predicted by the S-MDCPP model, which is the typical features of LDPE melts flow through the contraction channel. During the branched polymer melts flow into the contraction entry, the stronger shear stress occurs near the channel wall, while the higher extensional behavior appears along the centerline, which results in high extensional viscosity. These complex flow behaviors induce shear thinning and extension thickening (or strain hardening) of branched macromolecules occurring simultaneously in the contraction channel, which demonstrates the complicated butterfly pattern of PSD.
Figure 5 displays the horizontal velocity of the upstream channel in the cross section of x = 13 mm. It can be seen that the predicted velocities agree very well with the results measured by the PIV method except for the neighborhood of the channel wall. This indicates the pressure stabilized iterative fractional step algorithm adopted in this study is reliable, although the equal low-order triangle elements are applied to calculate the velocity field. The distributions of sectional velocities display a parabolic curve.
The pressure distribution in the viscoelastic contraction flow is shown in Figure 6. We can see the pressure field is stable, although the equal low-order triangle elements are applied to calculate the pressure, velocity, viscoelastic extra-stress tensor and discrete rate-of-strain tensor. This good stability feature attributes to the FIC approach introduced a pressure stabilization term.
The horizontal velocity and streamline are exhibited in Figure 7. We note that the horizontal velocity near the centerline of the downstream channel is larger than other locations. Moreover, we observe a large rotating vortex appears in the salient corner of the upstream channel (See Figure 7b), which is the typical feature of the LDPE melts flowing through a contraction channel. The corner vortex behavior is in agreement with the results of several researchers [1,13,15].
The components (i.e., τxx, τyy, τxy) of the viscoelastic extra-stress tensor are shown in Figure 8, Figure 9 and Figure 10, respectively. We may observe the viscoelastic extra-stress field is stable. The largest magnitude of the stresses is located in the sharp re-entrant corner. In addition, the stresses of τxx and τxy are larger in the vicinity of the solid wall, where the shear rates are higher than other regions of the channel.
The distribution of first normal stress difference N1 ( N 1 = τ x x τ y y ) is shown in Figure 11. It is observed that the distribution of N1 near the sharp corner is complex. Meanwhile, the largest first normal stress difference is located in this corner.
The backbone stretch distribution of the branched macromolecules is shown in Figure 12. We may observe the larger backbone stretch locates in the vicinity of the contraction sharp corner (See the inset in Figure 12), and the largest backbone stretch of 3.04 occurs in the neighborhood of the re-entrant corner near the downstream wall, where there is the maximum shear rate, as shown in Figure 13. The very high shear rate induces large stretch orientation of the macromolecules.
From Figure 14, we can see the backbone stretch closes approximately to a value of 1, when macromolecules reach far away from the contraction entrance wherever the upstream or downstream. However, when macromolecules flow through the contraction channel, because the larger axial velocity induces stronger orientation of the macromolecule, the larger backbone stretch appears in the vicinity of the contraction entrance. After that, the macromolecules begin to relax toward recovering the random coil status.

4. Conclusions

The viscoelastic flow of branched polymer melts (LDPE) in a 4:1 contraction flow channel was studied by using the FIB and PIV experiment devices, along with the pressured stabilized iterative fractional step algorithm. Firstly, the numerical results demonstrate that the S-MDCPP model is capable of reproducing the rheological behaviors of the branched polymers in the complex flow. Secondly, the results show that the predicted PSD pattern agrees well with the experimental results obtained by the FIB device, and the velocities calculated by the numerical scheme have good agreement with the PIV measured results. Consequently, the viscoelastic constitutive model (e.g., S-MDCPP model) originated from tube theory may capture the complex rheological behaviors of polymers. From the point of numerical result about the stretch behavior, the larger backbone stretch locates in the vicinity of the contract sharp corner, and the largest backbone stretch occurs in the neighborhood of the re-entrant corner near the downstream wall, where the maximum shear rate occurs, mainly because the high shear rate induces a large stretch orientation of the macromolecules. Based on the above results, we consider that the S-MDCPP model may be a good candidate for complex viscoelastic flow simulation in extrusion and injection molding in the future. The reliable and accurate numerical prediction may help us to reveal and deeply understand the complex rheological behaviors of polymers.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/polym13111876/s1, Table S1: Shear viscosities under different shear rates, Table S2: Uniaxial viscosities under different ex-tension rates.

Author Contributions

Conceptualization, W.W.; methodology, W.W.; software, W.W.; validation, W.W., and L.W.; investigation, W.W.; resources, W.W.; data curation, L.W.; writing—original draft preparation, W.W.; writing—review and editing, W.W.; visualization, L.W.; supervision, W.W.; project administration, W.W.; funding acquisition, W.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Shandong Provincial Natural Science Foundation, grant number ZR2018MEM022 and National Natural Science Foundation of China through Contract/Grant number 21274072.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hooshyar, S.; Germann, N. Shear banding in 4: 1 planar contraction. Polymers 2019, 11, 417. [Google Scholar] [CrossRef] [Green Version]
  2. Hooshyar, S.; Germann, N. The investigation of shear banding polymer solutions in die extrusion geometry. J. Non-Newton. Fluid Mech. 2019, 272, 104161. [Google Scholar] [CrossRef]
  3. Ganvir, V.; Gautham, B.P.; Pol, H.; Bhamla, M.S.; Sclesi, L.; Thaokar, R.; Lele, A.; Mackley, M. Extrudate swell of linear and branched polyethylenes: ALE simulations and comparison with experiments. J. Non-Newton. Fluid Mech. 2011, 166, 12–24. [Google Scholar] [CrossRef]
  4. Pol, H.V.; Joshi, Y.M.; Tapadia, P.S.; Lele, A.K.; Mashelkar, R.A. A Geometrical Solution to the Sharkskin Instability. Ind. Eng. Chem. Res. 2007, 46, 3048–3056. [Google Scholar] [CrossRef]
  5. 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]
  6. McLeish, T.C.B. Tube theory of entangled polymer dynamics. Adv. Phys. 2002, 51, 1379–1527. [Google Scholar] [CrossRef]
  7. 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–844. [Google Scholar] [CrossRef] [Green Version]
  8. Clemeur, N.; Rutgers, R.P.G.; Debbaut, B. On the evaluation of some differential formulations for the Pom-Pom constitutive model. Rheol. Acta. 2003, 42, 217–231. [Google Scholar] [CrossRef]
  9. 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]
  10. Wang, W.; Li, X.K.; Han, X.H. A numerical study of constitutive models endowed with Pom-Pom molecular attributes. J. Non-Newton. Fluid Mech. 2010, 165, 1480–1493. [Google Scholar] [CrossRef]
  11. Clemeur, N.; Debbaut, B. A pragmatic approach for deriving constitutive equations endowed with pom–pom attributes. Rheol. Acta. 2007, 46, 1187–1196. [Google Scholar] [CrossRef]
  12. Marín-Santibáñez, B.M.; Pérez-González, J.; Gómez-Herrera, G.; Rodríguez-González, F. Capillary extrusion of polypropylene/high-density polyethylene immiscible blends as studied by rheo-particle image velocimetry. Polym. Test. 2020, 84, 106390. [Google Scholar] [CrossRef]
  13. Clemeur, N.; Rutgers, R.P.G.; Debbaut, B. Numerical simulation of abrupt contraction flows using the Double Convected Pom–Pom model. J. Non-Newton. Fluid Mech. 2004, 117, 193–209. [Google Scholar] [CrossRef]
  14. Verbeeten, W.M.H.; Peters, G.W.M.; Baaijens, F.P.T. Viscoelastic Analysis of Complex Polymer Melt Flows using the eXtended Pom-Pom model. J. Non-Newton. Fluid Mech. 2002, 108, 301–326. [Google Scholar] [CrossRef]
  15. Verbeeten, W.M.H.; Peters, G.W.M.; Baaijens, F.P.T. Numerical simulations of the planar contraction flow for a polyethylene melt using the XPP model. J. Non-Newton. Fluid Mech. 2004, 117, 73–84. [Google Scholar] [CrossRef]
  16. Wang, X.L.; Chen, R.H.; Wang, M.M.; Jin, G. Validation of Double Convected Pom-Pom Model With Particle Image Velocimetry Technique. Polym. Eng. Sci. 2015, 55, 1897–1905. [Google Scholar] [CrossRef]
  17. Owens, R.G.; Phillips, T.N. Computational Rheology; Imperial College Press: London, UK, 2002; pp. 201–246. [Google Scholar]
  18. Aguayo, J.P.; Tamaddon-Jahromi, H.R.; Webster, M.F. Extensional response of the pom-pom model through planar contraction flows for branched polymer melts. J. Non-Newton. Fluid Mech. 2006, 134, 105–126. [Google Scholar] [CrossRef]
  19. Aguayo, J.P.; Phillips, P.M.; Phillips, T.N.; Tamaddon-Jahromi, H.R.; Snigerev, B.A.; Webster, M.F. The numerical prediction of planar viscoelastic contraction flows using the pom-pom model and higher-order finite volume schemes. J. Comput. Phys. 2007, 220, 586–611. [Google Scholar] [CrossRef]
  20. Wang, W.; Wang, X.P.; Hu, C.X. A comparative study of viscoelastic planar contraction flow for polymer melts using molecular constitutive models. Korea-Aust. Rheol. J. 2014, 26, 365–375. [Google Scholar] [CrossRef]
  21. Wang, W.; Hu, C.X.; Li, W.W. Time-dependent rheological behavior of branched polymer melts in extensional flows. Mech. Time-Depend. Mater. 2016, 20, 123–137. [Google Scholar] [CrossRef]
  22. Li, X.K.; Han, X.H. Numerical modeling of viscoelastic flows using equal low-order finite elements. Comput. Meth. Appl. Mech. Engrg. 2010, 199, 570–581. [Google Scholar] [CrossRef]
  23. Oñate, E. A stabilized finite element method for incompressible viscous flows using a finite increment calculus formulation. Comput. Meth. Appl. Mech. Engrg. 2000, 182, 355–370. [Google Scholar] [CrossRef] [Green Version]
  24. Li, X.K.; Han, X.H. An iterative stabilized fractional step algorithm for numerical solution of incompressible N–S equations. Int. J. Numer. Meth. Fluids. 2005, 49, 395–416. [Google Scholar] [CrossRef]
  25. Guénette, R.; Fortin, M. A new mixed finite element method for computing viscoelastic flows. J. Non-Newton. Fluid Mech. 1995, 60, 27–52. [Google Scholar] [CrossRef]
  26. Brooks, A.N.; Hughes, T.J.R. Streamline upwind/Petrov–Galerkin methods for convection dominated flows with particular emphasis on the incompressible Navier–Stokes equations. Comput. Meth. Appl. Mech. Engrg. 1982, 32, 199–259. [Google Scholar] [CrossRef]
  27. Verbeeten, W.M.H. Computational Polymer Melt Rheology. Ph.D. Thesis, Technische Universiteit Eindhoven, Eindhoven, The Netherlands, 2001. [Google Scholar]
Figure 1. Schematic illustration of the FIB system.
Figure 1. Schematic illustration of the FIB system.
Polymers 13 01876 g001
Figure 2. Flow channel set-up.
Figure 2. Flow channel set-up.
Polymers 13 01876 g002
Figure 3. Schematic and computational mesh around the re-entrant corner of 4:1 planar contraction.
Figure 3. Schematic and computational mesh around the re-entrant corner of 4:1 planar contraction.
Polymers 13 01876 g003
Figure 4. Comparison of the isochromatic fringe patterns measured by Verbeeten [27] and our experimental result (a), our experimental and computational results (b).
Figure 4. Comparison of the isochromatic fringe patterns measured by Verbeeten [27] and our experimental result (a), our experimental and computational results (b).
Polymers 13 01876 g004
Figure 5. Comparison of the PIV measured velocity and predicted velocity in the upstream cross section of x = 13 mm.
Figure 5. Comparison of the PIV measured velocity and predicted velocity in the upstream cross section of x = 13 mm.
Polymers 13 01876 g005
Figure 6. Distribution of pressure in the flow field.
Figure 6. Distribution of pressure in the flow field.
Polymers 13 01876 g006
Figure 7. Distributions of horizontal velocity (a) and streamline (b) in the flow field.
Figure 7. Distributions of horizontal velocity (a) and streamline (b) in the flow field.
Polymers 13 01876 g007
Figure 8. Distribution of stress τxx in the flow field near the contraction region.
Figure 8. Distribution of stress τxx in the flow field near the contraction region.
Polymers 13 01876 g008
Figure 9. Distribution of stress τyy in the flow field near the contraction region.
Figure 9. Distribution of stress τyy in the flow field near the contraction region.
Polymers 13 01876 g009
Figure 10. Distribution of stress τxy in the flow field near the contraction region.
Figure 10. Distribution of stress τxy in the flow field near the contraction region.
Polymers 13 01876 g010
Figure 11. Distribution of the first normal stress difference N1 in the flow field near the contraction region.
Figure 11. Distribution of the first normal stress difference N1 in the flow field near the contraction region.
Polymers 13 01876 g011
Figure 12. Distribution of backbone stretch Λ in the flow field near the contraction region.
Figure 12. Distribution of backbone stretch Λ in the flow field near the contraction region.
Polymers 13 01876 g012
Figure 13. Distribution of shear rate in the flow field near the contraction corner.
Figure 13. Distribution of shear rate in the flow field near the contraction corner.
Polymers 13 01876 g013
Figure 14. Distribution of backbone stretch along the axis of symmetry.
Figure 14. Distribution of backbone stretch along the axis of symmetry.
Polymers 13 01876 g014
Table 1. Materials properties of the experiment.
Table 1. Materials properties of the experiment.
Density
(kg/m3)
MwMw/MnMelt Index (dg/min)
LDPE (Stamylan LD 2008 XC43)920155,00011.928
Table 2. Maxwell and nonlinear parameters of DSM Stamylan LD 2008 XC43 LDPE melt [20].
Table 2. Maxwell and nonlinear parameters of DSM Stamylan LD 2008 XC43 LDPE melt [20].
S–MDCPP Model
ModeMaxwell ParametersNonlinear Parameters
G (Pa)λ0b (s)q = 2/νrξ
11.0377 × 1040.40101.20.02
Note that q represents the number of arms at the end of a backbone, r equals λ0bs, where λ0b and λs denote orientation and stretch relaxation time, respectively, and ξ denotes the slip coefficient, which introduces the non-zero second normal stress difference as adopted in the DCPP model. Moreover, ν reflects the influence of the neighborhood polymer chains on the backbone tube stretch. In Table 2, r and ξ are dimensionless parameters.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, W.; Wang, L. An Experimental Investigation of Viscoelastic Flow in a Contraction Channel. Polymers 2021, 13, 1876. https://doi.org/10.3390/polym13111876

AMA Style

Wang W, Wang L. An Experimental Investigation of Viscoelastic Flow in a Contraction Channel. Polymers. 2021; 13(11):1876. https://doi.org/10.3390/polym13111876

Chicago/Turabian Style

Wang, Wei, and Linlin Wang. 2021. "An Experimental Investigation of Viscoelastic Flow in a Contraction Channel" Polymers 13, no. 11: 1876. https://doi.org/10.3390/polym13111876

APA Style

Wang, W., & Wang, L. (2021). An Experimental Investigation of Viscoelastic Flow in a Contraction Channel. Polymers, 13(11), 1876. https://doi.org/10.3390/polym13111876

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