Journal of Photonics

Volume 2014 (2014), Article ID 645207, 11 pages

http://dx.doi.org/10.1155/2014/645207

Research Article

Heat Conduction Modeling of Fiber Fuse in Single-Mode Optical Fibers

Yoshito Shuto

Ofra Project, Iruma 358-0023, Japan

Received 16 September 2013; Revised 28 December 2013; Accepted 8 January 2014; Published 24 February 2014

Academic Editor: Jorge Diego Marconi

Copyright © 2014 Yoshito Shuto. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The unsteady-state thermal conduction process in single-mode optical fiber was studied theoretically with the explicit finite-difference method. We assumed that the vitreous silica optical fiber underwent pyrolysis at elevated temperatures to form (). We also proposed a model in which the optical absorption coefficient of the core layer increased with increasing molar concentration of . The core-center temperature changed suddenly and reached over  K when a 1.064-μm laser power of 2 W was input into a short core layer of 40 μm length, which was heated at 2923 K. This thermal wave, that is, a fiber fuse, increased in size and propagated toward the light source at a rate of about 0.54 m/s. The calculated propagation velocity of the fiber fuse was in agreement with the experimental value. Moreover, the average temperature of the radiated region of the core layer gradually approached a temperature of about 5700 K. It was found that the final average temperature was close to the experimentally reported values.

1. Introduction

By utilizing dense wavelength-division-multiplexed (DWDM) transmission systems using an optical-fiber amplifier, large amount of data can be exchanged at a rate of over 60 Tbit/s [1]. In connection with this achievement, the danger of the “fiber fuse phenomenon” occurring has been pointed out. This occurs when high-power (W order) optical signals are transmitted in an optical-fiber cable [2].

The fiber fuse phenomenon was first observed in 1987 by British investigators [3–6]. Several review articles [7–9] have been recently published, which cover many aspects of the current understanding of the fiber fuse. Most experimental results of fiber fuse generation have focused on an intensity of 0.35–25 MW/cm² [3–21]. This is many orders of magnitude below the intrinsic damage limit for silica, which exceeds 10 GW/cm² [3]. A fiber fuse can be generated by bringing the end of a fiber into contact with an absorbing material or melting a small region of a fiber by using an arc discharge of a fusion splice machine [3]. If a fiber fuse is generated, an intense blue-white flash occurs in the fiber core, and this flash propagates along the core in the direction of the optical power source at a velocity on the order of 1 m/s. Fuses are terminated by gradually reducing the laser power to provide a termination threshold at which the speed of the fuse is reduced to zero.

When a fiber fuse is generated, the core layer in which the fuse propagates is seriously damaged, and the damaged fiber cannot be used in an optical communication system. The damage is made manifest by periodic or nonperiodic bullet-shaped cavities left in the core [11, 20]. It was found that molecular oxygen was released and remained in the cavities while maintaining high pressure (about 4 atmospheres) at room temperature [4].

Dynamics of cavity formation have been investigated since the discovery of the fiber fuse phenomenon. Kashyap reported that the cavity shapes were dependent on the nature of the input laser light (CW or pulses) when the input power was kept at 2 W average [4]. When CW light was input, the cavities appeared to be oblong and cylindrically symmetric. On the other hand, short asymmetric cavities were formed by entering (mode-locked) pulses with 100 ps FWHM (full width at half maximum), while long bullet-shaped cavities were observed by entering pulses with 190 ps FWHM [4]. Dianov et al. found that the formation of periodic bullet-shaped cavities was observed after 20–70 μs from the passage of a plasma leading edge [22, 23]. Todoroki reported that fiber fuse propagation was classified into three modes (unstable, unimodal, and cylindrical) according to the plasma volume with respect to the pump beam size [24]. He found that pump power increments and decrements that occurred faster than the plastic flow of the surrounding glass melt led to a cavity-free segment and an irregular cavity pattern, respectively.

These cavities have been considered to be the result of either the classic Rayleigh instability caused by the capillary effect in the molten silica that surrounds a vaporized fiber core [17] or the electrostatic repulsion between negatively charged layers induced at the plasma-molten silica interface [25]. Recently, Todoroki proposed that the bullet-like shape of the cavity train resulted from the formation of an intrinsic bridge inside the hollow silica behind the traveling plasma and the successive compression of detached cavities under a steep temperature gradient along the fiber [8].

It is very important to detect and/or terminate fiber fuse in order to protect expensive equipments for optical communication system. The detection systems for fiber fuse have been recently proposed by two research institutes, which use FBG sensors [26] or an optical time-domain reflectometer [27], respectively. On the other hand, fiber fuse terminators have been investigated by several researchers, which contain tapered fibers [28, 29], thermally diffused expanded core (TEC) fibers [30], and a hole-assisted fiber (HAF) [31].

Several hypotheses have been put forward to explain the fiber fuse phenomenon. These include a chemical reaction involving the exothermal formation of germanium (Ge) defects [10], self-propelled self-focusing [3], thermal lensing of the light in the fiber via a solitary thermal shock wave [5], and the radiative collision of SiO and O complexes [32].

