Next Article in Journal
Precise Aperture-Dependent Motion Compensation with Frequency Domain Fast Back-Projection Algorithm
Next Article in Special Issue
Highly Sensitive Strain Sensor Based on a Novel Mach-Zehnder Interferometer with TCF-PCF Structure
Previous Article in Journal
Waveguide Bragg Gratings in Ormocer®s for Temperature Sensing
Previous Article in Special Issue
Gyroscope Technology and Applications: A Review in the Industrial Perspective
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling of Thermal Phase Noise in a Solid Core Photonic Crystal Fiber-Optic Gyroscope

Department of Opto-electronics Engineering, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(11), 2456; https://doi.org/10.3390/s17112456
Submission received: 5 September 2017 / Revised: 10 October 2017 / Accepted: 23 October 2017 / Published: 26 October 2017
(This article belongs to the Special Issue Integrated Photonic Technologies for Sensing Applications)

Abstract

:
A theoretical model of the thermal phase noise in a square-wave modulated solid core photonic crystal fiber-optic gyroscope has been established, and then verified by measurements. The results demonstrate a good agreement between theory and experiment. The contribution of the thermal phase noise to the random walk coefficient of the gyroscope is derived. A fiber coil with 2.8 km length is used in the experimental solid core photonic crystal fiber-optic gyroscope, showing a random walk coefficient of 9.25 × 10−5 deg / h .

1. Introduction

Interferometric fiber-optic gyroscopes (IFOGs) measure the angular rate around a fixed axis based on the Sagnac effect, and have been increasingly used due to their particular properties such as high performance, low power consumption and high reliability as a solid state structure with no moving parts [1]. The typical implementation example is the LN-200 in ‘Clementine’ in 1994, which was the first IFOG utilization in aerospace navigation. With the development of inertial guidance and navigation systems, the accuracy of IFOGs is improving by using different signal processing and noise suppression methods, such as square-wave modulation, relative intensity noise (RIN) subtraction and modulation phase optimization [2,3]. According to a recent report, the MARINS series IFOGs of iXBlue.co have even reached a bias stability of 2 × 10−5 degrees per hour [4]. When the gyros are designed for high-sensitivity applications, the thermal phase noise seems to be a fundamental limitation. The thermal phase noise is usually caused by the thermal fluctuations in the refractive index of the fiber [5,6] and was first observed in a conventional fiber-optic gyro in [7]. Additionally, for an open-loop dynamically sinusoidal biased IFOG, it was demonstrated theoretically [8] and then experimentally verified by Moeller and Burns in 1996 [9]. However, in the abovementioned studies, they adopted single mode fibers (SMFs) and polarization maintaining fibers (PMFs), respectively, and there have been no prior studies regarding the thermal phase noise in solid core photonic crystal fiber-optic gyroscopes (SC-PCFOGs). In fact, during long-term exposure in a space environment, the components of the space used IFOG are generally affected by radiation, leading to the performance degradation of the system. Compared with the conventional PMFs, solid core photonic crystal fibers (SC-PCFs) has much lower sensitivity to temperature and radiation in the same conditions [10,11], which shows great prospect in space applications and SC-PCFOGs. As the SC-PCFs are different from conventional PMFs in structure and material, and the thermal phase noise can be the limitation of high-precision SC-PCFOGs, it is necessary to analyze and measure the thermal phase noise in SC-PCFOGs.
In this paper, a theoretical model of thermal phase noise in a square-wave modulated SC-PCFOG was derived, and a thermal phase noise measurement setup was built to verify the model. During the measurement, a RIN subtraction method was adopted and the noise in the output signal decreased by 18 dB with this method. Additionally, the measured thermal phase noise voltage coincided with the theoretical thermal phase noise calculated from the derived model. Furthermore, contribution from thermal phase noise to the random walk coefficient of the gyroscope was calculated and showed a value of 9.25 × 10 5 deg / h using our SC-PCFOG parameters.

2. Theory of Thermal Phase Noise

As illustrated in Figure 1, a SC-PCFOG operating in a typical reciprocal configuration consists of an amplified spontaneous emission (ASE) source with high power and broad spectrum, a coupler, a multifunction integrated optical circuit device (MIOC) and a SC-PCF coil.
For a square-wave phase modulation applied on the MIOC, the interference output of the IFOG becomes [1]:
P d ( t ) = 1 2 P 0 { 1 + cos [ Δ ϕ R ( t ) + Δ ϕ N ( t ) + Δ ϕ B ( t ) ] }
where 1 / 2 P 0 is the input optical power, Pd is the optical power at the detector, Δ ϕ R is the rotation-induced phase difference, Δ ϕ B is the square-wave modulation bias phase and Δ ϕ N is the phase perturbation owing to noise. If we further assume that the gyro remains at rest, Δ ϕ R 0 ° , Equation (1) can be expressed as:
P d ( t ) 1 2 P 0 [ 1 + cos Δ ϕ B ( t ) Δ ϕ N ( t ) sin Δ ϕ B ( t ) ]
The MIOC is driven by an oscillator to generate a square-wave modulation Δ ϕ B = ± ϕ 0 , where ϕ 0 is the amplitude of the square-wave. We can expand Δ ϕ B by Fourier series at a frequency ω m = 2 π f m , f m is the modulation frequency:
Δ ϕ B ( t ) = n = 1 [ 4 ϕ 0 n π sin 2 ( n π 2 ) sin ( n ω m t ) ] = 4 ϕ 0 π sin ( ω m t ) + 4 ϕ 0 3 π sin ( 3 ω m t ) + + 4 ϕ 0 ( 2 n - 1 ) π sin [ ( 2 n - 1 ) ω m t ] , n = 1 , 2 , 3 . .
Only the first three orders of Equation (3) were retained and substituted into Equation (2). Then through Jn Bessel function expansion, we obtain:
I d ( t ) = 1 2 I 0 1 + [ J 0 ( A ) + 2 n = 1 J 2 n ( A ) cos ( 2 n ω m t ) ] × [ J 0 ( A 3 ) + 2 n = 1 J 2 n ( A 3 ) cos ( 6 n ω m t ) ] [ 2 n = 0 J 2 n + 1 ( A ) sin { ( 2 n + 1 ) ω m t } ] × [ 2 n = 0 J 2 n + 1 ( A 3 ) sin { ( 6 n + 3 ) ω m t } ] Δ ϕ N ( t ) × [ 2 n = 0 J 2 n + 1 ( A ) sin { ( 2 n + 1 ) ω m t } × { J 0 ( A 3 ) + 2 n = 1 J 2 n ( A 3 ) cos ( 6 n ω m t ) } + { J 0 ( A ) + 2 n = 1 J 2 n ( A ) cos ( 2 n ω m t ) } × 2 n = 0 J 2 n + 1 ( A 3 ) sin { ( 6 n + 3 ) ω m t } ] ,
where the modulation depth A = 4 ϕ 0 / π . I0 and Id are the input and output optical intensities and can be expressed as I = ηP, where η is the detector responsivity. By performing a Fourier transform to Equation (4), and only considering the noise terms, the noise spectrum is given by:
| I N ( ω ) | = 1 2 I 0 4 π B S
S = J 0 2 ( A / 3 ) n = 0 J 2 n + 1 2 ( A ) { Δ ϕ N , r m s 2 [ ω + ( 2 n + 1 ) ω m ] + Δ ϕ N , r m s 2 [ ω ( 2 n + 1 ) ω m ] } + J 0 2 ( A ) n = 0 J 2 n + 1 2 ( A / 3 ) { Δ ϕ N , r m s 2 [ ω + ( 6 n + 3 ) ω m ] + Δ ϕ N , r m s 2 [ ω ( 6 n + 3 ) ω m ] } + m = 1 n = 1 J 2 m 1 2 ( A ) J 2 n 2 ( A / 3 ) i = 1 2 j = 1 2 { Δ ϕ N , r m s 2 [ ω + ( 1 ) i ( 2 m 1 ) ω m + ( 1 ) j ( 6 n ) ω m ] } + m = 1 n = 1 J 2 m 2 ( A ) J 2 n 1 2 ( A / 3 ) i = 1 2 j = 1 2 { Δ ϕ N , r m s 2 [ ω + ( 1 ) i ( 2 m ) ω m + ( 1 ) j ( 6 n 3 ) ω m ] }
Δ ϕ N , r m s 2 ( ω ) k B T 2 L κ λ 2 ( d n e f f d T + n e f f α L ) 2 × ln { [ ( 2 W 0 ) 2 + ( ω v ) 2 ] 2 + ( ω D i ) 2 [ ( 4.81 d ) 2 + ( ω v ) 2 ] 2 + ( ω D i ) 2 } × [ 1 sin c ( ω L v ) ] ,
Equation (5) represents the single-sided thermal phase noise in the IFOG when a square-wave bias signal is applied, where B is the electrical bandwidth. Equation (7) is phase noise spectral density caused by temperature induced index fluctuations in the fiber of a Sagnac interferometer as [7], where the length of fiber coil is defined by L, ν is the effective speed of light in the fiber, λ is the operational wavelength, W0 is the mode field radius and d is the cladding diameter of fiber. αL is the linear thermal expansion coefficient, Di is the thermal diffusivity, κ is the thermal conductivity, kB is the Boltzmann constant, T is the temperature and dneff/dT is the temperature coefficient of the fiber effective refractive index in the fiber. By multiplying Equation (5) by the detector responsivity η and relating I0 to the optical power P0, we can get the thermal phase noise at the detector as:
i N ( ω ) = 1 2 η P 0 4 π B S
where P0 is the output optical power without modulation. For small rotation rates, the signal current is is given by the coefficient of the first harmonic terms of Equation (4), as in the usual noise calculation in the appendix of [12]:
i S = Δ ϕ R J 1 ( A ) J 0 ( A / 3 ) I 0 2
By setting the signal current at the fundamental frequency component ωm equal to the thermal phase noise current and introducing the Sagnac phase shift, we can obtain the frequency-dependent minimum detectable rotation rate as:
Ω r m s , t h e r m a l ( ω ) = λ c 2 π L D 1 J 1 ( A ) J 0 ( A / 3 ) 2 π B S
where D is the diameter of fiber coil and c is the speed of light.

