Abstract
The effects of polymer concentration and electrospinning parameters on the diameter of electrospun polymethyl methacrylate (PMMA) fibers were experimentally investigated. It was also studied how the controlled factors would affect the output with the intention of finding the optimal electrospinning settings in order to obtain the smallest PMMA fiber diameter. Subsequently the solution feed rate, needle gauge diameter, supply voltage, polymer concentration and tip-to-collector distance were considered as the control factors. To achieve these aims, Taguchi’s mixed-level parameter design (L18) was employed for the experimental design. Optimal electrospinning conditions were determined using the signal-to-noise (S/N) ratio that was calculated from the electrospun PMMA fiber diameter according to “the-smaller-the-better” approach. Accordingly, the smallest fiber diameter observed was 228 (±76) nm and it was yielded at 15 wt% polymer concentration, 20 kV of supply voltage, 1 ml/h feed rate, 15 cm tip-to-distance and 19 needle gauge. Moreover, the S/N ratio response showed that the polymer concentration was the most effective parameter on determination of fiber diameter followed by feed rate, tip-to distance, needle gauge and voltage, respectively. The Taguchi design of experiments method has been found to be an effective approach to statistically optimize the critical parameters used in electrospinning so as to effectively tailor the resulting electrospun fiber diameters and morphology.
Keywords
Introduction
The electrospinning process provided the opportunity to fabricate micrometer and nanometer size diameters of polymeric and inorganic fibers.1,2 The first step of the electrospinning process is carried out by putting a polymer solution into a syringe supplied through a needle connected to a high-power source. The high voltage difference between the needle and the collector, generally in the range of 5–30 kV, can be produced by a high-power source. In the second step of the process, the solvent induces evaporation by electrical charges after the polymer solution is ejected from the syringe. As a result, an elongated polymer fiber will deposit on the collector. In the past, a number of factors were investigated in order to comprehend as well as manage the electrospinning method.1,3,4
There are number of scholars who were interested in the fiber diameters and morphology and, indeed, they made several evaluations. However, one prominent phenomenon can be seen when browsing through these articles concerning the influence of parameters on the properties of fibers is the manifestation of the vast amount of contradicting observations and also the possible occurrences that these parameters are not affecting in the same manner for all types of polymers and solvents.1,3,5 Furthermore, it is shown that there is a considerable amount of electrospinning process outputs that is dependable on the adjustment of its parameters. For that reason, this work is conducted to fulfill this requirement as further investigations addressing possible cross-talk between different processing parameters and modifications coming from changes in the spinning configurations are definitely needed.6,7
An appropriate static model can be utilized for the electrospinning process to evaluate the effects of the parameters on the fiber diameter. This model is capable of displaying the highest effective factors on the fiber diameters. For this reason, in the empirical studies these parameters are the ones to be focused on. The experimental time could be effectively reduced by a theoretical prediction and a route to comprehend the factors influencing the fiber diameter and morphology.8,9 For instance, controlling the diameter of the fibers is still a methodological bottleneck. The experimental results show that the lowest fiber diameter can be achieved by the lowest flow rates, while, conversely, the fabrication rate and solution concentration affect the fiber diameter.5,10 Fortunately, using an appropriate predictive method would supply a route to perform many various options rapidly, without having the cost or time problems of a trial-and-error experiment. Recent research has not yet holistically investigated the optimization of the electrospinning process in a systematic manner. It is conventionally limited to the determination of the physical, thermal and structural properties of electrospun nanofibers/nanocomposites with the fixed processing and material parameters.1,3,4,11
The design of experiments (DoE) is a very robust experimental design system technique that can explore the significant factorial effects and optimum conditions in the modern manufacturing industry setting.12,13 Orthogonal arrays are utilized in the Taguchi method to organize the factors that influence the process and the levels where they should vary. Alternatively, the Taguchi method checks a pair of combinations rather than having to examine every possible combination, such as the full factorial method. In this method, the parameters that influence product quality could be indicated with a minimal number of trials. 14 The Taguchi method is giving the opportunity to be able to evaluate all parameters together or independently. The purpose of this paper is to optimize fiber diameter and understand what are the individual effects of parameters that are solely achievable through this method, not only for the electrospinning process but also for many other production processes.8,15–17
Several researchers have focused on the fabrication of polymethyl methacrylate (PMMA) via electrospinning; different achievements have been reported from previous works in terms of fiber diameter, PMMA fiber properties and the effect of individual parameters on the fiber properties.18–21 In fact, none of them have focused on the individual effects and interaction of the electrospinning parameters on the prediction of fiber diameter. Also, a large fiber diameter and a large number of beads are found to be the main problems with their results. 18 In order to fully recognize the potential of electrospun fibers, it is significant to produce different fiber diameters, which will definitely have an important influence on the performance. For example, filtration and composite materials have potential applications of small fibers.3,4 It should be noted that other applications, such as drug delivery and catalyst support, have attracted interest due to their high surface area as well as electrospinning of tissue engineering, which is being developed.1,3,4
This paper aims to study the effect of polymer concentration and electrospinning parameters on PMMA nanofiber morphology so as to find the optimized condition that yields the finest fibers. In order to do so, a robust statistical Taguchi DoE has been employed. In general, we considered the following questions for careful deliberations: What are the factors instigating these contradicting research findings found predominantly in the PMMA nanofiber production? And also, is there any possibility to have all-inclusive laws of dependencies concerning parameters and fiber morphology?
Materials and methods
Materials
In this work, analytically pure PMMA (−CH2C−(CH3))CO2CH3−) n , Mw = 120,000) from Aldrich and N,N-dimethylformamide (DMF), which was received from Labchem Sdn Bhd Co., Malaysia, were utilized for the working fluid. Correspondingly, polymer solution samples were obtained through dissolving 15, 20 and 25 wt% PMMA in a DMF solvent.
Experimental procedures
A precursor solution was prepared using DMF as the PMMA solvent. The typical process was carried out by dissolving an adjusted quantity of PMMA in DMF through vigorous stirring for 24 hours. After stirring, different weight ratios of PMMA solutions were achieved. Subsequently, a 5 ml plastic syringe with different needle gauges, 19, 23 and 27, were utilized to place each polymer solution, attributable to the significance of orifice size in the initiation of the jet.6,20
ES30P-30 W/SDPM (Gamma High Voltage Research, Ormond Beach, FL) was used as the source for the high voltage power supply. An appropriate Taylor cone was created for experimentation by using two different voltages – 15 and 20 kV – in which the high voltage was connected to a needle through an alligator clip. Then, a ground target covered with aluminum foil was used as the counter electrode, at three different tip-to-distances of 10, 15 and 20 cm. In addition, a syringe pump (NE-300, New Era Pump Systems, Inc.) was used in order to control the feed rate at three different rates: 1, 2 and 4 ml/h.
In view of the fact that the choice of the temperature and relative humidity (RH) are important, 19 the electrospinning process was performed at 25℃ and 32% RH. Modification of the polymer concentration, needle gauge, the collection distance and the feed rate was executed to obtain controllable nano-fiber polymer materials.
The morphology of electrospinning experiments was investigated by using a Field Emission-Scanning Electron Microscope (FE-SEM), Zeiss (Auriga). The diameters of 20 fibers were randomly measured using Digimizer 4.1 software.
Design and analysis of experiments
Five factors and their levels selected in the L18 design of experiments array
Experimental results of electrospun polymethyl methacrylate fibers diameter and signal-to-noise (S/N) ratios based on “smaller is better”
Also, optimization of the electrospinning factors has been done accurately. The performance characteristic function that the smaller is better was utilized for obtaining the optimal electrospinning factors for the diameter of the fibers. The S/N ratio η formulate is expressed as
Results and discussion
Surface morphology and distribution of fiber diameter
The FE-SEM micrographs of the composite fiber morphology in the DoE study and the corresponding average fiber diameters determined are illustrated in Figure 1 and Table 2. Thus, hereby, the discussion of the results is both quantitative and qualitative; the coverage of the qualitative justifications for fiber morphology includes uniformity, bead formation and fiber diameters and, then, quantitative justifications are established on fiber diameter, S/N ratio and p-value. On a global scale, a random fiber distribution is evident in most mat samples due to using a stand-alone mesh collector with little fiber alignment control. Electrospinning polymer solutions involve evaporation of the solvent when the jet is accelerating to the collector; the surface area of the jet increases significantly in a couple of milliseconds and some thermodynamically driven actions can be concluded.
Field emission scanning electron microscope images of electrospun polymethyl methacrylate fibers used in the design of experiments study: (a) ES#1; (b) ES#2; (c) ES#3; (d) ES#4; (e) ES#5; (f) ES#6; (g) ES#7; (h) ES#8; (i) ES#9; (j) ES#10; (k) ES#11; (l) ES#12; (m) ES#13; (n) ES#14; (o) ES#15; (p) ES#16; (q) ES#17; (r) ES#18 (magnification is 1000× and scale bar is 10 µm).
Analysis of variance for signal-to-noise ratios of electrospun polymethyl methacrylate fiber diameter
DF: degree of freedom; SS: sum of squares; MS: mean of squares; PCR: percentage contribution ratio.
In addition, large entangled fiber aggregates are shown in most of the experiments, resulting from the hindrance of the fiber structures from the solution jet. This typical phenomenon could be due to the combined effect of the higher feed rate and, consequently, more solution droplets being ejected from the needle tip per hour that do not have sufficient time to be elongated to individual fibers with relatively small fiber sizes received directly from the mesh collector. As a result, fiber aggregates or bundles become more evident due to the lack of fiber stretching/elongation.5,6,19 On the other hand, short fiber collections are noticed in ES#2, ES#8, ES#11 and ES#17 in the similarity of tip-to-distance value, equal to 15 cm. These shortened fibers can occur as the result of jet whipping instability during electrospinning, which leads to bending and stretching of the jet. 23
The small-bead defects are seen in ES#1, ES#2, ES#4, ES#7, ES#12 and ES#15 (Figure 1), which may arise from the clogging problem of the middle-concentration (wt% = 20) of PMMA at the needle tip, except in ES#7. Therefore, droplets were formed with the flow instability deteriorating the elongated fiber structures. The number of beads is found to be more in ES#6, ES#8, ES#11, ES#14 and ES#18 as compared to other experiments. The results of the mentioned tests illustrate that production and morphology of fibers are highly associated with solution viscosity.6,19 It has been reported that when the polymer concentration is low, the number of beads or microspheres created in electrospun fibers can change to electrospraying when the concentration decreases significantly.19,24,25
In contrast, the electrospun PMMA fibers become bead-free, as evidenced in ES#3, ES#5, ES#9, ES#10, ES#13, ES#16 and ES#17 (Figure 1). As a result of the high viscosity of the solution (high polymer concentration), bead-free uniform nanofibers were created while 25 wt% polymer concentration was appropriate for PMMA. This may reveal that the creation of the beads is mostly as a result of surface tension, which reduces the surface area. It has been demonstrated that the jet tends to break down into drops without surface tension and that beads are likely to form in low surface tension; also, it can even be justified that the significance of the p-value for concentration is prominent for controlling bead formation too. Therefore, fewer beads are created by increasing the polymer concentration and the surface tension. Based on the summarized average fiber diameter data (daverage < 300 nm), the minimum average diameter of 228.10 (±76) nm is found in ES#11, followed by ES#1, ES#14 and ES#18 in the average diameter range of 200–300 nm (Table 2). These average values of fiber diameter are the lowest fiber diameters that have been achieved for electrospun PMMA nanofibers in comparison to previous works.
In order to portray a clear view of the subject, histogram plots showing the percentage contribution against the electrospun PMMA fiber diameter range achieved by experimental tests are demonstrated in Figure 2. Upon analyzing the plots it was obvious that ES#11 is the only case study with a high percentage of frequencies distributed in the diameter range of 100–300 nm. Moreover, ES#1, ES#6, ES#14 and ES#18 show a contribution in the range of less than 500 nm. The contribution that was observed in the ES#8, ES#12, ES#15 and ES#16 experiments was in the range of less than 1 µm, along with the previously mentioned experiments. In contrast, ES#2, ES#3, ES#4, ES#5, ES#7, ES#9, ES#10, ES#13 and ES#17 have been shown as being able to fabricate fibers with a diameter higher than 1 µm, taking into consideration that ES#2, ES#4 and ES#7 only had contribution in the more than 1 µm range. It should be noted that ES#5, ES#9, ES#10, ES#13 and ES#17 contributions were distributed over all ranges from less than 1 µm to higher. Overall, ES#17 presented the worst set-up contributions in the range of fiber diameters less than 1 µm, with other contributions mostly allocated to higher than 1.5 µm.
Frequency contribution of fiber diameter range in the design of experiments (DoE) study. Figure 1 field emission scanning electron microscope images of electrospun polymethyl methacrylate fibers used in the DoE study: (a) ES#1; (b) ES#2; (c) ES#3; (d) ES#4; (e) ES#5; (f) ES#6; (g) ES#7; (h) ES#8; (i) ES#9; (j) ES#10; (k) ES#11; (l) ES#12; (m) ES#13; (n) ES#14; (o) ES#15; (p) ES#16; (q) ES#17; (r) ES#18. AD: average diameter.
Analyzing and evaluating the results of the experiments using the Taguchi method
Experimental results and signal-to-noise ratios
When all the experiments were carried out based on the L18 orthogonal array, the average fiber diameter that occurred during the electrospinning process of the PMMA material was measured. Furthermore, in order to control the fiber diameter, the voltage, feed rate, tip-to-distance, needle gauge and polymer concentration were considered in the development of the mathematical models. It should also be noted that the S/N ratio is the most important deciding factor in the Taguchi method to analyze the experimental data. In this study, the S/N ratio should be the maximum value to get the optimal electrospinning condition in agreement with the Taguchi method. Table 2 shows the experimental results and S/N ratios, which are calculated based on Taguchi’s “smaller is better” quality characteristic.
Table 2 shows the S/N ratios of the fiber diameter data achieved from the experimental test results that were calculated according to Equation (1), which will be utilized for establishing the optimum levels of each factor.
Analysis of variance
Table 3 shows the results of the ANOVA that was carried out to conclude the significant levels of the factors on the fiber diameter that occurred in the electrospinning process of PMMA. Also, the p-values, which show the significance level of each variable and its degree of freedom (DF), sum of squares (SS), mean of squares (MS), F-ratios and their percentage contribution ratio (PCR), for which the principle equation can be found in Antony 12 and Eriksson, 13 are indicated in Table 3. These are used for determining the significance levels of the variables in the ANOVA table. Table 3 indicates that the polymer concentration has the lowest p-value among other variables, which is also the sole significant factor. On the hand, other variables have shown roughly high p-values (more than 0.05), which shows those are not significant factors in the process.
Optimum electrospinning parameter condition and parameters interaction
Figure 3 illustrates the plots of S/N ratios that were calculated for the fiber diameter of electrospun PMMA. By using Table 2 and Figure 3, the variables that will be utilized to achieve an optimal fiber diameter value and their respective levels are concluded as A2, B1, C2, D3 and E1. In addition, Table 4 shows the different values of the S/N ratio between the maximum and the minimum main effects. The polymer concentration (E), the feed rate (B) and the tip-to-distance (C) are three factors that have the highest difference between values, 12.57, 0.97 and 0.68, respectively. As derived from the Taguchi prediction, the larger differences between the values of S/N ratio will have a more significant effect on the fiber diameter. Consequently, it can be concluded that increasing the polymer concentration will also significantly increase the fiber diameter.
Main effects plot for signal-to-noise (S/N) ratios of electrospun polymethyl methacrylate fiber diameter. Signal-to-noise response, Delta and ranking table for electrospun polymethyl methacrylate fiber diameter
Subsequently, the smallest fiber diameter and its S/N ratio of level were calculated according to Equations (2) and (3). Finally, the fiber diameter value was calculated as 283 (±69) nm and its S/N ratio was determined as −48.4863 dB.
Accordingly, for the experiments to be conducted, the optimum electrospinning conditions should be determined under the same conditions, which will be at 20 kV for the voltage, 1 ml/h for the feed rate, 15 cm for the tip-to-distance, 27 gauge for the needle and 15 wt% for the polymer concentration. The evaluation of the diameter of fibers formed in electrospinning of PMMA was done based on the average fiber diameter (nm). From this point, Figure 4 shows the graphs given to illustrate the relationships between the fiber diameter and the main effects of the electrospinning parameters.
Main effects plot for means of electrospun polymethyl methacrylate fiber diameter.
The plot (Figure 5) represents the interaction of the electrospinning parameters with each other and, ultimately, their effect on the fiber diameter. This illustrates the combined effect of the factors on a response variable. The Y-axis corresponds to the fiber diameter and the X-axis denotes the different electrospinning parameters.
Interaction plot for means of electrospun polymethyl methacrylate fiber diameter.
In order to better interpret the figures, firstly, a general conception of the curves was given and, ultimately, in order to reach a more detailed conclusion the study focused solely on the relationship between each two parameters, and once the common trend was explored the spotlight shifted to the next two parameters.
Initially it was observed that the only parameter that follows a similar pattern in all cases was concentration. In other words, when concentration increases, the fiber diameter rises dramatically with a high slope, regardless of the other variable it was matched up with. This behavior, as mentioned before, shows independency of the concentration effect on fiber diameters and its substantial role in deciding the outcome of the process. In other cases a more haphazard pattern is observable, implying that the other factors were mostly reliant on each other and that they did not possess the supremacy to influence the results independently.
An increase in the voltage appeared to have a neutralizing effect on the diameter except in the pair of “voltage-needle to collector distance”, for which a vivid deviation can be observed. This indicates that the effect of voltage was more prominent in very short or very long distances and it appears that an increase in voltage would lead to larger fibers when the needle-to-collector distance was too close and vice versa when the distance was too long. It is worth mentioning that an increase in voltage was able to increase the diameter when the feed rate was too high or when the needle was large.
The feed rate almost showed a random trend, which was mostly dependent on the factor with which it was being paired up with. In the case of the “Feed rate-needle to collector distance” curve, when the distance was set on an average scale it looked like the feed rate was able to dominate the equation and that the increase in the feed rate would lead to an increase in the fiber diameter. In addition, the feed rate was generally decreasing the diameter in short or long distances. Examination of the “Feed rate-needle” graph revealed that an increase in feed rate would only lead to an increase in diameter when the needle was too large, while when the needle was of normal or small size the fiber diameter decreased after an increase in the feed rate.
The tip-to-collector distance figures displayed a more random trend but when they were studied in more detail (distance-needle curve) it was found that an increase in distance caused a decrease in the fiber diameter unless the needle was too small.
Confirmation experiments
Taguchi mentioned that a confirmation test is a critical step in verifying the experimental conclusion. The confirmation experiment results that were carried out at their respective optimal factor levels (A2, B1, C2, D3 and E1), for establishing the fiber diameter were assessed by considering the confidence interval (CI) calculated using Equations (4) and (5):
15
Comparison between confirmatory experiment result and calculated values
S/N: signal-to-noise.
Table 5 provides a comparison of the findings of the confirmation tests that have been carried out depending on the optimal levels of the factors and the values calculated using Equations (2) and (3). This result was within a 95% confidence interval of the predicted optimal value of the selected parameters. The optimal setting of the electrospinning process parameters, as predicted in the analysis, was implemented. After three repetitions of the optimal condition the average fiber diameter was achieved as 270 (±58) nm. In addition, it should be mentioned that the experimental fiber diameter value was very close to the predicted value; therefore, it is suitable for predicting the electrospun PMMA fiber diameter.
Conclusion
PMMA fibers were electrospun based on five factors of the L18 Taguchi DoE. Polymer concentration and processing parameters on the fiber diameter were assessed in order to discover the optimal factor levels for a smaller fiber diameter during electrospinning. The results achieved from this investigation are as below.
The morphology of electrospun PMMA fibers displayed large entangled fiber aggregates in most of the experiments. Larger nanofiber diameters and fewer beads were created by increasing the polymer concentration. However, bead-free and uniform nanofibers were frequently formed in high polymer concentrations.
The smallest average fiber diameter values obtained during the electrospinning of PMMA were around 228 (±76) nm (it was yielded at 15 wt% polymer concentration, 20 kV of supply voltage, 1 ml/h feed rate, 15 cm tip-to-distance and 19 needle gauge), while the largest average fiber diameter was approximately 1686 (±455) nm. Different fiber diameter distributions were experienced in each electrospinning condition. PMMA nanofibers with a diameter range of 100–300 nm were achieved in the ES#11 experiment condition.
The effects of the factors on electrospun PMMA fiber diameter were determined using ANOVA. The most significant factor for electrospun PMMA fiber diameter was found to be the polymer concentration, followed by the feed rate, the tip-to-distance, needle and voltage supply. The optimum electrospinning conditions were determined to be as follows: 20 kV for the voltage, 1 ml/h for the feed rate, 15 cm for the tip-to-distance, 27 gauge for the needle and 15 wt% for the polymer concentration. In addition, the experimental fiber diameter value was very close to the predicted value.
Footnotes
Acknowledgement
This work was supported by the University of Malaya (grant numbers UMRG RP021-2012B and UMRG RP021-2012C) and the Ministry of Higher Education Malaysia through High Impact Research Grant MOHE-HIR D000008-16001.