The similarities between the fiber fuse propagation and the combustion flame propagation were pointed out by Facão et al. [33] and Todoroki [8]. A fast detonation-like mode of fiber fuse propagation with a velocity of 3.2 km/s was observed under intense laser radiation intensity of 4000 MW/cm² [34]. It is noteworthy here that a fiber fuse process itself is not an oxidation process such as the combustion and/or the detonation process. As described in this paper and also pointed out in [7], fiber fuse generation is related to a reduction (pyrolysis) process of SiO₂.

The optical absorption coefficient of a fiber core at high temperatures is closely related to the generation of the fiber fuse. Kashyap reported a remarkable increase in the value of a Ge-doped silica core above the critical temperature (~1323 K), while the value of about 0 dB/km at room temperature remained unchanged until the temperature () approached [4]. The value increased by nearly 1900 dB/km (~0.44 m⁻¹) at a wavelength of 1.064 μm when changed from to  K. Furthermore, Kashyap et al. reported that the best fit between the experimental and theoretical fiber fuse velocities was obtained when the value of the Ge-doped silica core at high temperatures of above was fixed to be 4.0 × 10⁴ m⁻¹ at 1.064 μm [13].

Hand and Russell found that this phenomenon was initiated by the generation of large numbers of Ge-related defects at high temperatures of above 1273 K, and the values at a wavelength of 0.5 μm obtained at temperatures of below 1873 K were modeled quite accurately using an Arrhenius equation [5, 6]. By contrast, they reported that the best fit between the experimental and theoretical fiber fuse velocities was obtained when the value of the Ge-doped silica core at 2293 K was assumed to be 5.6 × 10⁴ m⁻¹ at 0.5 μm [5]. This large value, however, could not be estimated using their Arrhenius equation [6]. As the focal length of thermal lens effect is inversely as the value [35], large value of 5.6 × 10⁴ m⁻¹ is necessary to obtain small value of 10 μm order, which is comparable with observed interval and/or large front size of the cavities (see the Appendix).

Furthermore, Hand and Russell reported that the electrical conductivity of the fiber core increased with the temperature and that the hot spot at the fiber fuse center was plasma-like [5]. Kashyap considered that the large values may be attributable to an increase in the value of the fiber core at high temperatures of above [4].

It is well known that silica glass is a good insulator at room temperature, and the electrical conductivity in silica glass below 1073 K is due to positively charged alkali ions moving under the influence of an applied field [36, 37]. The ionic conduction in the glass is not related solely to optical absorption.

We previously investigated the optical absorption mechanism causing the increase in the value and reported the relationship between and in silica glass at high temperatures of above 1273 K [18]. It was found that the increase in loss observed at 1.064 μm can be well explained by the electronic conductivity induced by the thermal ionization of a Ge-doped silica core, and it is not directly related to the absorption of Ge centers. However, the calculated values resulting from the electronic conductivity at 1.064 μm were of 10² m⁻¹ order at 2873 K, about two orders smaller than the values (1.0–4.0 × 10⁴ m⁻¹) reported by Kashyap et al. [4, 13]. Therefore, we need another mechanism to explain the increase in loss at high temperatures of above 2273 K to account for the large (10⁴ m⁻¹ order) values.

To satisfy this requirement, we proposed a thermochemical SiO production model in 2004 [18]. Using this model, we theoretically estimated large values of 10⁴ m⁻¹ order as a result of SiO absorption at high temperatures of 2800 K or above. This model was able to quantitatively explain the relation between the fiber fuse propagation velocity and the incident laser power intensity previously reported by other research institutions.

On the other hand, since the parameters (particularly the light-absorbing parameter) used for the numerical simulation were not optimized, the shortcoming that the maximum temperature of the core center obtained by calculation became unusually high (10⁵ K order) was observed. We have improved this model by optimizing several parameters required for the numerical computation.

In this paper, we describe the improved model in detail. That is, we explain the mechanism of the increased absorption in optical fibers at high temperatures due to the thermochemical production of and estimate high-temperature values at a laser wavelength of 1.064 μm. Then, using these values, we theoretically study the nonsteady-state thermal conduction process in a single-mode optical fiber using the explicit finite-difference technique.

2. High-Temperature Optical Absorption in Optical Fibers

2.1. Effect of Formation on Absorption

It has been reported that, at elevated temperatures, silica glass is thermally decomposed by the reaction [38]

Among the reductants of the silica generated by this pyrolysis reaction with the formula , the most thermally stable material is SiO ().

The internal core space heated at the elevated temperature is shielded from the external conditions. Thus, with increasing temperature, the internal pressure increases and the internal volume decreases. Dianov et al. reported that the internal core temperature is about 10000 K and the internal pressure is about 10000 atmospheres at the time of fiber fuse evolution [29]. It is thought that generated under the high-temperature and high-pressure conditions is densely packed into the internal core space because it is not allowed to expand, and it exists in a liquid-like form. For this reason, it is thought that the optical absorption spectrum of the high density in the core space is similar to that of solid .

Philipp reported that the optical absorption spectrum of a film is similar to that of amorphous Si because of the many Si–Si bonds in a film and that the absorption coefficient for SiO () near the threshold energy should be about one-twentieth of that for amorphous Si  [39]. Furthermore, Philipp estimated the theoretical concentration of Si–() when the five possible tetrahedral conformations centering on Si, Si–(Si₄), Si–(Si₃O), Si–(Si₂O₂), Si–(SiO₃), and Si–(O₄), were completely distributed at random within amorphous [40]. The relationship between and for is illustrated in Figure 1. According to Figure 1, for a SiO () film is 6.25% (=1/16). This is very close to the factor by which the value was reduced (about 1/20) as reported by Philipp [39].

Figure 1: Relative content of Si–(Si₄) in .

The optical absorption near the absorption edge of a film is dominated by that of amorphous Si with the Si–(Si₄) tetrahedral conformation, and the optical absorption coefficient near the absorption edge of a film can be calculated by multiplying the value of for amorphous Si by . That is, if the production rate of given by (1) is denoted by , then at temperature is given by

2.2. Production Rate of

The most thermally stable compound in generated via (1) is SiO (). Therefore, the production rate of in (2) is assumed to be almost equal to the production rate of SiO.

at temperature can be calculated as the ratio of the molar concentration of SiO at temperature to the maximum value . This is the molar concentration when changes to SiO via (1) with a yield of about 100%.

We assumed that the pyrolysis reaction system given by (1) reaches its equilibrium state during fiber fuse propagation. The equilibrium constant for (1) is denoted as . The value of at temperature was calculated using as described below.

First, the initial molar concentration (=0.0366 mol cm⁻³) is denoted by . We consider the case that the concentration decreases to via the pyrolysis reaction of (1). In this case, the molar concentration of SiO and the molar concentration of O₂ are expressed as and , respectively. is given in terms of and as follows:

Rearranging (3), we obtain a cubic equation in . The solution of this equation is given by

where

It is well known that the equilibrium constant is related to the standard Gibbs energy change for (1). The relationship between and is given by [41]

where is the gas constant. The value for (1) is given by

where is the standard production Gibbs energy of a reactant and/or a product.

Vitreous silica () is a solid at the standard temperature (298.15 K). It melts at high temperatures of above 1996 K and becomes a liquid. It also becomes a vapor at temperatures of above 3000 K.

The standard production Gibbs energies ,  , and in each phase have been published [42]. Thus, using these values, we first calculated the standard Gibbs energy change for (1). Next, we calculated by substituting the value into (6). In this way, we computed (=) at temperature by substituting and (=0.0366 mol ) into (4) and (5). The relationship between and is shown in Figure 2. increases with increasing and gradually approaches its maximum value ( mol cm⁻³) at of about 3200 K.

Figure 2: Molar concentration of SiO versus temperature.

Therefore, at temperature was estimated by dividing the value of calculated above by (0.0366 mol cm⁻³).

2.3. Absorption Coefficient of Amorphous Si

The optical absorption spectrum of amorphous Si was reported by Brodsky et al. [43]. The absorption coefficient of amorphous Si in (2) was estimated as follows.

First, we consider the absorption coefficient of an interband transition region, where the photon energy is larger than the energy gap (=1.26 eV) of amorphous Si. In this case, the values of (cm⁻¹ unit) near are given by

where the parameter is . The values calculated by (8) are in good agreement with the experimental values reported by Brodsky et al. [43] for the case of  eV.

Next, we consider the values of in the low-energy region, where . In this region, an absorption edge spectrum is broadened as a result of the interaction between optical phonons and electrons (or excitons). in this region exhibits exponential behavior (referred to as an “exponential tail” or “Urbach tail”) as follows [44]:

where and are parameters and is the Boltzmann constant. is the effective temperature. If the characteristic temperature of phonons is denoted by , then is given by [45]

For Si, = 600 K [46]. and are estimated to be and 0.14, respectively, by simulation using the experimental values reported by Brodsky et al. [43].

It is known that the energy gap of Si decreases linearly with increasing at > 200 K [47]. The temperature dependence of (eV unit) at > 300 K is given by

where is the temperature coefficient of the energy gap and takes a value of 4.2 × 10⁻⁴ eV/K [48].

2.4. Absorption Coefficient of

On the basis of the above results, the temperature dependence of at 1.064 μm ( eV) was calculated using (2) and the values shown in Figure 2. () was used in the calculation. The results are shown in Figure 3.

Figure 3: Absorption coefficient of at 1.064 μm versus temperature. The thick solid line was calculated using (2). The dotted line shows the values estimated by Kashyap et al.

The values (shown as the thick solid line in Figure 3) are about m⁻¹ at 2800 K and about  m⁻¹ at 2950 K. These values are close to  m⁻¹ (shown as the dotted line in Figure 3), which was estimated by simulation using the finite element method [13].

Moreover, it turns out that reaches its maximum value (about  m⁻¹) at about 3150 K then decreases gradually with increasing temperature.

Compared with the stepwise shape (dotted line in Figure 3) of the absorption coefficient assumed by Kashyap et al., the maximum value (about  m⁻¹) of is about 2.4 times the assumed value ( m⁻¹), and the graph has the shape of a somewhat distorted echelon form.

As mentioned above, is produced in the pyrolysis of vitreous silica, and it induces a very large amount of optical absorption with a large absorption coefficient of  m⁻¹ order at high temperatures of above 2800 K. It is thought that the large amount of optical absorption at high temperatures is the cause of the genesis of the fiber fuse phenomenon.

We investigated the thermal conduction behavior within a single-mode (SM) optical fiber by numerical computation using the thermochemical production model. In the next section, we describe the results of some numerical calculations related to the thermal conduction process in an SM optical fiber.

3. Simulation of Fiber Fuse in SM Optical Fiber

3.1. Heat Conduction in SM Optical Fiber

We assume the SM optical fiber to have a radius of and to be in an atmosphere with temperature . We also assume that part of the core layer of length is heated to a temperature of (see Figure 4). Such a region, called the “hot zone” in Figure 4, can be generated by heating the fiber end faces using the arc discharge of a fusion splice machine.

Figure 4: Hot zone in the core layer.

As explained above, the optical absorption coefficient of the core layer in an optical fiber is a function of temperature , and increases with increasing . In the hot zone in Figure 4, the values are larger than those in the remainder of the core region because of its elevated temperature. Thus, as light propagates along the positive direction (away from the light source) in this zone, considerable heat is produced by light absorption.

In the case of a heat source in part of the core layer, the nonsteady heat conduction equation for the temperature field in the SM optical fiber is given by [49]

where ,  , and are the density, specific heat, and thermal conductivity of the fiber, respectively. The values of ,  , and used for the calculation are described in the next subsection.

The last term in (12) represents the heat source resulting from light absorption, which is only required for the hot zone in the fiber core. can be expressed by

where is the optical power intensity in the core layer, which can be estimated by dividing the incident optical power by the effective area of the fiber.

3.2. Heat Conduction Parameters

In the heat conduction calculation for SM optical fibers, we used the following values of   (W ),   (kg ), and (J ) in each temperature range. The unit of is K (Kelvin).

Parameters in the temperature range from room temperature (298 K) to 1996 K [50] are as follows:

Parameters in the temperature range from 1996 to 3000 K [50] are as follows:

Parameters for  K are as follows:

In , the first term in the expression for was estimated by substituting the parameters for SiO in the following equation, which represents the thermal conductivity of diatomic molecules derived from the kinetic theory of gases [51]:

where is Avogadro's number,   (=44.0854) is the molecular weight of SiO, and (=1.5 Å) is half of the collision diameter.

Moreover, in , the second term in the expression for was estimated by multiplying the following thermal conductivity equation for weakly ionized gas plasma [52] by the correction factor of 1/20:

where is the electron mass, is the number density of electrons in the plasma, and is the collision frequency. If we assume that electrons mainly collide with ions with a charge of +1 in the plasma, then is given by [53]

where is the electron charge and is the permittivity of free space. is the so-called Coulomb logarithm, and it takes values of 4–34 in the electron temperature range of  K and the range of 1–  [54]. In our calculation, we used , which corresponds to the case of  K and  .

3.3. Boundary and Initial Conditions for Heat Conduction

We solved (12) using the explicit finite-difference method (FDM) [55] under the boundary and initial conditions described below.

The area in the numerical calculation had a length of (=4 cm) in the axial () direction and a width of (=125 μm) in the radial () direction. There were 24 and 2000 divisions in the and directions, respectively, and we set the calculation time interval to 1 μs. We assumed that the hot zone was located at the center of the fiber (length ) and that the length of the hot zone was 40 μm.

The boundary conditions are as follows.

Since the temperature distribution of the optical fiber is axisymmetric, the amount of heat conducted per unit area (heat flux) in the direction is set to at the fiber center () as follows:

At the outer fiber surface (), the amount of heat conducted per unit area (heat flux) is dissipated by radiative transfer or heat transfer to the open air and/or the jacketing layer () as follows:

where is the Stefan-Boltzmann constant and (0.9) is the emissivity of the surface. is the thickness of the thermal boundary layer and is assumed in our calculation, where is the step size along the -axis.

At the center () of the hot zone, the temperature of the fiber core center is . Also, there is heat inflow along the -axis only at the core center (), which is attributable to light absorption. Moreover, when and , there is neither heat inflow nor heat outflow along the -axis. These conditions are given by

At both fiber ends (), the amount of heat conducted per unit area (heat flux) is dissipated by radiation transfer to the open air () as follows:

In contrast, as initial conditions, in the optical fiber at , except in the hot zone, and the core-center temperature in the hot zone is equal to (>) as follows:

When light propagates through the fiber core () along the direction (away from the light source), the incident laser power decreases because of the nonzero optical absorption coefficient . When the laser light propagates from to along the -axis at , the value is given by

where is the initial laser power. The second term on the right-hand side expresses the optical absorption loss when the light propagates through a distance of .

The results described above assume that the laser light propagates through the fiber core along the positive direction (away from the light source).

When the core layer is heated to above the vaporization point of silica (3273 K), an enclosed hollow cavity is produced in the core center. This cavity contains oxygen, which is produced by the pyrolysis reaction of (1). The heat conductivity of the oxygen (0.03 W ) is two orders smaller than that of the silicate glass. Therefore, the heat transferred in the silica core is stopped at the cavity.

Moreover, as the cavity has a refractive index of , which is smaller than that of the silica core (), the light propagating in the core layer is reflected at the cavity wall. When the light direction is reversed at the cavity, the heat source term resulting from the optical absorption of the reflected light is added to in (13), where is the optical power intensity of the reflected laser light.

We consider the value at a position located near the cavity wall. The length of this position is assumed to be . The laser light reaches , propagates through a distance of , and then reaches the cavity wall, whose coordinate is . Next, the light is reflected at the cavity wall and propagates in the negative direction and then reaches again. In such a case, the value at is given by

where is the reflectivity at the boundary of the silica core and the cavity. is given by

The second integral on the right-hand side of (26), which is related to the reflected light, slightly affects the occurrence of the fiber fuse. However, by taking this term into consideration, the calculated fiber fuse velocities fit the experimental values. For this reason, in the present work, we decided to take into account the effect of the reflected light.

In the following section, we describe the calculated time () dependence of in an SM optical fiber.

3.4. Propagation of Fiber Fuse in SM Optical Fiber

In the calculation, we used  K and  K. It was assumed that laser light of wavelength  m and  W was incident to an SMF-28 optical fiber, which has a core diameter of 8.2 m, a refractive index difference = 0.36%, and  .

We calculated the values at , 11, and 21 ms after the incidence of the 2 W laser light. The calculated results are shown in Figures 5–7, respectively.

Figure 5: Temperature field in SMF-28 optical fiber after 1 ms when  W at  μm.

Figure 6: Temperature field in SMF-28 optical fiber after 11 ms when  W at  μm.

Figure 7: Temperature field in SMF-28 optical fiber after 21 ms when  W at  μm.

As shown in Figure 5, the core center temperature near the end of the hot zone ( mm) changes abruptly to a large value of about  K after 1 ms. This rapid rise in the temperature initiates the fiber fuse phenomenon as shown in Figures 6 and 7. After 11 and 21 ms, the high-temperature front in the core layer reached values of −6.1 and −11.5 mm, respectively. The average propagation velocity was estimated to be 0.54 m/s using these data.

When the laser light of = 1.064 m and  W is incident to the SMF-28 optical fiber, the optical power intensity is  MW/cm. The fiber fuse velocity at this value of is estimated to be about 0.55 m/s (see Figure 10 in [18]). This value is in very good agreement with the upper value obtained by calculation (0.54 m/s).

On the other hand, Hand and Russell measured the fiber fuse temperature to be 5400 K [5], and Dianov et al. obtained a temperature of 4700–10500 K [23, 56] by measurement. They estimated the temperatures from precisely measured spectral data in the 600–1400 nm [5] and 500–800 nm [23, 56] regions, while assuming blackbody radiation.

In our calculation, the temperature distribution of the fiber fuse in the core center is shown in Figure 8. Similar shapes of temperature distribution were reported by Kashyap et al. [13] and Facão et al. [33].

Figure 8: Temperature distribution of the SMF-28 core center versus length along the direction when  W at  μm. The time after the laser incidence is 10 ms.

As shown in Figure 8, it is clear that the sharp temperature peak is located near the light source, and a relatively high-temperature plateau of about 5000 K extends over about 1.5 mm behind the sharp peak. This region, called the “radiation zone” in Figure 8, exhibits high temperatures of above 323 K.

When gaseous SiO and/or molecules are heated to high temperatures above 5000 K, they decompose to form Si and O atoms and finally become and ions and electrons in the ionized gas plasma state. The numbers of electrons and ions in the plasma front, which exhibits sharp temperature peak, are larger than those in the plateau region. However, as the plasma tends to restore electrical neutrality, the motions of the electrons and ions will not produce any change on the initial temperature distribution shown in Figure 8, except for the peak temperature reduction of the plasma front due to energy loss induced by electron-ion collisions.

In the ionized gas plasma, electron-ion collisions generate electromagnetic radiation because of the deceleration during the collisions. This bremsstrahlung emission [57, 58] is a universal and irreducible process of energy loss. If we assume that electrons mainly collide with ions with a charge of +1 in the plasma, the spectral radiance function for the bremsstrahlung emission is given by [59]

where is the number density of electrons in the plasma, is Planck’s constant, and is the optical frequency. This function is directly proportionate to in the case of .

On the other hand, in the case of , the spectral radiance function for blackbody radiation is given by [60]

This is well known as the Rayleigh-Jeans formula, and this function is proportional to , too. Therefore, we assumed that the radiation zone, in which the bremsstrahlung emission of the plasma is liberated, can be treated as a blackbody because of its similar dependence on .

If we consider the radiation zone as a blackbody, that is, isolated from the surrounding nonheated regions, it is expected that the radiation zone will exhibit a radiation spectrum with a broad range of optical frequencies ranging from ultraviolet to infrared. The blackbody temerature of the zone is related to the frequency of the spectral peak as follows [60]:

If is defined as the average temperature of the radiation zone, the will be close to the .

The relationship between the value and the time after fiber fuse generation is shown in Figure 9. exceeds 10000 K immediately after fiber fuse generation but is less than 7000 K after 4 ms. Then, it gradually approaches about 5700 K. The value of 5760 K shown in Figure 9 is the average value from 4 to 30 ms. This temperature (5760 K) is close to the reported temperatures of 5400 K [5] and about 5800 K [56].

Figure 9: Average temperature of radiation zone versus time after fiber fuse generation.

Thus, it was found that if a fiber fuse, which exhibits a sharp temperature peak located near the light source, is approximated as a blackbody isolated from the surrounding nonheated regions, its average temperature from 4 ms after the generation of the fiber fuse approaches the experimentally estimated radiation temperatures.

4. Conclusion

We investigated the unsteady thermal conduction status in a single-mode optical fiber by numerical computation in order to visualize the mode of fiber fuse propagation. We assumed that the vitreous silica optical fiber underwent pyrolysis at elevated temperatures to form   (). We also proposed a model in which the optical absorption coefficient of the core layer increased with increasing molar concentration of . By using the model, we calculated the temperature distribution in the fiber with the explicit finite-difference method. It was found that when a short core with 40 m length was heated to 2923 K and a 2 W laser light (wavelength of 1.064 m) entered the core layer of an SMF-28 optical fiber, a thermal wave, that is, a fiber fuse, with a peak temperature of about 34000 K was generated at the boundary of the heating region near the light source. The fiber fuse was enlarged and propagated toward the light source at a rate of about 0.54 m/s. The calculated propagation velocity of the fiber fuse was in agreement with the experimental value. Moreover, the average temperature of the radiated region of the core layer was less than 7000 K at a time of 4 ms after the generation of the fiber fuse and gradually approached a temperature of about 5700 K. The final average temperature was close to the experimentally reported values.

Recently, Toratani et al. demonstrated femtosecond laser fabrication of submicrometer-sized voids in fused silica, and both reflector and optical waveguides were fabricated by femtosecond laser only [61]. In this fabrication technique, self-focusing plays an important role to generate the voids. This technique would be useful to develop 3-dimensional optical devices and fabricate precise sensors. Our study may contribute somewhat to analyzing or controlling such a fabrication process.

The cause of mechanism of the fiber fuse phenomenon has yet to be sufficiently clarified. It largely depends upon future multilateral studies.

Appendix

Focal Length of Thermal Lens Effect

If part of the core layer is heated by light absorption, the refractive index gradient in the core is induced by the thermal coefficient of the refractive index of the silica glass and the temperature distribution in the core layer. The refractive index of the heated core is expressed as follows:

where (=1.46) is the characteristic value of the core layer, is the radial distance from the center of the optical fiber, and is time, respectively.

If we assume that optical power in the optical fiber takes on Gaussian distribution, in the core layer is given by [35, 62]

where is the absorption coefficient, is the incident optical power, is the spot size radius of the laser beam, and is the thermal conductivity of the silica glass, respectively. Parameter , where and are the specific heat and density of the silica glass.

Substitution of (A.2) into (A.1) gives

When refractive index of the core layer takes a radial distribution as shown in (A.3), propagating laser beam will be focussed as a result of the thermal lens effect. In this case, the focal length is given by [35, 62]

where is the length of the heating core, where exhibits large value. In typical SM fiber,  W/mK,  J/kgK,  kg/m³ [5], and m. If we insert these values into (A.5), we obtain  s. This means that when is 10 s or above. If we assume  W,  m,   [5], and   [5], is given by

This value is of the same order of the observed interval and/or large front size of the cavities [17, 20].

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

References

A. Sano, T. Kobayashi, E. Yoshida, and Y. Miyamoto, “Ultra-high capacity optical transmission technologies for 100 Tbit/s optical transport networks,” IEICE Transactions on Communications, vol. 94, no. 2, pp. 400–408, 2011.

T. Morioka, “New generation optical infrastructure technologies: “EXAT initiative” towards 2020 and beyond,” in Proceedings of the 14th OptoElectronics and Communications Conference (OECC '09), July 2009.

R. Kashyap and K. J. Blow, “Observation of catastrophic self-propelled self-focusing in optical fibres,” Electronics Letters, vol. 24, no. 1, pp. 47–49, 1988.

R. Kashyap, “Self-propelled self-focusing damage in optical fibres,” in Proceedings of the 10th International Conference on Lasers, pp. 859–866, 1988.

D. P. Hand and P. S. J. Russell, “Solitary thermal shock waves and optical damage in optical fibers: the fiber fuse,” Optics Letters, vol. 13, no. 9, pp. 767–769, 1988.

D. P. Hand and P. S. J. Russell, “Soliton-like thermal shock-waves in optical fibres: origin of periodic damage tracks,” in Proceedings of the 14th European Conference on Optical Communication, pp. 111–114, September 1988.

R. Kashyap, “The fiber fuse—from a curious effect to a critical issue: a 25th year retrospective,” Optics Express, vol. 21, no. 5, pp. 6422–6441, 2013.

S. Todoroki, “Fiber fuse propagation behavior,” in Selected Topics on Optical Fiber Technology, Y. Moh, S. W. Harun, and H. Arof, Eds., pp. 551–570, InTech, Croatia, 2012.

P. André, A. Rocha, F. Domingues, and M. Facão, “Thermal effects in optical fibres,” in Developments in Heat Transfer, M. A. D. Bernardes, Ed., pp. 1–20, InTech, Croatia, 2011.

T. J. Driscoll, J. M. Calo, and N. M. Lawandy, “Explaining the optical fuse,” Optics Letters, vol. 16, pp. 1046–1048, 1991.

D. D. Davis Jr., S. C. Mettler, and D. J. DiGiovanni, “Experimental data on the fiber fuse,” in 27th Annual Boulder Damage Symposium: Laser-Induced Damage in Optical Materials, vol. 2714 of Proceedings of SPIE, pp. 202–210, October 1995.

D. D. Davis, S. C. Mettler, and D. J. DiGiovanni, “A comparative evaluation of fiber fuse models,” in Laser-Induced Damage in Optical Materials, vol. 2966 of Proceedings of SPIE, pp. 592–606, October 1996.

R. Kashyap, A. H. Sayles, and G. F. Cornwell, “Heat flow modeling and visualization of catastrophic self-propagating damage in single-mode optical fibers at low powers,” in Laser-Induced Damage in Optical Materials, vol. 2966 of Proceedings of SPIE, pp. 586–591, October 1996.

E. M. Dianov, I. A. Bufetov, A. A. Frolov et al., “Catastrophic destruction of optical fibres of various composition caused by laser radiation,” Quantum Electronics, vol. 32, no. 6, pp. 476–478, 2002.

E. M. Dianov, I. A. Bufetov, A. A. Frolov et al., “Catastrophic destruction of fluoride and chalcogenide optical fibres,” Electronics Letters, vol. 38, no. 15, pp. 783–784, 2002.

R. Kashyap, “High average power effects in optical fibres and devices,” in Reliability of Optical Fiber Components, Devices, Systems, and Networks, vol. 4940 of Proceedings of SPIE, pp. 108–117, October 2002.

R. M. Atkins, P. G. Simpkins, and A. D. Yablon, “Track of a fiber fuse: a Rayleigh instability in optical waveguides,” Optics Letters, vol. 28, no. 12, pp. 974–976, 2003.

Y. Shuto, S. Yanagi, S. Asakawa, M. Kobayashi, and R. Nagase, “Fiber fuse phenomenon in step-index single-mode optical fibers,” IEEE Journal of Quantum Electronics, vol. 40, no. 8, pp. 1113–1121, 2004.

S.-I. Todoroki, “In-situ observation of fiber-fuse propagation,” Japanese Journal of Applied Physics, vol. 44, no. 6, pp. 4022–4024, 2005.

S.-I. Todoroki, “Origin of periodic void formation during fiber fuse,” Optics Express, vol. 13, no. 17, pp. 6381–6389, 2005.

E. D. Bumarin and S. I. Yakovlenko, “Temperature distribution in the bright spot of the optical discharge in an optical fiber,” Laser Physics, vol. 16, no. 8, pp. 1235–1241, 2006.

L. A. Bufetov, A. A. Frolov, E. M. Dianov, V. E. Fortov, and V. P. Efremov, “Dynamics of fiber fuse propagation,” in Proceedings of the Optical Fiber Communication Conference (OFC/NFOEC '05), pp. 501–503, March 2005.

E. M. Dianov, V. E. Fortov, I. A. Bufetov et al., “High-speed photography, spectra, and temperature of optical discharge in silica-based fiberss,” IEEE Photonics Technology Letters, vol. 18, no. 6, pp. 752–754, 2006.

S. Todoroki, “Fiber fuse propagation modes in typical single-mode fibers,” in Proceedings of the Optical Fiber Communication Conference (OFC/NFOEC '13), 2013.

S. I. Yakovlenko, “Plasma behind the front of a damage wave and the mechanism of laser-induced production of a chain of caverns in an optical fibre,” Quantum Electronics, vol. 34, no. 8, pp. 765–770, 2004.

A. M. Rocha, P. F. D. C. Antunes, M. D. F. F. Domingues, M. Facão, and P. S. D. B. André, “Detection of fiber fuse effect using FBG sensors,” IEEE Sensors Journal, vol. 11, no. 6, pp. 1390–1394, 2011.

K. S. Abedin and M. Nakazawa, “Real time monitoring of a fiber fuse using an optical time-domain reflectometer,” Optics Express, vol. 18, no. 20, pp. 21315–21321, 2010.

D. P. Hand and T. A. Birks, “Single-mode tapers as “fibre fuse” damage circuit-breakers,” Electronics Letters, vol. 25, no. 1, pp. 33–34, 1989.

E. M. Dianov, I. A. Bufetov, and A. A. Frolov, “Destruction of silica fiber cladding by the fuse effect,” Optics Letters, vol. 29, no. 16, pp. 1852–1854, 2004.

S. Yanagi, S. Asakawa, M. Kobayashi, Y. Shuto, and R. Nagase, “Fiber fuse terminator,” in Proceedings of the Pacific Rim Conference on Lasers and Electro-Optics, vol. 1, p. 386, 2003.

K. Kurokawa and N. Hanzawa, “Suppression of fiber fuse propagation and its break in compact fiber fuse terminator,” in Proceedings of the OECC/PS, pp. 4–5, 2013.

S. I. Yakovlenko, “On reasons for strong absorption of light in an optical fibre at high temperature,” Quantum Electronics, vol. 34, no. 9, pp. 787–789, 2004.

M. Facão, A. M. Rocha, and P. S. De Brito André, “Traveling solutions of the fuse effect in optical fibers,” Journal of Lightwave Technology, vol. 29, no. 1, pp. 109–114, 2011.

E. M. Dianov, V. E. Fortov, I. A. Bufetov et al., “Detonation-like mode of the destruction of optical fibers under intense laser radiation,” JETP Letters, vol. 83, no. 2, pp. 75–78, 2006.

M. Katayama, Laser Chemistry; Nonlinear Spectroscopy and Laser Induced Chemical Processes, Syokabo, Tokyo, Japan, 1985.

A. E. Owen and R. W. Douglas, “The electrical properties of vitreous silica,” Journal of the Society of Glass Technology, vol. 43, pp. 159–178, 1959.

R. H. Doremus, “Ionic transport in amorphous oxides,” Journal of The Electrochemical Society, vol. 115, pp. 181–186, 1968.

H. L. Schick, “A thermodynamic analysis of the high-temperature vaporization properties of silica,” Chemical Reviews, vol. 60, no. 4, pp. 331–362, 1960.

H. R. Philipp, “Optical properties of non-crystalline Si, SiO, SiOx and SiO₂,” Journal of Physics and Chemistry of Solids, vol. 32, no. 8, pp. 1935–1945, 1971.

H. R. Philipp, “Optical and bonding model for non-crystalline SiOx and SiOxNy materials,” Journal of Non-Crystalline Solids, vol. 8–10, pp. 627–632, 1972.

W. J. Moore, Physical Chemistry, chapter 8, Prentice-Hall, 4th edition, 1972.

JANAF Thermochemical Tables, U.S. Department of Commerce and National Bureau of Standards, 2nd edition, 1971.

M. H. Brodsky, R. S. Title, K. Weiser, and G. D. Pettit, “Structural, optical, and electrical properties of amorphous silicon films,” Physical Review B, vol. 1, no. 6, pp. 2632–2641, 1970.

D. L. Dexter, “Interpretation of urbach's rule,” Physical Review Letters, vol. 19, no. 24, pp. 1383–1385, 1967.

H. Mahr, “Ultraviolet absorption of KI Diluted in KCl crystals,” Physical Review, vol. 125, no. 5, pp. 1510–1516, 1962.

G. G. Macfarlane and V. Roberts, “Infrared absorption of silicon near the lattice edge,” Physical Review, vol. 98, no. 6, pp. 1865–1866, 1955.

G. G. Macfarlane, T. P. McLean, J. E. Quarrington, and V. Roberts, “Fine structure in the absorption-edge spectrum of Si,” Physical Review, vol. 111, no. 5, pp. 1245–1254, 1958.

E. Burstein, G. Picus, and N. Sclar, “Optical and photoconductive properties of silicon and germanium,” in Proceedings of the Photoconductivity Conference, pp. 353–413, John Wiley & Sons, New York, NY, USA, 1956.

H. S. Carslaw and J. C. Jaeger, Conduction of Heat in Solids, chapter 13, Oxford University Press, Oxford, UK, 2nd edition, 1959.

I. Barin and O. Knacke, Thermochemical Properties of Inorganic Substances, Springer, Berlin, Germany, 1973.

M. Shoji, Heat Transfer Textbook, chapter 2, University Tokyo Press, Tokyo, Japan, 1995.

T. Okuda, Plasma Engineering, chapter 3, Corona-Sha, Tokyo, Japan, 1975.

T. Okuda, Plasma Engineering, chapter 1, Corona-Sha, Tokyo, Japan, 1975.

L. Spitzer Jr., Physics of Fully Ionized Gases, John Wiley & Sons, 2nd edition, 1962.

G. D. Smith, Numerical Solution of Partial Differential Equations: Finite Difference Methods, chapter 2, Clarendon Press, Oxford, UK, 3rd edition, 1985.

E. M. Dianov, V. E. Frotov, I. A. Bufetov, et al., “Temperature of optical discharge under action of laser radiation in silica-based fibers,” in Proceedings of the ECOC, 2005.

L. Spitzer, Physics of Fully Ionized Gases, chapter 5, John Wiley & Sons, 2nd edition, 1962.

Y. B. Zel’dovich and Y. P. Raizer, Physics of Shock Waves and High-Temperature Hydrodynamic Phenomena, chapter 5, Dover, New York, NY, USA, 2002.

T. Sekiguchi, Plasma Engineering, chapter 1, Ohm-Sha, Tokyo, Japan, 1997.

G. R. Fowles, Introduction to Modern Optics, chapter 7, Dover, New York, NY, USA, 2nd edition, 1975.

E. Toratani, M. Kamata, N. N. Nedyalkov, and M. Obara, “Femtosecond laser processing of subwavelength-sized voids for compact optical devices,” in Commercial and Biomedical Applications of Ultrafast Lasers, vol. 6108 of Proceedings of SPIE, pp. 1–8, January 2006.

J. P. Gordon, R. C. C. Leite, R. S. Moore, S. P. S. Porto, and J. R. Whinnery, “Long-transient effects in lasers with inserted liquid samples,” Journal of Applied Physics, vol. 36, no. 1, pp. 3–8, 1965.