3. Experimental Setup

The experimental set-up for measurement of thermal phase noise in the SC-PCFOG is shown schematically in Figure 2. The corresponding RIN subtraction method was adopted to produce a gyro signal with reduced RIN. Light from an ASE with a maximum output power of 4.8 mW was focused into the end of the input coupler, which delivered 5.85 μW at the gyro detector when we took into account the splice and optical components insertion losses. A MIOC was added on the unused lead of the 50/50 coupler in the reference arm to ensure that its transmission axis was parallel to the transmission of gyro arm. For maximizing the signal level, the square wave modulation phase was 1.37 rad at the proper frequency of the fiber coil. And the MIOC was operated in a push-pull mode, which led to a net phase bias of 7π/8. As the principle of RIN subtraction is to extract the reference signal from the unused lead of the couple, and to subtract it from gyro signal after latency through fiber, the reference signal was then delayed through a PMF coil equivalent to the length of the SC-PCF coil. An adjustable attenuator was then spliced to the output lead of the coupler. Both signals were detected by photo-detectors and the output signals were then amplified separately and subtracted in a differential amplifier. The output noise was measured by a spectrum analyzer.
All the parameters of experimental IFOG are listed in Table 1. As some fiber-optic parameters vary with fiber composition and geometry, the measured fiber structure parameters are also given in Table 1. The SC-PCF we used was fabricated in collaboration with Fiber Home Telecommunication Technologies CO., Ltd. The fiber cross section and the scanning electron microscope (SEM) photograph of the PCF are shown in Figure 3. Considering the applicability to the SC-PCF, we adopted the commercial software COMSOL Multiphysics 4.4 based on full vector finite element method to generate the numerical simulation of thermal conductivity κ, temperature coefficient of refractive index d n e f f / d T and linear thermal expansion coefficient α L . The emulation values as well as typical parameters used in Equations (5)–(7) are presented in Table 1.
The key of RIN subtraction is to ensure that the RIN of the reference signal and the gyro signal are coherent. Hence, prior to subtraction, the optical path difference (OPD) between gyro signal and reference signal was adjusted first to ensure the efficiency of RIN subtraction. The delay time between both signals can be obtained by estimating the position of the maximum correlation coefficient of the correlation array. In our experiment, when a 295.28 m of typical PM fiber was added to the reference arm, the OPD between gyro signal and reference signal was less than 8.28 cm. During the measurements, the optical power arriving at the detector was held constant by adjusting the attenuator, so that the reference signal and gyro signal had the same intensity but different optical paths. In Figure 4 the output signals are presented with and without noise subtraction respectively. Experimentally, we measured 18 dB reduction by subtracting the two outputs when the SC-PCFOG was without modulation, which quantified that the RIN subtraction scheme is available for our measurement.

4. Results and Discussion

The measured noise and the theoretical thermal phase noise calculated from our model are presented in Figure 5. Due to the Earth’s rotation, the measured signal contains a peak at the modulation frequency (36.337 kHz). Despite this, the agreement between the calculated and experimental results is quite good with only 1–2 dB difference. This difference could be induced by the error of parameters and the experimental uncertainties.
The thermal phase noise, shot noise, RIN and detection thermal noise are calculated and then plotted in Figure 6. And the theoretical calculation formulas of above mentioned noises can be expressed as [13,14]:
σ s h o t 2 = 2 e I p B
σ R I N 2 = I p 2 B / Δ v
σ T 2 = 4 k T B / R
where σ s h o t 2 , σ R I N 2 and σ T 2 represent shot noise power, RIN power and detection thermal noise power, respectively. Both the relative intensity noise and shot noise are calculated with the measured detector average current I p of 5.56 μA, while the bandwidth of the optical source Δ v is 12 nm. And other parameters for the calculation are listed in Table 1.
From Figure 6 we see that, for our configuration, the level of thermal phase noise is higher than shot noise and detection thermal noise, but lower than RIN at the modulation frequency. Clearly, just off the modulation frequency the wings of thermal phase noise are quite large. Based on the theoretical thermal phase noise model of Equation (7), the relationship of phase noise spectral density and frequency is represented in Figure 7. The phase noise spectral density at 1 Hz is 2.8 × 10−9 rad / Hz , and the contribution from thermal phase noise in SC-PCFOG to random walk coefficient (RWC) can be obtained by Ω r m s / B . According to Equation (10), the calculated result is 9.25 × 10−5 deg / h which indicates the detection limitation of SC-PCFOG.

5. Conclusions

A theoretical thermal phase noise model in a square-wave modulated SC-PCFOG was derived. The thermal phase noise in an experimental SC-PCFOG was measured, showing a good agreement between theory and experiment. The contribution of thermal phase noise to the RWC was also calculated, giving an estimate of the detection limit of SC-PCFOG as 9.25 × 10−5 deg / h . Compared with the RIN, the thermal phase noise is not the dominant noise source in SC-PCFOG, but it is of the same order of magnitude as the shot noise. For our experimental SC-PCFOG, RIN subtraction was necessary to measure the thermal phase noise and the results confirm that the thermal phase noise can be a fundamental limitation in the high-sensitivity SC-PCFOG when a RIN subtraction method is employed.

Acknowledgments

This work was supported by National Natural Science Foundation of China (NSFC) under grants No. 61575012 and 61575013, and National Key Scientific Instrument and Equipment Development Project of China under grant No. 2013YQ040877.

Author Contributions

Ningfang Song and Jing Jin promoted the experimental methods and setup; Kun Ma performed the experiments and wrote the paper; Kun Ma and Fei Teng analyzed the data; Wei Cai simulated the parameters of fiber used in the experiments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lefevre, H.C. The Fiber-Optic Gyroscope, 2nd ed.; Artech House: London, UK, 2014. [Google Scholar]
  2. Deppe, O.; Dorner, G.; König, S.; Martin, T.; Voigt, S.; Zimmermann, S. MEMS and FOG Technologies for Tactical and Navigation Grade Inertial Sensors-Recent Improvements and Comparison. Sensors 2017, 17, 567. [Google Scholar] [CrossRef] [PubMed]
  3. Pan, X.; Liu, P.; Zhang, S.; Jin, J.; Song, N. Novel Compensation Scheme for the Modulation Gain to Suppress the Quantization-Induced Bias in a Fiber Optic Gyroscope. Sensors 2017, 17, 823. [Google Scholar] [CrossRef] [PubMed]
  4. Paturel, Y.; Honthaas, J.; Lefèvre, H.; Napolitano, F. One nautical mile per month fog-based strapdown inertial navigation system: A dream already within reach? Gyroscopy Navig. 2014, 5, 1–8. [Google Scholar] [CrossRef]
  5. Krakenes, K.; Bløtekjær, K. Comparison of fiber-optic Sagnac and Mach-Zehnder interferometers with respect to thermal processes in the fiber. J. Lightwave Technol. 1995, 13, 682–686. [Google Scholar] [CrossRef]
  6. Bartolo, R.E.; Tveten, A.B.; Dandridge, A. Thermal phase noise measurements in optical fiber interferometers. IEEE J. Quantum Electron. 2012, 48, 720–727. [Google Scholar] [CrossRef]
  7. Knudsen, S.; Tveten, A.B.; Dandridge, A. Measurements of fundamental thermal induced phase fluctuations in the fiber of a Sagnac interferometer. IEEE Photonics Technol. Lett. 1995, 7, 90–92. [Google Scholar] [CrossRef]
  8. Knudsen, S.; Blotekjaer, K. Interferometric fiber-optic gyroscope performance owing to temperature-induced index fluctuations in the fiber: Effect on bias modulation. Opt. Lett. 1995, 20, 1432–1434. [Google Scholar] [CrossRef] [PubMed]
  9. Moeller, R.P.; Burns, W.K. Observation of thermal noise in a dynamically biased fiber-optic gyro. Opt. Lett. 1996, 21, 171–173. [Google Scholar] [CrossRef] [PubMed]
  10. Ma, P.; Song, N.; Jin, J.; Song, J.; Xu, X. Birefringence sensitivity to temperature of polarization maintaining photonic crystal fibers. Opt. Laser Technol. 2012, 44, 1829–1833. [Google Scholar] [CrossRef]
  11. Kim, D.H.; Kang, J.U. Analysis of temperature-dependent birefringence of a polarization-maintaining photonic crystal fiber. Opt. Eng. 2007, 46, 075003. [Google Scholar] [CrossRef]
  12. Moeller, R.P.; Burns, W.K.; Frigo, N.J. Open-Loop Output and Scale Factor Stability in a Fiber-Optic Gyroscope. J. Lightwave Technol. 1989, 7, 262–269. [Google Scholar] [CrossRef]
  13. Renato, C.R.; Ricardo, T.C.; James, B. SNR enhancement of intensity noise-limited FOGs. IEEE J. Lightwave Technol. 2000, 18, 2146–2150. [Google Scholar]
  14. Iwatsuki, K. Excess noise reduction in fiber gyroscope using broader spectrum linewidth Er-doped superfluorescent fiber laser. IEEE Photonics Technol. Lett. 1991, 3, 281–283. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of a solid core photonic crystal fiber-optic gyroscope.
Figure 1. Schematic diagram of a solid core photonic crystal fiber-optic gyroscope.
Sensors 17 02456 g001
Figure 2. Experimental set-up for the measurement of the thermal phase noise in the SC-PCFOG.
Figure 2. Experimental set-up for the measurement of the thermal phase noise in the SC-PCFOG.
Sensors 17 02456 g002
Figure 3. Solid core photonic crystal fiber used in the experiment (a) Drawing of the cross section (b) Scanning electron micrograph
Figure 3. Solid core photonic crystal fiber used in the experiment (a) Drawing of the cross section (b) Scanning electron micrograph
Sensors 17 02456 g003
Figure 4. The output signals without phase modulation from spectrum analyzer: without noise subtraction (curve a), with noise subtraction (curve b).
Figure 4. The output signals without phase modulation from spectrum analyzer: without noise subtraction (curve a), with noise subtraction (curve b).
Sensors 17 02456 g004
Figure 5. Experimental thermal phase noise voltages obtained with the spectrum analyzer with 1.37 rad square-wave modulation and the corresponding calculated result.
Figure 5. Experimental thermal phase noise voltages obtained with the spectrum analyzer with 1.37 rad square-wave modulation and the corresponding calculated result.
Sensors 17 02456 g005
Figure 6. Calculated noise voltages for the solid core photonic crystal fiber-optic gyroscope.
Figure 6. Calculated noise voltages for the solid core photonic crystal fiber-optic gyroscope.
Sensors 17 02456 g006
Figure 7. The phase noise spectral density as a function of frequency.
Figure 7. The phase noise spectral density as a function of frequency.
Sensors 17 02456 g007
Table 1. Parameters of the experimental IFOG used in the computations.
Table 1. Parameters of the experimental IFOG used in the computations.
ParameterPropertyValue
η (A/W)Detector responsivity0.95
R (kΩ)Load resistance288
λ (nm)Operational wavelength1550
B (Hz)Reference bandwidth30
T (K)Temperature293.15
P0 (μW)Average optical power5.85
L (km)Length of fiber coil2.8
ν (m/s)Effective speed of light in the fiber2.079 × 108
ϕ 0 (V)Modulation phase1.4
A (rad)Modulation depth1.8
d (μm)Fiber cladding diameter100
κ (W/(m K))Thermal conductivity1.02
Di (m2/s)Thermal diffusity0.82 × 10−6
αl (ppm/°C)Linear thermal expansion coefficient1.02
dneff/dT (/°C)Temperature coefficient of refractive index9.9 × 10−6
D (cm)Diameter of fiber coil16
neffEffective refractive index1.435
W0Mode field radius3.1
Δd (μm)Diameter of the air holes in the cladding3.4
ΔD (μm)Diameter of the two enlarged air holes in x-direction5.8
Λ (μm)Distance of two adjacent air hole centers5.9

Share and Cite

MDPI and ACS Style

Song, N.; Ma, K.; Jin, J.; Teng, F.; Cai, W. Modeling of Thermal Phase Noise in a Solid Core Photonic Crystal Fiber-Optic Gyroscope. Sensors 2017, 17, 2456. https://doi.org/10.3390/s17112456

AMA Style

Song N, Ma K, Jin J, Teng F, Cai W. Modeling of Thermal Phase Noise in a Solid Core Photonic Crystal Fiber-Optic Gyroscope. Sensors. 2017; 17(11):2456. https://doi.org/10.3390/s17112456

Chicago/Turabian Style

Song, Ningfang, Kun Ma, Jing Jin, Fei Teng, and Wei Cai. 2017. "Modeling of Thermal Phase Noise in a Solid Core Photonic Crystal Fiber-Optic Gyroscope" Sensors 17, no. 11: 2456. https://doi.org/10.3390/s17112456

APA Style

Song, N., Ma, K., Jin, J., Teng, F., & Cai, W. (2017). Modeling of Thermal Phase Noise in a Solid Core Photonic Crystal Fiber-Optic Gyroscope. Sensors, 17(11), 2456. https://doi.org/10.3390/s17112456

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop