검색
검색 팝업 닫기

Ex) Article Title, Author, Keywords

Article

Split Viewer

Article

Current Optics and Photonics 2020; 4(4): 293-301

Published online August 25, 2020 https://doi.org/10.3807/COPP.2020.4.4.293

Copyright © Optical Society of Korea.

Chip-scale Temperature-compensated Superstructured Waveguide Bragg Grating Based Multiparametric Sensor

Naik Parrikar Vishwaraj1,*, Chandrika Thondagere Nataraj2,**, Ravi Prasad Kogravalli Jagannath3, Prashanth Gurusiddappa1,***, and Srinivas Talabattula4

1Department of Electronics and Communication Engineering, National Institute of Technology Goa, Ponda 403401, India, 2Department of Electronics and Telecommunication Engineering, Siddaganga Institute of Technology, Tumakuru-Karnataka 572103, India, 3Department Applied Sciences, National Institute of Technology Goa, Ponda 403401, India, 4Department of Electrical Communication Engineering, Indian Institute of Science, Bangalore-Karnataka 560012, India

Corresponding author: rajparrikar@gmail.com

Received: December 26, 2019; Revised: February 27, 2020; Accepted: March 16, 2020

In this paper we propose and theoretically analyze a monolithic multiparametric sensor consisting of a superstructure of surface-relief waveguide Bragg gratings (WBGs), a micro-machined diaphragm, and a cantilever beam. Diaphragms of two different configurations, namely circular and square, are designed and analyzed separately for pressure measurement. The square diaphragm is then selected for further study, since it shows relatively higher sensitivity compared to the circular one, as it incurs more induced stress when any pressure is applied. The cantilever beam with a proof mass is designed to enhance the sensitivity for acceleration measurement. A unique mathematical method using coupled-mode theory and the transfer-matrix method is developed to design and analyze the shift in the Bragg wavelength of the superstructure configuration of the gratings, due to simultaneously applied pressure and acceleration. The effect of temperature on the wavelength shift is compensated by introducing another Bragg grating in the superstructure configuration. The measured sensitivities for pressure and acceleration are found to be 0.21 pm/Pa and 6.49 nm/g respectively.

Keywords: Waveguide Bragg gratings, Pressure sensor, Accelerometer, Optical sensors, MOEMS

In these days, fiber-Bragg-grating (FBG) sensors are used successfully for a broad spectrum of applications in automobiles, aerospace, and biosensing due to their multiplexed sensing capabilities, lightweight design, and compactness [1-5]. Measurements of pressure, temperature, and acceleration are critical in industrial applications, and subsequently there is high demand for integrated composite sensors [6, 7]. However, these are based on the piezoresistive principle of sensing, in which piezoresistive elements are placed over a diaphragm and cantilever beam to measure pressure and acceleration respectively. Compact, planar photonic components are being considered for such applications, due to their potential for large-scale integration and reduced fabrication cost.

Micro-opto-electro-mechanical systems (MOEMS), realized by combining silicon photonics with micro-electro-mechanical systems (MEMS), offer the advantages of a small device footprint, cheapness, batch-fabrication possibilities, and easier CMOS integration with less complexity. MOEMS pressure sensors using a Mach-Zehnder interferometer (MZI), ring resonator, directional coupler, and gratings have been reported [3, 8-10].

Waveguide Bragg gratings (WBGs) on the silicon-on-insulator (SOI) platform have attracted interest for optical communication and various sensing applications. WBGs were first demonstrated by Murphy et al. [1] in 2001. The Bragg gratings are realized by physically corrugating dielectric waveguides; their sensitivity is obtained by observing the shift in the reflection wavelength. This method of detection, using the shift in wavelength, is preferred over methods that utilize the change in intensity of propagated light, as the former is more immune to noise [11]. However, Bragg-grating-based sensors are equally sensitive to temperature as to the parameters being measured. Hence it remains a difficult task to discriminate between the shift in wavelength due to temperature change and that due to the parameter being measured. If the shift due to temperature is not compensated, a MOEMS sensor gives inaccurate results, so it is necessary to subtract the wavelength shift due to temperature effects, so that we can attribute the wavelength shift only to changes in the actual parameter being measured, such as pressure or acceleration. In earlier works a chirped grating and dual gratings have been proposed to compensate for the effect of temperature in pressure sensors [12, 13]. However, surface-relief waveguide Bragg gratings for monolithic multiparametric sensing with temperature compensation have not been addressed so far. In addition, so far there has been no unique mathematical model to analyze such systems theoretically. Such systems could be an affordable solution for automobile, aerospace, and health applications, since they could be easily manufactured using standard microfabrication technologies.

In this paper we propose a WBG-based monolithic multiparametric sensor that is immune to temperature aberrations, for both pressure and acceleration measurements. A unique mathematical model using coupled-mode theory and the transfer-matrix method is developed to design the configuration of the superstructure of surface-relief Bragg gratings. The proposed sensor is made up of a superstructure configuration of three identical surface-relief waveguide Bragg gratings, designed with distinct Bragg wavelengths. Two of these gratings are located respectively on a micromachined diaphragm structure and a cantilever beam structure, where pressure and acceleration are applied. The wavelength shifts of these gratings are proportional to both pressure and acceleration, along with temperature. The third grating is placed separately, and is sensitive only to temperature variations. We theoretically analyze the sensor for measuring pressure and acceleration simultaneously while providing temperature compensation. The analysis shows that there are distinct shifts in the Bragg wavelength when measuring pressure and acceleration.

Figure 1 gives the schematic representation of the proposed temperature-independent multiparametric sensor. It consists of three identical Bragg gratings (Grating 1, Grating 2, and Grating 3) in a superstructure configuration, with distinct Bragg wavelengths λ1, λ2, and λ3 respectively. Grating 1 is located on the square diaphragm, and Grating 2 is placed on the cantilever beam with a proof mass to obtain high sensitivity, whereas Grating 3 is located away from the diaphragm and cantilever beam. The distance between gratings is 500 μm.

Figure 1.Schematic diagram of a temperature-compensated, multiparametric WBG sensor.

When pressure and acceleration are respectively applied to the diaphragm and the cantilever beam, the mechanical displacement induces stress in the sensing structure. The effective refractive index (RI) and pitch of the grating varies as a result of this induced stress, which in turn leads to the shifts in the Bragg wavelengths of Grating 1 and Grating 2. Grating 3, however, is unaffected by the applied pressure and acceleration. Although the variation in temperature results in a Bragg-reflection-wavelength shift for all three gratings, this shift is equal, as the temperature change is the same for all the three gratings. While Grating 1 is affected by both pressure and temperature, Grating 3 is affected only by temperature. The difference between the two shifts eliminates the shift that is due purely to temperature. The same concept is applied for acceleration measurement also, wherein the difference of wavelength shifts between Grating 2 and Grating 3 is considered for temperature compensation.

The sensor design includes a superstructure grating, micro-machined diaphragm, cantilever beam, and optomechanical coupling between them.

3.1. Superstructure Surface-relief Waveguide Bragg Grating

The superstructure configuration consists of three surface-relief Bragg gratings connected by a waveguide between each pair, as shown in Fig. 2.

Figure 2.Superstructure pattern of a surface-relief WBG.

The surface-relief WBG consists of a silicon-slab wave-guide, where the top surface is etched to obtain the periodic structure. Light propagates in the z-direction. The core material is silicon, and the cladding is air. The height of the waveguide is t, the grating depth is tg, the length of the grating is L, and d is the ridge width, which is half of the grating pitch. Pitch can be calculated by using Bragg’s equation:

where λB, neff, and Λ are Bragg’s wavelength, the effective RI, and the pitch of the grating respectively. The transmitted and reflected fields are represented by a and b for light incident from the left, with the respective boundaries at z = 0 and z = L. Only the fundamental mode and first-order diffraction are considered. Coupled-mode theory is used to analyze the WBG [14, 15].

As light propagates along a grating, partial reflections occur at each of the etchings, due to the surface etching of the waveguide. At the Bragg condition, the forward wave (transmitted field) and backward wave (reflected field) are in phase. As a result, coupling occurs between the forward and backward traveling waves. The coupling between them is determined by utilizing the coupled-mode equations:

where a and b respectively represent the amplitudes of the forward and backward traveling waves, and Δβ is calculated as [14, 15]

where β is the propagation constant of the waveguide, k = (2π)/λ, and KC is the coupling coefficient, given by

The solution of the coupled-mode equations is given as [14, 15]

where A(0) is the input amplitude and α is given as

The transmittance T and reflectance R of the Bragg grating are

The superstructure configuration consists of five segments: Grating 1 of length L1 and Bragg wavelength 1550 nm, Grating 2 of length L3 and Bragg wavelength 1554 nm, Grating 3 of length L5 and Bragg wavelength 1557 nm, and two waveguides of length L2 and L4 respectively. It is analyzed by the transfer-matrix method for multiple segments [16]:

where the index i refers to the ith segment of the superstructure grating. In this case, the second and fourth segments are only a waveguide, so there is no coupling between forward and backward waves; only the forward wave is transmitted. Hence the coefficients M2 and M4 are adjusted to replicate a waveguide structure.

3.2. Micromachined Diaphragm

In this section, pressure measurements for two different diaphragm structures, circular and square, are analyzed.

3.2.1. Circular diaphragm

A circular diaphragm with radius f and thickness h is clamped at the edge where r = f, as shown in Fig. 3.

Figure 3.Schematic diagram of a circular diaphragm.

The displacement of the diaphragm under pressure p is given by [17]

D is the flexural rigidity of the diaphragm, given by

where E and v represent the Young’s modulus and Poisson’s ratio of the diaphragm respectively.

The stresses on the diaphragm are

where subscripts r and θ indicate the radial and tangential directions respectively.

The strain components are related to stress as

Figure 4 shows the stress distribution on the circular diaphragm, when simulated with the COMSOL Multiphysics software package.

Figure 4.Stress distribution due to applied pressure on the surface of the circular diaphragm, using COMSOL Multiphysics.

The maximum stress of 6.85 MPa is induced the along the perimeter of the circular diaphragm, indicating that this is the ideal location for the grating.

3.2.2. Square diaphragm

A square diaphragm with side length of 2a and thickness h under applied pressure is shown in Fig. 5.

Figure 5.Schematic representation of a square diaphragm.

The displacement undergone by the diaphragm under pressure p is [17]

Displacement of the diaphragm is linearly dependent on the applied pressure, and also varies inversely with the cube of the thickness, so ideally the thickness can be reduced to achieve increased sensitivity (though this also increases the chance of damage to the structure).

The stress distribution across the diaphragm’s surface along the z and y axes in Cartesian coordinates zy is

The corresponding strain components are

Stress variation along the length of the square diaphragm is shown in Fig. 6. The maximum stress is found at the edges of the diaphragm.

Figure 6.Stress distribution along the length of the square diaphragm.

Figure 7 shows the stress distribution on the square diaphragm, as simulated with COMSOL Multiphysics.

The maximum stress of 22 MPa is induced at the center of each edge of the diaphragm, as seen in Fig. 7. The significantly higher induced stress seen in the square diaphragm makes it a far more effective candidate for pressure sensing than the circular diaphragm, as a higher pressure sensitivity can be obtained. As a consequence, henceforth a square diaphragm is considered for pressure measurement, by placing a grating at the edge of the structure.

Figure 7.Stress distribution due to applied pressure on the surface of a diaphragm, using COMSOL Multiphysics.

With applied stress there is a proportional change in the effective RI and pitch of the grating. Thus Grating 1 in the superstructure shown in Fig. 1 is placed at the edge of the diaphragm, such that it undergoes a larger change in Bragg wavelength.

3.3. Cantilever Beam

A cantilever beam is shown in Fig. 8, with length, width, and thickness of a1, b1, and h1 respectively. The length, width, and thickness of the proof mass are (a2 - a1), b2, and h2 respectively.

Figure 8.Schematic diagram of a cantilever beam.

The acceleration is applied in the negative x direction, and the resulting stress is given by [17]

where m is the mass of the proof mass, A is the acceleration, and L is The maximum stress occurs at z = 0, so Grating 2 in the superstructure is placed at this position.

The strain component, which is related to stress, is

The variation in stress along the length of the cantilever beam is plotted in Fig. 9. The maximum stress on the beam is at z = 0.

Figure 9.Stress distribution for the cantilever beam along its length.

Figure 10 shows the stress distribution for the cantilever beam, as simulated using COMSOL Multiphysics.

As seen in Fig. 10, the maximum stress is at the anchor point, so this location is ideal for placement of the grating.

Figure 10.Stress distribution due to applied acceleration on the surface of the cantilever beam, using COMSOL Multiphysics.

Even though the acceleration is applied along the normal direction (x direction), its effect is also seen in the other directions y and z, and the resulting stresses are shown in Eqs. (26) and (27). This is referred to as cross-axis sensitivity.

When comparing the stress in the z direction to the maximum stress in the normal direction, we get Tz = (h1/(6L))T. This gives a very small value of Tz, due to low thickness and large proof mass length, and hence it can be neglected. On the contrary, the stress in the y direction is given by Ty = (h1/b1)T; this component cannot be neglected, and actually could be considered for some practical applications. This component can be reduced by using a beam of smaller thickness or larger width.

Based on the application, the fundamental frequency of acceleration can be chosen using Eq. (28):

The design parameters used for the structures are outlined in Table 1.

TABLE 1. Design Parameters for the diaphragm and cantilever structures


3.4. Optomechanical Coupling

Under applied pressure and acceleration, stress is induced in both the diaphragm and cantilever beam. This stress causes a change in the effective RI and pitch [18]. Eqs. (29) and (30) respectively quantify these changes in effective RI and pitch for TE-mode (Transverse Electric) analysis of the waveguide on the diaphragm.

where p11 and p12 are the stress-optic coefficients. The changes in effective RI and pitch of the grating in the case of the cantilever beam are given by

3.5. Temperature Compensation

In a Bragg grating, the wavelength shift depends on both strain and temperature. As a result, it is not possible to measure the pressure or acceleration directly by utilizing a single grating and observing the wavelength shift; it is necessary to compensate for the effect of temperature. In the proposed sensor, we have used another grating on the same waveguide, protected against pressure and acceleration yet at the same temperature as its neighbors. We have designed the superstructure configuration with three surface-relief gratings on the same waveguide and providing reflections at three distinct wavelengths. The wavelength shift in the reflection spectrum in Grating 1 depends on the applied pressure and temperature; similarly, for Grating 2 the shift depends on the acceleration as well as the temperature. On the other hand, the shift for Grating 3 depends on temperature only; as a result, Grating 3 can be used to compensate for the temperature-based shift.

The wavelength shift in the reflection of a grating when it undergoes a change in strain caused by the application of pressure and temperature is given by [19]

where Kp is the pressure sensitivity coefficient and KT is the temperature sensitivity coefficient. Δp and ΔT represent the variation in pressure and temperature respectively. Similarly, the wavelength shift of the grating incurred as a result of strain caused by any applied acceleration and temperature is expressed as

where Kacc is the acceleration sensitivity coefficient and ΔA is the change in acceleration.

In our case, we have considered a superstructure of three identical gratings on the same waveguide, but having reflections at three different Bragg wavelengths. The corresponding Bragg-wavelength shifts are given by

Since the material and the waveguide are common to all the three gratings, we can write

The difference between the shifts in wavelength for Grating 1 and Grating 3 cancels out any effect of temperature, and we obtain the shift due to the pressure applied on the diaphragm alone. Similarly, the difference between Grating 2 and Grating 3 wavelength shifts is taken to effectively cancel any temperature effect and capture the acceleration applied to the beam.

Simulations are performed using MATLAB and COMSOL Multiphysics for cases both with and without applied pressure and acceleration. The power-spectral analysis of reflected light for the superstructure WBG sensor without applying any pressure or acceleration is carried out and presented in Fig. 11. Three Bragg-wavelength reflections at 1550 nm, 1553 nm, and 1557 nm were observed, due to the three gratings of different pitch values 348.5 nm, 349.4 nm, and 350 nm respectively.

Figure 11.Reflection spectrum of the superstructure WBG in the steady state.

Figure 12 shows the shifts in the Bragg wavelengths Δλ1 and Δλ2 due to simultaneously applied pressure of 7 kPa and acceleration of 2 g on the sensor. Here 1 g equals 9.8 m/s2. Since there is no variation in temperature, there is no shift in Bragg wavelength of Grating 3; hence overlap of the spectra at its Bragg wavelength is observed.

Figure 12.Reflection spectrum of the superstructure WBG after application of pressure and acceleration.

Figures 13 and 14 show linear variation of Bragg wavelength with pressure and acceleration respectively. Additionally, the pressure sensitivity was obtained as 0.21 pm/Pa, and the sensitivity to acceleration was found to be 6.49 nm/g.

Figure 13.Change in Bragg wavelength with applied pressure.
Figure 14.Change in Bragg wavelength with applied acceleration.

We have theoretically analyzed and simulated a temperature-compensated monolithic WBG-based multiparametric sensor for simultaneous measurement of pressure and acceleration, using a unique mathematical model. Three Bragg gratings are used in a superstructure configuration, in which two gratings are used for pressure and acceleration measurement respectively and the remaining one for temperature compensation. A square diaphragm and cantilever beam are designed to obtain high sensitivity. A pressure sensitivity of 0.21 pm/Pa and an acceleration sensitivity of 6.49 nm/g are obtained. Since the fabrication protocols for the proposed superstructured WBG-based sensor fall in line with standard CMOS foundry processes, the device could be integrated at chip scale, with source and detectors on a single chip as a photonic integrated circuit. Thus such highly sensitive devices could be an affordable solution for problems where simultaneous measurement of multiple parameters is required.

  1. A. D. Kersey, M. A. Davis, H. J. Patrick, M. LeBlanc, K. P. Koo, C. G. Askins, M. A. Putnam, and E. J. Friebele, "Fiber grating sensors," J. Lightwave Technol. 15, 1442-1463 (1997).
    CrossRef
  2. Y. Zhang, D. Feng, Z. Liu, Z. Guo, X. Dong, K. S. Chiang, and B. C. B. Chu, "High-sensitivity pressure sensor using a shielded polymer-coated fiber Bragg grating," IEEE Photonics Technol. Lett. 13, 618-619 (2001).
    CrossRef
  3. H. J. Sheng, W. F. Liu, T. C. Chen, S. S. Bor, and M. Y. Fu, "A lateral pressure sensor using a fiber Bragg grating," in Proc. Pacific Rim Conference on Lasers and Electro-Optics (CLEO) (Taipei, Taiwan, Dec. 2003), 2, pp. 674.
  4. A. R. Sankar, J. G. Jency, and S. Das, "Design, fabrication and testing of a high-performance silicon piezoresistive Z-axis accelerometer with proof mass-edge-aligned-flexures," Microsyst. Technol. 18, 9-23 (2012).
    CrossRef
  5. D. Feng, X. Qiao, H. Yang, Q. Rong, R. Wang, Y. Du, M. Hu, and Z. Feng, "A fiber Bragg grating accelerometer based on a hybridization of cantilever beam," IEEE Sens. J. 15, 1532-1537 (2015).
    CrossRef
  6. J. Dong, Z. J. Long, H. Jiang, and L. Sun, "Monolithic-integrated piezoresistive MEMS accelerometer pressure sensor with glass-silicon-glass sandwich structure," Microsyst. Technol. 23, 1563-1574 (2017).
    CrossRef
  7. J. Xu, Y. Zhao, Z. Jiang, and J. Sun, "A monolithic silicon multi-sensor for measuring three-axis acceleration, pressure and temperature," J. Mech. Sci. Technol. 22, 731-739 (2008).
    CrossRef
  8. P. K. Pattnaik, B. Vijayaaditya, T. Srinivas, and A. Selvarajan, "Optical MEMS pressure and vibration sensors using integrated optical ring resonators," in Proc. SENSORS (Irvine, CA, USA, Nov. 2005), pp. 636-639.
  9. G. N. D. Brabander, J. T. Boyd, and G. Beheim, "Integrated optical ring resonator with micromechanical diaphragm for pressure sensing," IEEE Photonics Technol. Lett. 6, 671-673 (1994).
    CrossRef
  10. C. Thondagere, A. Kaushalram, T. Srinivas, and G. Hegde, "Mathematical modeling of optical MEMS differential pressure sensor using waveguide Bragg gratings embedded in Mach Zehnder interferometer," J. Opt. 20, 085802 (2018).
    CrossRef
  11. É. Pinet, "Pressure measurement with fiber-optic sensors: commercial technologies and applications," Proc. SPIE 7753, 775304 (2011).
  12. S. Kim, J. Kwon, S. Kim, and B. Lee, "Temperature-independent strain sensor using a chirped grating partially embedded in a glass tube," IEEE Photonics Technol. Lett. 12, 678-680 (2000).
    CrossRef
  13. V. Neeharika, and P. K. Pattnaik, "Optical MEMS pressure sensors incorporating dual waveguide Bragg gratings on diaphragms," IEEE Sens. J. 16, 681-687 (2016).
    CrossRef
  14. C. R. Pollock, and M. Lipson. Integrated Photonics (Springer US 2003) p. 241-269.
  15. A. Yariv, "Coupled-mode theory for guided-wave optics," IEEE J. Quantum Electron. 9, 919-933 (1973).
    CrossRef
  16. A. Yariv, and P. Yeh. Photonics: Optical Electronics in Modern Communications (Oxford Series in Electrical and Computer Engineering), 6th ed., Sedra A. S. (Oxford University Press, NY, USA 2007) p. 539-601.
  17. M. H. Bao. Micro Mechanical Transducers: Pressure Sensors, Accelerometers and Gyroscopes (Handbook of Sensors and Actuators Series), Middelhoek S. (Elsevier, The Netherlands 2000), 8 p. 23-88.
  18. Holmes C.. Direct UV written planar devices for sensing and telecommunication applications, Ph. D. Dissertation (University of Southampton, Southampton 2009).
  19. M. M. Werneck, R. C. S. B. Allil, B. A. Ribeiro, and F. V. B. D. Nazare. Current Trends in Short-and Long-period Fiber Gratings, Laborde C. C. (InTech, Rijeka, Croatia 2013) p. 2-24.

Article

Article

Current Optics and Photonics 2020; 4(4): 293-301

Published online August 25, 2020 https://doi.org/10.3807/COPP.2020.4.4.293

Copyright © Optical Society of Korea.

Chip-scale Temperature-compensated Superstructured Waveguide Bragg Grating Based Multiparametric Sensor

Naik Parrikar Vishwaraj1,*, Chandrika Thondagere Nataraj2,**, Ravi Prasad Kogravalli Jagannath3, Prashanth Gurusiddappa1,***, and Srinivas Talabattula4

1Department of Electronics and Communication Engineering, National Institute of Technology Goa, Ponda 403401, India, 2Department of Electronics and Telecommunication Engineering, Siddaganga Institute of Technology, Tumakuru-Karnataka 572103, India, 3Department Applied Sciences, National Institute of Technology Goa, Ponda 403401, India, 4Department of Electrical Communication Engineering, Indian Institute of Science, Bangalore-Karnataka 560012, India

Correspondence to:rajparrikar@gmail.com

Received: December 26, 2019; Revised: February 27, 2020; Accepted: March 16, 2020

Abstract

In this paper we propose and theoretically analyze a monolithic multiparametric sensor consisting of a superstructure of surface-relief waveguide Bragg gratings (WBGs), a micro-machined diaphragm, and a cantilever beam. Diaphragms of two different configurations, namely circular and square, are designed and analyzed separately for pressure measurement. The square diaphragm is then selected for further study, since it shows relatively higher sensitivity compared to the circular one, as it incurs more induced stress when any pressure is applied. The cantilever beam with a proof mass is designed to enhance the sensitivity for acceleration measurement. A unique mathematical method using coupled-mode theory and the transfer-matrix method is developed to design and analyze the shift in the Bragg wavelength of the superstructure configuration of the gratings, due to simultaneously applied pressure and acceleration. The effect of temperature on the wavelength shift is compensated by introducing another Bragg grating in the superstructure configuration. The measured sensitivities for pressure and acceleration are found to be 0.21 pm/Pa and 6.49 nm/g respectively.

Keywords: Waveguide Bragg gratings, Pressure sensor, Accelerometer, Optical sensors, MOEMS

I. INTRODUCTION

In these days, fiber-Bragg-grating (FBG) sensors are used successfully for a broad spectrum of applications in automobiles, aerospace, and biosensing due to their multiplexed sensing capabilities, lightweight design, and compactness [1-5]. Measurements of pressure, temperature, and acceleration are critical in industrial applications, and subsequently there is high demand for integrated composite sensors [6, 7]. However, these are based on the piezoresistive principle of sensing, in which piezoresistive elements are placed over a diaphragm and cantilever beam to measure pressure and acceleration respectively. Compact, planar photonic components are being considered for such applications, due to their potential for large-scale integration and reduced fabrication cost.

Micro-opto-electro-mechanical systems (MOEMS), realized by combining silicon photonics with micro-electro-mechanical systems (MEMS), offer the advantages of a small device footprint, cheapness, batch-fabrication possibilities, and easier CMOS integration with less complexity. MOEMS pressure sensors using a Mach-Zehnder interferometer (MZI), ring resonator, directional coupler, and gratings have been reported [3, 8-10].

Waveguide Bragg gratings (WBGs) on the silicon-on-insulator (SOI) platform have attracted interest for optical communication and various sensing applications. WBGs were first demonstrated by Murphy et al. [1] in 2001. The Bragg gratings are realized by physically corrugating dielectric waveguides; their sensitivity is obtained by observing the shift in the reflection wavelength. This method of detection, using the shift in wavelength, is preferred over methods that utilize the change in intensity of propagated light, as the former is more immune to noise [11]. However, Bragg-grating-based sensors are equally sensitive to temperature as to the parameters being measured. Hence it remains a difficult task to discriminate between the shift in wavelength due to temperature change and that due to the parameter being measured. If the shift due to temperature is not compensated, a MOEMS sensor gives inaccurate results, so it is necessary to subtract the wavelength shift due to temperature effects, so that we can attribute the wavelength shift only to changes in the actual parameter being measured, such as pressure or acceleration. In earlier works a chirped grating and dual gratings have been proposed to compensate for the effect of temperature in pressure sensors [12, 13]. However, surface-relief waveguide Bragg gratings for monolithic multiparametric sensing with temperature compensation have not been addressed so far. In addition, so far there has been no unique mathematical model to analyze such systems theoretically. Such systems could be an affordable solution for automobile, aerospace, and health applications, since they could be easily manufactured using standard microfabrication technologies.

In this paper we propose a WBG-based monolithic multiparametric sensor that is immune to temperature aberrations, for both pressure and acceleration measurements. A unique mathematical model using coupled-mode theory and the transfer-matrix method is developed to design the configuration of the superstructure of surface-relief Bragg gratings. The proposed sensor is made up of a superstructure configuration of three identical surface-relief waveguide Bragg gratings, designed with distinct Bragg wavelengths. Two of these gratings are located respectively on a micromachined diaphragm structure and a cantilever beam structure, where pressure and acceleration are applied. The wavelength shifts of these gratings are proportional to both pressure and acceleration, along with temperature. The third grating is placed separately, and is sensitive only to temperature variations. We theoretically analyze the sensor for measuring pressure and acceleration simultaneously while providing temperature compensation. The analysis shows that there are distinct shifts in the Bragg wavelength when measuring pressure and acceleration.

II. SENSOR CONFIGURATION

Figure 1 gives the schematic representation of the proposed temperature-independent multiparametric sensor. It consists of three identical Bragg gratings (Grating 1, Grating 2, and Grating 3) in a superstructure configuration, with distinct Bragg wavelengths λ1, λ2, and λ3 respectively. Grating 1 is located on the square diaphragm, and Grating 2 is placed on the cantilever beam with a proof mass to obtain high sensitivity, whereas Grating 3 is located away from the diaphragm and cantilever beam. The distance between gratings is 500 μm.

Figure 1. Schematic diagram of a temperature-compensated, multiparametric WBG sensor.

When pressure and acceleration are respectively applied to the diaphragm and the cantilever beam, the mechanical displacement induces stress in the sensing structure. The effective refractive index (RI) and pitch of the grating varies as a result of this induced stress, which in turn leads to the shifts in the Bragg wavelengths of Grating 1 and Grating 2. Grating 3, however, is unaffected by the applied pressure and acceleration. Although the variation in temperature results in a Bragg-reflection-wavelength shift for all three gratings, this shift is equal, as the temperature change is the same for all the three gratings. While Grating 1 is affected by both pressure and temperature, Grating 3 is affected only by temperature. The difference between the two shifts eliminates the shift that is due purely to temperature. The same concept is applied for acceleration measurement also, wherein the difference of wavelength shifts between Grating 2 and Grating 3 is considered for temperature compensation.

III. THEORY AND DESIGN

The sensor design includes a superstructure grating, micro-machined diaphragm, cantilever beam, and optomechanical coupling between them.

3.1. Superstructure Surface-relief Waveguide Bragg Grating

The superstructure configuration consists of three surface-relief Bragg gratings connected by a waveguide between each pair, as shown in Fig. 2.

Figure 2. Superstructure pattern of a surface-relief WBG.

The surface-relief WBG consists of a silicon-slab wave-guide, where the top surface is etched to obtain the periodic structure. Light propagates in the z-direction. The core material is silicon, and the cladding is air. The height of the waveguide is t, the grating depth is tg, the length of the grating is L, and d is the ridge width, which is half of the grating pitch. Pitch can be calculated by using Bragg’s equation:

where λB, neff, and Λ are Bragg’s wavelength, the effective RI, and the pitch of the grating respectively. The transmitted and reflected fields are represented by a and b for light incident from the left, with the respective boundaries at z = 0 and z = L. Only the fundamental mode and first-order diffraction are considered. Coupled-mode theory is used to analyze the WBG [14, 15].

As light propagates along a grating, partial reflections occur at each of the etchings, due to the surface etching of the waveguide. At the Bragg condition, the forward wave (transmitted field) and backward wave (reflected field) are in phase. As a result, coupling occurs between the forward and backward traveling waves. The coupling between them is determined by utilizing the coupled-mode equations:

where a and b respectively represent the amplitudes of the forward and backward traveling waves, and Δβ is calculated as [14, 15]

where β is the propagation constant of the waveguide, k = (2π)/λ, and KC is the coupling coefficient, given by

The solution of the coupled-mode equations is given as [14, 15]

where A(0) is the input amplitude and α is given as

The transmittance T and reflectance R of the Bragg grating are

The superstructure configuration consists of five segments: Grating 1 of length L1 and Bragg wavelength 1550 nm, Grating 2 of length L3 and Bragg wavelength 1554 nm, Grating 3 of length L5 and Bragg wavelength 1557 nm, and two waveguides of length L2 and L4 respectively. It is analyzed by the transfer-matrix method for multiple segments [16]:

where the index i refers to the ith segment of the superstructure grating. In this case, the second and fourth segments are only a waveguide, so there is no coupling between forward and backward waves; only the forward wave is transmitted. Hence the coefficients M2 and M4 are adjusted to replicate a waveguide structure.

3.2. Micromachined Diaphragm

In this section, pressure measurements for two different diaphragm structures, circular and square, are analyzed.

3.2.1. Circular diaphragm

A circular diaphragm with radius f and thickness h is clamped at the edge where r = f, as shown in Fig. 3.

Figure 3. Schematic diagram of a circular diaphragm.

The displacement of the diaphragm under pressure p is given by [17]

D is the flexural rigidity of the diaphragm, given by

where E and v represent the Young’s modulus and Poisson’s ratio of the diaphragm respectively.

The stresses on the diaphragm are

where subscripts r and θ indicate the radial and tangential directions respectively.

The strain components are related to stress as

Figure 4 shows the stress distribution on the circular diaphragm, when simulated with the COMSOL Multiphysics software package.

Figure 4. Stress distribution due to applied pressure on the surface of the circular diaphragm, using COMSOL Multiphysics.

The maximum stress of 6.85 MPa is induced the along the perimeter of the circular diaphragm, indicating that this is the ideal location for the grating.

3.2.2. Square diaphragm

A square diaphragm with side length of 2a and thickness h under applied pressure is shown in Fig. 5.

Figure 5. Schematic representation of a square diaphragm.

The displacement undergone by the diaphragm under pressure p is [17]

Displacement of the diaphragm is linearly dependent on the applied pressure, and also varies inversely with the cube of the thickness, so ideally the thickness can be reduced to achieve increased sensitivity (though this also increases the chance of damage to the structure).

The stress distribution across the diaphragm’s surface along the z and y axes in Cartesian coordinates zy is

The corresponding strain components are

Stress variation along the length of the square diaphragm is shown in Fig. 6. The maximum stress is found at the edges of the diaphragm.

Figure 6. Stress distribution along the length of the square diaphragm.

Figure 7 shows the stress distribution on the square diaphragm, as simulated with COMSOL Multiphysics.

The maximum stress of 22 MPa is induced at the center of each edge of the diaphragm, as seen in Fig. 7. The significantly higher induced stress seen in the square diaphragm makes it a far more effective candidate for pressure sensing than the circular diaphragm, as a higher pressure sensitivity can be obtained. As a consequence, henceforth a square diaphragm is considered for pressure measurement, by placing a grating at the edge of the structure.

Figure 7. Stress distribution due to applied pressure on the surface of a diaphragm, using COMSOL Multiphysics.

With applied stress there is a proportional change in the effective RI and pitch of the grating. Thus Grating 1 in the superstructure shown in Fig. 1 is placed at the edge of the diaphragm, such that it undergoes a larger change in Bragg wavelength.

3.3. Cantilever Beam

A cantilever beam is shown in Fig. 8, with length, width, and thickness of a1, b1, and h1 respectively. The length, width, and thickness of the proof mass are (a2 - a1), b2, and h2 respectively.

Figure 8. Schematic diagram of a cantilever beam.

The acceleration is applied in the negative x direction, and the resulting stress is given by [17]

where m is the mass of the proof mass, A is the acceleration, and L is The maximum stress occurs at z = 0, so Grating 2 in the superstructure is placed at this position.

The strain component, which is related to stress, is

The variation in stress along the length of the cantilever beam is plotted in Fig. 9. The maximum stress on the beam is at z = 0.

Figure 9. Stress distribution for the cantilever beam along its length.

Figure 10 shows the stress distribution for the cantilever beam, as simulated using COMSOL Multiphysics.

As seen in Fig. 10, the maximum stress is at the anchor point, so this location is ideal for placement of the grating.

Figure 10. Stress distribution due to applied acceleration on the surface of the cantilever beam, using COMSOL Multiphysics.

Even though the acceleration is applied along the normal direction (x direction), its effect is also seen in the other directions y and z, and the resulting stresses are shown in Eqs. (26) and (27). This is referred to as cross-axis sensitivity.

When comparing the stress in the z direction to the maximum stress in the normal direction, we get Tz = (h1/(6L))T. This gives a very small value of Tz, due to low thickness and large proof mass length, and hence it can be neglected. On the contrary, the stress in the y direction is given by Ty = (h1/b1)T; this component cannot be neglected, and actually could be considered for some practical applications. This component can be reduced by using a beam of smaller thickness or larger width.

Based on the application, the fundamental frequency of acceleration can be chosen using Eq. (28):

The design parameters used for the structures are outlined in Table 1.

TABLE 1.. Design Parameters for the diaphragm and cantilever structures.


3.4. Optomechanical Coupling

Under applied pressure and acceleration, stress is induced in both the diaphragm and cantilever beam. This stress causes a change in the effective RI and pitch [18]. Eqs. (29) and (30) respectively quantify these changes in effective RI and pitch for TE-mode (Transverse Electric) analysis of the waveguide on the diaphragm.

where p11 and p12 are the stress-optic coefficients. The changes in effective RI and pitch of the grating in the case of the cantilever beam are given by

3.5. Temperature Compensation

In a Bragg grating, the wavelength shift depends on both strain and temperature. As a result, it is not possible to measure the pressure or acceleration directly by utilizing a single grating and observing the wavelength shift; it is necessary to compensate for the effect of temperature. In the proposed sensor, we have used another grating on the same waveguide, protected against pressure and acceleration yet at the same temperature as its neighbors. We have designed the superstructure configuration with three surface-relief gratings on the same waveguide and providing reflections at three distinct wavelengths. The wavelength shift in the reflection spectrum in Grating 1 depends on the applied pressure and temperature; similarly, for Grating 2 the shift depends on the acceleration as well as the temperature. On the other hand, the shift for Grating 3 depends on temperature only; as a result, Grating 3 can be used to compensate for the temperature-based shift.

The wavelength shift in the reflection of a grating when it undergoes a change in strain caused by the application of pressure and temperature is given by [19]

where Kp is the pressure sensitivity coefficient and KT is the temperature sensitivity coefficient. Δp and ΔT represent the variation in pressure and temperature respectively. Similarly, the wavelength shift of the grating incurred as a result of strain caused by any applied acceleration and temperature is expressed as

where Kacc is the acceleration sensitivity coefficient and ΔA is the change in acceleration.

In our case, we have considered a superstructure of three identical gratings on the same waveguide, but having reflections at three different Bragg wavelengths. The corresponding Bragg-wavelength shifts are given by

Since the material and the waveguide are common to all the three gratings, we can write

The difference between the shifts in wavelength for Grating 1 and Grating 3 cancels out any effect of temperature, and we obtain the shift due to the pressure applied on the diaphragm alone. Similarly, the difference between Grating 2 and Grating 3 wavelength shifts is taken to effectively cancel any temperature effect and capture the acceleration applied to the beam.

IV. RESULTS AND DISCUSSION

Simulations are performed using MATLAB and COMSOL Multiphysics for cases both with and without applied pressure and acceleration. The power-spectral analysis of reflected light for the superstructure WBG sensor without applying any pressure or acceleration is carried out and presented in Fig. 11. Three Bragg-wavelength reflections at 1550 nm, 1553 nm, and 1557 nm were observed, due to the three gratings of different pitch values 348.5 nm, 349.4 nm, and 350 nm respectively.

Figure 11. Reflection spectrum of the superstructure WBG in the steady state.

Figure 12 shows the shifts in the Bragg wavelengths Δλ1 and Δλ2 due to simultaneously applied pressure of 7 kPa and acceleration of 2 g on the sensor. Here 1 g equals 9.8 m/s2. Since there is no variation in temperature, there is no shift in Bragg wavelength of Grating 3; hence overlap of the spectra at its Bragg wavelength is observed.

Figure 12. Reflection spectrum of the superstructure WBG after application of pressure and acceleration.

Figures 13 and 14 show linear variation of Bragg wavelength with pressure and acceleration respectively. Additionally, the pressure sensitivity was obtained as 0.21 pm/Pa, and the sensitivity to acceleration was found to be 6.49 nm/g.

Figure 13. Change in Bragg wavelength with applied pressure.
Figure 14. Change in Bragg wavelength with applied acceleration.

V. CONCLUSION

We have theoretically analyzed and simulated a temperature-compensated monolithic WBG-based multiparametric sensor for simultaneous measurement of pressure and acceleration, using a unique mathematical model. Three Bragg gratings are used in a superstructure configuration, in which two gratings are used for pressure and acceleration measurement respectively and the remaining one for temperature compensation. A square diaphragm and cantilever beam are designed to obtain high sensitivity. A pressure sensitivity of 0.21 pm/Pa and an acceleration sensitivity of 6.49 nm/g are obtained. Since the fabrication protocols for the proposed superstructured WBG-based sensor fall in line with standard CMOS foundry processes, the device could be integrated at chip scale, with source and detectors on a single chip as a photonic integrated circuit. Thus such highly sensitive devices could be an affordable solution for problems where simultaneous measurement of multiple parameters is required.

Fig 1.

Figure 1.Schematic diagram of a temperature-compensated, multiparametric WBG sensor.
Current Optics and Photonics 2020; 4: 293-301https://doi.org/10.3807/COPP.2020.4.4.293

Fig 2.

Figure 2.Superstructure pattern of a surface-relief WBG.
Current Optics and Photonics 2020; 4: 293-301https://doi.org/10.3807/COPP.2020.4.4.293

Fig 3.

Figure 3.Schematic diagram of a circular diaphragm.
Current Optics and Photonics 2020; 4: 293-301https://doi.org/10.3807/COPP.2020.4.4.293

Fig 4.

Figure 4.Stress distribution due to applied pressure on the surface of the circular diaphragm, using COMSOL Multiphysics.
Current Optics and Photonics 2020; 4: 293-301https://doi.org/10.3807/COPP.2020.4.4.293

Fig 5.

Figure 5.Schematic representation of a square diaphragm.
Current Optics and Photonics 2020; 4: 293-301https://doi.org/10.3807/COPP.2020.4.4.293

Fig 6.

Figure 6.Stress distribution along the length of the square diaphragm.
Current Optics and Photonics 2020; 4: 293-301https://doi.org/10.3807/COPP.2020.4.4.293

Fig 7.

Figure 7.Stress distribution due to applied pressure on the surface of a diaphragm, using COMSOL Multiphysics.
Current Optics and Photonics 2020; 4: 293-301https://doi.org/10.3807/COPP.2020.4.4.293

Fig 8.

Figure 8.Schematic diagram of a cantilever beam.
Current Optics and Photonics 2020; 4: 293-301https://doi.org/10.3807/COPP.2020.4.4.293

Fig 9.

Figure 9.Stress distribution for the cantilever beam along its length.
Current Optics and Photonics 2020; 4: 293-301https://doi.org/10.3807/COPP.2020.4.4.293

Fig 10.

Figure 10.Stress distribution due to applied acceleration on the surface of the cantilever beam, using COMSOL Multiphysics.
Current Optics and Photonics 2020; 4: 293-301https://doi.org/10.3807/COPP.2020.4.4.293

Fig 11.

Figure 11.Reflection spectrum of the superstructure WBG in the steady state.
Current Optics and Photonics 2020; 4: 293-301https://doi.org/10.3807/COPP.2020.4.4.293

Fig 12.

Figure 12.Reflection spectrum of the superstructure WBG after application of pressure and acceleration.
Current Optics and Photonics 2020; 4: 293-301https://doi.org/10.3807/COPP.2020.4.4.293

Fig 13.

Figure 13.Change in Bragg wavelength with applied pressure.
Current Optics and Photonics 2020; 4: 293-301https://doi.org/10.3807/COPP.2020.4.4.293

Fig 14.

Figure 14.Change in Bragg wavelength with applied acceleration.
Current Optics and Photonics 2020; 4: 293-301https://doi.org/10.3807/COPP.2020.4.4.293

TABLE 1. Design Parameters for the diaphragm and cantilever structures


References

  1. A. D. Kersey, M. A. Davis, H. J. Patrick, M. LeBlanc, K. P. Koo, C. G. Askins, M. A. Putnam, and E. J. Friebele, "Fiber grating sensors," J. Lightwave Technol. 15, 1442-1463 (1997).
    CrossRef
  2. Y. Zhang, D. Feng, Z. Liu, Z. Guo, X. Dong, K. S. Chiang, and B. C. B. Chu, "High-sensitivity pressure sensor using a shielded polymer-coated fiber Bragg grating," IEEE Photonics Technol. Lett. 13, 618-619 (2001).
    CrossRef
  3. H. J. Sheng, W. F. Liu, T. C. Chen, S. S. Bor, and M. Y. Fu, "A lateral pressure sensor using a fiber Bragg grating," in Proc. Pacific Rim Conference on Lasers and Electro-Optics (CLEO) (Taipei, Taiwan, Dec. 2003), 2, pp. 674.
  4. A. R. Sankar, J. G. Jency, and S. Das, "Design, fabrication and testing of a high-performance silicon piezoresistive Z-axis accelerometer with proof mass-edge-aligned-flexures," Microsyst. Technol. 18, 9-23 (2012).
    CrossRef
  5. D. Feng, X. Qiao, H. Yang, Q. Rong, R. Wang, Y. Du, M. Hu, and Z. Feng, "A fiber Bragg grating accelerometer based on a hybridization of cantilever beam," IEEE Sens. J. 15, 1532-1537 (2015).
    CrossRef
  6. J. Dong, Z. J. Long, H. Jiang, and L. Sun, "Monolithic-integrated piezoresistive MEMS accelerometer pressure sensor with glass-silicon-glass sandwich structure," Microsyst. Technol. 23, 1563-1574 (2017).
    CrossRef
  7. J. Xu, Y. Zhao, Z. Jiang, and J. Sun, "A monolithic silicon multi-sensor for measuring three-axis acceleration, pressure and temperature," J. Mech. Sci. Technol. 22, 731-739 (2008).
    CrossRef
  8. P. K. Pattnaik, B. Vijayaaditya, T. Srinivas, and A. Selvarajan, "Optical MEMS pressure and vibration sensors using integrated optical ring resonators," in Proc. SENSORS (Irvine, CA, USA, Nov. 2005), pp. 636-639.
  9. G. N. D. Brabander, J. T. Boyd, and G. Beheim, "Integrated optical ring resonator with micromechanical diaphragm for pressure sensing," IEEE Photonics Technol. Lett. 6, 671-673 (1994).
    CrossRef
  10. C. Thondagere, A. Kaushalram, T. Srinivas, and G. Hegde, "Mathematical modeling of optical MEMS differential pressure sensor using waveguide Bragg gratings embedded in Mach Zehnder interferometer," J. Opt. 20, 085802 (2018).
    CrossRef
  11. É. Pinet, "Pressure measurement with fiber-optic sensors: commercial technologies and applications," Proc. SPIE 7753, 775304 (2011).
  12. S. Kim, J. Kwon, S. Kim, and B. Lee, "Temperature-independent strain sensor using a chirped grating partially embedded in a glass tube," IEEE Photonics Technol. Lett. 12, 678-680 (2000).
    CrossRef
  13. V. Neeharika, and P. K. Pattnaik, "Optical MEMS pressure sensors incorporating dual waveguide Bragg gratings on diaphragms," IEEE Sens. J. 16, 681-687 (2016).
    CrossRef
  14. C. R. Pollock, and M. Lipson. Integrated Photonics (Springer US 2003) p. 241-269.
  15. A. Yariv, "Coupled-mode theory for guided-wave optics," IEEE J. Quantum Electron. 9, 919-933 (1973).
    CrossRef
  16. A. Yariv, and P. Yeh. Photonics: Optical Electronics in Modern Communications (Oxford Series in Electrical and Computer Engineering), 6th ed., Sedra A. S. (Oxford University Press, NY, USA 2007) p. 539-601.
  17. M. H. Bao. Micro Mechanical Transducers: Pressure Sensors, Accelerometers and Gyroscopes (Handbook of Sensors and Actuators Series), Middelhoek S. (Elsevier, The Netherlands 2000), 8 p. 23-88.
  18. Holmes C.. Direct UV written planar devices for sensing and telecommunication applications, Ph. D. Dissertation (University of Southampton, Southampton 2009).
  19. M. M. Werneck, R. C. S. B. Allil, B. A. Ribeiro, and F. V. B. D. Nazare. Current Trends in Short-and Long-period Fiber Gratings, Laborde C. C. (InTech, Rijeka, Croatia 2013) p. 2-24.