Drawing of the hollow all-polymer Bragg fibers

Drawing of the hollow all-polymer Bragg fibers Elio Pone, Charles Dubois, Ning Guo, Yan Gao, Alexandre Dupuis, Francis Boismenu, Suzanne Lacroix, Maksim Skorobogatiy École Polytechnique de Montréal, C.P. 6079 Succursale centre-ville, Montréal (Québec) H3C 3A7, Canada [email protected] www.photonics.phys.polymtl.ca Abstract: Drawing of the hollow all-polymer Bragg fibers based on PMMA/PS and PVDF/PC materials combinations are demonstrated. Hole collapse during drawing effects the uniformity of a photonic crystal reflector in the resultant fiber. We first investigate how the hole collapse effects fiber transmission properties. We then present modelling of fluid dynamics of hollow multilayer polymer fiber drawing. Particularly, hole collapse during drawing and layer thickness non-uniformity are investigated as a function of draw temperature, draw ratio, feeding speed, core pressurization and mismatch of material properties in a multilayer. Both the newtonian and generalized newtonian cases are considered assuming slender geometries. © 2006 Optical Society of America OCIS codes: (060.2280) Fiber design and fabrication; (060.2290) Fiber materials; (160.5470) Polymers; (220.4000) Microstructure fabrication; (230.1480) Bragg reflectors; (230.4170) Multilayers; (999.9999) Polymer microstructured fiber drawing References and links 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. T. Katsuyama and H. Matsumura. Infrared Optical Fibers, (Adam Hilger, Bristol, England, 1989). M. Saito and K. Kikuchi, “Infrared optical fiber sensors,” Opt. Rev. 4, 527-538 (1997). J. Sanghera and I. Aggarwal. Infrared Fiber Optics, (CRC, Boca Raton, USA, 1998). S. Martellucci (editor), A.N. Chester (editor), and A.G. Mignani (editor). Optical Sensors and Microsystems : New Concepts, Materials, Technologies, 1st ed. (Springer, New York, USA, 2000). J.A. Harrington, “A review of IR transmitting, hollow waveguides,” Fib. Integr. Opt. 19, 211 (2000). Y.W. Shi, K. Ito, Y. Matsuura, M. Miyagi, “Multiwavelength laser light transmission of hollow optical fiber from the visible to the mid-infrared,” Opt. Lett. 30, 2867-2869 (2005). P. Russell, “Photonic crystal fibers,” Science 299, 358-362 (2003). J. Canning, E. Buckley, and K. Lyytikainen “Propagation in air by field superposition of scattered light within a Fresnel fiber,” Opt. Lett. 28, 230-232 (2003). J. Canning, E. Buckley, K. Lyttikainen, and T. Ryan, “Wavelength dependent leakage in a Fresnel-based air-silica structured optical fibre,” Optics Communications 207, 35 (2002). C.M. Smith, N. Venkataraman, M.T. Gallagher, D. Muller, J.A. West, N.F. Borrelli, D.C. Allan, K.W. Koch, “Low-loss hollow-core silica/air photonic bandgap fibre,” Nature 424, 657-659 (2003). H. Han, H. Park, M. Cho and J. Kim “Terahertz pulse propagation in plastic photonic crystal fibers,” CLEO/Pac - Conference on Lasers and Electro-Optics, postdeadline paper Cat. No. 01TH8557, 22 (2001). J. Choi, K.Y. Kim, U.C. Paek, POF - Proceedings of Plastic Optical Fibres, 355 (2001). M.A. van Eijkelenborg, A. Argyros, G. Barton, I.M. Bassett, M. Fellew, G. Henry, N.A. Issa, M.C.J. Large, S. Manos, W. Padden, L. Poladian, J. Zagari, “Recent progress in microstructured polymer optical fibre fabrication and characterisation,” Opt. Fiber Techn. 9, 199-209 (2003). T.Katagiri, Y. Matsuura, M. Miyagi, “Photonic bandgap fiber with a silica core and multilayer dielectric cladding,” Opt. Lett. 29, 557-559 (2004). #69874 - $15.00 USD (C) 2006 OSA Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5838 15. B. Temelkuran, S.D. Hart, G. Benoit, J.D. Joannopoulos, Y. Fink “Wavelength-scalable hollow optical fibres with large photonic bandgaps for CO2 laser transmission,” Nature 420, 650 (2002). 16. G. Vienne, Y. Xu, C. Jakobsen, H.J. Deyerl, J. Jensen, T. Sorensen, T. Hansen, Y. Huang, M. Terrel, R. Lee, N. Mortensen, J. Broeng, H. Simonsen, A. Bjarklev, and A. Yariv, “Ultra-large bandwidth hollow-core guiding in all-silica Bragg fibers with nano-supports,” Opt. Express 12, 3500-3508 (2004). 17. Alexander Argyros, Martijn A. van Eijkelenborg, Maryanne C. J. Large, and Ian M. Bassett, “Hollow-core microstructured polymer optical fiber,” Opt. Lett. 31, 172-174 (2006). 18. A. Dupuis, Y. Gao, N. Guo, E. Pone, N. Godbout, S. Lacroix, C. Dubois and M. Skorobogatiy “Biodegradable, Double-Core, Porous Optical Fiber,” CLEO - Conference on Lasers and Electro-Optics, postdeadline paper Cat. No. CPDB5, (2006). 19. M. Skorobogatiy, “Efficient anti-guiding of TE and TM polarizations in low index core waveguides without the need of omnidirectional reflector,” Opt. Lett. 30, 2991 (2005). 20. Y. Gao, N. Guo, B. Gauvreau, M. Rajabian, O. Skorobogata, E. Pone, O. Zabeida, L. Martinu, C. Dubois, M. Skorobogatiy, “Consecutive Solvent Evaporation and Co-Rolling Techniques for Polymer Multilayer Hollow Fiber Preform Fabrication,” to appear in september issue of the J. Materials Research, 2006. 21. M.R. Matovich and J.R.A. Pearson, “Spinning a molten threadline - Steady-state isothermal viscous flows,” Ind. Eng. Chem. Fundam. 8, 512-520 (1969). 22. Y.T. Shah and J.R.A. Pearson, “On the stability of nonisothermal fiber spinning,” Ind. Eng. Chem. Fundam. 11, 145-149 (1972). 23. J.A. Burgman “Liquid glass jets in the forming of continuous fibers,” Glass Technol. 11, 110-116 (1970). 24. F.T. Geyling, “Basic fluid dynamic consideration in the drawing of optical fibers,” Bell Sys. Tech. J. 55, 10111056 (1976). 25. B.D. Freeman, M.M. Denn, R. Keunings, G.E. Molau and J. Ramos, “Profile development in drawn hollow tubes,” J. Polym. Eng. 6, 171-186 (1986). 26. A.D. Fitt, K. Furusawa, T.M. Monro, C.P. Please and D.J. Richardson, “The mathematical modelling of capillary drawing for holey fibre manufacture,” J. Eng. Math. 43, 201-227 (2002). 27. K. Lyytikainen, J. Canning, J. Digweed, J. Zagari, “Geometry control of air-silica structured optical fibres using pressurisation”, IMOC Proceedings, International Microwave and Optoelectronics Conference (Cat. No.03TH8678) 2, 1001-1005 (2003). 28. K. Lyytikainen. Control of complex structural geometry in optical fibre drawing, (Phd Thesis, University of Sydney, 2005). 29. S.C. Xue, R.I. Tanner, G.W. Barton, R. Lwin, M.C.J. Large and L. Poladian, “Fabrication of microstructured optical fibers Part I: Problem formulation and numerical modeling of transient draw process,” J. Lightwave Technol. 23, 2245-2254 (2005). 30. S.C. Xue, R.I. Tanner, G.W. Barton, R. Lwin, M.C.J. Large and L. Poladian, “Fabrication of microstructured optical fibers Part II: Numerical modeling of steady-state draw process,” J. Lightwave Technol. 23, 2255-2266 (2005). 31. S.C. Xue, M.C.J. Large, G.W. Barton, R.I. Tanner, L. Poladian and R. Lwin, “Role of material properties and drawing conditions in the fabrication of microstructured optical fibers,” J. Lightwave Technol. 24, 853-860 (2006). 32. Steven G. Johnson, Mihai Ibanescu, M. Skorobogatiy, Ori Weisberg, Torkel D. Engeness, Marin Soljačić, Steven A. Jacobs, J. D. Joannopoulos, and Yoel Fink, “Low-loss asymptotically single-mode propagation in large-core OmniGuide fibers,” Opt. Express 9, 748 (2001), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-9-13-748 33. L.J. Cummings, P.D. Howell, “On the evolution of non-axisymmetric viscous fibres with surface tension, inertia and gravity,” J. Fluid mech. 389, 361-389 (1999). 34. R.B. Bird, R.C. Amstrong, and O. Hassager. Dynamics of polymeric liquids, vol. 1, Fluid mechanics, (John Wiley & Sons, New York, 1987). 35. H.M. Reeve, A.M. Mescher and A.F. Emery, “Investigation of steady-state drawing force and heat transfer in polymer optical fiber manufacturing,” J. Heat Transfer 126, 236-243 (2004). 36. S. Wu, “Surface and interfacial tensions of polymer melts. II. Poly(methylmethacrylate), poly(nbutylmethacrylate), and polystyrene,” J. Phys. Chem. 74, 632-638 (1970). 37. T. Sedlacek, M. Zatloukal, P. Filip, A. Boldizar and P. Saha “On the effect of pressure on the shear and elongational viscosities of polymer melts,” Polymer eng. sci., 44, 1328-1337 (2004). 1. Introduction Hollow core multilayer and microstructured optical fibers (MOF) for radiation guiding in the near and mid-infrared (IR) [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17] have recently received close attention as they promise considerable advantage over their solid core counterparts in applications related to power guidance at almost any IR wavelength for military, #69874 - $15.00 USD (C) 2006 OSA Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5839 industry and medical applications, IR imaging and sensing, and even THz transmission. Due to its complexity, fabrication of such waveguides is an active field of research. Four main methods have been identified for hollow core fiber manufacturing, each offering its own advantages and challenges. First method is a deposition of metallo-dielectric films on the inside of a drawn capillary by liquid-phase coating [5, 6]; technical challenges in enforcing thickness uniformity in the resultant coatings limit fiber length to the distances of ∼ 10m. Second method is a capillary stacking or preform drilling [7, 10, 8, 9, 11, 12, 13] where glass or plastic capillaries are arranged in a predefined manner and then co-drawn, or the predefined array of holes is drilled in a plastic preform; so far such fibers have been mostly demonstrated to guide below 3µ m due to the non-transparency of silica and polymer materials used in the fabrication. Third method is a deposition of radially uniform thin films on a drawn substrate fiber by means of physical or chemical vapor deposition methods [14]; main challenge of this technology is presumably uniformity of the resultant coatings, and a throughput due to a relatively slow deposition process. Finally, film rolling process [15] starts with a deposition of a glass (chalcogenide) film on top of a polymer film with consecutive rolling around a mandrel tube, tube etching and drawing; potential challenges include fiber profile optimization which is somewhat nontrivial due to a strictly periodic reflector geometry imposed by the fabrication method, another potential challenge is controlling bio-compatibility of a resultant fiber. In our research group we study fabrication of all-polymer hollow multilayer Bragg fibers. We find that polymer material combination is very well suited for this purpose. The key issue for Bragg fiber fabrication is a choice of a pair of thermo-mechanically compatible polymer materials with sufficiently different refractive indexes that can be co-drawn into a fiber. Because of the commercial abundance of various moderately priced polymer materials the process of material selection is relatively short time and low cost. Moreover, polymer thermo-mechanical properties can be fine-tuned by simple variation of their molecular mass. Many polymers are available in various physical forms such as rods, tubes, films, granules and powders which enables multiple design decisions for preform fabrication. Many polymers are swelled or solvable in organic solvent which enables incorporation of active materials in their matrix such as dopants, laser dyes, nanoparticles, etc. Polymer materials are generally bio-friendly and some of them are biodegradable which is of potential interest for fiber based bio-medical technologies [18]. Finally, drawing of large core Bragg fibers does not typically require core pressurization which simplifies fabrication infrastructure. Although refractive index contrast between layers in an all-polymer Bragg fiber is relatively small (at most 1.3/1.7), as demonstrated in [19] liquid core all-polymer Bragg fibers can be designed to guide very well both TE and TM like modes, while gas filled all-polymer Bragg fibers can guide effectively a TE polarized mode. Moreover, in liquid filled core fibers high index-contrast omnidirectional designs typically hurt fiber performance making lower index-contrast systems superior in their transmission properties. Recently, our research group has succeeded in developing two methodologies for fabrication of multilayered all-polymer hollow preforms. One approach uses consecutive deposition of layers of two different polymers by solvent evaporation on the inside of a rotating polymer cladding tube [20]. Orthogonal solvents were found, and solvent evaporation process was developed for both PMMA(Polymethyl methylacrylate)/PS(Polystyrene) and PVDF(Polyvinylidene fluoride)/PC(Polycarbonate) material combinations. In the left of Fig. 1(a), a 30cm long allpolymer preform with 10 consecutive PMMA/PS layers deposited on the inside of a PMMA cladding tube is presented; on the right, preform cross section is shown. In the left of Fig. 1(b), again a 30cm long all-polymer preform with 15 consecutive PVDF/PC layers deposited on the inside of a PC cladding tube is presented; while on the right, preform cross section is shown. Alternative preform fabrication method uses a co-rolling of two dissimilar polymer films similarly to [15], where both commercial and home-made films were used. In the left of Fig. 1(c), #69874 - $15.00 USD (C) 2006 OSA Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5840 an end part of an all-polymer preform with 19 consecutive PVDF/PC layers is presented; in the middle of the figure a cross section of a drawn fiber is shown with a drawdown ration of 1:20; finally in the right of Fig. 1(c) another drawn fiber with 32 layers PMMA/PS layers is demonstrated with a similar drawdown ratio. PMMA/PS PVDF/PC 100µm 100µm PM PMMA tube tu PVDF PC PVDF PC PVDF PC PVDF PS PS PS PS PS PMMA PM PMMA PM PM PMMA PM PMMA PM PMMA a) c) PVDF/PC 10 µm PC PVDF PC PVDF c) PC PVDF PMMA/PS Fig. 1. a) Left - 30cm long all-polymer preform with 10 consecutive PMMA/PS layers deposited on the inside of a PMMA cladding tube. Right - PMMA/PS preform cross section. b) Left - 30cm long all-polymer preform with 15 consecutive PVDF/PC layers deposited on the inside of a PC cladding tube. Right - PVDF/PC preform cross section c) Left - end part of a rolled 19 layer PVDF/PC preform. Middle - cross section of a drawn PVDF/PC fiber with a 1:20 drawdown ratio. Right - cross section of a drawn PMMA/PS fiber with a 1:20 drawdown ratio. After preform is fabricated, hollow MOFs are manufactured by preform heating and drawing. Geometry of the final fiber can be significantly influenced by controlling various parameters in the drawing process such as temperature distribution in a furnace, fiber drawing and preform feed velocities, as well as pressurization of the hollow core. When preform is made of polymer materials non-Newtonian nature of polymer viscosity can be of importance. Moreover, if several materials are used in a single preform, drawing process can be influenced greatly by the mismatch in the viscosities of the constitutive materials. For the problem of hollow Bragg fiber drawing, geometries of all the resultant fibers can be simply parametrized by only two parameters - a drawdown ratio, and a hole collapse parameter. We further find that fiber transmission properties depend strongly on a hole collapse parameter which can only be controlled indirectly during drawing process. The purpose of this paper is to characterize hole collapse and multilayer non-uniformity during the draw of a polymer multilayer fiber as a function of standard control variables. We investigate importance of mismatch in the viscosities of two different materials during co-drawing and compare predictions of Newtonian and generalized Newtonian flow models. Previous studies on fiber drawing have focused mainly on spinning molten threadlines [21, 22] or drawing conventional solid optical fibers [23, 24]. Drawing of hollow fibers was first studied in [25] where the asymptotic “thin-filament“ equations were obtained but the effects of surface tension were neglected. A more complete analysis, although confined to Newtonian flow is given in [26, 27, 28, 29, 30, 31]. When preform radius is comparable to neck-down length (which means that the “thin-filament“ assumption is no more valid), the surface tension force is relatively small and the tube thickness is comparable to the hole radius, it is demonstrated in [29, 30, 31] than the hole experiences an expansion instead of a collapse. The hole enlargement is also observed in the fabrication of microstructured polymer fibers reported in [17], where a two-stage draw process was used. The analysis presented here is more pertinent #69874 - $15.00 USD (C) 2006 OSA Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5841 to a situation where the geometry is slender and the surface tension effects are important which could be the case of a second draw stage. 2. Effect of a hole collapse on the transmission properties of hollow Bragg fibers 102 5 rfi = Cr rfd i 10 4 3 rfd o rfd o 0.7 10 2 rfd i rfi 0.8 0.9 10 1 Cr=1.0 (design) Cr=1.0 0.4 0.6 0.8 1.0 Cr=0.6 1.2 1.4 λ (µm) a) 1.6 Water filled low index-contrast fiber 101 100 10-1 0.6 0.7 10-2 0.8 C =1.0 10-3 0.6 10 HE11 radiation loss (dB/m) TE 01 radiation loss (dB/m) 10 TE 01 radiation loss (dB/m) Air filled high index-contrast fiber 1.8 2.0 0.9 102 101 100 0.6 10-1 0.7 0.8 C =1.0 10-2 0.7 2.2 b) 0.8 0.9 1.0 0.9 1.1 1.2 1.3 1.4 1.5 1.6 1.7 λ (µm) Fig. 2. a) Radiation loss of the bandgap guided T E01 core modes for the high index-contrast (2.0/1.5) air filled fibers with different hole collapse ratios Cr , while the same outside radii ft Ro . Hole collapse leads to the shift of a bandgap center into the longer wavelength, as well as to a considerable increase in the modal radiation losses. b) Radiation losses of the bandgap guided T E01 and HE11 core modes for the low index-contrast (1.6/1.4) water filled fibers with different hole collapse ratios. Hole collapse and draw induced non-uniformity of a multilayer reflector have direct impact on the transmission properties of a hollow photonic crystal fiber. We define R if and Rof to be the inside and outside radii of a hollow Bragg fiber (insets in Fig. 2), while R ip and Rop are the corresponding radii of a hollow preform. The first parameter that relates preform and fiber dimensions is a drawdown ratio D which is defined as a ratio of the outside preform diameter to f that of a fiber D = R op /Ro . Drawdown ratio can be set during drawing process, and is typically well maintained by a feedback loop from a laser micrometer to a tractor assembly. The second parameter characterizes hole collapse in a fiber as compared to a preform, and is defined as C r = (Rif /Rof )/(Rip /Rop ). This parameter is important when preform core is made of a compressible material such as gas, which is the case during hollow core Bragg fiber fabrication. Thus, C r = 0 signifies that during drawing, the hole collapses completely resulting in a solid core fiber; while Cr = 1 signifies that there is no hole collapse and all the fiber dimensions can be calculated from the corresponding preform dimensions by a simple division by a drawdown ratio. Given these two drawing parameters and assuming that all the materials in a melted state are incompressible fluids, then a circular contour of radius r p in a preform translates into a circular contour of radius r f in a fiber, and they are related by (the ratio between the areas outside and inside the contour is conserved): R p (Rop )2 − (r p )2 rp . (1) (r f )2 = ( )2 − (1 − Cr2)( i )2 p 2 D D (Ro ) − (Rip )2 To understand the effect of a hole collapse on the transmission properties of the resultant fibers, in Fig. 2(a) we first present a set of theoretical curves showing radiation losses of the T E 01 core modes for the high index-contrast air filled fibers drawn with different values of C r , while featuring the same outside diameter R of . In this example, Cr = 1 corresponds to a target hollow core fiber n c = 1 with a strictly periodic 15 layer quarter-wave reflector having material #69874 - $15.00 USD (C) 2006 OSA Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5842  refractive indices n h = 2.0, nl = 1.5 and layer thicknesses d ht = 0.25λct / n2h − n2c = 144nm,  dlt = 0.25λct / n2l − n2c = 224nm, where λ ct = 1µ m (for more details on design of such fibers see [32]). Target fiber inside and outside radii are chosen to be R if t = 5µ m, Rof t = 12.72µ m. By design, such a fiber has a large band gap centered at λ ct . In the presence of a hole collapse C r < 1 (assuming the same value of a drawdown ratio D) two major changes in the the fiber geometry happen. First, while the outside fiber radius is ft f ft fixed Ro , the fiber core radius is reduced R i = Ri Cr . Second, from Eq. (1) it can be shown that thicknesses of the reflector layers become non-uniform, increasing towards the fiber core, ¯ ∼ d t /Cr . These geometrical changes can while, on average, layer thicknesses increase as dh,l h,l significantly modify fiber transmission spectra. Thus, as the center wavelength λ c of a photonic bandgap is proportional to the average reflector layer thickness, then, in the presence of a hole collapse, center of a bandgap is expected to shift to the longer wavelengths λ c ∼ λct /Cr ( Fig. 2(a)). We find, however, that in the presence of a hole collapse the ratio of a band gap to a mid-gap (relative bandgap) stays almost uneffected. Another prominent effect of a hole collapse is on the core mode radiation losses. From ([32]), f radiation losses of the bandgap guided core modes scale as (λ c ) p−1 /(Ri ) p , where exponent p equals 3 for the T E 0n modes, while for the HE, EH, and TM modes, exponent p is in the range [1, 3] and depends strongly upon the fiber core size. Thus, in the presence of a hole collapse, due to reduction of the core radius, and due to shift in the center of a bandgap, we expect core −(2p−1) mode radiation loss to increase as Loss ∼ Cr , which for T E 01 mode gives Loss ∼ Cr−5 . From more detailed simulations we find that for T E 01 mode, actual scaling exponent varies from −5 when Cr ≃ 1 to almost −7 when Cr ≃ 0.5 signifying additional degradation of modal confinement due to nonuniformity in the reflector layer thicknesses. In Fig. 2(b) we present another example of the effect of a core collapse on the radiation losses of the T E01 and HE11 core modes of the low index-contrast water filled fibers. In this example, a target water filled core n c = 1.332 fiber features a strictly periodic 16 layer quarterwave reflector with refractive indices n h = 1.6, nl = 1.4 and layer thicknesses d ht = 282nm, dlt = 580nm, where band gap center frequency is λ ct = 1µ m; inside and outside radii are chosen ft ft to be Ri = 10µ m, Ro = 21.90µ m. Behavior of modal losses as a function of a core collapse parameter in this case is similar to that of a high index-contrast system presented in Fig. 2(a). However, in this particular case due to accidental degeneracy of the studied modes with higher order ones, for some values of C r modal losses can exhibit a double dip profile instead of a single dip profile corresponding to a classical band gap. From the analysis above it follows that hole collapse mainly leads to the linear shift in the bandgap frequency and a super-linear increase in the radiation losses of the core guided modes. In what follows we quantify hole collapse and layer non-uniformity as a function of the standard control parameters during drawing. 3. Basic equations Schematic of a hollow multilayer preform profile during drawing is shown in Fig. 3. As the density of polymers varies little in the range of temperatures and pressures considered here we consider an incompressible axisymmetric steady flow. The equations for conservation of mass and momentum in cylindrical coordinates are as follows: 1 ∂ (rvr ) ∂ vz + =0 r ∂r ∂z #69874 - $15.00 USD (C) 2006 OSA (2) Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5843   ∂ vr ∂ vr ∂ p 1 ∂ (rτrr ) τθ θ ∂ τrz + vz + − + ρ vr =− ∂r ∂z ∂r r ∂r r ∂z   ∂ vz ∂ vz ∂ p 1 ∂ (rτrz ) ∂ τzz + vz + + + ρg ρ vr =− ∂r ∂z ∂z r ∂r ∂z (3) where r and z are the radial and axial coordinates, v r and vz are the r and z components of the velocity vector v, ρ is a constant density, p is a pressure, τ i j is an extra-stress and g is a gravitational acceleration. The components of a total stress tensor σ̄¯ are σi j = −pδi j + τi j (4) The definition of τ i j depends upon the polymer model, and is discussed in details later. For these equations, we need to specify the boundary conditions. At the interfaces between different layers the kinematic conditions are vr = R′j vz at r = Rj (5) where R j = R j (z) denote the interfaces between layers and the index j = 1, 2 . . . N is used to number them starting from the inner one. The primes denote the derivative with respect to z. Since the first and the N-th interfaces are external interfaces, we will distinguish them by denoting R i ≡ R1 and Ro ≡ RN for the inner and outer boundaries respectively. 0 Ri(z) Vf Ro(z) r Hollow core p = Pi to ni no ti Vd z Fig. 3. Schematic of a hollow multilayer preform during drawing. Different colors correspond to different materials in a multilayer. Hollow core can be pressurized in order to control its collapse under the action of a surface tension. In this case, at the inner interface the dynamic boundary conditions are σ̄¯ · ni = (γκi − Pi ) ni σ̄¯ · ti = 0 #69874 - $15.00 USD (C) 2006 OSA (6) Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5844 where γ denotes the surface tension coefficient, and κi = 1 Ri (1 + R′i2 )1/2 − R′′i (1 + R′i2 )3/2 (7) is the curvature, while Pi is the hole overpressure (constant ambient pressure has no effect on the flow). Outward-pointing normal at the inner boundary n i is defined as: ⎞ ⎛ ′ R 1 ⎠, , 0,  i ni T = (nr , nθ , nz ) = ⎝−  (8) 2 ′ 1 + Ri 1 + R′i2 while ti T = (nz , 0, −nr ) (9) is the unit tangent vector. In a similar way, the dynamic boundary conditions at the outer boundary are σ̄¯ · no = −γκo no σ̄¯ · to = 0 (10) where no , to and κo satisfy the same equations as −n i , −ti and κi respectively with R i replaced by Ro . At the interfaces between the internal layers, we will consider a continuous stress and velocity. In the axial direction the boundary conditions are the known values of the drawing (V d ) and the feeding (V f ) velocities. Furthermore, as an initial condition, the values R j (0) are known. 4. Thin filament equations One of the basic dimensionless parameters in the problem is the ratio between the preform radius and the length of the neck down region defined as ε . In the case when ε ≪ 1 a so called thin filament approximation can be used. There are two different approaches for simplifying the equations in this case. In the first approach [26, 33], the variables are expanded as power series in ε 2 and only the dominant terms are retained in the equations. In the second approach [25], which we also follow in this paper, the equations are averaged over the cross-section at each value of z. The average ϕ (z) of a variable ϕ (z) is defined as ϕ (z) = 1 2 π Ro − Ri 2 Ro 2π rϕ (r, z)dr (11) Ri For the axial velocity, the assumption v z = vz is made explicitly. We note first that for a slow varying thin filament; R ′j ≪ 1, and by neglecting terms of the order R ′ 2j the boundary conditions Eqs. (6) and (10) take the following form: #69874 - $15.00 USD (C) 2006 OSA Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5845 γ − Pi R  i γ σrz = − + Pi R′i Ri   γ 2 − Pi R′ i σzz = Ri σrr = r = Ri at (12) σrr = − γ Ro γ ′ R Ro o γ 2 σzz = − R′ o Ro σrz = at r = Ro Multiplying the r-component of the momentum equation (3) by 2π r 2 , integrating from R i to d Ro 2 Ro , considering dz Ri 2π r τrz dr ≈ 0 (see Ref. [34] p.382), neglecting the inertial term because of the small value of the radial velocity and using the boundary values of σ rr given by Eq. (12), we obtain τ rr + τ θ θ γ (Ro + Ri ) − R2i Pi (13) p= + 2 R2o − R2i Multiplying the z-component of the momentum equation (3) by 2π r, integrating from R i to Ro , using the boundary values of τ rz = σrz given by Eq. (12) as well as the Eq. (13) and neglecting terms of relative order R ′ 2j , we obtain ρ Qv′z ′   Q τ rr + τ θ θ Q = τ zz − + γ (Ro + Ri ) + ρ g vz 2 vz (14) where π Q = π R2o − R2i vz is the constant volumetric flow rate. This is the axial force balance equation. 5. Newtonian flow The constitutive equation for the Newtonian fluid is τ = η (r, z) ∇v + ∇vT (15) where η is the viscosity which in our case depends on r and z. In this paper we consider a uniform temperature in any given cross-section, thus T is only a function of z. We also assume an axial variation of the temperature dependent viscosity for each constituent material. From the continuity equation one finds: rv′ A (16) vr = − z + 2 r where A = A(z) is a function to be determined later. The extra-stress tensor takes the form: ⎛  ⎞  0 0 −η v′z + 2A r2 ⎟ ⎜   ⎟ τ =⎜ (17) 0 −η v′z − 2A 0 ⎝ ⎠ 2 r ′ 0 0 2η v z #69874 - $15.00 USD (C) 2006 OSA Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5846 where the component τ rz = η ∂∂vzr has been neglected. From Eq. (17) the following relation holds for the averaged components of an extra-stress tensor τ zz − τ rr + τ θ θ = 3η v′z 2 and the axial force balance Eq. (14) can be written in the form ′ v′ Q ρ Qv′z = 3η Q z + γ (Ro + Ri ) + ρ g vz vz (18) (19) The expression for A is obtained by directly integrating the r-component of the momentum equation (3) from R i to Ro and using the boundary values of σ rr given by Eq. (12)   1 1 + Pi − γ Ri Ro (20) A= Ro η (r) 4 dr r3 Ri Combining the kinematic boundary conditions Eq. (5) with relation Eq. (16), we obtain ′ R2j vz = 2A (21) 6. Numerical solution for the Newtonian case The last three equations may be considered as a system of coupled differential equations for vz (z), A(z) and R j (z). This system of equations can be solved with an iterative method. Starting from an arbitrary initial distribution of v z (z), say a linear distribution between the feeding and drawing speed, the initial value equations (21) can be integrated in order to obtain R j (z) with the value of A given by Eq. (20). This functions are then used to solve the boundary value problem Eq. (19) to obtain a new function v z (z) passing so at the next iteration. For the examples given later in this paper we have tested this procedure and it converges very fast (less than 200 iterations). In most cases of practical importance, inertial, gravitational and capillary terms in Eq. (19) can be neglected and it takes a simple form η v′z =C vz (22) where C is a constant. This equation, which is now uncoupled from the other two, can be easily integrated to give ⎛ ⎞ z dz ⎜ Vd ⎟ 0 η (z) ⎟ vz (z) = exp ⎜ (23) ⎝lnV f + L dz ln V f ⎠ 0 η (z) where L is a furnace length. Once the axial velocity is known, the initial value equations (21) can be easily integrated to obtain the profile of a drawn structure. One of the key aspects of hollow preform drawing is the partial or even complete collapse of a compressible core as a result of the surface tension forces acting at the free boundaries. As introduced earlier, we characterize the hole collapse by the parameter C r . Hole collapse typically results in a faster reduction of a smaller core radius compared to the larger outer radius. Thus, starting with identical thicknesses of the same material layers in a preform, in #69874 - $15.00 USD (C) 2006 OSA Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5847 z (m) 0.3 0.25 0.2 0.15 0.1 0.05 0 120 140 160 180 T ( oC) Fig. 4. Temperature distribution in the furnace. a drawn fiber the inner layers will become thicker than the outer ones. We will characterize the thickness non-uniformity by the ratio hhoi between the thickness of the outer layer and the thickness of the inner one, assuming they were equal in a preform. In the following we investigate how the hole collapse and layer thickness non-uniformity are influenced by various control parameters. As an example, we consider drawing of a multilayer hollow Bragg fiber preform where cladding tube and one of the two materials of a multilayer is PMMA, while the other material is a different polymer. Materials in consecutive polymer layers are alternated to create a periodic multilayer structure. 6.1. Effects of draw ratio, temperature and viscosity mismatch In our calculations we assume a uniaxial temperature distribution with a maximum at the furnace center (Fig. 4).In the following when we vary the maximum value of the temperature, we simply rescale the whole profile. We assume that Newtonian viscosity of PMMA obeys an Arrhenius type dependency:   1 1 − η (T ) = η0 exp α (24) T T0 where T (◦ C) is the temperature, η 0 (Pa) and α (◦ C) are constant coefficients and T0 is a reference temperature. For PMMA their values are given in Ref. [35]; η 0 = 1.506 × 10 5Pa, α = 2935◦C and T0 = 170◦ C. In our first calculation viscosity of the other polymer is assumed to be simply two times higher than that of PMMA. Although we recognize that to model correctly the flow of a particular polymer we need to use its proper temperature dependent viscosity, in our first simulation we rather want to highlight the major effect of adding another material into a preform. Particularly, we want to investigate how the hole collapse is affected by the viscosity of the second material despite of its small volume fraction in the preform. We consider drawing of a preform with external and internal diameter 31.75 and 25.4mm respectively. PMMA tube is coated on the inside with 25 alternated layers of PMMA and another polymer with a viscosity two times higher, each one of them having a thickness of 50µ m. The value of surface tension coefficient is considered constant for exterior interfaces γ = 0.032N/m #69874 - $15.00 USD (C) 2006 OSA Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5848 [36] and both densities are considered to be 1195kg/m 3 [35]. We also assume a furnace length L = 30cm, a constant preform feeding velocity V f = 2.5µ m/s and a zero pressurization Pi = 0. First, we consider the effects of varying the drawing ratio defined as D r = Vd /V f (not to be confused with a drawdown ratio defined as D = R op /Rof ) and the maximum temperature in a furnace. In Fig. 5(a), solid lines represent the parameter of a hole collapse C r as a function of a draw ratio D r = Vd /V f for different values of the maximal temperature. Dashed lines represent the parameter curves resulting in a constant outside diameter D o ≡ 2Ro = 125µ m and Do = 250µ m after the draw. For comparison, in dotted curves we present the hole collapse if no other polymer is present in the preform (drawing of a simple PMMA tube of the same inner and outer radii as a multilayer preform). o o 1 160 C 0.8 180 C 160 C 1 o 170 C o 170 C 0.6 o 180 C 0.4 Do = 250 µ m Do = 125 µ m 0.2 0 o 190 C 0 1 2 a) 3 4 Dr =Vd /Vf 5 6 h o /h i Cr 0.8 0.6 Do = 250 µ m Do = 150 µ m 0.4 o 190 C 0.2 0 7 x10 4 o b) 0 1 2 3 4 D r =Vd /Vf 5 6 7 4 x10 Fig. 5. (a) Hole collapse parameter Cr as a function of the draw ratio Dr for different temperatures. Solid lines correspond to multilayer preform. Dotted lines correspond to a simple tube with the same thickness as the multilayer preform. Dashed lines represent the curves of a constant outside diameter. (b) Ratio ho /hi between the inner and outer layer thicknesses as a function of the draw ratio for different temperatures. We see that to prevent hole collapse higher draw ratios and lower temperatures have to be used. Both cases demand higher draw force which might lead to fiber breakage. We also observe that the collapse of the tube with the same thickness as the multilayer preform is more pronounced because of a lower average viscosity. Effect of the temperature on a hole collapse is easy to understand as viscosity decreases rapidly with reduction of the temperature, thus hindering the hole collapse. Effect of the draw ratio is more subtle. Starting with preforms of the same diameter, draw ratio increase leads to the reduction of a resultant fiber cross-section. As a consequence, the forces of surface tension become more pronounced, thus favoring the hole collapse. This is overcompensated by the fact that increasing the draw ratio leads to the higher axial velocities, thus the time a cross-section spends in a melted zone diminishes which works against the hole collapse. In Fig. 5(b), thickness non-uniformity parameter is presented as a function of the draw ratio for different values of the maximum temperature. The curves are similar to those for the hole collapse in Fig. 5(a). From mass conservation, the ratio between the cross-section areas of different layers remains constant from which it follows that the thickness non-uniformity parameter is proportional to the hole collapse. We now describe in more details the effect of mismatch in the polymer viscosities when two different materials are used in the same preform. As seen in Fig. 5(a) the hole collapse is considerably less pronounced for the multilayer structure despite the fact that the higher viscosity material occupies only a very small fraction of the total volume and its viscosity is only two times higher than that of PMMA. In what follows we assume that multilayer preform is made of PMMA and another polymer. For the viscosity of a second material a similar Arrhenius law as for PMMA is assumed. #69874 - $15.00 USD (C) 2006 OSA Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5849 First, we investigate the effects of changing η 0 while keeping the other parameters unchanged, which corresponds to the case of using the same polymer, but with a different molecular mass. It should be mentioned here that the polymer’s molecular mass determines whether the fiber is drawable in the first place. Second, we investigate the effects of changing T 0 which corresponds to the case of varying the polymer material. In Fig. 6 we consider drawing of preforms of various compositions at a fixed maximum furnace temperature of 190 ◦ C and a draw rate Dr = 30000. Multilayer preform geometry is the same as described above; mismatch in the polymer viscosities is described in terms of the ratios of the material parameters η 0 /η0,PMMA and (T0 − T0,PMMA )/T0,PMMA . From Fig. 6 we see that the hole collapse depends significatively on the viscosity of a second material and can be prevented by choosing a polymer with an appropriate viscosity. 1 1 0.8 0.8 0.6 0.6 ho/hi 0.4 Cr Cr 0.4 0.2 0.2 0 ho/hi 0 20 40 a) 60 0 / 80 0 100 0 -5 b) 0,PMMA 5 10 15 20 (T0 - T0,PMMA)/T0,PMMA 25 x10-2 Fig. 6. Effect of mismatch in the viscosities of materials in a multilayer on hole collapse (solid lines), and layer non-uniformity (dotted lines).Maximum furnace temperature is T = 190◦ C, draw ratio Dr = 30000. (a) Effects of η0 . (b) Effects of (T0 − T0,PMMA )/T0,PMMA . 6.2. Effects of the pressurization and preform feeding velocity Other parameters that influence hole collapse are the hole overpressure and the feeding speed. By increasing the hole pressure we expect to reduce the hole collapse. Also, for a given draw ratio, by increasing the preform feeding speed we expect reduction of the hole collapse as fiber cross-section would spend less time in the melted zone. 1 1 0.8 ho/hi 0.6 Cr 0.6 0.4 0.4 Cr 0.2 0 ho/hi 0.8 0 1 2 a) 3 4 Pi (Pa) 5 0.2 6 7 0 0 b) 5 10 15 20 25 Vf ( m/s) 30 35 40 Fig. 7. The hole collapse and thickness non-uniformity as a function of the hole overpressure and feeding speed. Maximal furnace temperature is T = 190◦ C and draw ratio Dr = 5000. (a) Effects of hole overpressure Pi . (b) Effects of feeding speed V f . We consider drawing of the same preform as in a previous section using the draw ratio Dr = 5000, while the other draw parameters remain unchanged. In Fig. 7 hole collapse and layer thickness non-uniformity are presented as functions of the hole overpressure P i and feeding speed V f . As expected, the hole collapse is very sensitive to pressurization, and in principle, can #69874 - $15.00 USD (C) 2006 OSA Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5850 be reduced by increasing the pressure. Time evolving transient draw simulations, not discussed in this paper, also show that above a certain critical value for an overpressure, which in our case is less than 7Pa, even if the fiber does not blow up immediately, the drawing process never reaches its steady state. A more subtle way of controlling the hole collapse is by changing the preform feeding speed, although, for a given draw ratio, this could be limited by the maximal draw velocity. 7. Generalized Newtonian model In the previous sections we considered only the temperature dependence of the viscosity, but for polymers such as PMMA and PS the viscosity depends also on the kinematics of the flow. In this case we refer to η as a non-Newtonian viscosity or the generalized viscosity. Polymers in general are also viscoelastic, however we will neglect elasticity effects in our analysis which may be considered as a first step toward a complete modelling of polymeric flow. For the non-Newtonian viscosity we will use the Carreau-Yasuda model [34] which was also considered in [37] and the values of different coefficients for some important polymers are given in that paper. According to this model the viscosity is given in the form η (T, IID ) =  η0 f √ a  1−n a 1 + K1 f 2IID where η0 , K1 , a,and n  are constant coefficients, D = 1 2 (25) ∇v + ∇vT is the rate of the deformation tensor, IID = tr D · D is its second invariant and the term f gives the temperature dependence α  1 1 T − T0  . Considering the Eqs. (17) and (15) and of the viscosity in the Arrhenius form f = e the fact that A is negligible compared to v ′z , one obtains IID = 23 v′z 2 . The numerical procedure is the same as in the Newtonian case. The system of coupled partial differential equations (19-21) can be solved again by the same iterative procedure. At each step of the iteration we must solve the boundary value problem Eq. (19) which, in this case, becomes more complex because the viscosity depends on the rate of deformation. Thus, this equation must be solved itself by an iterative procedure (so in overall there are two nested iterations). However, since the inertial and gravitational terms are usually negligible, the axial velocity vz (z) can be evaluated separately by using repeatedly Eq. (23) with a viscosity update according to Eq. (25), until the convergence is attained. Once v z (z) is known, Eq. (21) can be easily integrated. As an example, we consider drawing of the same preform as for the Newtonian case, which is a PMMA tube with OD and ID diameters 31.75 and 25.4mm respectively, coated on the inside with 25 alternated layers of PMMA. Viscosity of PS is assumed to be two times higher than that of a PMMA. Layer thickness is assumed to be 50µ m for all the layers. We consider γ = 0.032N/m, ρ =1195kg/m 3, L = 30cm, V f = 2.5µ m/s and Pi = 0. Parameters for the nonNewtonian viscosity Eq. (25) are K 1 = 0.0861 n = 0.1401 and a = 0.7347 for PMMA, and K1 = 0.3891 n = 0.2194 and a = 0.6097 for PS [37]. Finally we assume α = 2935, T 0 = 170◦C, and the value of the maximum temperature in the furnace being T = 180 ◦C. Numerical results are presented in Fig.8(a) where hole collapse parameter is plotted as a function of the draw rate and maximum furnace temperature for Newtonian and non-Newtonian models. We observe that hole collapse is more pronounced in the non-Newtonian case because, as seen from Eq. (25), the non-Newtonian viscosity is lower than the Newtonian one, making it easier for the surface tension to reduce the hole. As expected, non-Newtonian nature of polymer viscosity becomes more apparent with increase in D r . This can be clearly seen in Fig. 8(b) where the viscosity distribution is plotted as a function of z for different draw ratios. For #69874 - $15.00 USD (C) 2006 OSA Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5851 comparison, we have also included the Newtonian viscosity at the same temperature. 1 o 0.3 o 0.25 o 0.2 0.15 0.1 160 C 0.8 170 C 0.4 180 C 0.2 190 C 0 z (m) Cr 0.6 o 0 1 a) 2 3 4 Dr =Vd /Vf 5 6 0.05 0 x10 4 b) Newtonian D = 5000 r Dr= 35000 Dr= 65000 11 12 13 14 log(η(Pa.s)) 15 Fig. 8. a) Comparison between Newtonian and generalized Newtonian model. Solid lines correspond to generalized Newtonian model and dotted lines to Newtonian model. b) Viscosity distribution in the furnace for different draw ratios at T = 180◦ C. Solid lines correspond to generalized Newtonian model and dotted line to Newtonian model. 8. Conclusions Drawing of multilayer hollow polymer fibers is studied using the thin filament approximation. Both, the Newtonian and generalized Newtonian models of polymer flow are considered. Hole collapse is identified as a key parameter effecting transmission properties of the resultant hollow Bragg fiber. The hole collapse caused by surface tension is characterized as well as the closely related layer thickness non-uniformity. It is demonstrated that by varying various control parameters such as furnace temperature, feeding speed and pressurization it is possible to reduce hole collapse. While hole pressurization provides a very effective way of compensating for the hole collapse, it is found that the final fiber dimensions are very sensitive to the value of an overpressure. Moreover, the draw process could not reach a steady state if the overpressure was larger than a critical value. Finally, under the same draw conditions, the hole collapse is more pronounced when non-Newtonian viscosity is taken into account. Comparison between the experimental results and theoretical modelling is currently underway. In general, qualitative agreement between theory and experiment of the dependence of a core collapse on the draw parameters is readily observable. As the results of our numerical analysis correspond to the equilibrium drawing conditions, to make a quantitative comparison with experiments we have to ensure that our experimental drawing has reached equilibrium and is not dominated by transients. Acknowledgment Funding for this work comes partially from the contribution of Canada Research Chairs, NSERC, FQRNT, Canada Institute for Photonic Innovations, projects BP5 and FP3. #69874 - $15.00 USD (C) 2006 OSA Received 10 April 2006; revised 7 June 2006; accepted 9 June 2006 26 June 2006 / Vol. 14, No. 13 / OPTICS EXPRESS 5852