Analysis of Temperature Dependence of Thermally Induced Transient Effect in Interferometric Fiber-optic Gyroscopes
- Author: Choi Woo-Seok
- Organization: Choi Woo-Seok
- Publish: Journal of the Optical Society of Korea Volume 15, Issue3, p237~243, 25 Sep 2011
Thermal characteristics, such as diffusivity and temperature induced change in the fiber mode index of rotation sensing fiber coil are critical factors which determine the time varying, thermo-optically induced bias drift of interferometric fiber-optic gyroscopes (IFOGs). In this study, temperature dependence of the transient effect is analyzed in terms of the thermal characteristics of the fiber coil at three different temperatures. By applying an analytic model to the measured bias in the experiments, comprehensive thermal factors of the fiber coil could be extracted effectively. The validity of the model was confirmed by the fact that the extracted values are reasonable results in comparison with well known properties of the materials of the fiber coil. Temperature induced changes in the critical factors were confirmed to be essential in compensating the transient effect over a wide temperature range.
Fiber-optic gyroscope , Shupe effect , Thermal characteristics of fiber coil , (060.2800) Gyroscopes , (060.2370) Fiber optic sensors , (120.5790) Sagnac effect , (120.6810) Thermal effect
The interferometric fiber-optic gyroscope (IFOG) is a representative fiber-optic sensor based on the Sagnac effect to measure absolute rotation rates in inertial reference frames [1-2]. IFOGs have many advantages over conventional gyroscopes, such as long life, high reliability, light weight,low cost, and so on. Many useful research results related with IFOGs have been achieved so that the IFOG is now a good prospective candidate for future inertial navigation systems. However, in spite of such positive advances, there are still some difficulties to be overcome to improve IFOG performance.
The thermally induced non-reciprocity of IFOGs has been one of main obstacles to applying them in practical applications.This phenomenon, the so called Shupe effect, yields transient bias drift when a temporarily changing thermal gradient passes through the rotation-sensing fiber coil .With the aid of invaluable earlier works, it could be proved that thermal transient non-reciprocity can be considerably reduced by applying a quadrupolar winding pattern to the fiber coil [4-5]. Nonetheless, this does not sufficiently address the problem because there can still be residual drift much larger than steady-state bias stability depending on the degree of imperfection in the thermal design and fabrication of the fiber coil. Quantitative analysis based on a thermal model considering heat conduction and physical stress is essential to identify proper methodologies which can address the imperfection [5-6]. Because thermal characteristics, such as diffusivity and temperature induced change in the effective index of the fiber coil, are critical factors determining the Shupe effect, it is necessary to examine the behavior of the thermal characteristics in detail to further improve the temporal thermal performance. However, all the established studies in analytic modeling of the Shupe effect have been performed with an assumption that thermal characteristics are invariant over temperature changes even though this is not true in the real world. For an accurate analysis of the thermal transient over a wide temperature range, it is indispensible to examine changes in thermal characteristics as a function of the average temperature measured during experiments.
In this study, temperature dependencies of the transient responses are analyzed focusing on changes in the thermal characteristics of the fiber coil at several fiber coil temperatures.The related theory with analytic model is presented in Section 2, experimental methods, analyzed results, and a summary are given in Sections 3, 4, and 5, respectively.
Ais the amplitude of the Shupe effect, ΔT&is the time derivative of ΔT, Land l are the length of the fiber and the length of the curvilinear abscissa along the fiber,respectively. Eq. 1 is obvious under the assumption that the thermally induced temporal phase error is linearly proportional to the accumulated amount of non-uniform temperature change difference along the whole fiber. The amplitude, Ahas the important role of describing phase imbalance between two counter-propagating optical beams along the Sagnac loop in terms of effective index and length variation of the fiber by thermal perturbations that can be expressed as
kis the wave vector in vacuum, nis the mode index of the fiber, c is the speed of light in vacuum, and α is the linear thermal expansion coefficient of the fiber, respectively.Generally, as α is much smaller than one order of magnitude, the α nterm is ignored in the following .The critical factor, dn/dT, is not easy to know accurately in this model even though there is well-known information for the silica fiber itself. This is because the fiber experiences more severe physical circumstance when it is utilized in the fabrication of a fiber coil. As seen in Fig. 1, the fiber coil is usually comprised of various materials. Long coated fiber needs to be wound in a metallic spool and to be encapsulated by an adhesive potting material to make it a suitable coil with a homogeneous solid state. When some transient thermal perturbations are applied to the fiber coil,secondary effects such as stress or strain can occur due to elastic deformation caused by complex thermal interaction between the constituting materials [8-9]. The effective optical path along the fiber is also susceptible to such physical effects. Therefore, in this study, dn/dTis not a fixed constant,but it is considered a critical factor that can change according to environmental conditions and effects resulting from the fabrication of the fiber coil.
For the thermal perturbation applied to the radial direction of the fiber coil with quadrupolar winding pattern,Eq. (1) can be rewritten as a discrete form:
where N is the total layer number,
Δ？(n)and lnare temperature change and the length in the nth fiber layer, respectively.To calculate the amount of accumulated transit phase error in Eq. (3), it is necessary to determine the ΔT&
for a fixed coil geometry. There have been two main methodologies to determine the Δ？ ; one is based on the finite element method (FEM) , and the other utilizes a numerical calculation with an analytic heat conduction equation . Even though the latter methodology may be less adaptive for correct estimation of
ΔT&, it is very advantageous to take analytic insight into the thermal characteristics in the process of solving the problem. In this study, a one-dimensional analytic heat conduction model based on the following equation is selected to catch the thermal diffusivities of the fiber coil directly :
Kis the thermal diffusivity of the material. Here, Kis a critical constant expressed by thermal-related properties as
k, ρ, and cpare thermal conductivity, density, and specific heat capacity, respectively. Because the coil is composed of various materials that can be largely divided by two different parts, namely, the metallic spool and arranged fiber encapsulated by potting material, it is necessary to simplify the calculation by treating them as a homogeneous body with a rescaled coordinate, □ for the nth layer:
lsand dcare the length of metallic spool and the average diameter of fiber composite per layer, and Ksand Kcare the thermal diffusivity of spool and fiber composite,respectively. For the given initial condition ΔT(x,0)= 0 and the boundary condition ΔT(0,t), the rescaled coordinate Δ？ (x,t )can be calculated numerically. If ΔT& (x,t )is determined,the rotation rate error due to the thermal transient perturbation can be also determined by applying the Eq.(1) to following equation :
where λ and D are the central wavelength of the light source and the average loop diameter of the fiber, respectively.
Note that the transient rate error due to the thermal perturbation is apparently dependent on
dn/dT, Ks, and Kcfor the given ΔT(x,t)and fixed coil geometry. This implies that even small changes in these three critical factors can influence the Shupe effect of IFOG. Also, as confirmed in, strictly speaking, the Shupe effect is not only a function of the temperature rate but also of temperature itself.As Eq. (1) seems to be solely related to the relative temperature gradient, it is necessary to consider the critical factors as variables on temperature. Therefore, for precise analysis of the Shupe effect, it is indispensable to include the variant effects of dn/dT, Ksand Kcin terms of temperature change. In this study, temperature dependence on these factors is investigated to identify the variant susceptibility of the Shupe effect on the temperature change.
The configuration of IFOG tested in this work is depicted in Fig. 2. The light from a broadband erbium-doped fiber source (EDFS) with a 1550 nm center wavelength is split by a 3-dB optical coupler and then coupled to a fiber which is pigtailed to a LiNbO3 integrated optic chip (IOC)which functions as a TM-mode-rejecting polarizer, a Y-branch 3-dB beam splitter, and an electro-optical phase modulator.A polarization maintaining fiber (PMF) with 80 ㎛ cladding and 155 ㎛ coating diameter, and 1.45 core index is wound in a 40 mm diameter spool for realization of the fiber coil.The total length of the PMF in the coil is about 900 m,and the fiber coil is composed of 60 fiber layers with a quadrupolar winding pattern. After passing through the fiber coil, two counter propagating light beams are interfered at the common input port of the Sagnac loop, and an interference signal is detected at the photo diode. The rate signal is obtained by a closed-loop approach which is based on the serrodyne modulation and demodulation technique. This ensures a linear and stable scale factor .
It is very important to check whether other non-reciprocal error sources in addition to the Shupe effect are present in the measured rate bias of IFOG. The presence of any additional non-reciprocal errors causes overestimation in the analysis because they usually cannot be distinguished from the Shupe effect induced errors. The broad spectrum of the utilized EDFS, about 20 nm, reduces not only the nonreciprocal phase error due to back reflections and backscattering in the Sagnac interferometer  but also the Kerr
effect . The IOC with a high polarization extinction ratio over 60 dB and birefractive PMF with an h-parameter of less than 10-5 utilized in the fiber coil can suppress the phase error due to the effect of polarization cross-coupling to a considerably lower level [14-15]. The Faraday bias error which is sensitive to magnetic fields is also eliminated by using a magnetic shield case with highly permeable metal. Before the experiment to analyze the Shupe effect, it was confirmed that this configuration has a very small phase error of less than 10-7 radian by testing the IFOG in a static environment without thermal transient perturbation.Because the peak level of the expected phase error due to transient thermal perturbation in this study is larger than an order of 10-6, the configuration can be considered suitable for analyzing the Shupe effect.
After the preliminary investigation, experiments with the Shupe effect were performed at several surrounding temperatures,namely, -30℃, 25℃, and 50℃, to identify temperature dependence on the critical factors,
dn/dT, Ks, and Kc. The IFOG was mounted in the temperature chamber,and the initial surrounding temperature condition was achieved as a state of thermal equilibrium. Then, an electrically induced local heat source was turned on so that thermal transient perturbation was applied asymmetrically to the fiber coil by exposing one of its sides to thermal flux coming from the heat source. The rate signal of the IFOG and temperature information were measured for 3 hours. The two thermocouples were located at the center of the fiber coil and outside of the fiber coil. In this work, thermal perturbation is required to have adequate temperature variation over a range of several degrees without affecting significant changes in the critical factors while confirming the Shupe effect in the measured results. The relatively long measurement duration of 3 hours is an important condition. The long-term stability of the rate signal including the steadystate are identified for the given thermal perturbation; such necessary information is very helpful in analyzing the effects more accurately.
For the given coil geometry and
ΔT(x,t), dn/dT, Ks, and Kcwere extracted by applying the analytic model described in Section 3 to the measured rate error of the IFOG during the experiment. The measured boundary conditions at the two thermocouples are shown in Fig. 3, and the measured transient rate error due to the perturbation is shown in Fig.4. The Shupe effect induced rate error could be obtained after eliminating the noise in the measured output by using the 100-second adjacent average and compensating the fixed bias by assuming that there is no Shupe effect induced error in the steady- state. The fitted values in Fig. 4 are analyzed results with the extracted three critical factors,and they seem to be almost the same as the measured results. This implies that the analysis was accurate, and the analytic model proposed in this study can be applied in compensation of the transient rate error. In Fig. 5, the compensated outputs using the fitted results are shown.
Actually, there may be some overestimation in identifying
the correct values for the thermal characteristics. This is because the experiments and analytic model in this study are not perfect to characterize correct values for
dn/dT, Ks,and Kc. For example, non-ideal quadrupolar winding in the fabrication of the fiber coil can cause overestimation in analyzing the dn/dT. In addition, the one-dimensional simple analytic model cannot treat various thermal features of the non-uniform radiant heat flux utilized during the experiment which are also important variables determining the Shupe effect. However, despite such imperfections of the experiments and the analytic model, they are still acceptable in that the temperature-induced changes in the Shupe effect can be analyzed effectively with the extracted results.
All the extracted critical factors from each experiment are summarized in Figs. 6 and 7. It is possible to confirm two noteworthy tendencies in Figs. 6 and 7. One is that
Ksand Kcdecreased, and another is that dn/dTincreased with increase in the surrounding temperature. Considering that all the constituting materials of the fiber coil presented in Fig. 1 (b) have almost linearly decreased features in thermal diffusivity, while increased features in thermal expansion
coefficient which can enhance the stress induced index change as the temperature of the fiber coil is increased , these results seem reasonable. Although it is difficult to estimate the correct values for the critical factors by using the thermal characteristics of the constituting material themselves, this fact is sufficient to conclude that the comprehensive thermal characteristics of the fiber coil have the same tendencies as the constituting materials on the temperature changes.The fractional changes in the extracted values can be evaluated effectively by assuming that they comply with linear temperature increase as seen in Figs. 6 and 7. By applying the 1st-order least-square fitting technique to the extracted values, the constant change rates of
Ks, Kc, and dn/dTon temperature increase were analyzed to be -41 ppm?cm2/s/℃,-0.11 ppm?cm2/s/℃, and 53 ppm/℃2, respectively.
The analyzed simple linear equations from Figs. 6 and 7 are very advantageous when they are utilized in the estimation of the Shupe effect based on the analytic model because they can generate reasonable values on the critical factors at arbitrary temperatures within the considered range.The calculated Shupe effect induced biases for five different temperatures utilizing the extracted linear equations and the analytic model are presented in Fig. 8. The boundary conditions which are necessary for the calculation of the Shupe effect induced biases are assumed to be same as that measured at 25℃. The calculation results in Fig. 8 show that the peak-to-peak quantity and the zero-crossing time of the Shupe effect induced bias increase with increase in the surrounding temperature. These calculation results have a tendency similar to that of the measured results seen in Fig. 4 in the aspect of the broadening of the transience.However, the variant features of the peak-to-peak quantity seen in Fig. 8 do not seem similar to that of Fig.4. The peak-to-peak quantity of the Shupe effect induced bias measured at -30℃ is obviously larger than that of the results measured at higher temperatures of 25℃ and 50℃.These differences between Figs. 4 and 8 are due to the different boundary conditions in the process of estimating
the Shupe effect induced error. Considering that the nominal values of the boundary condition measured at -30℃ are larger than those of the boundary condition measured at 25℃ and 50℃, as seen in Fig. 3, the variant features in the peak-to-peak quantity of the Shupe effect induced bias in Fig. 4 are essentially identical to the features of Fig. 8.
From the results of Figs. 4 and 8, it can be easily confirmed that the variation of the Shupe effect induced error per unit of temperature change is very small relative to the intrinsic outputs in steady-state, which drift and swing within the bounded range of ±0.1 deg/hr. Because the Shupe effect is analyzed based on the measured outputs of IFOG,any measurements for numerous surrounding temperature that do not exhibit distinct variations in the Shupe effect induced error more than the intrinsic bias uncertainty of IFOG can be supernumerary. Because the transient-state is usually inclined to have poorer performance than the steady-state section even when the Shupe effect induced error is not considered, for an accurate analysis, it is more important to achieve good bias repeatability in the transient-state. In this study, it was confirmed by taking repeated measurements at the predetermined surrounding temperatures that bias repeatability in the transient-state is as good as that in the steady-state section. The confirmation of bias performance and the corresponding three temperature points for measurement and analysis in this study are the minimum suitable conditions for analyzing the Shupe effect with the extracted variant equations over the considered range of temperatures.For 25℃, the calculated Shupe effect induced error with the variant equations is presented together with the fitted results based on the analytic model and the measured data in Fig. 9. It can be easily found out that there are negligible differences of less than 0.03 deg/hr between the fitted and calculated Shupe effect induced bias error.
In Table 1, the results of Shupe effect compensation using the linear equations are summarized. The compensation is carried out with only the transient data up to initial 1
hour of the measured results for 3 hours and the presented values in Table 1 are worst cases of the compensated results. Even though the compensation is not as precise as the results with the extracted values seen in Fig 5, this compensation still considerably improves the bias performance of the IFOG over a wide range of temperatures.
Temperature dependence on the Shupe effect of IFOGs was analyzed in terms of the thermal characteristics of the fiber coil. It was confirmed that thermal characteristics, such as diffusivity and temperature-induced change in the fiber mode index of the fiber coil, were critical factors that determined the Shupe effect of IFOGs. By applying an analytic model with 1-dimensional heat conduction and the geometry of the fiber coil to the measured transient responses of an IFOG, the critical factors were accurately identified. Temperature-induced changes in the critical factors were also analyzed by experiments at three different temperatures,and these factors exhibited variant tendencies to temperature changes. They were confirmed to be essential in compensating the transient effect over a wide range of temperatures.
[FIG. 1.] Schematic view of (a) fiber coil and (b) its1-dimensional cross-section with quadrupolar windingpattern (one half coil length is represented by + and the otherby -).
[FIG. 2.] Schematic view of IFOG configuration.
[FIG. 3.] Measured boundary conditions.
[FIG. 4.] Comparison of the measured and fitted rate errors dueto the Shupe effect with the radial boundary conditions seenin Fig. 3.
[FIG. 5.] Compensated output obtained by subtracting thefitted values from the measured results presented in Fig. 4.
[FIG. 6.] Extracted thermal diffusivities of spool and fibercomposite from the measured results.
[FIG. 7.] Extracted temperature induced changes of fiber modeindex from the measured results.
[FIG. 8.] Calculated rate error with the three linear equations ofFigs. 6 and 7 and the measured boundary condition at 25℃.
[FIG. 9.] Comparison of the measured fitted and calculatedrate error at 25℃.
[TABLE 1.] Shupe effect compensation results using linearequations presented in Fig. 6 and 7