Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Hybrid method for monitoring large Fabry-Pérot cavity displacements with nanometer precision

Open Access Open Access

Abstract

The change in length of an optical fiber-based Fabry-Pérot cavity (FPC) can be precisely measured using phase tracking, but the displacement range is limited by phase ambiguity. Period tracking techniques determine the absolute FPC length, but with larger uncertainties from tracking the spacing between multiple peaks. A hybrid method is demonstrated that identifies appropriate peaks for phase tracking using a coarse estimate obtained from the free spectral range to effectively maintain the high precision (∼1 nm) of phase tracking techniques to measure ∼24 µm displacements, well beyond the range limitations (typically <1 µm) of phase tracking methods.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

Optical fiber-based sensors are routinely used for extremely high-precision measurements of various parameters, including length, displacement [1], and various engineering process parameters that can be derived from changes in length, such as temperature [2], pressure [3], fluid flow rate [4,5], liquid level [6,7], and vibration [8]. Their immunity to electromagnetic interference, small footprint, remote interrogation, and ability to survive harsh environmental conditions also make them candidates for process monitoring in industrial applications, such as chemical processing or energy production applications [9,10]. For example, fused silica [1113] and single-crystal sapphire [14,15] optical fibers transmit light at temperatures approaching 1,000 and 1,400$^{\circ }$C, respectively. Their ability to survive relatively high doses of neutron [1621] and gamma [2224] radiation also makes them attractive candidates for nuclear reactor applications.

Recent work has focused on leveraging the high-precision displacement monitoring capabilities of fiber-optic sensors to develop a sensor to measure corrosion in nuclear reactors and in experiments performed in research reactors [25]. The sensor is fabricated using ultrasonic additive manufacturing (UAM) [26,27] to directly embed an optical fiber within a metal matrix [2831]. The surrounding metal fixes the terminal end of the fiber within a metal housing. The metal housing terminates with a thin metal diaphragm that deflects when subjected to a differential pressure [32]. In this way, a Fabry-Pérot cavity (FPC) is established between the terminal end of the embedded optical fiber and at the inner surface of the metal diaphragm. This approach protects the optical fiber from the surrounding corrosive fluid and enables the sensor to be actuated via internal pneumatic pressurization. Corrosion can be monitored by calculating the change in the diaphragm’s thickness based on the analytical expression that relates the measured diaphragm displacement to the applied internal pressure [25]. With an appropriate interrogation technique, the same sensor could also measure higher frequency dynamic fluid interactions such as flow-induced turbulence or vibration.

The ability to accurately resolve micrometer-scale changes in thickness requires measuring diaphragm displacements on the order of micrometers to tens of micrometers with an accuracy on the order of tens of nanometers. Prior works [33,34] have applied spectral transforms (fast Fourier transforms and cepstral analyses), up-sampling, and curve fitting to the measured interference spectra in an effort to achieve nanometer-scale precision. These spectral techniques come at the cost of increased processing complexity. Wavelet transform-based techniques have also been utilized to demodulate change in FPC length. For example, a Fourier transform was used to generate a plausible range of FPC deflections, which was searched using wavelet transforms to identify the FPC deflection with greater precision [35,36]. Similar to other spectral transforms, the complexity of the wavelet transform analysis makes it challenging to implement and quickly process data, particularly for high frequency applications.

More recently, a custom-designed low-coherence interferometer (LCI) was developed and used novel phase-demodulation scheme to resolve changes in FPC length [32,37,38]. Although this system achieved the required resolution, the range of displacements that could be measured were limited to $\pm$ $\sim$15 $\mu$m based on the coherence length of the light source. Moreover, because this system measured total displacement by tracking the cumulative sum of incremental changes in FPC length (i.e., phase tracking), it was susceptible to errors that could be introduced during rapid changes in the FPC length caused by changes in temperature or differential pressure, which could propagate through future measurements [39]. LCI systems can also suffer from signal attenuation, even with compensation techniques [40,41].

This motivated interest in developing a simple method for monitoring larger FPC length changes with a higher dynamic range (e.g., systems that use tunable lasers) and a resolution on the order of tens of nanometers or better. Such a method would have a much broader applicability to any FPC-based sensor that would benefit from the high precision of phase tracking methods over a wider measurement range. Optical coherence tomography techniques have been developed using spectral multiplexing [42] and phase gradients [43] to introduce secondary features into the interference pattern to extend the dynamic range of the measurements while maintaining the precision of phase tracking. Other interrogation methods have achieved similar results with free-space optics by using the optical frequency comb of a femtosecond laser pulse [44] and a pseudo-dithering scheme based on laser feedback [45]. Because the ultimate application for the present work is for sensors deployed in harsh environments, free-space optics would be challenging to implement. Therefore, this work focuses on optical fiber-guided signals.

2. Methods

2.1 Measurement theory

An FPC is an optical resonator that comprises of two parallel reflective surfaces separated by a medium with group index $n_g$ and cavity length $l$. Light that is incident on the first surface—usually launched via an optical fiber—is partially reflected. The transmitted portion of the light passes through the medium and is reflected at the second surface before it recombines with the light that was reflected from the first surface to generate an interference pattern. The intensity of the reflected interference pattern (reflectance, $I_R$) generated by an ideal lossless FPC is given by

$$I_R(\lambda) = \bigg(\dfrac{\big( \sqrt{R_1} - \sqrt{R_2}\big)^2 + 4 \sqrt{R_1 R_2}\sin^2 (\phi/2) }{\big(1-\sqrt{R_1R_2}\big)^2 + 4\sqrt{R_1R_2}\sin^2 (\phi/2)}\bigg) I_0,$$
where $R_1$ and $R_2$ are the reflectivities of the first and second reflections, respectively; $I_0$ is the incident light intensity; $\phi =4 \pi n_g l/\lambda$ is the phase difference between the two reflected beams; and $\lambda$ is the wavelength of the incident light. This expression is consistent with previous derivations when $R_1 \neq R_2$ [46] and $R_1 = R_2$[47]. Interrogating an FPC with a broadband light source produces a periodic spectral interference pattern. The free spectral range can be calculated based on the peaks or valleys in $I_R$. Henceforth, only the term peaks will be used in this paper. For consecutive interference peaks $m$ and $m-1$ with wavelengths $[\lambda ]_m$ and $[\lambda ]_{m-1}$, the free spectral range (FSR) is equal to $[\Delta \lambda _{fsr}]_{m} = [\lambda ]_m - [\lambda ]_{m-1}$. Assuming a normal incidence, the FSR is related to $l$ at time $t$ according to Eq. (2):
$$[\Delta \lambda_{fsr}]_{m}(t) = \dfrac{[\lambda]_m(t)^2}{2 n_g l(t) + [\lambda]_m(t)} .$$
Based on peak $m$, the cavity length can be calculated per Eq. (3). The true time-dependent cavity length $l(t)$ can be determined by using various statistical methods, and one such one method is described in more detail in Section 2.2.
$$l_m(t) = \dfrac{[\lambda]_m(t)\big([\lambda]_m(t) - [\Delta \lambda_{fsr}]_{m}(t)\big)}{2 n_g [\Delta \lambda_{fsr}]_{m}(t)} .$$
The method for determining the cavity length is generally called period tracking. The main advantage of period tracking techniques is that they offer an absolute measurement of the cavity length. Period tracking methods can also use a high-intensity tunable laser source (TLS) to provide a high signal-to-noise ratio. The disadvantage of period tracking techniques is their relatively low displacement resolution.

Relative displacements can also be measured from the same spectral interference pattern using phase tracking methods. With this approach, one or multiple peak wavelengths in the interference pattern are tracked over time, and the relative change in cavity length is proportional to the relative change in the peak wavelengths. Equation (4) describes the relationship between the relative change in cavity length ($\Delta l_m(t) = l_m(t) - l_m(t_0)$) and the change in wavelength $[\Delta \lambda ]_m$ for peak $m$ relative to the initial wavelength at time $t_0$:

$$\dfrac{\Delta l_m(t) }{l_m(t_0)} = \dfrac{[\Delta \lambda]_m(t) }{[\lambda]_m(t_0)}.$$
Because phase tracking methods depend only on the change in one peak, they can resolve much smaller displacements [48]—nanometer resolution is possible—but the maximum resolvable change in cavity length is limited because of phase ambiguities. Per the definition of $\phi$, the maximum resolvable change in cavity length for an air cavity ($n_g=1$) is equal to $\lambda /2$, or less than 1 $\mu$m for typical wavelengths on the order of 1,550 nm. This range limitation can be overcome by performing a cumulative sum of incremental phase shifts, but aliasing errors can still occur if the incremental phase shift exceeds 2$\pi$ [32].

2.2 Approach

Drawing inspiration from [35], a predictor-corrector scheme is described herein to realize the high precision of phase tracking techniques while also measuring absolute deflections greater than the limitations imposed by the FSR (i.e., phase ambiguity). Figure 1 illustrates the scheme with an example interference pattern, and Fig. 2 provides a flow diagram of the proposed scheme. First, Eq. (3) is used to provide a coarse measurement of the initial cavity length, $l_m(t_0)$, and peak locations $[\lambda ]_m(t_0)$. A mean operation, including a removal of the upper and lower 25% of the length estimates $l_m(t_0)$ to mitigate the effect of outliers, is used to determine $l(t_0)$ from $l_m(t_0)$. This process is then repeated for all $t$ to determine the time-dependent cavity lengths $l(t)$. The predicted peak wavelengths ($[\lambda _P]_m(t)$) for $t>0$ are calculated using the change in $l$ and the initial peak locations ($[\lambda ]_m(t_0)$) as

$$[\lambda_{P}]_m(t) = \bigg( \dfrac{l(t)-l(t_0)}{l(t_0)} + 1\bigg) [\lambda]_m(t_0).$$
From this coarse estimate, the corrected peak locations ($[\lambda _C]_m(t)$) are identified by searching for the closest peaks to $[\lambda _P]_m(t)$. The corrected peak locations are then used to precisely determine the change in cavity length $\Delta l_m(t)$ (Eq. [6]), or $\Delta l(t)$ after taking the median over all $m$ for which $\lambda _{min}<[\lambda _P]_m<\lambda _{max}$:
$$\Delta l_m(t) = \bigg(\dfrac{[\lambda_C]_m(t) - [\lambda]_m(t_0)}{[\lambda]_m(t_0)} \bigg) l(t_0).$$
By defining $[\lambda _C]_m(t)$ as the peak $[\lambda ]_m(t_0)$ nearest to $[\lambda _P]_m(t)$, the predicted peak can move closer to a peak that is higher or lower than the true peak. When this occurs, $[\lambda _C]_m(t)$ will differ from $[\lambda ]_m(t)$ by $[\Delta \lambda _{fsr}]_m$. A median filter ensures that these extrema do not significantly affect the determination of $\Delta l(t)$. Finally, a median filter may be applied in the time domain over a fixed-length sample window. This operation further reduces the effects of misidentifying the true peak, but it also could filter out shorter transient displacements.

 figure: Fig. 1.

Fig. 1. Schematic of (a) an FPC and (b) the associated interference patterns showing the reflected intensity $I_R$ from an idealized FPC and the peaks shifting from blue to red as the cavity length changes. The true wavelength shift $\Delta \lambda _m$ is shown, along with the coarse estimate $\lambda _P$. The upper plot shows the probability density $\rho$ for the predicted wavelength $\lambda _P$ determined using period tracking. The peak nearest $\lambda _P$ is then identified as $\lambda _C$.

Download Full Size | PDF

 figure: Fig. 2.

Fig. 2. Schematic flow diagram for the hybrid phase-period demodulation scheme where $k$ is the sample window used for the temporal median filter.

Download Full Size | PDF

2.3 Hardware and data acquisition

The predictor-corrector scheme was evaluated using a sensor that was described previously for measuring the external pressure [32] or corrosion [25] of a thin metal diaphragm. The FPC was formed by first by using UAM to embed an optical fiber in a nickel-200 block, which was then welded within a cylindrical nickel-200 housing. The housing was terminated with a mirror-polished thin diaphragm with the fiber terminus separated from the diaphragm’s surface by $404 \pm 5$ $\mu$m (measured using period tracking). Even after polishing, a thin film could remain on the inner surface of the diaphragm, which could lead to the creation of a dual-cavity FPC, as described in [49,50]. However, no such artifacts were observed during previous testing of similar sensors [25,32]. Diaphragm deflections were initiated either by manually depressing the outer surface of the diaphragm or via internal pneumatic pressurization of the sensor. The FPC was interrogated using a Hyperion si155 optical sensing instrument (Luna Innovations Inc., Blacksburg, Virginia, USA). The sensing instrument contains a TLS with a wavelength range of 1,460 to 1,620 nm (Fig. 3). A photodetector (PD1) records light launched by the TLS, which passes through a reference fiber located on the lower-amplitude arm of a 90/10 coupler. The other arm guides the light from the TLS into the FPC through a 50/50 coupler. The reflected light from the FPC passes through the other input arm of the 50/50 coupler and is recorded by a second photodetector (PD2). The optical spectra recorded by PD1 and PD2 are combined to produce the interference pattern resulting from light reflected from the FPC. The optical interrogator includes a peak-finding algorithm that returns the peak locations at a sampling rate of 1 kHz. These peaks were first analyzed using phase tracking, period tracking, and hybrid algorithms implemented in Python 3.8.11. These same algorithms were later integrated into the data acquisition software in LabVIEW 2019 (version 19.0f2; National Instruments, Austin, Texas, USA).

 figure: Fig. 3.

Fig. 3. Schematic showing the optical network in the Hyperion si155 optical interrogator used to measure the FPC length.

Download Full Size | PDF

3. Results

3.1 Simulated deflections

Although the proposed demodulation scheme is not subject to phase ambiguity like phase tracking techniques, it is limited by the range of the TLS. Rearranging Eq. (4) such that $[\lambda ]_m(t)$ spans the TLS range $[\lambda _{min}, \lambda _{max}]$, the FPC length at time $t$, relative to the initial length, is bounded by

$$\dfrac{\lambda_{min}}{\lambda_{max}} < \dfrac{l(t)}{l(t_0)} < \dfrac{\lambda_{max}}{\lambda_{min}}.$$
These bounds consider the FPC length changes required to shift a peak across the wavelength range of the TLS: from $\lambda _{min}$ to $\lambda _{max}$ or from $\lambda _{max}$ to $\lambda _{min}$. Using a 0.4 mm FPC length and a TLS wavelength range of $160$ nm centered at 1,540 nm, the maximum detectable FPC deflection is approximately $\pm 0.04$ mm (10%). Beyond these limits, the peak wavelength shifts outside of the TLS range.

To characterize the theoretical performance of the proposed hybrid demodulation scheme, FPC length changes were simulated over 99% of the allowable deflection range (Fig. 4). Figure 4(a) shows the imposed FPC deflections used to generate the simulated interference patterns and the FPC deflections that were determined using the hybrid scheme with the addition of random errors in the peak identification (wavelength domain). These random errors are meant to simulate the effects of false peak identification resulting from the coarse estimate obtained from period tracking. As the cavity length was varied vs. time (Fig. 4(a)), the added error was randomly sampled across a uniform distribution between $\pm 1$ nm ($<[\Delta \lambda _{fsr}]/2 \approx 1.5$ nm per Eq. (2)) and $\pm 3$ nm ($>[\Delta \lambda _{fsr}]/2$). Figure 4(b) shows the error between the imposed vs. demodulated changes in the FPC length as a function of time (left side) and as a histogram (right side). The apparent discrete nature of the error distribution results from false peak identification during the coarse estimate (predictor) step in the algorithm, which is evident when the peak wavelength error exceeds $[\Delta \lambda _{fsr}]/2$. This phenomenon is discussed in more detail in Section 3.4.

 figure: Fig. 4.

Fig. 4. Simulated FPC deflections demodulated using the proposed hybrid scheme. The calculated deflections with $\pm 1$ nm and $\pm 3$ nm error applied to the peak wavelength selection, along with the imposed deflections (a), and the difference between the calculated vs. imposed deflections are also shown (b), with a histogram with 40 nm bins showing the symmetric discrete error distribution that results from identifying false peaks in the demodulation scheme.

Download Full Size | PDF

3.2 Repeated large negative deflections

To demonstrate the proposed hybrid technique and compare the measured displacements with those obtained using conventional period and phase tracking techniques, the diaphragm was repeatedly manually depressed to produce deflections on the order of 15–20 $\mu$m (Fig. 5). As expected, the period tracking and hybrid methods were able to resolve deflections beyond the limitations of the phase-tracking method. The displacements determined by using the hybrid method matched those determined by using the period tracking method but with a noticeably higher precision.

 figure: Fig. 5.

Fig. 5. Experimental measurements of FPC deflections demodulated using period tracking (gray), phase tracking (red), and the proposed hybrid technique (blue) as the diaphragm was manually depressed many times over the course of 10 s. The upper plot shows the phase tracking data with the initial phase measurement used as a reference without performing the cumulative summation. This illustrates how phase tracking is unsuitable for tracking such large FPC deflections. The cumulative phase tracking method is shown in the lower plot to illustrate how cumulative sums of incremental displacements determined using phase tracking can lead to error propagation.

Download Full Size | PDF

3.3 Step-wise positive deflections

More controlled diaphragm deflections were initiated by pneumatically pressurizing the interior of the sensor. A split-rail proportional–integral–derivative (PID) controller, described previously [25], was used to produce a step-wise pressure waveform. An independent pressure transducer was used to record the internal pressure simultaneously with the FPC displacements. Both sets of measurements were recorded at a rate of 1 kHz as the pressure set point was increased from 100 to 150 kPa in steps of 12.5 kPa, each lasting 30 s (Fig. 6). These data were analyzed with and without using median time-domain filtering with a 25 sample window (i.e., affecting $k$ in the gray box in Fig. 2). Figure 6(a) demonstrates that both the period tracking and hybrid demodulation methods produce deflections that closely follow the increasing pressure steps at steady state, which was considered to be the last $\sim$3 s of each hold-step before the PID set point was changed. To characterize the variation of the data around each steady-state point, distributions were created for all data, both hybrid and period tracking methods, within each $\sim$3 s steady-state window after subtracting the corresponding median deflection obtained using the unfiltered hybrid data at each pressure set point. For both the hybrid and period tracking techniques, results are shown with and without a time-domain median filter with a 25 sample window. Without time-domain median filtering, the period tracking method yields a similar steady-state standard deviation of $\pm 0.07$ to $0.09$ $\mu$m, whereas the hybrid method produces a standard deviation of $\pm 0.002$ to $0.2$ $\mu$m (Fig. 6(b), Table 1). However, the primary source of error when using the hybrid method originates when $\lambda _C$ is incorrectly identified as a peak to the left or right ($\lambda _{m-1}$ or $\lambda _{m+1}$ in Fig. 1), resulting in a discrete phase shift error of $\pm 2\pi$, or a discrete error of $\pm \lambda /2$ in the FPC length. These errors collectively manifest as impulse or salt-and-pepper noise and are visible near $\Delta l - \bar {l} = -0.7 \mu m \approx -\lambda /2$ in Fig. 6(b). Therefore, this collection of incorrect peaks inflates the standard deviation associated with deflection measurements made with the hybrid method. Applying median time-domain filters (Fig. 6(c)) is a simple way to mitigate these errors [51]. The median operation, rather than the mean, is more appropriate for removing errors that occur when the incorrect peak selection is skewed to one side (Fig. 6(b)), especially because the error distribution is discrete and appears in $<10\%$ of the measurements. This filter improved the precision of the period tracking measurement by approximately two- to three-fold and had relatively little effect on the steady-state deflection measurements (Table 1). However, the same filter applied to deflection measurements produced using the hybrid method resulted in up to a hundred-fold increase in precision because it mitigated the effects of incorrect peak selection.

 figure: Fig. 6.

Fig. 6. The internal FPC pressure was increased between 100 and $\sim$150 kPa and held for $\sim$30 s at each step. (a) The pressure measurements (gray, right axis) are plotted along with the deflection measurements (left axis) determined using the unfiltered hybrid (transparent blue) and period tracking (transparent red) methods, as well as the filtered hybrid (opaque blue) and period tracking (opaque red) methods. Histograms with a bin size of 20 nm were generated from the steady-state measurements comprising the last $\sim$3 s of each $\sim$30 s hold period and were compared for the (b) unfiltered and (c) filtered hybrid and period tracking methods. The histogram data are represented as the percent of total counts in each bin. The error is defined as the difference between the demodulated FPC length change and the value obtained using the mean of the unfiltered hybrid data at each pressure set point.

Download Full Size | PDF

Tables Icon

Table 1. Steady-state internal FPC pressures and deflections measured using the period tracking and hybrid methods during the last 3 s of each 30 s constant-pressure period. Data shown as mean(median) $\pm$ standard deviation.

3.4 Propagation of uncertainty

The propagation of uncertainty for each FPC interrogation method was analyzed using the variance formula, as in [48]. The error in $\Delta l$ when using phase tracking ($[\sigma _{phase}]_m$) for peak $m$ is given by

$$\dfrac{[\sigma_{phase}]_m}{l(t_0)} = \dfrac{\sigma_\lambda}{[\lambda]_{m}(t_0)}.$$
Relative to the initial FPC length, this error is independent of time and only depends on the wavelength uncertainty $\sigma _\lambda$ relative to the initial peak wavelength $[\lambda ]_{m}(t_0)$. The time-dependent error in in $\Delta l$ when using period tracking ($[\sigma _{period}]_m(t)$) is given by
$$\dfrac{[\sigma_{period}]_m(t)}{l(t_0)} = \bigg(\dfrac{\sqrt{[\lambda]_m(t)^4 + [\lambda]_{m+1}(t)^4}}{2 n_g l(t_0)\big([\lambda]_m(t) - [\lambda]_{m+1}(t)\big)^2}\bigg) \sigma_\lambda.$$
If $[\lambda ]_{m+1}(t) \approx [\lambda ]_m(t)$, the numerator in Eq. (4) is simplified to $([\lambda ]_m(t))^2$, which allows the equation to be rewritten as [48]
$$\dfrac{[\sigma_{period}]_m(t)}{l(t_0)} \approx \dfrac{1}{2n_g l(t_0)} \bigg( \dfrac{[\lambda]_m(t)}{[\Delta \lambda_{fsr}]_{m}(t)} \bigg)^2 \sigma_\lambda.$$
Similar to phase tracking uncertainties, this uncertainty is proportional to the wavelength resolution but it is also inversely proportional to the square of $[\Delta \lambda _{fsr}]_m$. This explains why period tracking techniques have poorer uncertainty, particularly for larger cavity lengths (i.e., smaller $[\Delta \lambda _{fsr}]_m$).

The proposed hybrid method is essentially identical to phase tracking as long as the coarse estimate obtained using period tracking (Eq. (5)) properly identifies the true peak. If $[\lambda _P]_m(t)$ provides a sufficiently good approximation for $[\lambda ]_m(t)$ (i.e., within $-[\Delta \lambda _{fsr}]_m(t)/2 < [\lambda _C]_m(t) - [\lambda ]_m(t) < [\Delta \lambda _{fsr}]_m(t)/2$), then the correct corresponding peak is selected so that $[\lambda _C]_m(t) = [\lambda ]_m(t)$. However, if $[\lambda _P]_m(t)$ is outside this window, an incorrect peak will be selected for $[\lambda _C]_m(t)$, which differs from $[\lambda ]_m(t)$ by approximately an integer multiple of $[\lambda _{fsr}]_m(t)$ in addition to the small correction required to identify $[\lambda ]_m(t)$. Therefore, the error in the hybrid method $[\sigma _{hybrid}]_m$ can be represented as

$$[\sigma_{hybrid}]_m = [\sigma_{phase}]_m + k_m l(t_0)\bigg( \dfrac{[\Delta \lambda_{fsr}]_m(t)}{[ \lambda]_m(t_0)} \bigg),$$
where $k_m$ is an integer. If more than half of the peaks are properly selected ($k_m = 0$), the median filters that are applied in both the $m$ and $t$ domains could essentially eliminate the contribution due to incorrect peak selection. This would be the case even if the distribution of $k_m$ is not perfectly symmetric about $k_m = 0$ (Fig. 6(b-c)). In this hypothetical scenario, as the window length of the time-domain filter increases, the probability of incorrect peak selection decreases and the error in the hybrid method approaches the error for phase tracking. The cost of increasing the time-domain window length is the inability to accurately resolve transient displacements that occur within the filter window.

The data in Table 1 demonstrates that standard deviation of the filtered measurements made using the hybrid method are on the order of 1 nm. If the distribution of $k_m$ is truly symmetric around $k_m = 0$ such that $[\sigma _{hybrid}]_m = [\sigma _{phase}]_m \approx 1$ nm, then Eq. (3) can be used with $l(t_0)=404$ $\mu$m and $[\lambda ]_{m}(t_0) \approx 1540$ nm to give $\sigma _\lambda = 4$ pm. The manufacturer of the optical interrogator quotes a wavelength accuracy and repeatability of 1 pm, which is a similar order of magnitude. This suggests that the error for the hybrid does indeed approach that of phase tracking.

4. Discussion

The hybrid method proposed herein for demodulating FPC lengths from their optical spectra demonstrably improves precision compared with period tracking methods without relying on cumulative sums of incremental phase tracking measurements, which suffer from error propagation and drift. Rather than relying on introducing secondary features into the interference pattern by manipulating the incident light or affecting the optical signal returning to the interferometer, as in previous work [32,4245], the hybrid method demonstrated in this work uses the primary interference pattern from the FPC alone via a filtering scheme. This is particularly attractive because it enables the hybrid method to be implemented as a simple drop-in replacement for demodulating data from FPC-based sensors. Although the hybrid method overcomes the challenges associated with conventional peak-tracking algorithms, it is still limited by the wavelength range of the tunable laser used to interrogate the FPC. Ideally, the maximum wavelength shift observable using the hybrid method is equivalent to the tuning range of the TLS, although in practice this is more limited. In the present study, >50 identifiable peaks initially appeared in the interference pattern recorded by the optical interrogation system. As the FPC length changed, many of those initial peaks were shifted outside the range of the TLS. When the number of remaining peaks that fell within the range of the TLS decreased below approximately 10–15, the algorithm failed to reliably demodulate the FPC length. This limitation in the maximum resolvable FPC displacement could be overcome by updating the reference spectrum at intermediate steps similar to the adaptive methods developed for optical frequency domain reflectometry [52,53].

The application of a time-domain median filter to the FPC deflection measurements produced using the hybrid method effectively removes the impulse noise from the measurement but also could introduce a delay and smooths over short-duration physical FPC length changes. The effect of this filtering depends on the size of the sample window. Depending on the desired accuracy and measurement frequencies of interest, the time-domain filtering could be adjusted, as needed. For the original application mentioned herein (i.e., corrosion monitoring via step-wise, internal static pressurization of an FPC), the frequency limitations imposed by the time-domain filtering would not affect the ability to perform these measurements.

5. Conclusion

Optical fiber-based FPCs are a popular sensor for measuring fine displacements because of their high precision. Spectral interrogation methods for FPCs are particularly popular because they produce measurements with precision on the order of nanometers using phase tracking methods. However, the maximum resolvable displacements using methods are limited by phase ambiguities. Period tracking methods overcome this limitation and enable the calculation of an absolute FPC length but do not offer the same precision. To overcome the range limitation of phase tracking methods and the relatively poor precision of period tracking methods, a hybrid approach was demonstrated. This method used a period tracking method to predict a coarse wavelength shift, which was then fine-tuned to select the nearest peak in the interference spectrum to determine a much more precise change in the FPC length. Using a median filter to remove impulse noise, this hybrid method was demonstrated to produce consistent measurements with standard deviations on the order of nanometers, a hundred-fold lower than the precision obtained using similarly filtered period tracking methods. Furthermore, this precision is maintained while measuring cavity deflections >20 $\mu$m (>5% of the initial FPC length).

Funding

Office of Nuclear Energy.

Acknowledgments

This manuscript has been authored by UT-Battelle, LLC, under contract DE-AC05-00OR22725 with the US Department of Energy (DOE). The US government retains and the publisher, by accepting the article for publication, acknowledges that the US government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for US government purposes. DOE will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan). The authors would like to thank Dr. Brandon Wilson and Dr. Adrian Schrell for their valuable technical feedback regarding this work. This work was supported by the Versatile Test Reactor Program of the US Department of Energy, Office of Nuclear Energy.

Disclosures

DCS, AB, and CMP have submitted a patent application related to this work.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. G. Berkovic and E. Shafir, “Optical methods for distance and displacement measurements,” Adv. Opt. Photonics 4(4), 441–471 (2012). [CrossRef]  

2. A. Ukil, H. Braendle, and P. Krippner, “Distributed temperature sensing: Review of technology and applications,” IEEE Sens. J. 12(5), 885–892 (2012). [CrossRef]  

3. Q. Yu and X. Zhou, “Pressure sensor based on the fiber-optic extrinsic fabry-perot interferometer,” Photonic Sens. 1(1), 72–83 (2011). [CrossRef]  

4. J. X. Fang, H. F. Taylor, and H. S. Choi, “Fiber-optic fabry–perot flow sensor,” Microw. Opt. Technol. Lett. 18(3), 209–211 (1998). [CrossRef]  

5. J.-N. Wang and C.-Y. Luo, “Long-period fiber grating sensors for the measurement of liquid level and fluid-flow velocity,” Sensors 12(4), 4578–4593 (2012). [CrossRef]  

6. S. Khaliq, S. W. James, and R. P. Tatam, “Fiber-optic liquid-level sensor using a long-period grating,” Opt. Lett. 26(16), 1224–1226 (2001). [CrossRef]  

7. C. M. Petrie and J. L. McDuffee, “Liquid level sensing for harsh environment applications using distributed fiber optic temperature measurements,” Sens. Actuators, A 282, 114–123 (2018). [CrossRef]  

8. X. Liu, B. Jin, Q. Bai, Y. Wang, D. Wang, and Y. Wang, “Distributed fiber-optic sensors for vibration detection,” Sensors 16(8), 1164 (2016). [CrossRef]  

9. J. M. López-Higuera, L. R. Cobo, A. Q. Incera, and A. Cobo, “Fiber optic sensors in structural health monitoring,” J. Lightwave Technol. 29(4), 587–608 (2011). [CrossRef]  

10. S. D. Rountree, O. J. Ohanian, A. Boulanger, D. Kominsky, M. Davis, M. Wang, K. Chen, A. Leong, J. Zhang, F. Scurti, N. Oliver, J. Schwartz, T. E. Blue, and J. Jones, “Multi-parameter fiber optic sensing for harsh nuclear environments,” in Fiber Optic Sensors and Applications XVII, vol. 11739 (International Society for Optics and Photonics, 2021), p. 117390H.

11. A. Rose, “Devitrification in annealed optical fiber,” J. Lightwave Technol. 15(5), 808–814 (1997). [CrossRef]  

12. T. W. Wood, B. Blake, T. E. Blue, C. M. Petrie, and D. Hawn, “Evaluation of the performance of distributed temperature measurements with single-mode fiber using rayleigh backscatter up to 1000 °c,” IEEE Sens. J. 14(1), 124–128 (2014). [CrossRef]  

13. J. T. Jones, D. C. Sweeney, A. Birri, C. M. Petrie, and T. E. Blue, “Calibration of distributed temperature sensors using commercially available smf-28 optical fiber from 22 °c to 1000 °c,” IEEE Sens. J. 22(5), 4144–4151 (2022). [CrossRef]  

14. B. A. Wilson, C. M. Petrie, and T. E. Blue, “High-temperature effects on the light transmission through sapphire optical fiber,” J. Am. Ceram. Soc. 101(8), 3452–3459 (2018). [CrossRef]  

15. C. M. Petrie and T. E. Blue, “In situ thermally induced attenuation in sapphire optical fibers heated to 1400 °c,” J. Am. Ceram. Soc. 98(2), 483–489 (2015). [CrossRef]  

16. G. Cheymol, H. Long, J. F. Villard, and B. Brichard, “High level gamma and neutron irradiation of silica optical fibers in cea osiris nuclear reactor,” IEEE Trans. Nucl. Sci. 55(4), 2252–2258 (2008). [CrossRef]  

17. C. M. Petrie, A. Birri, and T. E. Blue, “Optical transmission and dimensional stability of single-crystal sapphire after high-dose neutron irradiation at various temperatures up to 688 °c,” J. Nucl. Mater. 559, 153432 (2022). [CrossRef]  

18. C. M. Petrie, D. P. Hawn, W. Windl, and T. E. Blue, “Reactor radiation-induced attenuation in fused silica optical fibers heated up to 1000°c,” J. Non-Cryst. Solids 409, 88–94 (2015). [CrossRef]  

19. C. M. Petrie, A. Birri, and T. E. Blue, “High-dose temperature-dependent neutron irradiation effects on the optical transmission and dimensional stability of amorphous fused silica,” J. Non-Cryst. Solids 525, 119668 (2019). [CrossRef]  

20. C. M. Petrie and T. E. Blue, “In situ reactor radiation-induced attenuation in sapphire optical fibers heated up to 1000°c,” Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 342, 91–97 (2015). [CrossRef]  

21. C. M. Petrie, W. Windl, and T. E. Blue, “In-situ reactor radiation-induced attenuation in sapphire optical fibers,” J. Am. Ceram. Soc. 97(12), 3883–3889 (2014). [CrossRef]  

22. A. Gusarov, F. Berghmans, A. Fernandez, O. Deparis, Y. Defosse, D. Starodubov, M. Decreton, P. Megret, and M. Bondel, “Behavior of fibre bragg gratings under high total dose gamma radiation,” IEEE Trans. Nucl. Sci. 47(3), 688–692 (2000). [CrossRef]  

23. D. P. Hawn, C. M. Petrie, T. E. Blue, and W. Windl, “In-situ gamma radiation induced attenuation in silica optical fibers heated up to 600°c,” J. Non-Cryst. Solids 379, 192–200 (2013). [CrossRef]  

24. C. M. Petrie, B. Wilson, and T. E. Blue, “In situ gamma radiation-induced attenuation in sapphire optical fibers heated to 1000°c,” J. Am. Ceram. Soc. 97(10), 3150–3156 (2014). [CrossRef]  

25. D. C. Sweeney, A. M. Schrell, and C. M. Petrie, “Pressure-driven fiber-optic sensor for online corrosion monitoring,” IEEE Trans. Instrum. Meas. 70, 1–10 (2021). [CrossRef]  

26. A. Hehr and M. Norfolk, “A comprehensive review of ultrasonic additive manufacturing,” Rapid Prototyping Journal (2019).

27. R. Friel and R. Harris, “Ultrasonic additive manufacturing – a hybrid production process for novel functional products,” Procedia CIRP 6, 35–40 (2013). [CrossRef]  

28. A. Hehr, M. Norfolk, J. Wenning, J. Sheridan, P. Leser, P. Leser, and J. A. Newman, “Integrating fiber optic strain sensors into metal using ultrasonic additive manufacturing,” JOM 70(3), 315–320 (2018). [CrossRef]  

29. C. M. Petrie, N. Sridharan, M. Subramanian, A. Hehr, M. Norfolk, and J. Sheridan, “Embedded metallized optical fibers for high temperature applications,” Smart Mater. Struct. 28(5), 055012 (2019). [CrossRef]  

30. C. M. Petrie, N. Sridharan, A. Hehr, M. Norfolk, and J. Sheridan, “High-temperature strain monitoring of stainless steel using fiber optics embedded in ultrasonically consolidated nickel layers,” Smart Mater. Struct. 28(8), 085041 (2019). [CrossRef]  

31. H. C. Hyer, D. C. Sweeney, and C. M. Petrie, “Functional fiber-optic sensors embedded in stainless steel components using ultrasonic additive manufacturing for distributed temperature and strain measurements,” Addit. Manuf. 52, 102681 (2022). [CrossRef]  

32. D. C. Sweeney, A. M. Schrell, Y. Liu, and C. M. Petrie, “Metal-embedded fiber optic sensor packaging and signal demodulation scheme towards high-frequency dynamic measurements in harsh environments,” Sens. Actuators, A 312, 112075 (2020). [CrossRef]  

33. D. Tosi, “Simultaneous detection of multiple fiber-optic fabry–perot interferometry sensors with cepstrum-division multiplexing,” J. Lightwave Technol. 34(15), 3622–3627 (2016). [CrossRef]  

34. H. Bae, D. Yun, H. Liu, D. A. Olson, and M. Yu, “Hybrid miniature fabry–perot sensor with dual optical cavities for simultaneous pressure and temperature measurements,” J. Lightwave Technol. 32(8), 1585–1593 (2014). [CrossRef]  

35. B. Zhang, X. Tong, P. Hu, Q. Guo, Z. Zheng, and C. Zhou, “Wavelet phase extracting demodulation algorithm based on scale factor for optical fiber fabry-perot sensing,” Opt. Express 24(26), 29506–29511 (2016). [CrossRef]  

36. Y. Wu, L. Xia, N. Cai, and L. Zhu, “A highly precise demodulation method for fiber fabry-perot cavity through spectrum reconstruction,” IEEE Photonics Technol. Lett. 30(5), 435–438 (2017). [CrossRef]  

37. B. Qi, D. E. Winder, and Y. Liu, “Quadrature phase-shifted optical demodulator for low-coherence fiber-optic fabry-perot interferometric sensors,” Opt. Express 27(5), 7319–7329 (2019). [CrossRef]  

38. Y. Liu, R. Strum, D. Stiles, C. Long, A. Rakhman, W. Blokland, D. Winder, B. Riemer, and M. Wendel, “Digital phase demodulation for low-coherence interferometry-based fiber-optic sensors,” Opt. Commun. 411, 27–32 (2018). [CrossRef]  

39. D. Sweeney, A. Schrell, and C. Petrie, “The transient thermal response of a pressure-driven fabry-pérot cavity,” in 12th Nuclear Plant Instrumentation, Control and Human-Machine Interface Technologies (NPIC & HMIT) (American Nuclear Society, 2021).

40. D. Sweeney, A. Schrell, and C. Petrie, “Compensation scheme for radiation-induced attenuation in optical fibers interrogated using low-coherence interferometry,” Transactions of the American Nuclear Society 122, 291–294 (2020). [CrossRef]  

41. Y. Liu, B. Qi, and D. E. Winder, “Faraday michelson interferometers for signal demodulation of fiber-optic sensors,” J. Lightwave Technol. 39(8), 2552–2558 (2021). [CrossRef]  

42. B. Dong, Y. Zhang, W. Zhang, Z. He, S. Xie, and Y. Zhou, “Highly sensitive, wide dynamic range displacement sensor combining chromatic confocal system and phase-sensitive spectral optical coherence tomography,” Opt. Express 25(5), 5426–5430 (2017). [CrossRef]  

43. V. Y. Zaitsev, A. L. Matveyev, L. A. Matveev, G. V. Gelikonov, E. V. Gubarkova, N. D. Gladkova, and A. Vitkin, “Hybrid method of strain estimation in optical coherence elastography using combined sub-wavelength phase measurements and supra-pixel displacement tracking,” J. Biophotonics 9(5), 499–509 (2016). [CrossRef]  

44. J. Park, J. Jin, J.-A. Kim, and J. Kim, “Absolute distance measurement method without a non-measurable range and directional ambiguity based on the spectral-domain interferometer using the optical comb of the femtosecond pulse laser,” Appl. Phys. Lett. 109(24), 244103 (2016). [CrossRef]  

45. S. Donati, G. Giuliani, and S. Merlo, “Laser diode feedback interferometer for measurement of displacements without ambiguity,” IEEE J. Quantum Electron. 31(1), 113–119 (1995). [CrossRef]  

46. Y.-J. Rao, “Recent progress in fiber-optic extrinsic fabry–perot interferometric sensors,” Opt. Fiber Technol. 12(3), 227–237 (2006). [CrossRef]  

47. J. Santos, A. Leite, and D. A. Jackson, “Optical fiber sensing with a low-finesse fabry–perot cavity,” Appl. Opt. 31(34), 7361–7366 (1992). [CrossRef]  

48. Y. Huang, T. Wei, Z. Zhou, Y. Zhang, G. Chen, and H. Xiao, “An extrinsic fabry–perot interferometer-based large strain sensor with high resolution,” Meas. Sci. Technol. 21(10), 105308 (2010). [CrossRef]  

49. J. Monzón, L. Sánchez-Soto, and E. Bernabeu, “Influence of coating thickness on the performance of a fabry–perot interferometer,” Appl. Opt. 30(28), 4126–4132 (1991). [CrossRef]  

50. Q. Ren, P. Jia, G. An, J. Liu, G. Fang, W. Liu, and J. Xiong, “Dual-wavelength demodulation technique for interrogating a shortest cavity in multi-cavity fiber-optic Fabry – Pérot sensors,” Opt. Express 29(20), 32658–32669 (2021). [CrossRef]  

51. R. H. Chan, C.-W. Ho, and M. Nikolova, “Salt-and-pepper noise removal by median-type noise detectors and detail-preserving regularization,” IEEE Transactions on Image Processing 14(10), 1479–1485 (2005). [CrossRef]  

52. D. C. Sweeney, D. M. Sweeney, and C. M. Petrie, “Graphical optimization of spectral shift reconstructions for optical backscatter reflectometry,” Sensors 21(18), 6154 (2021). [CrossRef]  

53. D. C. Sweeney, A. M. Schrell, and C. M. Petrie, “An adaptive reference scheme to extend the functional range of optical backscatter reflectometry in extreme environments,” IEEE Sens. J. 21(1), 498–509 (2020). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1.
Fig. 1. Schematic of (a) an FPC and (b) the associated interference patterns showing the reflected intensity $I_R$ from an idealized FPC and the peaks shifting from blue to red as the cavity length changes. The true wavelength shift $\Delta \lambda _m$ is shown, along with the coarse estimate $\lambda _P$. The upper plot shows the probability density $\rho$ for the predicted wavelength $\lambda _P$ determined using period tracking. The peak nearest $\lambda _P$ is then identified as $\lambda _C$.
Fig. 2.
Fig. 2. Schematic flow diagram for the hybrid phase-period demodulation scheme where $k$ is the sample window used for the temporal median filter.
Fig. 3.
Fig. 3. Schematic showing the optical network in the Hyperion si155 optical interrogator used to measure the FPC length.
Fig. 4.
Fig. 4. Simulated FPC deflections demodulated using the proposed hybrid scheme. The calculated deflections with $\pm 1$ nm and $\pm 3$ nm error applied to the peak wavelength selection, along with the imposed deflections (a), and the difference between the calculated vs. imposed deflections are also shown (b), with a histogram with 40 nm bins showing the symmetric discrete error distribution that results from identifying false peaks in the demodulation scheme.
Fig. 5.
Fig. 5. Experimental measurements of FPC deflections demodulated using period tracking (gray), phase tracking (red), and the proposed hybrid technique (blue) as the diaphragm was manually depressed many times over the course of 10 s. The upper plot shows the phase tracking data with the initial phase measurement used as a reference without performing the cumulative summation. This illustrates how phase tracking is unsuitable for tracking such large FPC deflections. The cumulative phase tracking method is shown in the lower plot to illustrate how cumulative sums of incremental displacements determined using phase tracking can lead to error propagation.
Fig. 6.
Fig. 6. The internal FPC pressure was increased between 100 and $\sim$150 kPa and held for $\sim$30 s at each step. (a) The pressure measurements (gray, right axis) are plotted along with the deflection measurements (left axis) determined using the unfiltered hybrid (transparent blue) and period tracking (transparent red) methods, as well as the filtered hybrid (opaque blue) and period tracking (opaque red) methods. Histograms with a bin size of 20 nm were generated from the steady-state measurements comprising the last $\sim$3 s of each $\sim$30 s hold period and were compared for the (b) unfiltered and (c) filtered hybrid and period tracking methods. The histogram data are represented as the percent of total counts in each bin. The error is defined as the difference between the demodulated FPC length change and the value obtained using the mean of the unfiltered hybrid data at each pressure set point.

Tables (1)

Tables Icon

Table 1. Steady-state internal FPC pressures and deflections measured using the period tracking and hybrid methods during the last 3 s of each 30 s constant-pressure period. Data shown as mean(median) ± standard deviation.

Equations (11)

Equations on this page are rendered with MathJax. Learn more.

I R ( λ ) = ( ( R 1 R 2 ) 2 + 4 R 1 R 2 sin 2 ( ϕ / 2 ) ( 1 R 1 R 2 ) 2 + 4 R 1 R 2 sin 2 ( ϕ / 2 ) ) I 0 ,
[ Δ λ f s r ] m ( t ) = [ λ ] m ( t ) 2 2 n g l ( t ) + [ λ ] m ( t ) .
l m ( t ) = [ λ ] m ( t ) ( [ λ ] m ( t ) [ Δ λ f s r ] m ( t ) ) 2 n g [ Δ λ f s r ] m ( t ) .
Δ l m ( t ) l m ( t 0 ) = [ Δ λ ] m ( t ) [ λ ] m ( t 0 ) .
[ λ P ] m ( t ) = ( l ( t ) l ( t 0 ) l ( t 0 ) + 1 ) [ λ ] m ( t 0 ) .
Δ l m ( t ) = ( [ λ C ] m ( t ) [ λ ] m ( t 0 ) [ λ ] m ( t 0 ) ) l ( t 0 ) .
λ m i n λ m a x < l ( t ) l ( t 0 ) < λ m a x λ m i n .
[ σ p h a s e ] m l ( t 0 ) = σ λ [ λ ] m ( t 0 ) .
[ σ p e r i o d ] m ( t ) l ( t 0 ) = ( [ λ ] m ( t ) 4 + [ λ ] m + 1 ( t ) 4 2 n g l ( t 0 ) ( [ λ ] m ( t ) [ λ ] m + 1 ( t ) ) 2 ) σ λ .
[ σ p e r i o d ] m ( t ) l ( t 0 ) 1 2 n g l ( t 0 ) ( [ λ ] m ( t ) [ Δ λ f s r ] m ( t ) ) 2 σ λ .
[ σ h y b r i d ] m = [ σ p h a s e ] m + k m l ( t 0 ) ( [ Δ λ f s r ] m ( t ) [ λ ] m ( t 0 ) ) ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.