Next Article in Journal
Metamorphic Testing of Relation Extraction Models
Next Article in Special Issue
Enhancing the Distributed Acoustic Sensors’ (DAS) Performance by the Simple Noise Reduction Algorithms Sequential Application
Previous Article in Journal
Assembly and Production Line Designing, Balancing and Scheduling with Inaccurate Data: A Survey and Perspectives
Previous Article in Special Issue
Mathematical Modeling of Capillary Drawing Stability for Hollow Optical Fibers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Analysis of the Methods for Fiber Bragg Structures Spectrum Modeling

1
Department of Radiophotonics and Microwave Technologies, Kazan National Research Technical University Named after A.N. Tupolev-KAI, 10 K. Marx St., 420111 Kazan, Russia
2
Department of Radio Communication and Radio Technology Systems, Institute of Telecommunications and Informatics of Turkmenistan, 68 Magtymguly St., Ashgabat 744000, Turkmenistan
*
Author to whom correspondence should be addressed.
Algorithms 2023, 16(2), 101; https://doi.org/10.3390/a16020101
Submission received: 17 January 2023 / Revised: 7 February 2023 / Accepted: 8 February 2023 / Published: 10 February 2023
(This article belongs to the Special Issue Algorithms and Calculations in Fiber Optics and Photonics)

Abstract

:
The work is dedicated to a comparative analysis of the following methods for fiber Bragg grating (FBG) spectral response modeling. The Layer Sweep (LS) method, which is similar to the common layer peeling algorithm, is based on the reflectance and transmittance determination for the plane waves propagating through layered structures, which results in the solution of a system of linear equations for the transmittance and reflectance of each layer using the sweep method. Another considered method is based on the determination of transfer matrices (TM) for the FBG as a whole. Firstly, a homogeneous FBG was modeled using both methods, and the resulting reflectance spectra were compared to the one obtained via a specialized commercial software package. Secondly, modeling results of a π-phase-shifted FBG were presented and discussed. For both FBG models, the influence of the partition interval of the LS method on the simulated spectrum was studied. Based on the analysis of the simulation data, additional required modeling conditions for phase-shifted FBGs were established, which enhanced the modeling performance of the LS method.

1. Introduction

Mathematical modeling of the reflection and/or transmission spectra of passive and active elements of fiber-optic circuits is one of the main approaches used in the design of both fiber communication lines and fiber-optic sensor systems. Mathematical modeling allows not only to predict the results of a certain impact on a fiber-optic system, but also to optimize individual elements of such systems. One of the most commonly used elements of fiber-optic sensor systems, as well as data transmission networks, is a periodic or quasi-periodic structure formed in the core of an optical fiber by changing the refractive index. In particular, a wide class of passive fiber-optic structures, such as fiber Bragg gratings (FBGs), belongs to this type of elements. Over the past decade, in addition to classical FBGs, numerous fiber-optic structures have become widespread: FBGs with one or more phase shifts, apodized and chirped FBGs, superstructured FBGs, FBGs combined with open or closed Fabry–Perot resonators (FPRs), etc.
The basic components of the abovementioned complex or combined fiber-optic structures are FBGs and FPRs. Combining these basic elements in different configurations makes it possible to obtain various sensor designs with specific reflection spectra and different properties and characteristics (Figure 1). An internal (or external) FPR (Figure 1a) can be used as a high-precision point temperature sensor [1]. A conventional uniform FBG (Figure 1b) is also used as a sensing element for temperature, deformation, vibration, etc. measurements [2]. Two FBGs recorded sequentially (Figure 1c) or one over the other [3,4] (Figure 1d) have found applications as addressed fiber Bragg structures [5]. An FBG-based structure with a discrete phase shift (Figure 1e) is used as a high-precision Bragg-type sensor due to the presence of an ultra-narrow transparency window in its reflectance spectrum, the wavelength position of which can be measured with higher accuracy than the shift of the entire FBG spectrum [6]. The FBG structures with two or more phase shifts can also be used as addressed [5] (Figure 1f) or even multi-addressed [7] structures. The sensing elements based on the combination of FBG and FPR (Figure 1g) represent a promising solution for gas concentration measurement systems [8,9,10].
To enhance the performance of the designed elements of the fiber-optic data transmission or sensor systems, it is necessary to have reliable methods and tools that enable spectral characteristics modeling of fiber-optic structures with high accuracy.
The basic element of such tools, which is of particular interest, is the mathematical apparatus that makes it possible to reconstruct the FBG reflection spectrum. The current paper investigated and compared various methods of FBG spectral response modeling. The first method is based on determination of the reflectance and the transmittance of plane waves propagating through a layered structure similarly to the layer peeling algorithm [11], which results in the solution of a system of linear equations for the transmittance and reflectance of each layer using the sweep method (further referred to as the Layer Sweep (LS) method). The second technique is based on the formulation of transfer matrices for the entire FBG structure [12,13] (Transfer Matrix (TM) method). The FBG transfer matrix is often reduced to a matrix with Chebyshev polynomials of the second kind [14], which requires their calculation for each wavelength. Modeling techniques based on recurrence relations, including the ones using Chebyshev polynomials, were deliberately excluded from consideration, since they do not provide performance benefits in comparison with the LS method.
In this paper, firstly, a homogeneous non-apodized FBG was modeled using both methods, and the resulting reflectance spectra were compared to the one obtained via a specialized commercial software. Secondly, modeling results of a π-phase-shifted FBG are presented. For both FBG models, the influence of the partition interval of the LS method on the simulated spectrum was studied. It was established that the asymmetry of the refractive index profile partition near the phase shift relative to its center causes the asymmetry of the spectral response of the simulated grating, resulting in the increased simulation error. Therefore, in the case of phase-shifted FBG, in order to increase the simulation accuracy, it is not enough to decrease the partition interval of the LS method, but it is also required to fulfill certain conditions of the refractive index definition in the model, which are discussed in the manuscript.

2. Modeling Methods

The interference scheme for recording FBGs, as one of the most commonly used schemes for recording fiber Bragg gratings, implements a change in the refractive index along the length of the grating according to the harmonic law [15,16]:
n ( z ) = n 0 + Δ n 2 + Δ n 2 sin ( 2 π Λ z )   ,
where Λ is the FBG period, n0 is the refractive index of the optical fiber core, and Δn is the value of the induced refractive index, which generally can be represented as a function of the grating length, for example, as the Gaussian apodization function [15,17]:
Δ n ( z ) = Δ n exp ( 1 2 ( z z M σ ) 2 ) ,
where zM is the FBG midpoint and σ is the apodization parameter defining the ratio p of the induced refractive index at the FBG edges to its value at the midpoint:
σ 2 = 1 8 H 2 ln p   ,
where H is the FBG length. The FBG is formed in the optical fiber; therefore, the refractive index of the core of the optical fiber is the minimum value of the refractive index in the FBG interval, and the maximum will be the sum of the refractive indices of the fiber and the induced refractive index.
Thus, the initial parameters for the FBG modeling are: the length of the grating H, the refractive index of the optical fiber core n0, the induced refractive index Δn, the FBG period Λ, and the apodization parameter σ.

2.1. Layer Sweep Method

Any medium inside the optical fiber with a known variation of the refractive index can be represented as a set of thin homogeneous “films”, in each of which the refractive index can be considered constant [18]. Thus, the entire FBG is divided by a one-dimensional coordinate grid {zi}, i = 1, N (Figure 2) into homogeneous layers, in which the refractive index is determined as an integral average over the layer thickness:
n i + 1 = 1 z i + 1 z i z i z i + 1 n ( z ) d z     .
It must be noted that the length of the partition interval Δz = zi+1zi may significantly affect the performance of the modeled structure, which is discussed in Section 3. Substituting (1) and (2) into (4), we obtain a function that describes the layer-by-layer dependence of the refractive index on the coordinate z:
n i + 1 = 1 z i + 1 z i z i z i + 1 [ n 0 + 1 2 Δ n exp ( 1 2 ( z z M σ ) 2 ) ( 1 + sin ( 2 π Λ z ) ) ] d z     .
Equation (5) is not an integral-elementary function, but integration in (5) is not required in this case, since the refractive index of the medium must be constant inside each layer, and expression (5) does not take this requirement into account, since the exponential multiplier that normalizes the amplitude of the refractive index variation implies that its amplitude within the layer changes along the coordinate z. The constancy of the induced refractive index within each layer implies the constancy of the exponent in (5) over the layer thickness. Therefore, the exponential multiplier in (5) must be replaced by the average value of the refractive index, which, on the interval from [zi, zi+1] can be taken as its value in the middle of the interval:
exp ( 1 2 ( z ( z M ) k σ k ) 2 )     exp ( 1 2 ( z i + 1 + z i 2 ( z M ) k σ k ) 2 )   + o ( Δ z i ) ,
where Δzi = zi+1zi, and o(x) is the error function, which denotes that the specified function decreases faster than its argument. Substituting (6) into (5), we obtain:
n i + 1 = 1 z i + 1 z i z i z i + 1 [ n 0 + 1 2 Δ n ( exp ( 1 2 ( z i + 1 + z i 2 z M 2 σ ) 2 ) + o ( Δ z i ) ) ( 1 + sin ( 2 π Λ z ) ) ] d z       .
After that, the multiplier of the amplitude becomes independent of the integration variable and the integral can be taken analytically; thus, the dependence of the refractive index over all the layers of an arbitrary layered structure can be expressed as follows:
n i + 1 = n 0 + Δ n 2 exp ( 1 2 ( z i + 1 + z i 2 z M 2 σ ) 2 )   ( 1 Λ π ( z i + 1 z i ) sin ( π ( z i + 1 + z i ) Λ ) sin ( π ( z i + 1 z i ) Λ ) ) + o ( Δ z i )     .
In order to model the reflection spectrum of such a structure, it is necessary to initialize the parameters of the layers based on relations (8), with further application of the classical algorithms of backward and forward sweeps for the layered structure calculation. The electric and magnetic components of the electromagnetic field at each point of each layer are the superposition of two waves (Figure 3) moving in different directions—the direct and the reflected ones:
E i ( z ) = t i exp ( j ( k i z + ω τ ) ) + r i exp ( j ( k i z ω τ ) )   , H i ( z ) = t i w i exp ( j ( k i z + ω τ ) ) r i w i exp ( j ( k i z ω τ ) ) ,     i = 0 , N + 1 ¯     ,
where τ is the time, ti and ri are the transmittance and reflectance of the i-th layer; ki is the wavenumber and wi is the impedance for the i-th layer. The first summand in (9) for both the electric and magnetic components describes the wave moving in the direction of propagation, and the second term describes the waves reflected from the interfaces of the media and moving against the direction of propagation.
The wave number and wave impedance depend on the radiation wavelength, permittivity, and permeability of the medium:
k i = 2 π f ε μ = 2 π c ε μ λ ε i μ i = 2 π λ ε i μ i ,         w i = μ i μ / ε i ε   ,
where λ = c/f is the wavelength, εi is the permittivity and μi is permeability of the i-th layer, respectively, ε and μ are the permittivity and permeability of the vacuum, respectively, and c is the speed of light in vacuum.
The condition for the continuity of electromagnetic waves in space presupposes the equality of electric and magnetic fields at every point in space, including the interfaces between the layers:
{ t i exp ( j k i z i ) + r i exp ( j k i z i ) = t i + 1 exp ( j k i + 1 z i ) + r i + 1 exp ( j k i + 1 z i ) t i exp ( j k i z i ) r i exp ( j k i z i ) w i = t i + 1 exp ( j k i + 1 z i ) r i + 1 exp ( j k i + 1 z i ) w i + 1 ,     i = 0 , N ¯ .
where λ = c/f is the wavelength, εi is the permittivity and μi is permeability of the i-th layer.
It must be noted that the common multiplier exp(−jωτ), which describes the variations of the electric and magnetic waves over time, is reduced in the left and right parts of Equation (11). As a result, the relations (11) give a system of 2(N + 1) linear equations for determination of 2(N + 2) unknown transmission coefficients ti and the reflection coefficients ri in each of the (N + 2) layers.
An additional analysis of the system of Equation (11) makes it possible to exclude from them the permittivity ε and permeability μ of the vacuum, using the following relation:
c ε μ 1   ,
where c is the speed of light. In the system of Equation (11), the left and right parts are reduced by the square root of the product of the absolute permittivity and permeability, and the wave number is reformulated in terms of the wavelength.
Permittivity ε and permeability μ are related to the refractive index as follows:
n = ε μ .
In order to close the resulting system of equations, it is necessary to supplement the system (11) with two more relations. As we postulated above, all radiation coming from the source to the zeroth layer (“Layer 0”) passes through it without loss, and the reflection from the second boundary of the last layer (“Layer N + 1”) is absent due to the infinite remoteness of the boundary. These conditions make it possible to supplement the system of equations with two missing relations for the transmission coefficients of the zeroth layer and the reflection of the last layer, respectively:
t 0 = 1 ,     r N + 1 = 0   .
Combining (11) and (14) we obtain a closed system of 2(N + 1) linear equations for 2(N + 1) unknown transmission and reflection coefficients of each of the layers, which can be written in matrix form:
[ 1 1 1 0 0 0 0 0 0 w 1 w 0 w 0 0 0 0 0 0 0 0 a 1 b 1 c 1 d 1 0 0 0 0 0 e 1 f 1 g 1 h 1 0 0 0 0 0 0 0 a i b i c i d i 0 0 0 0 0 e i f i g i h i 0 0 0 0 0 0 0 0 0 a N b N c N 0 0 0 0 0 0 0 e N f N g N ] × [ r 0 t 1 r 1 t 2 r 2 t i r i t N + 1 ] = [ 1 w 1 0 0 0 0 0 0 0 0 ]   ,  
where the following denotations are introduced:
a i = exp ( j κ i z i ) b i = exp ( j κ i z i ) c i = exp ( j κ i + 1 z i ) d i = exp ( j κ i + 1 z i ) e i = w i + 1 exp ( j κ i z i ) f i = w i + 1 exp ( j κ i z i ) g i = w i exp ( j κ i + 1 z i ) h i = w i exp ( j κ i + 1 z i )     .
A system of linear equations formulated in the matrix form can be solved using the Gaussian method by reducing the matrix of coefficients to a triangular form. Given the specific shape of the matrix elements, each internal row contains only four elements; moreover, rows containing elements in the same positions are grouped by two. Thus, the coefficient matrix is a five-diagonal matrix containing, in addition to the main diagonal, two over-diagonals, and two under-diagonals. Besides that, the two over- and under-diagonals further from the central diagonal are sparse. For systems of linear equations with a matrix of coefficients of three-, four-, and five-diagonal form, the sweep method can be applied, which involves the successive transformation of the matrix rows, thus bringing it to a triangular form. Reduction the coefficient matrix to a diagonal form immediately allows one to determine the value of the unknown standing at the top of the triangle. The matrix of coefficients in (15) can be reduced with equal ease both to the upper triangular form and to the lower triangular form. In the case of reducing the matrix to the upper triangular form, the value of the transmittance of the last layer is automatically calculated. In the case of the lower triangular form, it becomes possible to calculate the reflection coefficient of the zeroth layer.
The reduction of the matrix of coefficients in (15) to the lower triangular form can be performed using the backward sweep formulas for the recurrence relation, setting the initial value of γ in the following form:
γ N = 1 + w N w N 1 w N + w N + 1 exp ( j 2 k N ( z N z N 1 ) )   ,
then, in the reverse cycle, the remaining values of γi are determined using the recursive relation:
γ i = 1 + ( w i w i + 1 ) γ i + 1 + 2 w i + 1 ( w i + w i + 1 ) γ i + 1 2 w i + 1 exp ( j 2 k i ( z i z i 1 ) ) ,     i = N 1 , 1 , 1 ¯   ,
which ultimately allows calculating the transmittance and reflectance for the zeroth layer:
{ t 0 A = 1 r 0 A = ( w 1 + w 0 ) γ 1 2 w 1 ( w 1 w 0 ) γ 1 2 w 1   .

2.2. Transfer Matrix Method

The reflectance spectrum of a homogeneous non-apodized FBG is modeled by the formulation of the resulting transmission matrix in the following form [12,13,19]:
T ( λ ) = [ cosh ( γ H ) j Δ β γ sinh ( γ H ) j κ γ sinh ( γ H ) j κ γ sinh ( γ H ) cosh ( γ H ) + j Δ β γ sinh ( γ H ) ]   ,
where H is the length of the FBG, and the parameters γ, κ, and Δβ are determined by the following formulas:
Δ n eff = Δ n 2 ,     n eff = n 0 + Δ n eff ,     κ = π Δ n eff λ ,   Δ β = 2 π n eff λ π Λ ,     γ = | κ | 2 Δ β 2   .  
The dependence of T(λ) on the wavelength λ in (20) is implicit and expressed in the terms of κ, Δβ, and γ, which are the functions of the wavelength.
The transfer matrix of the FBG with one discrete phase shift (Figure 1e) is formulated as the product of three transfer matrices, two of which (the first and the third ones) are the transfer matrices of homogeneous FBGs with the length of H1 and H2, respectively, and the second matrix describes the section of the optical fiber with the length h forming the following phase shift:
T π - FBG ( λ ) = T ( λ , H 1 ) × T φ ( λ , h ) × T ( λ , H 2 )   ,
where Tφ is the transfer matrix of the π-phase shift:
T φ ( λ , h ) = [ exp ( j 2 π n e f f λ h ) 0 0 exp ( j 2 π n e f f λ h ) ]   .
In (23), the length h corresponds to the phase shift of the grating: h = Λ/2.
The reflection coefficient at each wavelength is defined as the square of the ratio modulus of the resulting transfer matrix elements:
R ( λ ) = | T 1 , 0 T 0 , 0 | 2   .

3. Results and Discussion

3.1. Modeling Parameters of Homogeneous FBG

In order to verify the implementation of the layer sweep (LS) and transfer matrix (TM) methods, a homogeneous non-apodized FBG was modeled and the results were compared with the spectrum obtained using specialized software package Optiwave OptiGrating 4.2.2. The spectra simulated via the OptiGrating were taken as reference, since it is a commercial software widely used in industry.
To compare the simulation results, we chose a homogeneous non-apodized FBG 6 × 10−3 m long with a Bragg wavelength of 1550 nm, a refractive index n0 = 1.4682, and an induced refractive index Δn = 2 × 10−4. The chosen Bragg wavelength corresponds to the commonly used C-band wavelength range [20]. The refractive index of fiber core n0 is chosen according to the specifications of the common optical fiber SMF-28 at 1550 nm [21], while the induced refractive index lies within the typical range for the FBG [12,22,23]. The effective refractive index of the FBG model is defined as the average value of the refractive index over the fiber section containing the grating:
n eff = n 0 + Δ n 2   ,
which gives the average value of the effective refractive index of the FBG equal to neff ≈ 1.4683. The FBG period is determined by the Wulff–Bragg relation:
Λ = λ Br 2 n eff = λ Br 2 ( n 0 + Δ n / 2 )   ,
which is equal to Λ = 5.278213 × 10−7 m. The defined FBG parameters are listed in Table 1.

3.2. Results of Homogeneous FBG Modeling

An important note should be made regarding the choice of the partitioning interval used in the Layer Sweep method, since the modeling results using the LS method may significantly vary depending on the chosen partition interval. In order to establish the optimal interval, the following approach can be taken. A decreasing sequence of partition intervals {Δzn} is assigned so that
Δ z n = 1 2 n Δ z 0   .
We require that two spectra obtained with different partition intervals Δzn and Δzn+1 differ from each other by an amount not exceeding a small arbitrary pre-defined value ε. The spectrum obtained with the partition interval Δzn is denoted as a function depending on the wavelength and the partition interval Φ(λ, Δzn). Then, there is an integer number N such that for all n greater than N, and for any partition {λi} of an arbitrary interval [λBr − Δλ, λBr + Δλ], where i takes the values from 1 to M, the norm of the relative differences of the two spectra obtained for partitions zn and zn+1 differ by an amount not exceeding the given value ε.
N , M ; ε , λ Br , Δ λ : Δ z n = 2 n Δ z , n N ; λ i = λ Br Δ λ + i 2 Δ λ M , i = 0 , M ¯ i = 0 M Φ ( λ i , Δ z n + 1 ) Φ ( λ i , Δ z n ) Φ ( λ i , Δ z n )   ε   .
In other words, the partition interval is successively halved each time until the difference between the spectra obtained for two successive intervals becomes less than the pre-defined small value ε. However, as it will be shown further, the difference between the modeled and the reference spectra does not always monotonically decrease with the decrease of the partition interval, especially in non-homogeneous FBG structures, such as π-phase-shifted FBGs. In such cases, it is required to take into account additional criteria of the structure partitioning.
Figure 4 presents the results of the homogeneous FBG modeling with the parameters listed in the Section 3.1 and the partition interval Δz = 0.04 × Λ.
As it can be seen from Figure 4, the results of both methods were very close to the spectrum obtained using the proprietary software, although having slightly lower reflectance values in general. For further illustration, the deviations of spectra obtained using each method were plotted against the wavelength, which is presented in Figure 5. The root-mean-square deviations (RMSD) were calculated for both methods according to the equation
RMSD = 1 N i = 1 N ( R ( λ i ) R ref ( λ i ) ) 2 ,
where N = 201 is the number of wavelength values for which the reflectance R is estimated using LS or TM methods, and Rref is the reference values of FBG reflectance derived from the OptiGrating software. In this case, RMSDLS = 0.0043 and RMSDTM = 0.0024 for the LS and TM methods, respectively.
It should be noted that the deviation peaks of the LS and TM methods did not coincide besides the central wavelength, which is the evidence of differences in the mathematical modeling approaches. Nevertheless, the deviations of the modeled FBG spectra across the whole relevant wavelength range can be considered negligibly small for most applications, since the deviations do not exceed −40 dB, while the typical signal-to-noise ratio of the FBG interrogation device is between 30 to 40 dB [24]. Therefore, the results confirm the correctness of the modeling methods implementation.
In addition, the RMSD values of the LS method were obtained for various partition intervals Δz, and the corresponding dependence is presented in Figure 6.
The results shown in Figure 6 indicate that the RMSDLS values in the case of homogeneous FBG modeling decrease monotonically with the decrease of partition interval Δz. The data points were approximated by the second order polynomial (blue line in Figure 6) with the following coefficients: 1.184, 1.311 × 10−3, 2.366 × 10−3 and the norm of residuals equal to 3.087 × 10−5. Thus, with the vanishing partition interval of the homogeneous FBG model, the RMSDLS converges to the value slightly lower than the RMSDTM.

3.3. Results of a π-Phase-Shifted FBG Modeling

As it was mentioned before, the fiber Bragg structures comprising a phase shift (Figure 1e) are of particular interest for sensing applications since their spectral response includes an ultra-narrow transparency window, the wavelength position of which can be measured with higher accuracy than the shift of the entire FBG spectrum. Such a fiber Bragg structure can be represented as two identical FBGs located sequentially with the distance between them equal to the phase π-shift of their periodic refractive index variation. The modeling parameters of the sub-gratings are listed in Table 2. The simulated spectra of the π-phase-shifted FBG are shown in Figure 7, while the wavelength dependences of the spectra deviations from the OptiGrating data are presented in Figure 8. The partition interval was equal to Δz = 0.04 × Λ. In this case, RMSDLS = 3.89 × 10−3 and RMSDTM = 3.68 × 10−3 for the LS and TM methods, respectively.
As can be seen from the Figure 7 and Figure 8, the results of both methods are also very close to the spectrum obtained using the specialized software. The RMSD values were obtained for various partition intervals as well, which is demonstrated in Figure 9.
Unlike the previously considered example, the dependence of the RMSDLS values on the partition interval Δz in the case of π-phase-shifted FBG modeling is not monotonical and cannot be approximated with an acceptable precision. There are certain values of the partition interval, at which the RMSDLS significantly increases.
In order to establish the cause of such phenomenon, the modeled FBG refractive index profiles near the phase shift with different partition intervals Δz and the resulting spectra were compared to each other, as it is illustrated in Figure 10, where the intervals Δz = 0.13 × Λ and Δz = 0.12 × Λ are considered as an example.
As it can be seen from Figure 10, the spectrum obtained using the LS method with the partition interval Δz = 0.12 × Λ (Figure 10a, red solid line) has an asymmetric shape, which results in the increased RMSD in comparison with the spectrum acquired with the interval Δz = 0.13 × Λ despite the lower value of the partition interval. In this case, the displacement of transparency window of the phase-shifted FBG is comparable to its FWHM, which is not acceptable, since such structures are often used for high-resolution sensing [6]. If we compare the refractive index profiles of the grating (Figure 10b,c), we can see that in the case of the interval Δz = 0.12 × Λ (Figure 10b), the partitioning of the refractive index profile near the phase shift is asymmetrical relative to the center of the shift, while the partitioning with the interval Δz = 0.13 × Λ gives symmetrical profile relative to the phase shift center. The same effect was observed at all of the partition intervals indicated in Figure 9, at which the RMSDLS exceeded 0.02.
Therefore, based on the analysis of the refractive index profile of the grating model, one can conclude that the asymmetry of the refractive index profile partition near the phase shift relative to its center causes the asymmetry of the spectral response of the simulated grating, which results in the significant increase of the RMSDLS values at certain partition intervals, as it was reported in Figure 9. In the next sub-section of the paper, the authors propose a solution to eliminate the discussed negative effect in phase-shifted FBG simulation.

3.4. Amendment to Layer Sweep Method for Phase-Shifted FBG Modeling

In order to eliminate the spectral response asymmetry of the phase-shifted FBG modeled using the Layer Sweep method, it is necessary to ensure the symmetry of the partitioned refractive index profile of the FBG relative to the center of the phase shift. For that, the authors introduce an additional condition to the choice of the partition interval, according to which each FBG period must contain only integer and odd number of partition points. The simulation results obtained using the proposed approach are presented in Figure 11, where the value of Δz = (1/7) × Λ is taken as the partition interval.
The spectrum obtained using the LS method shown in Figure 11a closely replicates the spectrum derived from the OptiGrating software despite rather coarse partition of the refractive index profile depicted in Figure 11b. In the considered case, the RMSDLS = 0.0151, which is significantly lower even than the deviation demonstrated at Δz = 0.03 × Λ (see Figure 9). The RMSD values for various partition intervals corresponding to the introduced condition are shown in Figure 12.
The data points can be approximated by the fourth-order polynomial (blue line in Figure 12) with the following coefficients: −27.929, 9.476, −0.2491, 3.734 × 10−3, 3.667 × 10−3 and the norm of residuals equal to 2.123 × 10−5. Hence, the RMSD of the phase-shifted FBG spectrum modeled using the LS method with integer number of partition points in the FBG period decrease monotonically with the reduction of the partition interval. With the vanishing partition interval, the RMSDLS of the phase-shifted FBG model converges to the value slightly lower than the RMSDTM similarly to the homogeneous FBG modeling.

4. Conclusions

The results of the FBG modeling using both the Layer Sweep (LS) and Transfer Matrix (TM) methods are generally close to the ones obtained using the specialized software. However, in order to mitigate the deviation of the LS method spectrum from the reference one, it is necessary to decrease the partition interval of the refractive index profile.
It was also established that at certain values of the partition interval of the LS method, the deviations of the modeled phase-shifted FBG spectrum significantly increase due to its asymmetrical distortion. Based on the analysis of the refractive index profile of the grating model, it was concluded that the asymmetry of the refractive index profile partition near the phase shift relative to its center causes the asymmetry of the spectral response of the simulated grating.
A solution to eliminate the discussed negative effect in phase-shifted FBG simulation was proposed, which presupposes an additional condition to the choice of the partition interval, according to which each FBG period must contain only integer (odd) number of partition points.
It must be noted that the Layer Sweep method offers more flexibility in terms of the modeled fiber optic structure configuration, since it can simulate any refractive index profile of the fiber medium, such as two or more FBGs recorded one over the other (also known as the moiré recording) [4], which is not possible to model using the transfer matrix (TM) approach, while the latter offers benefits in terms of computational performance, especially when small partition interval is required to provide the desired accuracy of the LS method implementation.
Further research will be dedicated to the modeling of other types of fiber-optic structures, including the combined structures incorporating FBGs and Fabry–Perot interferometers, moiré-recorded FBGs, etc., as well as to the enhancement of the LS method computational performance. The latter can be achieved, for instance, by applying the variable partition interval.

Author Contributions

Conceptualization, A.S. and V.A.; methodology, A.S. and T.A.; software, A.S. and B.V.; validation, T.A., B.V., A.S., A.N. and Y.G.; formal analysis, V.A.; investigation, T.A. and B.V.; resources, O.M.; data curation, A.S.; writing—original draft preparation, T.A.; writing—review and editing, A.S., V.A., T.A. and A.N.; visualization, A.S., T.A. and Y.G.; supervision, O.M..; project administration, O.M. and A.S.; funding acquisition, O.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Ministry of Science and Higher Education as part of the “Priority 2030” program.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tseng, S.-M.; Chen, C.-L. Optical Fiber Fabry-Perot Sensors. Appl. Opt. 1988, 27, 547–551. [Google Scholar] [CrossRef] [PubMed]
  2. Hegde, G.; Asokan, S.; Hegde, G. Fiber Bragg Grating Sensors for Aerospace Applications: A Review. ISSS J. Micro Smart Syst. 2022, 11, 257–275. [Google Scholar] [CrossRef]
  3. Legoubin, S.; Fertein, E.; Douay, M.; Bernage, P.; Niay, P.; Bayon, F.; Georges, T. Formation of Moiré Grating in Core of Germanosilicate Fibre by Transverse Holographic Double Exposure Method. Electron. Lett. 1991, 27, 1945–1947. [Google Scholar] [CrossRef]
  4. Gillooly, A.M.; Dobb, H.; Zhang, L.; Bennion, I. Distributed Load Sensor by Use of a Chirped Moiré Fiber Bragg Grating. Appl. Opt. 2004, 43, 6454–6457. [Google Scholar] [CrossRef] [PubMed]
  5. Morozov, O.G.; Sakhabutdinov, A.J. Addressed fiber bragg structures in quasi-distributed microwave-photonic sensor systems. Comput. Opt. 2019, 43, 535–543. [Google Scholar] [CrossRef]
  6. Deepa, S.; Das, B. Interrogation Techniques for π-Phase-Shifted Fiber Bragg Grating Sensor: A Review. Sens. Actuators A Phys. 2020, 315, 112215. [Google Scholar] [CrossRef]
  7. Morozov, O.; Sakhabutdinov, A.; Anfinogentov, V.; Misbakhov, R.; Kuznetsov, A.; Agliullin, T. Multi-Addressed Fiber Bragg Structures for Microwave-Photonic Sensor Systems. Sensors 2020, 20, 2693. [Google Scholar] [CrossRef] [PubMed]
  8. Ma, W.; Wang, R.; Rong, Q.; Shao, Z.; Zhang, W.; Guo, T.; Wang, J.; Qiao, X. CO2 Gas Sensing Using Optical Fiber Fabry–Perot Interferometer Based on Polyethyleneimine/Poly(Vinyl Alcohol) Coating. IEEE Photonics J. 2017, 9, 2700053. [Google Scholar] [CrossRef]
  9. Bremer, K.; Lewis, E.; Leen, G.; Moss, B.; Lochmann, S.; Mueller, I.A.R. Feedback Stabilized Interrogation Technique for EFPI/FBG Hybrid Fiber-Optic Pressure and Temperature Sensors. IEEE Sens. J. 2012, 12, 133–138. [Google Scholar] [CrossRef]
  10. Tosi, D.; Poeggel, S.; Iordachita, I.; Schena, E. Fiber Optic Sensors for Biomedical Applications. In Opto-Mechanical Fiber Optic Sensors; Alemohammad, H., Ed.; Butterworth-Heinemann: Oxford, UK, 2018; pp. 301–333. ISBN 978-0-12-803131-5. [Google Scholar]
  11. Skaar, J.; Wang, L.; Erdogan, T. On the Synthesis of Fiber Bragg Gratings by Layer Peeling. IEEE J. Quantum Electron. 2001, 37, 165–173. [Google Scholar] [CrossRef]
  12. Erdogan, T. Fiber Grating Spectra. J. Light. Technol. 1997, 15, 1277–1294. [Google Scholar] [CrossRef]
  13. Ikhlef, A.; Hedara, R.; Chikh-Bled, M. Uniform Fiber Bragg Grating Modeling and Simulation Used Matrix Transfer Method. Int. J. Comput. Sci. Issues (IJCSI) 2012, 9, 7. [Google Scholar]
  14. Tai, H. Theory of Fiber Optical Bragg Grating: Revisited. In Proceedings of the Optical Modeling and Performance Predictions, San Diego, CA, USA, 22 January 2004; SPIE: Bellingham, WA, USA, 2004; Volume 5178, pp. 131–138. [Google Scholar]
  15. Khalid, K.S.; Zafrullah, M.; Bilal, S.M.; Mirza, M.A. Simulation and Analysis of Gaussian Apodized Fiber Bragg Grating Strain Sensor. J. Opt. Technol. 2012, 79, 667–673. [Google Scholar] [CrossRef]
  16. Hamarsheh, M.N.; Abdullah, M.K. Analysis of Fiber Bragg Gratings Apodized with Linearly Approximated Segmented Gaussian Function. In Proceedings of the 9th Asia-Pacific Conference on Communications (IEEE Cat. No.03EX732), Penang, Malaysia, 21–24 September 2003; Volume 1, pp. 163–167. [Google Scholar]
  17. Halder, P.; Deyasi, A. Calculating Reflectivity of Fiber Bragg Grating for Different Apodization Techniques. MMC_A 2017, 90, 116–124. [Google Scholar] [CrossRef]
  18. Capmany, J.; Sales, S.; Muriel, M.A.; Rubio, J.J. Novel Layer Peeling Algorithm for the Synthesis of Fiber Bragg Gratings Yielding Smoother Fabrication Profiles. In Proceedings of the Optical Fiber Communication Conference, 2004, OFC 2004, Los Angeles, CA, USA, 22 February 2004; Volume 1, p. 77. [Google Scholar]
  19. Fazzi, L.; Klyukin, D.; Groves, R.M. Transfer Matrix Method for Fundamental LP01 Core Mode Coupling in a Tilted FBG Sensor. AIP Conf. Proc. 2020, 2293, 200010. [Google Scholar] [CrossRef]
  20. Tosi, D. Review and Analysis of Peak Tracking Techniques for Fiber Bragg Grating Sensors. Sensors 2017, 17, 2368. [Google Scholar] [CrossRef] [PubMed]
  21. SMF-28® Ultra Optical Fibers|Single Mode Optical Fiber|Corning. Available online: https://www.corning.com/optical-communications/worldwide/en/home/products/fiber/optical-fiber-products/smf-28-ultra.html (accessed on 30 January 2023).
  22. Ashry, I.; Elrashidi, A.; Mahros, A.; Alhaddad, M.; Elleithy, K. Investigating the Performance of Apodized Fiber Bragg Gratings for Sensing Applications. In Proceedings of the 2014 Zone 1 Conference of the American Society for Engineering Education, Bridgeport, CT, USA, 3–5 April 2014; pp. 1–5. [Google Scholar]
  23. Prashar, S.; Engles, D.; Malik, S.S.; Vohra, R. Investigative Study of Transmission Spectra of FBG at Varying Induced Index & Grating Length. Procedia Eng. 2012, 38, 3031–3036. [Google Scholar] [CrossRef]
  24. Pustakhod, D.; Kleijn, E.; Williams, K.; Leijtens, X. High-Resolution AWG-Based Fiber Bragg Grating Interrogator. IEEE Photonics Technol. Lett. 2016, 28, 2203–2206. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Several types of layered structures formed in the optical fiber core: (a) Fabry–Perot resonator (FPR); (b) fiber Bragg grating (FBG); (c) two sequential FBGs; (d) two FBGs recorded one over the other; (e) FBG with π phase shift; (f) FBG with two π phase shifts; (g) FBG and FPR combination.
Figure 1. Several types of layered structures formed in the optical fiber core: (a) Fabry–Perot resonator (FPR); (b) fiber Bragg grating (FBG); (c) two sequential FBGs; (d) two FBGs recorded one over the other; (e) FBG with π phase shift; (f) FBG with two π phase shifts; (g) FBG and FPR combination.
Algorithms 16 00101 g001
Figure 2. Scheme of a layered structure of a substance with an arbitrary dependence of the refractive index over its thickness.
Figure 2. Scheme of a layered structure of a substance with an arbitrary dependence of the refractive index over its thickness.
Algorithms 16 00101 g002
Figure 3. One-dimensional model of a layered structure.
Figure 3. One-dimensional model of a layered structure.
Algorithms 16 00101 g003
Figure 4. Simulated spectra of the homogeneous FBG obtained using the Layer Sweep method (red solid line), Transfer Matrix method (green solid line), OptiGrating software (blue dotted line), and amplified spectra near the reflectance peak (insert).
Figure 4. Simulated spectra of the homogeneous FBG obtained using the Layer Sweep method (red solid line), Transfer Matrix method (green solid line), OptiGrating software (blue dotted line), and amplified spectra near the reflectance peak (insert).
Algorithms 16 00101 g004
Figure 5. Deviations of the homogeneous FBG spectra from the reference (OptiGrating) spectrum (left axis): deviation of the Layer Sweep method (red line), deviation of the Transfer Matrix method (green line), reference OptiGrating spectrum (black line, right axis).
Figure 5. Deviations of the homogeneous FBG spectra from the reference (OptiGrating) spectrum (left axis): deviation of the Layer Sweep method (red line), deviation of the Transfer Matrix method (green line), reference OptiGrating spectrum (black line, right axis).
Algorithms 16 00101 g005
Figure 6. Values of the root-mean-square deviation (RMSD) of the homogeneous FBG spectrum modeled using the LS method (blue dots) depending on the partition interval Δz, and their approximation (blue solid line).
Figure 6. Values of the root-mean-square deviation (RMSD) of the homogeneous FBG spectrum modeled using the LS method (blue dots) depending on the partition interval Δz, and their approximation (blue solid line).
Algorithms 16 00101 g006
Figure 7. Simulated spectra of the π-phase-shifted FBG obtained using the Layer Sweep method (red solid line), Transfer Matrix method (green solid line), OptiGrating software (blue dotted line), and amplified spectra near the reflectance peak (insert).
Figure 7. Simulated spectra of the π-phase-shifted FBG obtained using the Layer Sweep method (red solid line), Transfer Matrix method (green solid line), OptiGrating software (blue dotted line), and amplified spectra near the reflectance peak (insert).
Algorithms 16 00101 g007
Figure 8. Deviations of the π-phase-shifted FBG spectra from the reference (OptiGrating) spectrum (left axis): deviation of the Layer Sweep method (red line), deviation of the Transfer Matrix method (green line); and reference OptiGrating spectrum (black line, right axis).
Figure 8. Deviations of the π-phase-shifted FBG spectra from the reference (OptiGrating) spectrum (left axis): deviation of the Layer Sweep method (red line), deviation of the Transfer Matrix method (green line); and reference OptiGrating spectrum (black line, right axis).
Algorithms 16 00101 g008
Figure 9. Values of the root-mean-square deviation (RMSD) of the π-phase-shifted FBG spectrum modeled using the LS method (blue dots) depending on the partition interval Δz.
Figure 9. Values of the root-mean-square deviation (RMSD) of the π-phase-shifted FBG spectrum modeled using the LS method (blue dots) depending on the partition interval Δz.
Algorithms 16 00101 g009
Figure 10. Comparison of modeling results of the π-phase-shifted FBG using the Layer Sweep method at different partition intervals: (a) resulting spectra: at Δz = 0.12 × Λ (red solid line), at Δz = 0.13 × Λ (green solid line), OptiGrating spectrum (blue dotted line); (b) refractive index profile of the grating at Δz = 0.12 × Λ; and (c) refractive index profile of the grating at Δz = 0.13 × Λ.
Figure 10. Comparison of modeling results of the π-phase-shifted FBG using the Layer Sweep method at different partition intervals: (a) resulting spectra: at Δz = 0.12 × Λ (red solid line), at Δz = 0.13 × Λ (green solid line), OptiGrating spectrum (blue dotted line); (b) refractive index profile of the grating at Δz = 0.12 × Λ; and (c) refractive index profile of the grating at Δz = 0.13 × Λ.
Algorithms 16 00101 g010
Figure 11. Modeling results of the π-phase-shifted FBG using the Layer Sweep method with the partition interval of Δz = (1/7) × Λ: (a) resulting spectra: of LS method (red solid line), OptiGrating spectrum (blue dotted line); (b) and refractive index profile of the modeled grating.
Figure 11. Modeling results of the π-phase-shifted FBG using the Layer Sweep method with the partition interval of Δz = (1/7) × Λ: (a) resulting spectra: of LS method (red solid line), OptiGrating spectrum (blue dotted line); (b) and refractive index profile of the modeled grating.
Algorithms 16 00101 g011
Figure 12. Values of the root-mean-square deviation (RMSD) of the phase-shifted FBG spectrum modeled using the LS method with odd number of partition points in the FBG period (blue dots) depending on the partition interval Δz, and their approximation (blue solid line).
Figure 12. Values of the root-mean-square deviation (RMSD) of the phase-shifted FBG spectrum modeled using the LS method with odd number of partition points in the FBG period (blue dots) depending on the partition interval Δz, and their approximation (blue solid line).
Algorithms 16 00101 g012
Table 1. Modeling parameters of the homogeneous non-apodized FBG.
Table 1. Modeling parameters of the homogeneous non-apodized FBG.
ParameterValue
Grating length H, m6 × 10−3
Refractive index of the fiber n0, RIU1.4682
Induced refractive index Δn, RIU2 × 10−4
Grating period Λ, m5.278213 × 10−7
Table 2. Modeling parameters of the π-phase-shifted FBG.
Table 2. Modeling parameters of the π-phase-shifted FBG.
ParameterValue
Sub-grating length H, m5 × 10−3
Refractive index of the fiber n0, RIU1.4682
Induced refractive index Δn, RIU2 × 10−4
Sub-grating period Λ, m5.278213 × 10−7
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Agliullin, T.; Anfinogentov, V.; Morozov, O.; Sakhabutdinov, A.; Valeev, B.; Niyazgulyeva, A.; Garovov, Y. Comparative Analysis of the Methods for Fiber Bragg Structures Spectrum Modeling. Algorithms 2023, 16, 101. https://doi.org/10.3390/a16020101

AMA Style

Agliullin T, Anfinogentov V, Morozov O, Sakhabutdinov A, Valeev B, Niyazgulyeva A, Garovov Y. Comparative Analysis of the Methods for Fiber Bragg Structures Spectrum Modeling. Algorithms. 2023; 16(2):101. https://doi.org/10.3390/a16020101

Chicago/Turabian Style

Agliullin, Timur, Vladimir Anfinogentov, Oleg Morozov, Airat Sakhabutdinov, Bulat Valeev, Ayna Niyazgulyeva, and Yagmyrguly Garovov. 2023. "Comparative Analysis of the Methods for Fiber Bragg Structures Spectrum Modeling" Algorithms 16, no. 2: 101. https://doi.org/10.3390/a16020101

APA Style

Agliullin, T., Anfinogentov, V., Morozov, O., Sakhabutdinov, A., Valeev, B., Niyazgulyeva, A., & Garovov, Y. (2023). Comparative Analysis of the Methods for Fiber Bragg Structures Spectrum Modeling. Algorithms, 16(2), 101. https://doi.org/10.3390/a16020101

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