# Investigation on the Effect of Underwater Acoustic Pressure on the Fundamental Mode of Hollow-Core Photonic Bandgap Fibers

Authors: Adel Abdallah; Zhang Chaozhu; Zhong Zhi
Journal: International Journal of Optics
 (2015)

---

## Abstract

Recently, microstructured optical fibers have become the subject of extensive research as they can be employed in many civilian and military applications. One of the recent areas of research is to enhance the normalized responsivity (NR) to acoustic pressure of the optical fiber hydrophones by replacing the conventional single mode fibers (SMFs) with hollow-core photonic bandgap fibers (HC-PBFs). However, this needs further investigation. In order to fully understand the feasibility of using HC-PBFs as acoustic pressure sensors and in underwater communication systems, it is important to study their modal properties in this environment. In this paper, the finite element solver (FES) COMSOL Multiphysics is used to study the effect of underwater acoustic pressure on the effective refractive indexneff of the fundamental mode and discuss its contribution to NR. Besides, we investigate, for the first time to our knowledge, the effect of underwater acoustic pressure on the effective area Aeff and the numerical aperture (NA) of the HC-PBF.


---

## Body

## 1. Introduction

Optical fiber hydrophone is a significant area of research. For many years, researchers have shown the prospects of using the conventional SMF interferometric hydrophones as alternatives to the conventional sound navigation and ranging (SONAR) systems [1]. However, the conventional SMF is made of glass that has high Young’s modulus (E); in addition, its change of effective refractive index of the fundamental mode (neff) due to applied acoustic pressure has opposite sign with respect to length change, and hence both affect NR [2]. As a result, researchers began to search for alternatives to enhance NR. One possible solution is to test the microstructured optical fibers that can be classified into two classes: solid-core PCF (SC-PCF) and HC-PBF. It was shown experimentally that SC-PCF has about the same phase sensitivity to axial strain as SMF [3]. As a result, using SC-PCF as alternatives to SMF is infeasible. However, because SC-PCFs are bend-insensitive, they can be used to reduce the hydrophone size [4]. It was reported that HC-PBFs have many advantages over SMFs that provide better responsivity to measurands and eligibility for many sensing applications [2, 5–7]. Some of these advantages are as follows. (1) The design and manufacturing flexibility of HC-PBFs reduces the effective Young’s modulus of the fiber and enhances the NR of the HC-PBF to acoustic pressure [2], (2) matching between the mode indexes of the HC-PBF with the ambient air helps to reduce the back-reflected light at the fiber end faces which is useful for many applications [7], (3) holes of the HC-PBF can be filled with a substance with opposite thermal expansion to make the material insensitive to temperature [5], (4) in the HC-PBF, because the optical mode propagates in air, it has smaller Faraday, Kerr, and thermal constants than solid-silica cores, and this reduces the dependencies on temperature, magnetic field, and power fluctuations [5], and (5) HC-PBFs are almost entirely bend-insensitive and can be bent to very small diameters (<1 cm) with minimal loss, and this makes it suitable for small-size hydrophone systems [6, 8]. However, using HC-PBFs as underwater acoustic sensors needs further investigation to be feasible alternatives to their counterparts of SMFs. Effective mode area (Aeff) is a key factor in designing PCFs [9]. It determines how tightly the mode is confined to the core of the fiber. The effective area can be used to study the nonlinearities, mode-field diameter (MFD), confinement losses, bending losses, splicing losses, and numerical aperture (NA) of the optical fiber [9, 10]. For this reason, studying the effect of acoustic pressure on Aeff allows getting information and important relations about how the acoustic pressure affects other important quantities. For accurate results, it is important to accurately model Aeff of the HC-PBF. It was proposed in [11] that the electrical field distribution of the PCFs is non-Gaussian and cannot be determined by assuming a conventional step-index distribution. Also, it was found that the ITU-T Petermann II definition is the most suitable for describing Aeff and MFD of PCFs with non-Gaussian distribution; as a result, it is adopted in this paper. ITU-T stands for the International Telecommunication Union-Telecommunication Standardization Sector [12]. To our knowledge, the effect of underwater acoustic pressure on Aeff of HC-PBF has not been studied. In this paper, we study the interferometric optical fiber hydrophone, which is based on acoustic pressure to phase transduction mechanism, in which the investigated HC-PBF represents its measuring arm. The effect of underwater acoustic pressure on the fundamental mode of hollow-core photonic bandgap fibers is investigated by coupling between the acoustic-solid interaction (ASI) and the electromagnetic waves (EMW) modules in the FES COMSOL Multiphysics. The ASI module is used to apply acoustic pressures of different amplitudes and frequencies that cause structural deformation and the induced stresses and strains for the investigated HC-PBF are calculated. The EMW module in the FES is used to calculate neff for undeformed HC-PBF; then the coupled ASI and EMW modules are used to calculate neff for the deformed fiber. This enables us to study and analyse the effect of acoustic pressure on neff, Aeff, the MFD, and NA of the investigated HC-PBF. The investigated HC-PBF is the commercial HC-1550 with air-filling ratio (η) of 92% that represents the ratio of the air hole diameter (dh) in the microstructured cladding to the pitch (Λ) which is the central distance between two adjacent air holes. The used parameters of the investigated fiber are listed in Table 1.Table 1
Physical parameters of the investigated optical fiber [2, 6].
Region
HC-PBF
Core
Material: aird1 = 10.1 μm
Clad
Air-silica honeycomb
d
2 = 70 μm, ν2=0.17,  Λ=3.8 μmdh = 3.5 μm, η = 92%Er=Eθ=3/21-η3Eo Ez=1-ηEo
Solid-silica clad
d
3 = 120 μm, E3 = 72 GPa, ν3 = 0.17
Polymer coating
Material: acrylatet = 93.5 μm, E4 = 0.75 GPa, ν4 = 0.45
E is Young’s modulus of each layer, Eo is Young’s modulus of the silica, ν is Poisson’s ratio, d is the diameter of each region, dh is the air hole diameter in the microstructured region, t is the thickness of the acrylate coating, and Er, Eθ, and Ez are used to describe the anisotropy of the microstructured area.


## 2. Mathematical Model

The cross-section of the investigated HC-PBF is shown in Figure1. The HC-PBF is modeled as four circular regions, an air-core, an air-silica microstructured inner cladding consisting of array of cylindrical air holes, a solid silica outer cladding, and an acrylate layer. The parameters of each region of the HC-PBF are denoted by a superscript (i=1,2,3, and 4). The HC-PBF is with 8 rings of arrays of air holes arranged in a triangular lattice with Λ=3.8 μm and η=92%. The air-core of diameter (d1) is formed by removing seven central air holes.Figure 1
Cross-section of the investigated HC-PBF with an air-core, a microstructured air-silica inner cladding, a solid silica outer cladding, and a polymer coating, where (a–d) represent the radius of each region.The equations that describe the sound propagation in fluids are derived from the governing equations of fluid flow such as the mass conservation described by the continuity equation; the conservation of momentum, which is known as the Navier-Stokes equation; an energy conservation equation; the constitutive equations; and an equation of state that describes the relation between thermodynamic variables [13]. The acoustic pressure (p) is governed by the wave equation and is given by(1)1ρoco2∂2p∂t2+∇·-1ρo∇p-q=Q,where t is the time, ρo is the density of the fluid, q and Q are the acoustic dipole and monopole source, respectively, and co is the velocity of the acoustic wave in the medium. The wave equation can be solved in the frequency domain to expand the acoustic signal into harmonic components via its Fourier series. A harmonic solution has the form(2)px,y=pxeiωt,where ω is the angular frequency and the actual physical value of the acoustic pressure is the real part of (2); consequently, the time-dependent wave equation reduces to the Helmholtz equation given by(3)∇·-1ρo∇p-q-ω2ρoco2p=Q.In the homogenous case where the two source terms q and Q are zero, the solution to the Helmholtz equation is the plane wave given by(4)px,y=poxei(ωt-k·x),where po is the amplitude of the wave and it is moving in the k direction with angular frequency ω and wave number k=k. The acoustic pressure (p) acting on the microstructured fiber induces a stress distribution in the fiber’s cross-section and structural deformation. Both factors affect the phase represented by a change in the fiber length and the effective refractive index and is given by(5)φ=2πλneffL,where λ is the wavelength of the propagating light, and L is the length of the fiber. The applied acoustic pressure induces stress distribution in the optical fiber material. As a result of the stress-optic effect, stress induces anisotropic change of the index of refraction within the optical fiber. In this study, this process is achieved by coupling between the ASI and the EMW module in COMSOL Multiphysics. Acoustic pressures with different amplitudes and frequencies are applied by the ASI to obtain the induced strain and stress vectors and finally the index of refraction due to the stress-optic effect is calculated by the general linear stress-optical relation as follows [14, 15]. The general linear stress-optical relation is given by the tensor notation as(6)Δnij=-Bijklσkl,where Δnij=nij-noIij, nij is the index of refraction tensor, no is the index of refraction for a stress-free material, Iij is the identity tensor, Bijkl is the stress-optical tensor, and σkl is the stress tensor. The number of independent parameters in the stress-optical tensor that characterizes this relation is reduced by symmetry. Because nij and σkl are both symmetric matrices, Bijkl=Bjikl and Bijkl=Bijlk. For the fiber material, the number of independent parameters is reduced to two independent parameters, B1=-6.9×10-13 m2/N and B2=-41.9×10-13 m2/N, that represent the first and second stress-optical coefficients, respectively [15]. In this case, the stress-optical relation simplifies to(7)ΔnxΔnyΔnz=-B1B2B2B2B1B2B2B2B1σxσyσz,where nx=n11, ny=n22, and nz=n33 are the refractive indices along the x-, y-, and z-axis, respectively, and σx=σ11, σy=σ22, and σz=σ33 are the principal components of the induced stresses in the three directions. Using the two parameters B1 and B2, the model assumes that the nondiagonal parts of nij are negligible. As a result, the used general linear stress-optical relation, which describes the relation between the induced refractive indices and the principle stresses in x, y, and z directions, is reduced to(8)nx=no-B1σx-B2σy+σz,ny=no-B1σy-B2σx+σz,nz=no-B1σz-B2σx+σy.For accurate calculations of modal characteristics of the investigated optical fiber, the full vectorial wave equations need to be solved. The FES COMSOL Multiphysics is used to solve the vectorial electric field wave equation as an eigenvalue problem and is given by [14, 16] (9)∇×μr-1∇×E-λeE=0,in which E is the electric field vector, μr is the relative permeability, and λe is the eigenvalue given by(10)λe=ko2εr,in which ko=2π/λ is the free-space wave number and εr is the relative permittivity. In this model the electromagnetic wave, frequency domain interface, in the FES COMSOL Multiphysics is used for the mode analysis. The simulation is set up with the electric field components Ex, Ey, and Ez which are the dependent variables. The EMW can be described by the form [17](11)E=Ex,yej-αz=Exx,y,Eyx,y,Ezx,yej-αz,where the parameter α=γ+jβ is the complex propagation constant, γ represents the attenuation along the propagation direction, and β is the propagation constant. The effective mode index of a confined mode is given by(12)neff=βko.The acoustic pressure primarily affects the L and neff terms in (5). The normalized responsivity (NR) is a figure of merit independent of wavelength and optical fiber dimensions are commonly used to compare between different hydrophone designs. NR of the HC-PBF is given by(13)NR=dφφ(dp)=1LdLdp+1neffdneffdp=εz2dp+1neffdneffdp,where εz2=dL/L is the axial strain of the microstructured region of the HC-PBF, and the superscript (2) is used to denote the microstructured region of the HC-PBF. In this paper, we are interested in calculating the index term of (13).We investigate another area of research which is the effect of the acoustic pressure on the effective areaAeff, MFD, and NA of the HC-PBF. Aeff depends on fiber’s index of refraction and the propagating wavelength, and accurate description of Aeff is based on ITU-T Petermann II definition and is given by [12, 17](14)Aeff=∬-∞∞E(x,y)2dxdy2∬-∞∞E(x,y)4dxdy,where E(x,y) is the modal field distribution inside the fiber. A general relation between the MFD and Aeff is given by(15)Aeff=kλπ4MFD2,where k(λ) is wavelength dependent mapping value. If the electrical field distribution were Gaussian as in the case of conventional SMFs, the k(λ) factor would be equal to one [11]. Furthermore, as d/Λ becomes larger, the correction factor k(λ) approaches 1.2 for PCFs [11]. In this paper, k(λ)=1.2 is considered.The NA of the SC-PCFs is investigated in [18]. However, the calculation of NA in HC-PBF is different where the refractive index of the cladding is larger than that in the core [19]. A formula for the NA of HC-PBF was proposed in [11] and is given by (16)NA=1-βUβo2,where βU is the propagation constant of the upper edge of the bandgap and is related to the effective refractive index by βU=2πneff/λU, where λU is the wavelength corresponding to the upper edge of the bandgap. The NA of the PBF is wavelength dependent [11]. In this study, the effect of the acoustic pressure on the NA is calculated using (16) where the calculations are based on the calculated bandgaps of HC-1550 proposed in [7] and the propagation wavelength λo=1550 nm and λU=1650 nm.As a conclusion, these calculations are performed by coupling between the ASI and the EMW frequency domain interface modules in the FES as follows: in the simulation setup of the EMW module,nx, ny, and nz are inserted as the main diagonal elements of the anisotropic index of refraction tensor. Acoustic pressures with different amplitudes and frequencies that act on the investigated optical fiber are applied by the ASI. Once the induced stress distribution in the optical fiber is calculated, the anisotropic refractive index elements are calculated using (8), and the EMW module is used for mode analysis. Mode analysis study allows us to calculate neff, Aeff, MFD, and NA. By this way it is easy to obtain the relation between the applied acoustic pressure and the investigated parameters, and the contribution of the acoustic pressure induced neff change on φ and NR is calculated.


## 3. Simulation Results and Analysis

In this section, simulation results are introduced. The material’s physical parameters of the investigated HC-PBF shown in Table1 are imported into the FES. The microstructured region of the HC-PBF is modeled as anisotropic material while the silica outer cladding and the acrylate regions are modeled as isotropic materials [2]. The investigated HC-PBF is based on two-dimensional triangular structure with the air-core being formed by removing seven central holes. The index of refraction of silica cladding is 1.444, and the number of the air hole rings is eight. Perfectly matched layer (PML) is used to remove spurious reflections. We study the response of the HC-PBF to acoustic pressure by coupling between the ASI and the EMW modules in the FES. This allows easy transfer of the required data between the acoustics and electromagnetic waves modules and provides accurate calculations. After setting up the model geometry and inserting the required equations and data to the FES, the first step of calculations is to perform mode analysis for undeformed fiber. The calculated fundamental mode’s intensity profile of the undeformed HC-PBF is shown in Figure 2, where neff=0.996128. The second step is to use the acoustics module to apply acoustic pressure of different amplitudes and frequencies that cause structural deformation and the induced stresses and strains of the investigated fiber are calculated. The EMW module exchanges data with the ASI module by coupling between them and then by performing mode analysis; neff and Aeff corresponding to the deformed structure are calculated. This allows calculating the phase change and NR given by (5) and (13), respectively.Figure 2
Fundamental mode confinement of the investigated optical fiber.Figure3 shows neff of the investigated HC-PBF as a function of the acoustic pressure at acoustic frequency 10 kHz. Due to applied acoustic pressure, the propagation constant changes because of two factors: the geometrical deformation and the change of the refractive index due to the stress-optic effect. From this figure, it can be seen that the calculated neff for the undeformed HC-PBF is 0.996128 and it increases as the acoustic pressure increases, where dneff/dp is 1.2×10-12 Pa−1, where the acoustic pressure is swept from 0 to 100 MPa at frequency 10 kHz. Figure 4 shows the calculated NR of the HC-1550 and SMF as a function of acoustic frequency obtained from a previous study [20]. It can be seen that the calculated NR of the investigated HC-PBF is -2.5577×10-11 Pa−1; consequently, from (13) the acoustic pressure-induced HC-PBF length change is -2.6777×10-11 Pa−1. Simulation results showed that the contribution of the index change term to the total sensitivity of the HC-PBFs to acoustic pressure is minor with respect to the fiber length change. However, generally, for accurate design and simulations of the HC-PBFs, the index change term should be taken into account.Figure 3
The calculated fundamental mode refractive index of the investigated HC-PBF as a function of the acoustic pressure at acoustic frequency 10 kHz.Figure 4
NR of the HC-1550 and SMF as a function of acoustic frequency.Figures5 and 6 show the calculated effective area and MFD of the investigated HC-PBF as a function of the acoustic pressure at acoustic frequency 10 kHz. The calculated Aeff for undeformed fiber is 57.1 μm2 while the MFD is 7.758×10-6 μm. It can be seen that the change in Aeff is inversely proportional to the acoustic pressure, where dAeff/dp is -2×10-8 μm2 Pa−1, while the calculated MFD change with pressure is dMFD/dp=-1.35×10-15 μm Pa−1.Figure 5
The calculated effective area of the investigated HC-PBF as a function of the acoustic pressure at acoustic frequency 10 kHz.Figure 6
The calculated mode-field diameter of the investigated HC-PBF as a function of the acoustic pressure at acoustic frequency 10 kHz.Figure7 shows the fundamental mode index corresponding to wavelengths 1550 nm and 1650 nm used to calculate NA of HC-PBF as a function of the acoustic pressure at acoustic frequency 10 kHz obtained by coupling the ASI and EMW modules in the FES, where the pressure is swept from 0 to 1 GPa. It can be seen that neff for the undeformed fiber of λ=1550 nm and λ=1650 nm is 0.996128 and 0.995313, respectively. Using (16) and the calculations proposed in this figure, we can preliminarily investigate the effect of acoustic pressure on the NA of the investigated HC-PBF. The calculated NA of the investigated HC-PBF as a function of the acoustic pressure at acoustic frequency 10 kHz is shown in Figure 8. It can be seen that NA for undeformed fiber is 0.3449383, and it is inversely proportional to the acoustic pressure where the change of NA with the acoustic pressure is dNA/dp=-5.43×10-13 Pa−1.Figure 7
Fundamental mode index corresponding to wavelengths 1550 nm and 1650 nm as a function of the acoustic pressure at acoustic frequency 10 kHz.Figure 8
The calculated NA of the investigated HC-PBF as a function of the acoustic pressure at acoustic frequency 10 kHz.


## 4. Conclusion

In this paper, the FES COMSOL Multiphysics is used to study the effect of acoustic pressure on the fundamental mode of a HC-PBF, mainly onneff, Aeff, MFD, and NA. The proposed simulation results showed that the acoustic pressure-induced change of neff is 1.2×10-12 Pa−1, where the acoustic pressure is swept from 0 to 100 MPa at frequency 10 kHz. Our previous study showed that the calculated NR of the investigated HC-PBF is -2.5577×10-11 Pa−1; consequently, the acoustic pressure-induced HC-PBF length change is -2.6777×10-11 Pa−1. It can be concluded that the contribution of the index change term to the total sensitivity of the HC-PBFs to acoustic pressure is minor with respect to the fiber length change but it should be taken into account for accurate design and simulations of the HC-PBFs. Another area of research which is the effect of the acoustic pressure on the effective area Aeff, MFD, and NA of the HC-PBF is investigated. It is shown that the change in Aeff is inversely proportional to the pressure, where dAeff/dp is -2×10-8 μm2 Pa−1, while the calculated MFD change with pressure is dMFD/dp=-1.35×10-15 μm Pa−1. Finally, we preliminarily investigated the effect of acoustic pressure on the NA of the investigated HC-PBF. It is shown that NA for undeformed fiber is 0.3449383, and it is inversely proportional to the acoustic pressure, where the change of NA with the acoustic pressure is dNA/dp=-5.43×10-13 Pa−1.



---

*Source: 815967-2015-03-29.xml*