Next Article in Journal
Intelligence-Aware Batch Processing for TMA with Bearings-Only Measurements
Previous Article in Journal
A Smart Capacitive Sensor Skin with Embedded Data Quality Indication for Enhanced Safety in Human–Robot Interaction
Previous Article in Special Issue
Wavelength-Tunable, Ultra-Broadband, Biconical, Long-Period Fiber Grating Mode Converter Based on the Dual-Resonance Effect
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sensors
Volume 21
Issue 21
10.3390/s21217212
sensors-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation of the Transmission Spectrum of Long-Period Fiber Gratings Structures with a Propagating Acoustic Shock Front

by
Oleg V. Ivanov
1,2,*,
Paulo Caldas
3,4 and
Gaspar Rego
3,4
1
Ulyanovsk Branch of Kotel’nikov Institute of Radio Engineering and Electronics of Russian Academy of Sciences, Ulitsa Goncharova 48, 432071 Ulyanovsk, Russia
2
Ulyanovsk State University, S.P. Kapitsa Research Institute of Technology, Ulitsa L. Tolstogo 42, 432017 Ulyanovsk, Russia
3
ProMetheus-Research Unit in Materials, Energy and Environment for Sustainability, Instituto Politécnico de Viana do Castelo, 4900-347 Viana do Castelo, Portugal
4
Center for Applied Photonics, Institute for Systems and Computer Engineering, Technology and Science (INESC TEC), Rua Dr. Roberto Frias, 4200-465 Porto, Portugal
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(21), 7212; https://doi.org/10.3390/s21217212
Submission received: 28 September 2021 / Revised: 22 October 2021 / Accepted: 26 October 2021 / Published: 29 October 2021
(This article belongs to the Special Issue 25 Years of Long-Period Fiber Gratings)
Browse Figures
Versions Notes

Abstract

:
In this paper, we investigate modification of transmission spectra of long-period fiber grating structures with an acoustic shock front propagating along the fiber. We simulate transmission through inhomogeneous long-period fiber gratings, π-shift and reflective π-shift gratings deformed by an acoustic shock front. Coupled mode equations describing interaction of co-propagating modes in a long-period fiber grating structures with inhomogeneous deformation are used for the simulation. Two types of apodization are considered for the grating modulation amplitude, such as uniform and raised-cosine. We demonstrate how the transmission spectrum is produced by interference between the core and cladding modes coupled at several parts of the gratings having different periods. For the π-shift long-period fiber grating having split spectral notch, the gap between the two dips becomes several times wider in the grating with the acoustic wave front than the gap in the unstrained grating. The behavior of reflective long-period fiber gratings depends on the magnitude of the phase shift near the reflective surface: an additional dip is formed in the 0-shift grating and the short-wavelength dip disappears in the π-shift grating.

1. Introduction

In the past two to three decades, fiber gratings have been intensively studied for applications in optical communications and sensing [1,2,3]. In the later domain, fiber Bragg gratings have also been used for detection of acoustic and ultrasound waves [4,5,6]. More recently, the scientific community focused their attention on the development of sensors able to detect deformations caused by shock waves working at frequencies from hundred kHz up to several GHz; examples of such acoustic shock waves are explosive detonation waves. In this context, fiber Bragg gratings have already proved their potential and several papers have been published in literature [7,8,9,10,11], while the interaction of a shock wave front with a fiber Bragg grating was also studied [12]. It would be of interest to investigate if similar methods could be used for long-period fiber gratings (LPFGs) and to analyze the difference between acoustic effects in short and long-period gratings. However, the use LPFGs for ultrafast applications is discussed only in scarce works related to the detection of ultrasound waves [13,14,15]. In general, the mechanism of detection relies on the effect of acoustic pressure on the grating that can be further amplified if the LPFG is bent. In fact, different configurations, included interferometric ones, can be employed for detection. So far, no theoretical work has been published regarding the effect of shock waves on LPFG transmission spectra.
Conversely, the interaction of ultrasound with grating sensors can be used to identify leakage of medical gas delivery systems in health care facilities. This issue has a tremendous financial impact, which is even more demanding in the current pandemic situation.
For the case of step-changed acoustic fronts, calculating transmission spectra is similar to the problem of inhomogeneous step-changed LPFGs where some part of the grating is strained, while the strain sensitivity of LPFGs has been analyzed in Ref. [16]; meanwhile, the properties of inhomogeneous LPFGs with apodization, chirp, and phase shift have been investigated in Refs. [17,18,19,20,21], and the step can be produced by concatenating two LPFGs with different periods [22,23], modulation amplitudes [24,25], and some phase shift between them [26]. Inhomogeneous gratings can be formed by heating part of the fiber grating and are used as gain-flattening filters [27].
The spectral response of LPFGs can be controlled by adding a phase shift in the middle of the grating [19,28,29], while the phase shift is formed by a section of fiber without grating, mechanically, or by the changes in its refractive index [30,31]. In contrast to homogeneous LPFGs, a pass band is created in the resonance notches of the phase-shifted LPFGs as a result of interference between two parts of the grating. The phase-shifted LPFGs have been used for sensing, in particular, for the measurements of temperature and refractive index [32], twist [33], strain [34], bending [23], as gain-flattening filters in erbium-doped fiber amplifiers [35,36], and optical differentiators [37]. Unique response characteristics of π-phase-shifted fiber Bragg gratings for ultrasonic detection have been described [38]. Higher sensitivity of the phase-shifted LPFG may be also advantageous when it is employed for the detection of acoustics waves.
The standard LPFG works in transmission and requires access to both end-sides of the fiber: one brings the light from a source and the other returns the remaining light to a photo-detector. In terms of sensing, it is more convenient to create a sensing head having access to a single end of the fiber. This can be attained by using a reflective surface at the end of the optical fiber with the LPFG in a self-interference Michelson configuration [39,40,41,42,43,44,45,46], and the size of the sensing head is halved in this case. An arbitrary phase shift can be introduced in the center of the grating, if a gap is added between the LPFG and the mirror. We can also expect different behavior of such a structure under the influence of an acoustic wave front.
In this work, we investigate modification of transmission spectra of long-period fiber grating structures with an acoustic shock front propagating along the fiber. We simulate transmission through homogeneous long-period fiber gratings, π-shift, and reflective gratings with propagating acoustic shock front producing non-uniform deformation along the optical fiber. Coupled mode equations describing interaction of co-propagating modes in a long-period fiber grating structures with inhomogeneous deformation are written and used for the simulation. We consider the influence of smoothness of the acoustic front and the magnitude of the phase shift on the spectra of various types of LPFG structures.

2. Long-Period Fiber Grating Structures

The long-period fiber grating has periodical modulation of its refractive index with period between 100 and 1000 μm (Figure 1a). The transmission spectrum of a LPFG is a function with several dips centered at resonance wavelengths, which are determined by the grating period and propagation constants of the core and cladding modes. The homogeneous LPFGs have symmetric dips with some sidelobes.
The phase-shifted LPFG has a section with modified fiber structure, where the grating has an additional space between the grating lines or changed propagation constants of the core and cladding modes (Figure 1b). Each spectral notch of the phase-shifted LPFG is split into two dips. The influence of the phase shift on the spectrum is strongest when it is in the center of grating. The relation between the amplitudes of the dips is determined by the magnitude of the phase shift. The amplitudes are equal for π phase shift; here, we consider only π-phase-shifted LPFGs in the center of the grating.
The reflective LPFG has a reflective coating at the end face of the optical fiber near the grating. If the grating ends exactly at the end face, no phase shift is introduced for the propagating modes (Figure 1c). Such a grating works as a single homogeneous LPFG with doubled length. However, its behavior deformed by acoustic wave differs from the homogeneous LPFG, because the acoustic wave covers the unwrapped grating starting from its center, not from one side.
If there is some space between the grating and the end face, a phase shift is formed (Figure 1d). We obtain a grating that is equivalent in unwrapped form to a phase-shifted grating or a cascaded grating, and again, when an acoustic wave propagates through the grating, it starts from the center of the unwrapped grating.
We assume that a longitudinal compressive shock front of acoustic shock wave propagates along the fiber and deforms the grating starting from z = L (Figure 1). The speed of the acoustic shock wave is defined by the speed of sound in the silica fiber. When the shock front propagates along the fiber with the LPFG, the grating is divided into two parts; one part is deformed and its period is decreased and the refractive index is changed due to strain-optic effect, while, in a homogeneous grating, this would shift the resonance wavelengths of the spectral peaks. The shift can be positive or negative for different cladding modes depending on specific fiber structure parameters [47,48]. The other part is unaffected and its period is not changed.
The front of the shock wave may have different profiles, which is determined by the time profile of the initial impact and by diffusion of the shock front due to dispersion of the speed of sound in the fiber and surrounding material. Assuming that the initial shock front is described by the step function, it can be represented as a Fourier transform of waves with different frequencies. These waves propagating in a dispersive medium have different speeds and arrive to the grating in different moments. The resulting sum of the waves is some smoothed stepwise function. The magnitude of smoothing is determined by the dispersion of the medium and the distance from the source to the grating. We consider only the effect of the leading edge of compression of the shock wave, which is much sharper compared to the trailing edge of expansion. Slow expansion has space distribution larger than the grating length and covers the whole grating, which acts as a homogeneous structure in this case.
The strain can be strong enough so that the wavelength shift is large and resonances of cladding mode with adjacent mode numbers can overlap. This is an alternative kind of interference, because cladding modes with different mode numbers do not interfere themselves; indeed, it would take place only interference between them when they couple back to the core mode. Another situation is possible when the LPFG excites modes of different types [49]: the overlap of resonances of different types of modes. We will not consider this situation in this paper.
The refractive index modulation in a fiber with a LPFG can be written as follows:
n ( r , z ) = n 0 ( r ) 1 + g ( r ) σ ( z ) 1 + cos ( 2 π z / Λ 0 )
where n 0 ( r ) is the refractive index of the fiber without the grating, σ ( z ) is a slowly varying envelope of modulation amplitude of the induced refractive index, g ( r ) is the function describing radial dependence refractive index modulation (in photoinduced gratings, g ( r ) = g co for r ≤ r co and g ( r ) = 0 for r > r co ), and Λ 0 is the grating period.
To simulate the process of propagation of acoustic shock front through the LPFG, we assume that the amplitude of the front profile is described by the following function, which changes the strain from one level to another with controllable smoothness defined by the parameter η :
s ( z ) = s 0 2 1 + tanh z − z 0 η
Here, z 0 is the position of the shock front. Figure 2 illustrates the dependence of the amplitude of the shock front on the longitudinal coordinate. The shock front compresses the fiber, and each point of the fiber is displaced by u ( z ) in − z direction:
u ( z ) = ∫ − ∞ z s ( z )   d z
The refractive index modulation in the deformed fiber takes the form
n ( r , z ) = n 0 ( r ) 1 + g ( r ) σ z + u ( z ) 1 + cos 2 π ( z + u ( z ) ) / Λ 0
The period of the deformed grating is determined by the function s ( z ) and decreases a result of strain: Λ ( z ) = Λ 0 ( 1 − s ( z ) ) .
In addition to displacement of the fiber grating, the refractive index of the fiber is modified due to strain-optic effect by
Δ n s ( r , z ) = n 0 ( r ) p a ( r ) s ( z )
where p a is the effective strain-optic coefficient p a = n 2 2 p 12 − ν ( p 11 − p 12 ) , p 12 and p 11 are the strain-optic coefficients, and ν is the Poisson’s ratio. In contrast to fiber Bragg gratings, the dependence of the refractive index and the strain-optic coefficient in the right part of the Equation (5) on radius cannot be neglected. Otherwise, the core and cladding modes would have identical sensitivity to strain, and their contributions in the wavelength shift, as we show below, would cancel each other.
The dielectric permittivity and its change in the deformed and strained fiber can be expressed from Equation (4) as
ε ( r , z ) = n 0 ( r ) 2 1 + 2 p a s ( z ) 1 + 2 g ( r ) σ ( z ) 1 + cos 2 π ( z + u ( z ) ) / Λ 0 Δ ε ( r , z ) = 2 n 0 ( r ) 2 p a s ( z ) + g ( r ) σ ( z ) 1 + cos 2 π ( z + u ( z ) ) / Λ 0
Here we have neglected the terms of the second order of smallness in the change of dielectric permittivity.

3. Coupled Mode Equations

The long-period fiber gratings couple modes propagating in the same direction. We consider the case of one core mode that is coupled to the forward propagating cladding modes. The interaction with backward propagating modes may occur only at the endface of the fiber for reflective LPFGs. The coupling coefficient between two co-propagating modes is defined as the following overlap integral over the fiber cross-section:
K μ ν ( z ) = ω ε 0 4 ∫ ∞ E μ * ( r , φ ) Δ ε ( r , z ) E ν ( r , φ ) d S
where Δ ε ( r , z ) is the perturbation of dielectric permittivity in the fiber as a result of grating inscription and propagating acoustic wave, E μ ( r , φ ) and E ν ( r , φ ) are the transverse distributions of electric field of HE1μ and HE1ν modes, respectively, ω is the frequency, and ε 0 is the dielectric constant of vacuum.
The coupled mode equations that describe the amplitudes of modes can be written as
d A μ d z =   i ∑ ν K μ ν A ν exp i ( β ν − β μ ) z
where β μ and β μ are the propagation constants of the μ -th and ν -th modes. We represent the coupling coefficient, taking into account that the perturbation of dielectric permittivity depends on coordinate z , as
K μ ν ( z ) = γ μ ν ( z ) + κ μ ν ( z ) 1 + cos 2 π ( z + u ( z ) ) / Λ 0
The two coupling constants γ μ ν ( z ) and κ μ ν ( z ) describe mode coupling due to strain and modulation in the grating, respectively:
γ μ ν ( z ) = ω ε 0 2 p a s ( z ) ∫ 0 r cl n 0 ( r ) 2 E μ * ( r ) E ν ( r ) 2 π r d r
κ μ ν ( z ) = ω ε 0 2 σ ( z ) ∫ 0 r cl n 0 ( r ) 2 g ( r ) E μ * ( r ) E ν ( r ) 2 π r d r
We have assumed here that only HE1μ modes without azimuthal structure interact, because the fiber is straight and has no birefringence and the grating is radially symmetric. The coupled mode equations are obtained by ignoring quickly oscillating terms and describe slowly varying amplitudes of two co-propagation modes:
d A d z =   i γ co - co ( z ) + κ co - co ( z ) A + i 2   κ co - ν ( z ) B exp ( − i 2 δ z ) d B d z =   i γ ν - ν ( z ) + κ ν - ν ( z ) B + i 2   κ co - ν ( z ) A exp ( i 2 δ z )
where A and B denote the amplitudes of co-propagating modes HE11 and HE1ν, respectively, and the detuning parameter is
δ ( z ) = 1 2 β co − β ν − 1 + u ( z ) z π Λ 0
For obtaining Equation (12), we use the equality of coupling constants κ co - ν = κ ν − co .
The resonance condition for mode coupling is
δ + γ co - co − γ ν - ν + κ co - co − κ ν - ν / 2 = 0
The effective refractive indices n ( eff ) = β / k 0 depend on wavelength ( k 0 = ω / c is the wavenumber in vacuum) and can be approximately written in a range of tens of nanometer around one of the cladding mode dips as n ( eff ) ( λ 0 ) + Δ λ d n ( eff ) / d λ . Then, we write for the detuning parameter
δ = π Δ λ λ 0 − n co ( eff ) − n ν ( eff ) λ 0 + d n co ( eff ) d λ − d n ν ( eff ) d λ − u ( z ) z
For simulation, we use parameters of SMF-28 fiber at wavelength λ 0 = 1550 nm: core radius r co = 4.2 μm, cladding refractive index n cl = 1.4440 , core refractive index n co = 1.4494 , and the effective refractive index of the core mode n co ( eff ) = 1.4465 . Integrating over the fiber cross-section, we obtain the following approximations for the coupling constants:
γ co - co ( z ) = 1 + α co - co n cl p a s ( z ) k 0 γ ν - ν ( z ) = 1 + α ν - ν n cl p a s ( z ) k 0 κ co - co ( z ) = 0.78   n cl σ ( z ) k 0   κ ν - ν ( z ) = 0 κ co - ν ( z ) = χ ν   n cl σ ( z ) k 0
The coefficient 0.78 in the self-coupling constant κ co - co is obtained as the self-overlap intergal of the core mode over the fiber core. The overlap integrals γ co - co and γ ν - ν are equal to n cl p a s ( z ) k 0 with some small corrections α co - co and α ν - ν , which are caused by several effects induced by deformation and strain in the fiber [12]. The cross-coupling coefficient κ co - ν depends on mode number via χ ν : χ ν = 0.030 , 0.055, 0.076, 0.092 for ν = 2 , 3 , 4 , 5 , respectively (see Appendix A).
As follows from Equation (14), when the fiber is strained, the wavelength shift of a spectral dip is determined by the difference of self-coupling constants depending on strain:
γ co - co − γ ν - ν = α co - co − α   ν - ν n cl p a s ( z ) k 0
We take α co - co − α   ν - ν = − 0.01 , which is close to the first and largest overlap integral for the compression of the core refractive index profile under the deformation (see Equations 21 and 35 in [12]). As compared to fiber Bragg gratings, the effect of changing phase of the strained grating is smaller, because the period is smaller by several orders of magnitude. The mechanical displacement of the long-period grating end under strain is of the same order as the grating period.

4. Simulation of Uniform LPFG

We begin simulation with a uniform grating having constant modulation amplitude of induced refractive index along the length of the fiber. Let the grating has a period of 430 μm. Then the resonance of the fifth mode HE15 is centered at 1550 nm. At this wavelength, we can find the effective refractive index of the fifth cladding mode n co ( eff ) = 1.44306 . We choose such amplitude of refractive index modulation σ = 1.48 ⋅ 10 − 4 that the transmission in the center of resonance is equal to zero. This corresponds to the coupling constant between the core the fifth cladding mode κ co - 5 = 1.6 ⋅ 10 − 4 μm–1. All parameters that we used in the simulation are listed in Table 1. The numerical simulation was made in Matlab (R2012b).
The transmission spectra of a regular LPFG are shown in Figure 3 for different values of strain in stepwise acoustic front ( η = 0 ) located in the center of the grating z 0 = L / 2 . The spectrum of the uniform grating without strain is shown in the lower plot. This spectral shape is a dip with minimum transmission in the center and represents the standard LPFG notch.
The acoustic front propagating through the grating divides the grating into two parts with different parameters. Such a grating is equivalent to a step-changed grating [22], and these parts have different resonance wavelengths. The core and cladding modes coupled at the two sections interfere and form an interference spectrum in transmission. It is symmetric, because the step is in the center of the grating. The spectrum shifts to shorter wavelengths with an increasing amplitude of the acoustic front due to decreasing average period of the grating. We can also see that the presence of the shock front in the fiber leads, first, to formation of a broader notch and, second, to splitting of the dip. For high strains, the spectrum completely splits into two dips. Each dip corresponds to its own homogeneous part of the grating.
Figure 4 shows the influence of smoothness parameter η on grating spectrum. The plots are presented for three magnitudes of the parameter η = 0 , 12.5 to 25 mm for the case of acoustic front in the center of the grating and strain s 0 = 18 mstrain. It can be seen that two separate peaks merge in one, when the smoothness becomes comparable with the grating length.
Figure 5 shows evolution of grating spectra with z 0 moving through the LPFG for three values of smoothness η : 0, 10 mm, and 20 mm ( s 0 = 18 mstrain). The dip transits continuously and smoothly from one wavelength to another, when the smoothness is 20 mm. In this case, interference effects are not observed. For η ≈ 10 mm, the shifting process in wavelength is accompanied by the appearance of a three-prong fork shape and splitting of the spectrum, when the front is in the middle of the grating. For step-changed shock front with η = 0 , one dip splits into several smaller disappearing dips, and a new dip appears at a shifted wavelength with multiple interference fringes around the main dips.
Let us consider the dependence of spectral characteristics of the LPFG with an acoustic front in its center. When the length increases, the transmission loss grows overall. The spectrum is wide for short grating lengths and forms a band centered at 1530 nm (as illustrated in Figure 6). Then, for lengths above 40 mm, two dips are formed at the wavelength corresponding to the strained and non-strained parts of the grating. Further increasing length results in more narrow dips and the appearance of higher sidelobes. This behavior is similar to the behavior of spectra of uniform long-period gratings. Something like a spectral standing wave with increasing number of sidelobes is formed between two main dips.
In practice, some apodization can be applied to gratings to produce a smoother spectrum. This is described in our case by the function σ ( z ) . As an example, we calculated the transmission spectrum of the LPFG with raised-cosine apodization (Figure 7). It is similar to the spectrum of the uniform grating but has no spectral sidelobes on long-wavelength side and dipper sidelobes on short-wavelength side. For the stepwise wave front, the transition of the dip from the initial wavelength to the shifted wavelength is asymmetric and irregular, as compared to uniform grating. There is also an effect of loss concentration at z 0 / L = 0.23 and wavelength 1497 nm. For η = 25 mm, one main dip is formed closer to the center wavelength, and it is accompanied by smaller dips on the short-wavelength side. Compared to Figure 5, no interference fringes are seen on the long-wavelength side of the picture.

5. Simulation of π-Shift LPFG

The structure of π-shift LPFG is shown in Figure 1b. We assume that the phase shift is in the center of the grating. For the simulation of transmission, we divide the whole structure into two parts. Then, we obtain solution of the coupled mode equations in the first part, which gives us the amplitudes of the core and cladding modes. Next, we add a phase shift to one mode and solve the coupled mode equations in the second part of the grating. Finally, we find transmission coefficient from the amplitudes of the modes at the end of the grating.
The spectrum of a π-shift LPFG has a transmission maximum in the center of the grating with two rejection bands on the sides, while the maximum is the result of destructive interference is caused by the introduction of the phase shift at the center. The same phase shift produces two rejection bands that are formed by constructive interference at the side wavelengths. The maximum between the two rejection bands is largest when the phase shift is at the center of the grating; if the amplitude of modulation in the grating is the same, the two rejection bands are shallower, compared to the single band the uniform grating.
When the acoustic front propagates through the phase-shifted LPFG, the fiber grating is divided at first into three sections: unstrained before phase shift, unstrained after phase, strained after phase. After the wave passes the phase shift, the three sections are the following: unstrained before phase shift, strained before phase, and strained after phase shift.
Evolution of the LPFG’s spectrum with the front moving along the grating is presented in Figure 8. We can see from the figure that the acoustic front shifts the two dips of the π-shift LPFG to the short-wavelength side of the spectrum. For small strain, the shift is accompanied by flow of intensity of one dip into another. With increasing strain, the gap between the dips becomes several times wider than for unstrained grating, while for high strain, the transition of the dips becomes discontinuous; they dissolve at the initial wavelengths and appear at the new wavelengths. The sensitivity of the π-shift LPFGs to the acoustic wave is similar to that of the uniform LPFGs; at the same time, the regular spectrum with a continuous shift of the dips is destroyed at lower strains due to narrower dips.
Figure 9 shows how the spectra are changed as a result of smoothing of the acoustic wave front. The smoothness parameter in this simulation is 10 mm. It can be seen that the rejection bands become smoother and wider when compared to the stepwise front when the front covers half of the grating. At the same time, the gap between the bands is almost unchanged and is not widened.

6. Simulation of Reflective LPFG

The reflective LPFG is formed by making a mirror near the end of the grating. Depending on the space between the end of the grating and the mirror, we can obtain the structure that works as either uniform, phase-shifted, or cascaded LPFG. In Figure 10, an equivalent scheme of reflective LPFG is shown with the right side produced by reflecting the left side in the mirror, which is in the center of the figure. Thus, when an acoustic wave propagates through the fiber, the resulting structure represents five sections: non-deformed, deformed, phase shift, deformed, non-deformed. Thus, the simulation was carried out for the LPFG with the account of these five sections. In the equivalent scheme, the acoustic wave propagates from the center to the sides.
When the phase shift is zero, the initial grating is uniform, and its spectrum contains one notch at 1550 nm (Figure 11a); the incoming acoustic wave introduces a phase shift in the middle of the grating and compresses the fiber. As a result, the notch is shifted to shorter wavelengths and split into two when the new notch appears at the initial wavelength, and additional interference dips appear on both sides of the spectrum. When the acoustic front covers the whole grating, the LPFG becomes uniform again and the single notch is reconstituted in a shifted wavelength.
For the reflective π-shift LPFG, the initial spectrum coincides with the spectrum of transmissive π-shift LPFG with two symmetric dips (see Figure 11b). When the acoustic front propagates through the grating, the short-wavelength dip quickly disappears and appears back when the larger part of the grating is strained by the front. The long-wavelength dip shifts almost unaffected. Thus, the reflective π-shift LPFG has a more stable spectral dip that can be used as a signal for sensing applications. In terms of wavelength versus strain sensitivity, different kinds of LPFGS are similar.

7. Conclusions

We have studied the spectral behavior of LPFGs deformed by propagating acoustic shock front. Transmission of light through uniform, π-shift, and reflective LPFGs has been simulated using coupled mode theory for co-propagating core and cladding modes with inhomogeneous coupling coefficients. We have shown that, in general, a compressive acoustic shock front propagating through various LPFG structures shifts the resonances of cladding modes to shorter wavelength due to strain in the fiber. However, there are more complex spectral effects, which depend on parameters of the acoustic wave and the grating, such as front amplitude and smoothness, location of the front in the grating, magnitude of the phase shift, and the grating length. The spectrum represents a result of interference of the core and cladding modes interacting differently in different parts of the structure: strained and unstrained, before and after phase shift.
When the smoothness of the acoustic front is compared to the grating length, the notch of a uniform LPFG shifts continuously to shorter wavelengths. Otherwise, when front is steep with smoothness parameter that is significantly less than the grating length, the shifting notch splits into two dips with interference fringes around; indeed, LPFGs with the raised-cosine profile have interference fringes at shorter wavelength side. The spectral behavior of uniform LPFGs strained by a propagating shock front is similar to Bragg gratings, but the spectrum of LPFG cannot be smoothed by using cosine apodization at grating edges due to periodic dependence of transmission on the length of LPFG. The period and length of LPFGs is larger than that of fiber Bragg grating, and characteristics of the structure related to the smoothness of the shock front are changed proportionally.
For the π-shift LPFG having split spectral notch, the gap between the dips becomes several times wider than the gap in the unstrained grating. Increasing smoothness of the acoustic front results in continuous transition of the grating spectra with reduced interference effects. The rejection bands become wider and the gap between the bands is not widened.
The notch of the reflective LPFG without additional phase is shifted to shorter wavelengths with appearance of a new notch at the initial wavelength under the moving acoustic front, while the long-wavelength dip of the reflective π-shift LPFG remains almost unaffected by the acoustic front. The reflective π-shift LPFG has advantages for application in sensors due to its more stable spectral characteristics and one fiber connection needed for the sensor. In terms of sensitivity, different types of LPFG structures have similar parameters.
The results of this investigation can be used for creating LPFG sensors based on the measurement of transmission spectra around one of the cladding mode resonances for detecting shock waves. Furthermore, this analysis is crucial to assess the possibility to implement long-period fiber gratings-based sensors in order to detect leakage of medical gases in health care unities (pressure > 10 MPa), by detection of ultrasound waves (frequency < 100 kHz and 40–60 dB) produced in pipeline systems exhibiting different kinds of anomalies.

Author Contributions

Conceptualization, methodology and software, O.V.I.; writing—review and editing, O.V.I., P.C. and G.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work was developed within the scope of the project proMetheus—Research Unit on Materials, Energy and Environment for Sustainability, FCT Ref. UID/05975/2020, financed by national funds through the FCT/MCTES. This work was also financed by National Funds through the Portuguese funding agency, FCT-Fundação para a Ciência e a Tecnologia, within project UIDB/50014/2020. O.V. Ivanov acknowledges support from the Russian Ministry of Higher Education and Science (project 075-15-2021-581).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The processing codes and data segments can be obtained by contacting the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Calculation of the Cross-Coupling Coefficient

The coupling between the core mode and different cladding modes in the LPFG is determined by the cross-coupling coefficient κ co - ν (see Equation (16)). We represent this coefficient as function of dimensionless constant χ ν expressed through the overlap integral:
κ co - ν ( z ) = χ ν   n cl σ ( z ) k 0
χ ν = ω ε 0   n cl k 0 ∫ 0 r cl n 0 ( r ) 2 g ( r ) E co * ( r ) E ν ( r ) π r d r
Assuming that the grating is written in the fiber core, we have g ( r ) = 1 for r ≤ r co and g ( r ) = 0 for r > r co . The refractive index distribution n 0 ( r ) in the core is, therefore, equal to n co and can be approximated as n cl . Then, we obtain
χ ν = ω ε 0 n cl   k 0 ∫ 0 r co E co * ( r ) E ν ( r ) π r d r
The electric fields of the modes are normalized to 1 W of power carried by each mode. In this case, the mode fields satisfy the following orthogonality relation:
∫ 0 ∞ E μ * ( r ) E ν ( r ) 2 π r d r = 2 δ μ ν ω μ 0 β ν
where μ 0 is the vacuum permeability, and δ μ ν is the Kronecker delta. In Figure A1, the electric fields of the core mode and four cladding modes are plotted.
Sensors 21 07212 g0a1 550
Figure A1. Electric field distributions of five fiber modes. The gray rectangle is the integration area for calculating the overlap integral.
Figure A1. Electric field distributions of five fiber modes. The gray rectangle is the integration area for calculating the overlap integral.
Sensors 21 07212 g0a1
By integration of the product of two electric fields over the fiber core, which is depicted a gray rectangle in the figure, we can find the overlap constant χ ν for coupling of the core mode with four cladding modes: χ ν = 0.030 , 0.055, 0.076, 0.092 for ν = 2 , 3 , 4 , 5 , respectively. The field amplitude in the core is highest for the core mode HE11 and grows with the mode number for the cladding modes HE1ν. The overlap constant χ ν changes correspondingly.

References

  1. Kashyap, R. Fiber Bragg Gratings, 2nd ed.; Academic Press: San Diego, CA, USA, 2009; 632p. [Google Scholar]
  2. Rego, G. Arc-induced long period fiber gratings. J. Sens. 2016, 2016, 3598634. [Google Scholar] [CrossRef] [Green Version]
  3. Heck, M.; Gauthier, J.-C.; Tünnermann, A.; Vallée, R.; Nolte, S.; Bernier, M. Long period fiber gratings for the mitigation of parasitic laser effects in mid-infrared fiber amplifiers. Opt. Express 2019, 27, 21347–21357. [Google Scholar] [CrossRef] [PubMed]
  4. Takahashi, N.; Yoshimura, K.; Takahashi, S.; Imamura, K. Development of an optical fiber hydrophone with fiber Bragg grating. Ultrasonics 2000, 38, 581–585. [Google Scholar] [CrossRef]
  5. Fisher, N.E.; Webb, D.J.; Pannell, C.N.; Jackson, D.A.; Gavrilov, L.R.; Hand, J.W.; Zhang, L.; Bennion, I. Medical ultrasound detection using fiber Bragg gratings. Proc. SPIE 1999, 3541, 27–32. [Google Scholar] [CrossRef]
  6. Silva, R.E.; Franco, M.A.R.; Neves, P.T.; Bartelt, H.; Pohl, A.A.P. Numerical and experimental analysis of the modulation of fiber Bragg gratings by low-frequency complex acoustic waves. Opt. Fiber Technol. 2016, 30, 17–22. [Google Scholar] [CrossRef]
  7. Rodriguez, G.; Jaime, M.; Balakirev, F.; Mielke, C.H.; Azad, A.; Marshall, B.; La Lone, B.M.; Henson, B.; Smilowitz, L. Coherent pulse interrogation system for fiber Bragg grating sensing of strain and pressure in dynamic extremes of materials. Opt. Express 2015, 23, 14219–14233. [Google Scholar] [CrossRef]
  8. Rodriguez, G.; Gilbertson, S.M. Ultrafast fiber Bragg grating interrogation for sensing in detonation and shock wave experiments. Sensors 2017, 17, 248. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Barbarin, Y.; Lefrancois, A.; Zaniolo, G.; Chuzeville, V.; Jacquet, L.; Magne, S.; Luc, J.; Osmont, A. Optimization of detonation velocity measurements using a chirped fiber Bragg grating. Proc. SPIE 2015, 9480, 94800S. [Google Scholar] [CrossRef]
  10. Udd, E.; Benterou, J.; May, C.; Mihailov, S.J.; Lu, P. Review of high-speed fiber optic grating sensor systems. Proc. SPIE 2010, 7677, 76770B. [Google Scholar]
  11. Wei, P.; Lang, H.; Liu, T.; Xia, D. Detonation velocity measurement with chirped fiber Bragg grating. Sensors 2017, 17, 2552. [Google Scholar] [CrossRef] [Green Version]
  12. Ivanov, O.V. Simulation of reflection from Bragg gratings deformed by acoustic shock front. Opt. Commun. 2021, 481, 126537. [Google Scholar] [CrossRef]
  13. Hu, X.; Yang, Z.; Wang, Y.; Ma, X. Investigation of detecting ultrasound using a long period grating. DEStech Trans. Eng. Technol. Res. 2017. [Google Scholar] [CrossRef] [Green Version]
  14. JGaudron, J.-O.; Surre, F.; Sun, T.; Grattan, K.T.V. Long period grating-based optical fibre sensor for the underwater detection of acoustic waves. Sens. Actuators A Phys. 2013, 201, 289–293. [Google Scholar] [CrossRef]
  15. Gaudron, J.-O.; Surre, F.; Sun, T.; Grattan, K.T.V. LPG-based optical fibre sensor for acoustic wave detection. Sens. Actuators A Phys. 2012, 173, 97–101. [Google Scholar] [CrossRef]
  16. Wang, G.; Wang, Y.; Li, N. Axial strain sensitivity analysis of long period fiber grating by new transfer matrix method. Front. Optoelectron. China 2011, 4, 430–433. [Google Scholar] [CrossRef]
  17. Gu, Y.; Chiang, K.S. Effects of average index variation in apodized long-period fiber gratings. Photonic Sens. 2013, 3, 102–111. [Google Scholar] [CrossRef] [Green Version]
  18. Chan, F.Y.M.; Chiang, K.S. Analysis of apodized phase-shifted long-period fiber gratings. Opt. Commun. 2005, 244, 233–243. [Google Scholar] [CrossRef]
  19. Ke, H.; Chiang, K.S.; Peng, J.H. Analysis of phase-shifted long-period fiber gratings. IEEE Photonics Technol. Lett. 1998, 10, 1596–1598. [Google Scholar] [CrossRef]
  20. Wu, Z.; Dong, X.; Xu, Z. Chirp characteristics of long-period fiber gratings. Proc. SPIE 2001, 4595, 89–93. [Google Scholar] [CrossRef]
  21. Stegall, D.B.; Erdogan, T. Dispersion control with use of long-period fiber gratings. J. Opt. Soc. Am. A 2000, 17, 304–312. [Google Scholar] [CrossRef]
  22. Rego, G.; Falate, R.; Ivanov, O.V.; Santos, J.L. Simultaneous temperature and strain measurements performed by a step-changed arc-induced long-period fiber grating. Appl. Opt. 2007, 46, 1392–1396. [Google Scholar] [CrossRef]
  23. Falate, R.; Frazão, O.; Rego, G.; Ivanov, O.; Kalinowski, H.J.; Fabris, J.L.; Santos, J.L. Bending sensitivity dependent on the phase-shift imprinted in the long-period fibre gratings. Meas. Sci. Technol. 2007, 18, 3123–3130. [Google Scholar] [CrossRef]
  24. Guan, B.-O.; Zhang, A.-P.; Tam, H.-Y.; Chan, H.L.W.; Choy, C.-L.; Tao, X.-M.; Demokan, M.S. Step-changed long-period fiber gratings. IEEE Photonics Technol. Lett. 2002, 14, 657–659. [Google Scholar] [CrossRef] [Green Version]
  25. Chen, X.; Nie, J.; Cui, Z.; Chen, W.; Liu, W. Analysis of transmission spectrum of step-changed cascaded LPG and its applications in multi-channel filter. In Advances in Optoelectronics and Micro/Nano-Optics; IEEE: New York, NY, USA, 2010; pp. 1–3. [Google Scholar] [CrossRef]
  26. Yan, M.; Luo, S.; Zhan, L.; Wang, Y.; Xia, Y.; Zhang, Z. Step-changed period chirped long-period fiber gratings fabricated by CO2 laser. Opt. Commun. 2008, 281, 2784–2788. [Google Scholar] [CrossRef]
  27. Bae, J.K.; Bae, J.; Kim, S.H.; Park, N.; Lee, S.B. Dynamic EDFA gain-flattening filter using two LPFGs with divided coil heaters. IEEE Photonics Technol. Lett. 2005, 17, 1226–1228. [Google Scholar] [CrossRef]
  28. Liu, Y.; Williams, J.A.R.; Zhang, L.; Bennion, I. Phase shifted and cascaded long-period fiber gratings. Opt. Commun. 1999, 164, 27–31. [Google Scholar] [CrossRef]
  29. Chiavaioli, F.; Baldini, F.; Trono, C. Manufacturing and spectral features of different types of long period fiber gratings: Phase-shifted, turn-around point, internally tilted, and pseudo-random. Fibers 2017, 5, 29. [Google Scholar] [CrossRef] [Green Version]
  30. Sakata, H.; Komori, K.; Ono, Y. Tunable bandpass filters based on π-shifted mechanical long-period fibre gratings in two-mode fibres. Electron. Lett. 2015, 51, 843–845. [Google Scholar] [CrossRef]
  31. Zhou, X.; Shi, S.; Zhang, Z.; Zou, J.; Liu, Y. Mechanically-induced π-shifted long-period fiber gratings. Opt. Express 2011, 19, 6253–6259. [Google Scholar] [CrossRef] [PubMed]
  32. Hishiki, K.; Li, H. Phase-shift formed in a long period fiber grating and its application to the measurements of temperature and refractive index. Opt. Express 2013, 21, 11901–11912. [Google Scholar] [CrossRef]
  33. Gao, R.; Jiang, Y.; Jiang, L. Multi-phase-shifted helical long period fiber grating based temperature-insensitive optical twist sensor. Opt. Express 2014, 22, 15697–15709. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, P.; Xian, L.; Li, H. Fabrication of phase-shifted long-period fiber grating and its application to strain measurement. IEEE Photonics Technol. Lett. 2015, 27, 557–560. [Google Scholar] [CrossRef]
  35. Qian, J.R.; Chen, H.F. Gain flattening fibre filters using phase-shifted long period fibre gratings. Electron. Lett. 1998, 34, 1132–1133. [Google Scholar] [CrossRef]
  36. Deparis, O.; Kiyan, R.; Pottiez, O.; Blondel, M.; Korolev, I.G.; Vasiliev, S.A.; Dianov, E.M. Bandpass filters based on π-shifted long-period fiber gratings for actively mode-locked erbium fiber lasers. Opt. Lett. 2001, 26, 1239–1241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Azana, J.; Kulishov, M. All-fibre ultrafast optical differentiator based on p phase-shifted long-period grating. Electron. Lett. 2005, 41, 1368–1369. [Google Scholar] [CrossRef]
  38. Liu, T.; Han, M. Analysis of π-phase-shifted fiber Bragg gratings for ultrasonic detection. IEEE Sens. J. 2012, 12, 2368–2373. [Google Scholar] [CrossRef]
  39. Swart, P.L. Long-period grating Michelson refractometric sensor. Meas. Sci. Technol. 2004, 15, 1576–1580. [Google Scholar] [CrossRef]
  40. Falate, R.; Frazão, O.; Rego, G.; Fabris, J.L.; Santos, J.L. Refractometric sensor based on a phase-shifted long-period fiber grating. Appl. Opt. 2006, 45, 5066–5072. [Google Scholar] [CrossRef] [PubMed]
  41. Xue, H.; Xu, Z.; Chen, H.; Yang, Y.; You, J.; Yan, J.; Fu, H.; Zhang, D. Continuous liquid level sensor based on a reflective long period fiber grating interferometer. Meas. Sci. Technol. 2015, 26, 037001. [Google Scholar] [CrossRef]
  42. Jiang, M.; Zhang, A.P.; Wang, Y.C.; Tam, H.Y.; He, S. Fabrication of a compact reflective long-period grating sensor with a cladding-mode-selective fiber end-face mirror. Opt. Express 2009, 17, 17976. [Google Scholar] [CrossRef] [Green Version]
  43. Kim, D.W.; Zhang, Y.; Cooper, K.L.; Wang, A. In-fiber reflection mode interferometer based on a long-period grating for external refractive-index measurement. Appl. Opt. 2005, 44, 5368–5373. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Quero, G.; Zuppolini, S.; Consales, M.; Diodato, L.; Vaiano, P.; Venturelli, A.; Santucci, M.; Spyrakis, F.; Costi, M.P.; Giordano, M.; et al. Long period fiber grating working in reflection mode as valuable biosensing platform for the detection of drug resistant bacteria. Sens. Actuat. B Chem. 2016, 230, 510–520. [Google Scholar] [CrossRef]
  45. Dey, T.K.; Biswas, P.; Basumallick, N.; Bandyopadhyay, S.; Bandyopadhyay, S. Realization of long period fiber grating in reflection mode operating near turn around point. IEEE Sens. J. 2017, 17, 4100–4106. [Google Scholar] [CrossRef]
  46. Rana, S.; Kandadai, N.; Subbaraman, H. Towards the design of a wideband reflective long period grating distributed sensor. J. Phys. Commun. 2020, 4, 065015. [Google Scholar] [CrossRef]
  47. Bhatia, V.; Vengsarkar, A.M. Optical fiber long-period grating sensors. Opt. Lett. 1996, 21, 692–694. [Google Scholar] [CrossRef]
  48. Ivanov, O. Coupling of hybrid modes in strained and heated fibers. Opt. Commun. 2004, 239, 311–321. [Google Scholar] [CrossRef]
  49. Rego, G.; Ivanov, O.V. Two types of resonances in long-period gratings induced by arc discharges in boron/germanium co-doped fibers. Opt. Lett. 2007, 32, 2984–2986. [Google Scholar] [CrossRef] [PubMed]
Sensors 21 07212 g001 550
Figure 1. Scheme of LPFG structures deformed by an acoustic shock front: (a) homogeneous, (b) π-sift, (c) reflective, (d) reflective π-shift.
Figure 1. Scheme of LPFG structures deformed by an acoustic shock front: (a) homogeneous, (b) π-sift, (c) reflective, (d) reflective π-shift.
Sensors 21 07212 g001
Sensors 21 07212 g002 550
Figure 2. Longitudinal coordinate dependence of amplitude of the shock front.
Figure 2. Longitudinal coordinate dependence of amplitude of the shock front.
Sensors 21 07212 g002
Sensors 21 07212 g003 550
Figure 3. Transmission coefficient of grating for LP05 cladding mode, when the acoustic front is in its center z 0 = L / 2 , for different strains s 0 = 0 , 6, 12, and 18 mstrain.
Figure 3. Transmission coefficient of grating for LP05 cladding mode, when the acoustic front is in its center z 0 = L / 2 , for different strains s 0 = 0 , 6, 12, and 18 mstrain.
Sensors 21 07212 g003
Sensors 21 07212 g004 550
Figure 4. Transmission of grating for LP05 cladding mode, when the acoustic front is in its center z 0 = L / 2 , for η = 0 ,   12.5 ,   and   25 mm ( s 0 = 18 mstrain).
Figure 4. Transmission of grating for LP05 cladding mode, when the acoustic front is in its center z 0 = L / 2 , for η = 0 ,   12.5 ,   and   25 mm ( s 0 = 18 mstrain).
Sensors 21 07212 g004
Sensors 21 07212 g005 550
Figure 5. Evolution of transmission spectra with transition point z 0 moving through the grating for three magnitudes of smoothness η : (a) 0, (b) 10 mm, and (c) 20 mm ( s 0 = 18 mstrain).
Figure 5. Evolution of transmission spectra with transition point z 0 moving through the grating for three magnitudes of smoothness η : (a) 0, (b) 10 mm, and (c) 20 mm ( s 0 = 18 mstrain).
Sensors 21 07212 g005
Sensors 21 07212 g006 550
Figure 6. Evolution of transmission spectra with increasing grating length, when the acoustic front is in its center z 0 = L / 2 , for η = 0 mm and s 0 = 18 mstrain.
Figure 6. Evolution of transmission spectra with increasing grating length, when the acoustic front is in its center z 0 = L / 2 , for η = 0 mm and s 0 = 18 mstrain.
Sensors 21 07212 g006
Sensors 21 07212 g007 550
Figure 7. Evolution of transmission spectra for raised-cosine grating and two magnitudes of smoothness (a) η = 0   and (b) 25 mm ( s 0 = 24 mstrain).
Figure 7. Evolution of transmission spectra for raised-cosine grating and two magnitudes of smoothness (a) η = 0   and (b) 25 mm ( s 0 = 24 mstrain).
Sensors 21 07212 g007
Sensors 21 07212 g008 550
Figure 8. Evolution of transmission spectra with transition point z 0 moving through the π-shift LPFG for different strains: (a) s 0 = 9 , (b) 18, and (c) 24 mstrain ( η = 0 ).
Figure 8. Evolution of transmission spectra with transition point z 0 moving through the π-shift LPFG for different strains: (a) s 0 = 9 , (b) 18, and (c) 24 mstrain ( η = 0 ).
Sensors 21 07212 g008
Sensors 21 07212 g009 550
Figure 9. Evolution of transmission spectra with transition point z 0 moving through the π-shift LPFG for smoothness η = 10 mm (strain s 0 = 18 ).
Figure 9. Evolution of transmission spectra with transition point z 0 moving through the π-shift LPFG for smoothness η = 10 mm (strain s 0 = 18 ).
Sensors 21 07212 g009
Sensors 21 07212 g010 550
Figure 10. Equivalent scheme of reflective LPFG structure with π-shift.
Figure 10. Equivalent scheme of reflective LPFG structure with π-shift.
Sensors 21 07212 g010
Sensors 21 07212 g011 550
Figure 11. Evolution of transmission spectra with transition point z 0 moving through the reflective LPFG for strain s 0 = 18 with (a) zero phase shift and (b) π-shift.
Figure 11. Evolution of transmission spectra with transition point z 0 moving through the reflective LPFG for strain s 0 = 18 with (a) zero phase shift and (b) π-shift.
Sensors 21 07212 g011
Table
Table 1. Parameters used in simulation.
Table 1. Parameters used in simulation.
ParameterDesignationNumerical Value
Wavelength λ 0 1550 nm
Core radius r co 4.2 μm
Core refractive index n co 1.4494
Cladding refractive index n cl 1.4440
Strain-optic coefficient p a 0.22
Grating length L 40 mm
Grating period Λ 430 μm
Refractive index modulation σ 1.48 ⋅ 10 − 4
Effective refractive indices n co ( eff ) 1.44649
n 5 ( eff ) 1.44306
Coupling coefficient κ co - 5 1.6 ⋅ 10 − 4 μm–1
Strain correction α co - co − α   ν - ν –0.01
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ivanov, O.V.; Caldas, P.; Rego, G. Simulation of the Transmission Spectrum of Long-Period Fiber Gratings Structures with a Propagating Acoustic Shock Front. Sensors 2021, 21, 7212. https://doi.org/10.3390/s21217212

AMA Style

Ivanov OV, Caldas P, Rego G. Simulation of the Transmission Spectrum of Long-Period Fiber Gratings Structures with a Propagating Acoustic Shock Front. Sensors. 2021; 21(21):7212. https://doi.org/10.3390/s21217212

Chicago/Turabian Style

Ivanov, Oleg V., Paulo Caldas, and Gaspar Rego. 2021. "Simulation of the Transmission Spectrum of Long-Period Fiber Gratings Structures with a Propagating Acoustic Shock Front" Sensors 21, no. 21: 7212. https://doi.org/10.3390/s21217212

APA Style

Ivanov, O. V., Caldas, P., & Rego, G. (2021). Simulation of the Transmission Spectrum of Long-Period Fiber Gratings Structures with a Propagating Acoustic Shock Front. Sensors, 21(21), 7212. https://doi.org/10.3390/s21217212

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 pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy