The elongational flow behavior of polyethylene, polypropylene, polystyrene, poly(methyl methacrylate), and polycarbonate, temperatures from 70 to 290 °C and pressures up to 70 MPa, is examined with the Yahsi-Dinc-Tav (YDT) model and its particular case known as the Cross model. The viscosity data employed in the range of 3-405 elongational rates were acquired from the literature at ambient and elevated pressures. The predictions and the fitting results of the proposed YDT model with the same measurement data are compared with the Cross model. The average absolute deviations of the viscosities predicted by the YDT model range from 0.54% to 9.44% at ambient and 1.95% to 6.28% at high pressures. Additionally, the linear formulations derived from the YDT model are employed to relate the viscosity with temperature and hole fraction (“thermooccupancy” function) at zero level of elongational rate and constant elongational rate along with constant elongational stress. The effects of the four viscosity parameters (such as transmission and activation energy coefficients in these equations) on the elongational viscosity are analyzed in detail and some conclusions on the structural differences for the polymers are discussed.

1. Introduction

In rheology, flow is mostly a mixture of elongation and shear. Elongational dominated, shear dominated, in rare cases pure shear, or purely elongational flow is the nature of this phenomenon. Shear flow is typically used to characterize flow behavior since it is easily produced under laboratory conditions. Elongational viscosity which is known to generate strong stress responses also plays a crucial role for material characterization [1]. In some polymer melt processing only elongational properties can identify the polymeric materials having the same shear viscosity properties. In various manufacturing processes such as film blowing, fiber spinning, film casting, etc., elongational viscosity gives more accurate knowledge for polymer processibility as well [2]. The important attempt to reveal the underlying mechanisms behind these applications is modeling.

Several Newtonian and non-Newtonian models are present in the literature. However, none of the existing Newtonian equations is able to predict well the elongational flow. Unlike Newtonian equations, non-Newtonian models have a capability to more correctly represent nonlinear behavior of the polymer melts. Those most frequently applied to elongational viscosity is Sarkar-Gupta elongational viscosity model [3] and a few-parametric modified White–Metzner model (mWM) [4, 5]. Gupta et al. [3] fitted their shear viscosity values in terms of the Carreau model for the shear-rate dependence. For elongational data, they introduced an equation which is a modificated form of the Carreau model to capture the strain hardening effect as it is in LDPE polymer [1, 3]. Several authors, however, have attempted to develop composite models combining the shear and elongational viscosities data fitting applications [4, 5]. Zatloukal employed the proposed modified Leonov, the eXtended Pom-Pom and the modification of White–Metzner equations for both shear and elongational flows of LDPE, mLLDPE, and PVB melts. He tested a new modification of the Leonov model with a proposed dissipation term and showed that the model results are very good fitting for all flow situations. In his study, while the XPP model accurately describes the elongational viscosity up to elongational rate about 10 , at higher elongational rates, poor estimation is accomplished. In a similar manner, the elongational flow, in this paper, has been represented by the YDT (Yahsi-Dinc-Tav) model which has been successfully employed for the shear viscosity previously [6]. The model works successfully on both temperature- and pressure-dependent zero level of viscosities and critical stress parameter with a single fitting parameter. Analyses were carried out on HDPE, LDPE, LLDPE, PP, PC, PMMA, and PS melts at both low (3-50 ) and high (up to about 405 ) elongational rates in a wide range of T and P data.

The aim of the present paper as the first part of the work is to evaluate the suitability of the YDT equation in modeling of elongational flow. The fitting/predicting performance of the proposed YDT model is compared with the Cross-like model and better results were obtained for the former one. The viscosity data and the parameters determined from elongational viscosity were analyzed and correlated in terms of temperature and pressure by the YDT model.

In recent years there has been a considerable interest in free volume concept as a means to describe the mechanical behavior of polymer materials, particularly viscosity [614]. The term “free volume”, which is expressed as an additional space needed for the thermal motion, is first introduced in Van der Waals’ thesis [15]. More than a century ago, Batschinski [16] suggested the fact that the free volume is more decisive than the temperature for determining the viscosity of low-molecular weight liquids. This approach was further elaborated on polymers in the 1950s by Doolittle [17, 18]. In the late 1960s, Simha and Somcynsky derived an equation of state in their lattice-hole theory of liquids and introduced the hole fraction quantity as a measure of free volume [19]. Utracki presented an inverse relation of the logarithm of zero shear viscosity with hole fraction using a modified functional form of Doolittle’s equation [710]. Sedlacek et al. [11] applied their modified Utracki’s hole fraction model to linearize the zero shear and constant shear-stress viscosities data. In recent articles, the YDT model by Dinc et al. [6, 13, 14] offered a linear relation between shear-rate/shear-stress viscosities and temperature- and hole fraction-dependent thermooccupancy function for various polymers. Although to date several derived models have been developed to account for the relation of shear viscosity to the free volume [614], almost no studies have been carried out for the elongational viscosity.

In this work, as a second part, we focused on the effect of hole fraction on elongational flow. The dependence of elongational viscosity on hole fraction is analyzed in terms of linear equations derived from the YDT model at zero level of elongational rate, constant elongational rate, and constant elongational stress. From these derivations we also obtained the viscosity parameters, e.g., transmission coefficient and a measure of activation energy coefficient, and related them with constant elongational rate and constant elongational stress.

2. Theories

2.1. The Simha–Somcynsky (SS) Lattice-Hole Theory

The Simha and Somcynsky (SS) [19] lattice-hole theory developed in 1969 describes the pressure, volume, and temperature (PVT) behavior of polymer melts and particularly introduces the term occupied site fraction, . The temperature- and pressure-dependent occupied lattice site fraction parameter, y, expressed as the ratio of the occupied sites on all of lattice sites, and the complementary hole fraction, h, are given by the following equation:where is the number of segments in a molecule and and are the number of molecules and holes, respectively.

The model includes P, V, T variables in a scaled form (, , and ), namely, , , and with the characteristic scaling parameters, , , and :where is the Avogadro’s number, is the molar mass of a segment, and k is the Boltzmann’s constant. and are, respectively, the characteristic molar volume and attractive energy parameter per segment. is the number of the first neighbor intermolecular pairs of the s-mer (chain length of a polymer) in a lattice of the coordination number . The theory incorporates a measure of chain flexibility expressed as the ratio of the external degree of freedom of the molecule, 3c, to the number of segments per molecule, s.

Using reduced variables , , and hole fraction parameter, h, the SS provides two coupled equations:where .

The characteristic scaling parameters, , , as well as temperature- and pressure-dependent hole fraction h(P,T) parameter can be determined by superimposing experimental PVT data on the theoretical surface with the use of (3) and (4). The model is advantageous in that it provides the complementary hole fraction, h, which has been correlated with transport properties and explains the mechanical behavior of polymers such as viscosity [6, 10], ionic conductivity [20], etc.

2.2. The Yahsi-Dinc-Tav (YDT) Model

We have shown previously that the Yahsi-Dinc-Tav (YDT) non-Newtonian viscosity model fitted well the shear viscosity data particularly at high pressure [6]. For the bulk structured polymers like PMMA, PC, and PS of which the viscosity values are high comparing the remaining polymers, the outcomes were found reasonable as well. This result led us to test the YDT model for the high viscosity values as in elongational flow. Hence, one of the aims in the present study is primarily not to observe the results, but to prove our foresight about the suitability of the YDT model for predicting elongational viscosity.

Here, elongational behavior of the same polymers is going to be modeled with the same equation which mathematically relates Newtonian viscosity (zero level of elongational viscosity), , and elongational rate, as :Neglecting 2q power term in the denominator of (5), we obtain the Cross-like model [6, 13, 14];where the nonlinearity parameter is the order of the fractional derivative of the exponential function in the Eyring’s strain rate [6] and denotes the gamma function. , is a measure of the characteristic ratio, involved the angle, , between the deformation stress and the displacement vector of the segment jumping into the neighboring th lattice position. n denotes the sequential order of the fractional derivative and the fractional integral in the fractional Taylor series expansion of the first exponential product in the Eyring’s rate of strain in (3) of [6]. is the critical stress parameter in which and are gas constant and the molar volume of a segment, respectively.

2.3. The Viscosity Model in Terms of Thermooccupancy Function

A linear relation between the viscosities and thermooccupancy function, , is [6, 1214]where is zero level of elongational viscosity (Newtonian viscosity) (), constant elongational rate viscosity (), or constant elongational stress viscosity (). depends on the structural occupancy-nonoccupancy ratio () and inversely related to the temperature. The viscosity-intercept (for the zero level of elongational rate, constant elongational rate, and constant elongational stress; , , or , respectively) represents the extrapolated viscosity values at zero (i.e., when the system is dominated by large hole fraction at high temperature and low pressure). The slope refers to a measure of activation energy at zero level of elongational rate (), constant elongational rate (), or constant elongational stress (). The parameters and are given bywhere and are Planck and Boltzmann’s constants, respectively. κ is the transmission coefficient which is a measure of possibility of a molecule to move from one site to a vacant site. is the proportionality constant of activation energy and (assumed to be constant [6]) is the interaction potential energy between a pair of segments of a polymer.

3. Calculations and Discussion

3.1. The SS Theory Parameters

We analyzed the PVT behavior of some commercially available polymers of high-density polyethylene (HDPE), low-density polyethylene (LDPE60 and LDPE65), linear low-density polyethylene (LLDPE), polypropylene (PP), polystyrene (PS), poly(methyl methacrylate) (PMMA), and polycarbonate (PC). Their experimental PVT data reported by Sedlacek et al. [11] were accomplished on a fully automated PVT instrument (pvT 100 SWO, Germany) in the temperature range 70-290 °C and the pressure range 15-70 MPa with the increments of 5 MPa.

In the present paper we used the values of the characteristic scaling parameters , , and the structural flexibility parameter 3c/s with computed from the experimental data by the use of the coupled equations (3) and (4) from our previous work [6] in which LDPE 65 was excluded (LDPE refers to LDPE 60 in ref [6]). For the corresponding polymer the scaling parameters and the segment number are computed as , , , , and with the same method explained before [6]. These parameters are used for the determination of hole fraction at each temperature and pressure and later used in viscosity calculations.

3.2. The Elongational Viscosity Data Fit

The experimental rheological data for the polymers studied were reported by Sedlacek et al. [4] for three tested temperatures and six pressures (0.1, 10, 20, 35, 50, and 70 MPa). The elongational rate dependent viscosity was measured in a modified capillary rheometer (Göttfert 2001, Germany). A back-pressure device in order to raise the pressure and two capillaries with length-to-diameter ratio of 0.12 and 20 were used. The measurements were complemented by employing the Cogswell analysis together with pressure-dependent entrance pressure drops. The effect of viscous heating was omitted during the measurements because the magnitudes of the Nahme number were found to be less than 1. A closer description of the experimental technique and instrumentation can be found in the preceding article [4].

Now, these sets of data are fitted for each polymer using (5) and (6) with the double and single nonvanishing term(s) in the denominator, respectively, to obtain the zero level of elongational viscosities, , and the stress parameters, , along with the fraction as a nonlinearity parameter in the power of elongational rate. Our motivation is to obtain and for each temperature and pressure data set while is a sole parameter for each polymer. To obtain and , we employ nonlinear fit scheme assuming that is taken as a disposable quantity. The best value is adjusted with the minimum average absolute deviation in viscosity defined as (relative mean absolute percentage error)

The viscosity data are given in the range of 3 to 405 elongational rate. These datasets can be used to obtain and as well as parameter explained above scheme. For high pressures, we obtain all and by fixing value at ambient pressure. The and values at ambient pressure are presented in Table 1 alongside with mean percentage relative errors in viscosity, , and correlation coefficient squared, . The latter, also called the coefficient of determination, is a measure of how much of the variance in experimental value of is explained by the model. The percentages relative errors are computed by using (9). The values in parentheses are calculated using (6) while the rest with (5). To reduce figure crowding, we demonstrated HDPE at Figures 1, 2, and 4 as well as PMMA at Figure 3. Logarithmic elongational viscosity versus logarithmic elongational rate of HDPE is plotted at various temperatures at ambient pressure in Figure 1 and for various pressures only at 210 °C in Figure 2. Similar plots are obtained for the other polymers. The solid line is drawn by (5) with the best fit parameters. These parameters at ambient pressure are collected in Table 1 with and . In a similar manner, the values of the parameters in (6) are also tabulated in parenthesis. When we compare the result we can say that does not differ according to both models for HDPE. Only for PP the Cross-like model has a better capability to describe the viscosity values, while for the remaining polymers the values are estimated better with the YDT model. Specifically, from Table 1, it can be observed that, with the YDT model, good estimation is achieved for the linear or nearly linear polymers (HDPE; LLDPE) and PP with simple structure. The fitting is reasonable for linear long chain branched (LDPEs) and large group attached polymers like PMMA, PC, and PS. The mean percentage errors in viscosity for all the polymers considered are less than 6.30% at ambient and 6.28% at higher pressures. However, using (6), the mean percentage errors are less than 7.58% at ambient pressure and 6.59% at higher pressures. Overall, according to (5) decreases almost 0.5% at ambient and 0.8% at high pressures comparing to (6). Hence (5) better fits the viscosity data than (6) both at ambient and at high pressures. In the previous findings, the observed non-Newtonian behavior is better described by the Cross model only at ambient pressure (see [6]).

Moreover, in Figure 3, the zero level of elongational rate viscosity () and the stress parameter () are plotted as a function of pressure at various temperatures for PMMA. For the zero level of elongational rate viscosity the graph includes the best fit line. The zero elongational viscosity increases linearly with increasing pressure, but the stress parameter, depending on pressure in terms of , increases but not consistently in temperature. On the other hand, the zero elongational viscosity decreases steadily with increasing temperature, but we may not draw an exact temperature dependency on the stress parameter. The behavior of and with respect to pressure at elongational flow is similar to the one at shear flow (see [6]).

3.3. Elongational Viscosity-Hole Fraction Dependence under Constant Elongational Rate and Stress

The relationship between hole fraction and viscosities at zero and constant elongational rates and constant elongational stress given in (7) is tested. The zero level of elongational viscosity values is estimated from the experimental data by using both the YDT and the Cross-like models (see (5) and (6), respectively), for the studied polymers. With regard to the constant elongational viscosity, the data are evaluated properly through the polynomial function fit within the measured elongational rate and pressure ranges as follows: the viscosity data are fitted with respect to elongational rate by a third-order (for LDPE 65, LLDPE, PMMA, PP, and PS), second-order (for LDPE 60 and HDPE), and first-order polynomial equations (for PC). The interpolating values are calculated at constant elongational rate values, , in Table 3, at the three temperatures and various pressure levels discussed earlier in explaining the computation of elongational viscosity. Similarly, the constant elongational stress viscosities, , are calculated at certain constant elongational stress, , in Table 4, at the same temperature and pressure levels. The obtained zero elongational, constant elongational rates and constant elongational stress viscosity data are fitted to (7), in which the hole fraction is calculated by the SS theory. The regression parameters derived from the equation, , , and , , , are presented in Tables 24 with the correlation coefficient and obtained from (9) for the polymers. Using the best fit parameters, the values of in zero elongational, constant elongational rates and constant elongational stress are ranging from 1.35% to 3.32%, 0.23% to 2.44% and 0.47% to 4.97%, respectively. in average for zero elongational viscosity is found to be 1.87% via the YDT model with a better accuracy compared to the Cross model of which deviation is 1.97% according to the results listed in Table 2.

In Figure 4, the logarithms of zero elongational rate, constant elongational rate, and constant elongational stress viscosities are plotted with respect to for HDPE in which the solid lines are drawn through the data with the best fit line. A good linearization, for all conditions above, is achieved for the three tested temperatures and for all given pressures. In the figure, the pressure increases from left (the lower part) to right (the higher part) along each line and at each temperature. As the hole fraction decreases, increases for each temperature, so the viscosity increases. Each line slope gives the value of α, in Kelvin unit, evident from (7) and (8). As expected, in Figure 4, the logarithm of viscosities, , decreases as the elongational rate and elongational stress increase for HDPE.

For comparison, the special cases of Figure 4 at 90 elongational rate and certain elongation stresses are delineated in Figures 5 and 6, respectively. Figure 5 shows a plot of versus at 90  and Figure 6 displays versus graph at certain elongational stresses (the elongational stresses used in Figure 6 indicated with italic font in Table 4). Overall, we conclude that, for the elongational rate and all elongational stress values, polymers with simpler structures, PEs and PP, show slowly varying viscous behavior with respect to the thermooccupancy function, , so their slopes are smaller than PMMA, PC, and PS. PMMA has the highest thermooccupancy functional dependency, which seems to result from two pendant groups (CH3 and CO2CH3) on each of its monomer units. Accordingly, the measure of the activation energy coefficient, , positively correlates with activation energy; thus this indicates that PMMA, PC, and PS polymers require more activation energy, , than polyolefins (PP and PEs). The differences in and values from Tables 3 and 4 can relatively be attributed to the large benzene ring on the former structures and related to the equilibrium shape of the chain [21]. Figures 7 and 8 graphs and with respect to elongational rate and stress and provide a clearer picture to this pattern.

The computed and parameters are analyzed at various constant elongational rate and constant elongational stress values including those in non-Newtonian region. At constant elongational rates of 0, 60, 90, and 120 , the slopes of the lines for HDPE are 126.93, 35.23, 35.15, and 34.58, respectively. As can be seen the slopes of the lines are inversely correlated to elongational rates. This can be explained by the fact that, with the increasing , the flowability increases, so decreases. The same pattern is observed for the other polymer types except PC ( and directly related to elongational rate) in Figure 7(a). A quite slight decrement is observed since the values of the strain rates and strain rate ranges are smaller in the elongational flow (3-405 ) than in the shear flow (0.1-8500 ) [6]. The activation energy coefficient, , is calculated for each polymer given in Table 3 and plotted in Figure 7(b). The behavior of is similar to . Figure 8 shows the change of and parameters with elongational stress for the polymers given in Table 4. These parameters slightly increase with elongational stress. However, in our previous work where the stress values are much lower than the elongational stress values as expected, and stay almost constant with shear stress. Figure 8 is drawn with the exception of PMMA since and inversely related to elongational stress. PC is also excluded in Figure 8(b) because of its largest values. In comparison, PMMA, PC, and PS have higher and values than polyolefins. As noticed above, this can be arising from the fact that, bulky pendant groups attached to the backbone of these polymers make them more elongated and increase the intermolecular friction (viscosity). As the viscosity increases, the activation energy also increases. Consequently, the more activation energy requires the higher values of and .

Further, the transmission coefficient, , in (8) is also calculated and plotted in Figures 9(a) and 10 with respect to the elongational rate and the elongational stress, respectively. It has been found that the transmission coefficient, , increases linearly along with the increasing elongational rate and elongational stress. In Figure 9(a), there is a slight increment because of the low and small range of elongational rates. This rise enables the molecules to jump into the holes much more easily. On the other hand, the tendency of the molecules to the flowability decreases the intercepts, , at the viscosity axis as shown in Figures 9(b) and 11. represents the extrapolated viscosity values when the system dominates sufficiently large hole fraction expectedly at high temperature and low pressure in (7). In other words, the thermooccupancy function, , goes to zero. This contrary behavior of to is obvious from Tables 3 and 4. Graphical results of this observation are depicted in Figures 911. Increasing the constant elongational rate given in Table 3 and the constant elongational stress given in Table 4, the intercepts for HDPE decrease steadily as 9.26, 8.98, 8.83 Pa·s and 9.25, 8.53, 8.11 Pa·s, respectively. The decrements for the elongational rate and the elongational stress are also observed for the other polymers shown in Figures 9(b) and 11. LDPE 65 in Figure 9 and PMMA in Figure 11 are excluded, because of the increments with elongational rate and stress, respectively. From Figures 911 we understand that the calculated values of of the constant elongational rate and stress of PC, PS, and PMMA are larger than of the polyolefins while they have smaller values of . This can be interpreted, such that the polymers, PMMA, PC, and PS, have bulky and larger pendant groups as allyl, methyl, and phenyl. Since the polymers with the side groups are very sensitive to produce more free volume and less viscous behavior as going to higher temperature or lower pressure, they take the higher values of and the smaller values of . The calculated values of of the constant elongational stress order of the materials are as follows: HDPE>LLDPE>LDPE>PP>PC>PS>PMMA. The transmission coefficient, , order for the constant elongational stress is vice versa. This result is the same for the shear viscosity we obtained before (see [6]).

4. Conclusions

The Yahsi-Dinc-Tav (YDT) non-Newtonian viscosity equation and the Cross-like model were employed to describe the correlation among the elongational viscosity and elongational rate as well as the rheological parameters of polymer materials. The results from the YDT model showed that the estimations were in good agreement with the experimental data both at ambient and at high pressures. In contrast to the findings for shear rate in our previous work (where the Cross model better fitted the data at ambient pressure), in this paper, the YDT model provides 0.5% less deviation than the Cross model at ambient pressure. The YDT equation gives better results than the Cross-like model not only at ambient pressure but also at elevated pressure. The mean percentage error decreases almost 0.8% at higher pressure by the use of the YDT model. Analyses also show that good estimation is observed for the polymers with high elongational rates (up to about 405 and 308 for PC and LLDPE, respectively). Moreover, the viscosity parameters determined from shear (from our earlier paper) and elongational flow maintained the similar temperature and pressure dependence. In brief, with this model, correct nonlinear behavior in both shear and elongation was achieved for the linear, branched, and large group attached polymers.

Regarding the dependence of viscosities on thermooccupancy function (), the YDT-derived linear equations, at zero level of elongational and constant elongational rates as well as at constant elongational stress, provide a good description of dependency of hole fraction on viscosity. It was found that, for individual polymers, a linear relationship exists between the logarithm of zero elongational, constant elongational rate, and constant elongational stress viscosities (at several temperatures and pressures) and . Similar to the results we obtained in our earlier study for the shear viscosity, we show that an increase in thermooccupancy function results in higher viscosity values at various elongational rate and elongational stress conditions, since decrease in hole fraction causes stronger interactions between molecules. Polymers with the simpler structures, (HDPE, LLDPE, LDPE, and PP) appear to have the lower thermooccupancy function dependence. Larger pendant groups linked to the backbone chains of the remaining polymers (PC, PS, and PMMA) result in a higher dependency on temperature- and pressure-dependent thermooccupancy function.

Further, the accessibility of rheological data over a broad range of pressures and temperatures and the hole fraction computed from the SS theory allow the viscosity parameters such as activation energy and transmission coefficients to be computed from such linear formulations. The changes in these viscosity parameters are discussed in terms of constant elongation rate and constant elongational stress. Estimates show that activation energy coefficient is negatively correlated with elongational rate. With regard to the elongational stress, it has a direct relationship contrary to what we have found for shear viscosity. Confirming the results, we have had for shear flow, transmission coefficient is positively correlated with constant elongational rate and constant elongational stress, while the viscosity-intercept has an inverse relationship.

In summary, the model predictions for the nonlinear behavior in both elongation and shear viscosities in this and our preceding work have been revealed to be in very good agreement with the reported experimental data. This viscosity equation in experimental studies is reliable and functional with only a fitting parameter as well as convenient not only for linear polymers, but also for the polymers with long chain branching and large attached groups. An intriguing forthcoming study may be on predicting other physical properties of different classes of neat polymers and their compositions with a similar methodology followed in the current work.


:Proportionality constant of activation energy
:Total external degrees of freedom
:Activation energy (J)
:Planck constant
:Hole fraction
:Boltzmann’s constant
:Segmental molar mass of molecules (kg)
:Avogadro’s number
:Reduced pressure, volume, and temperature
:Characteristic pressure, volume, and temperature
:Hole fraction and reduced volume dependent quantity in (3) and (4)
:Nonlinearity parameter in the power of elongational rate
:Number of interchain nearest neighbor pairs in a lattice of coordination number
:Gas constant (8.314 J/mol K)
:Number of segments of molecules
:Temperature (°C or K)
:Molar segmental volume
:Characteristic molar volume of a segment
:Occupied site fraction
:“Thermooccupancy function”, in (7)
:Coordination number.
Greek Letters
α:Slope of (7)
:Elongational rate ()
:Constant elongational rate
:Zero level of elongational rate
:Average percentage error in viscosity
:Attractive interaction parameter of a segment of the potential minimum
:Intercept of (7) and (8)
, :Zero level of elongational viscosity, constant elongational rate viscosity (Pa·s)
:The angle between the deformation stress and the displacement vector of th contiguity of a segment
:Transmission coefficient
:Interaction potential energy between a pair of segments (.K)
:Elongational stress (N/m2)2
:Critical stress parameter.

Data Availability

The non-Newtonian viscosity and the PVT data supporting this theoretical analysis are from previously reported studies and datasets which were published by the one of the authors Sedlacek et al. The data have been cited in the text as [4, 11] for the elongational non-Newtonian viscosity and the PVT data, respectively. The processed data are available from the corresponding author upon request.


This research has not been financially funded by any special and governmental project and support. And also this work has not performed as part of the employment of the authors.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.