Abstract

This article proposes a novel microfluidic sensor designed with a highly accurate Q-factor for liquid detection. The proposed sensor is developed and implemented with a gap waveguide cavity resonator (GWCR) approach. The GWCR approach is formed from the two metallic plates denoted as upper and lower plates. These plates are separated by an array of metallic pins attached to the lower plate, leading to high electric field concentration. A microfluidic channel is created at the midpoint of each plate to place the holder of liquid under test (LUT). The GWCR provides a high electric field, which increases Q-factor and is shown to exhibit a significant improvement in sensitivity and linearity. To characterise and evaluate the dielectric properties of the fluid, the LUT is placed inside a hairlike glass, which passes through the microfluidic channels. The LUT perturbs the electric field distribution inside the GWCR, known as the perturbation principle. The relation between the LUT and the electric field changes the electric field behaviours in terms of resonant frequency, Q-factor, and transmission coefficient. The analysis of these changes in the electric field behaviours leads to identifying the dielectric properties of the LUT. The anonymous dielectric characteristics of LUT, permittivity, and loss tangent formulas are derived utilising the polynomial fitting approach. The measurement outcomes reveal that the stated sensor can measure the permittivity and loss tangent for both LUT samples, such as ethanol and methanol, at 6.1 GHz and 23.4°C.

1. Introduction

Developments in liquid characterisation techniques have been grown sharply in the last few years. Microwave sensors have been used to identify the dielectric properties of liquids, which is essential in industries such as those in quality control in food, chemical, pharmaceutical, biosensor, and biomedical industries [13].

Material characterisation using microwave resonant sensors is considered one of the most accurate ways to obtain dielectric properties. Liquid characterisation has been realised using a conventional waveguide, a dielectric and coaxial probe that provides high Q-factor and accuracy [47]. However, these traditional resonators are usually complex and expensive to manufacture [8, 9]. Thus, different methods have been used for liquid characterisation, such as planar resonators [1013]. Planar resonator for petroleum oils characterisation was proposed in [14]. Substrate integrated waveguide (SIW) is another method for liquid detection [1517]. These resonators are simple to design and cheap to fabricate. However, usually, these resonators have low Q-factor, which limit their use in many critical applications. Some researchers tend to improve the Q-factor by using the metal waveguide method [18]. Most researchers used a cavity waveguide which is suitable for liquid sensing [19]. Some waveguide structures were designed at higher frequencies to miniaturise the size and characterise a nanoliter liquid sample [2022]. However, these sensors required high frequency to operate. On the other hand, the complete microfluidic system of the waveguide sensors designed at lower frequencies is complex, large, and expensive to build [23]. A food security sensing system was proposed using a waveguide microwave imaging antenna applied for solid egg detection using a transmitting waveguide antenna (TWA) [24]. The suggested microwave sensors had a dimension of 30 × 30 cm2 and operated at frequency bands 7–13 GHz. Another study was conducted using waveguide antenna sensors for crack detection applications [25].

The gap waveguide (GW) terminology was initially launched in 2009 by Per Simon Kildal [26]. In the GW, the waves propagate between the parallel metal plates. The advantage of the GW compared to planar structures, coplanar waveguide (CPW), and substrate integrated waveguide (SIW) is that it has lower losses. The principle of operation of GW is based on two parallel plates’ structure. Thus one plate acts as a perfect electric conductor (PEC), while the other plate acts as a perfect magnetic conductor (PMC) [26, 27]. A full electromagnetic wave propagates within the pins array boundaries. Hence the electric field is restricted along the pins array area, and this leads to avoiding the leakage of the electric field, which affects traditional technologies at higher frequencies [28]. There are four different configurations of the gap waveguide compositions, as displayed in Figure 1 [29, 30].

In the ridge gap waveguide presented in Figure 1(a), the field propagates through the metal ridge, while in the groove gap waveguide, shown in Figure 1(b), the field propagates through the groove between the parallel plates [28, 31]. The advantage of these two geometries is that they do not need any dielectric material between the plates and can be constructed by creating an array of metal pins under the plate [32]. In the microstrip gap waveguide displayed in Figure 1(c), the field propagates between the upper metal plate and the stripline [33]. In the microstrip gap waveguide exhibited in Figure 1(d), the field propagates along the strip with visas. Gap waveguide technology has gained researchers interest where it has been practised for antennas configuration as proposed in [34]. Moreover, the gap waveguide terminology has been used for millimetre wave applications such as couplers and filters [35]. In this paper, the pothole gap waveguide is selected to design the proposed sensor due to its advantages such as the simplicity, high electric field in the array pins area (pothole), and high Q-factor. These advantages make the groove gap waveguide suitable for liquid characterisation.

A new design of microwave sensors for liquid characterisation is introduced in this study. The novel design of the resonator sensor is formed and rested on the GWCR technique. The new sensor structure is constructed of two parallel plates and a pins array on the lower plate. It is worth mentioning that the materials of the proposed pins are all made from aluminium. The waves are propagating between the two parallel plates and restricted by the pins array. Thus, the new sensor affords a great electric field concentration in the pins array area, leading to a high Q-factor. In order to test the liquid, the LUT samples are installed in between the GWCR where the electric field will be as its maximum concentration utilising a glass capillary; therefore, a sensation phenomenon within the electric field and the LUT is acknowledged as perturbation theory. The stated resonator is verified through laboratory measurement and compared with other related literature studies. The corresponding equation between the permittivity and loss tangent is obtained based on the polynomial fitting approach. Lastly, the measurement outcomes indicated that the recommended sensor could detect ethanol and methanol efficiently.

2. Sensor Design and Configuration

The GWCR can be represented by two parallel plates and RLC equivalent circuit, as presented in Figure 2.

The variables of the equivalent circuit can be calculated according to the length of the plate (l), and the interval separating the two parallel plates (z) can be estimated using equations (1)–(4) [36]:

The air filled the space between the two parallel plates; thus and .

The modelled sensor is proposed based on GW technology as suggested in [26]. The working mechanism of the GWCR depends on the upper and lower plates, and it is presented as follows: the upper plate serves as a perfect conductor to induce the magnetic conductor, which can also be defined as a perfect dielectric, whereas the lower plate works as a perfect conductor in terms of electric conductor. The artificial magnetic conductor (AMC) is produced by the surface of the pins array on the lower plate, which restricts the waves propagation within the pins array area. These waves are larger for more petite holes among the pins and the upper plate [28]. In the GWCR, distribution takes place in the region enclosed by the pins on all sides. The GWCR case supports a TE10 mode of the rectangular waveguide, as presented in Figure 3. The GWCR has less loss and high distribution of the electric field. Thus the GWCR is suitable for liquid characterisation.

Electric (E) and magnetic (H) fields of the GWCR can be described by Maxwell’s equations [37, 38]:where ω, B, and M are the radian frequency, magnetic flux density, and magnetic current density, respectively. The magnetic field can be expressed bywhere D and J are the electric flux density and the electric current density, respectively; the electric flux density and magnetic flux density can be expressed bywhere ρe and ρm are the electric and the magnetic charge density.

For TE modes, a natural boundary condition is described by the expressionwhere n stands for a mode.

For harmonic TE mode wave is described by the expression:where Zo is the wave impedance; at the same time, (5) can be expressed by

In the microwave sensor design, the electric field is significantly considered to ensure a sufficient interaction between the LUT and the electric field of the resonator; thus, the sensor will be able to detect the dielectric characteristics of the material. The electronic concentration in the GWCR can be measured by the Q-factor. Therefore, the high Q-factor reflects the high resolution [39]. The loaded Q-factor is equal to half unloaded Q-factor, where the unloaded Qu-factor (Qu) is calculated according to the 3 dB bandwidth (BW) and resonant frequency (f) using the following equation [40]:

The GWCR is studied and examined at 6.1 GHz, employing CST software. Figure 4 shows that the proposed structure consists of two parallel plates fed by two SMA connectors attached to the upper plate and a pins array on the lower plate. The pin sizes are standard of 1 mm2 [31], and the gap between the upper cover of the pins and the top plate is found to be 1 mm, whereas the proportional thickness of the pin (h) is determined by [26]where corresponds to the wavelength, which can be defined using the equation below:where c denotes the light speed, which is equivalent to 3 × 108 m/s, the estimated gaping among the pins is greater than 1.5 mm, and the height of each plate (thickness) is larger than 5 mm regarding the CNC device abilities. Table 1 presents the progressive values of the intended GWCR.

The waves are propagated between the two parallel plates and confined by the array pins from all directions to avoid the leakage of the electric field. Thus the electric field increases in the array pins area, as displayed in Figure 5. A glass capillary with a thickness of 0.21 mm is positioned in the centre of the cavity to characterise the liquid specimen, as illustrated in Figure 6. The electric field of the GWCR interacts with the LUT in the sensing area, and this manages to change the behaviour of the electric field due to the dielectric properties of the LUT. Hence it can be assumed that the relationship between the electric field and the dielectric properties is linear. The shifting of the resonance frequency depends on the reaction between the LUT and the electric field distribution of the sensor. The changes in the behaviour of the electric field are investigated and analysed to identify the dielectric properties of the LUT.

3. Results and Discussion

The proposed sensor’s resonant frequency is controlled by changing the height of the pins on the lower plate, as illustrated in Figure 7. It can be remarked that the resonant frequency can be increased by reducing the height of the pins and vice versa. Therefore, when the pins’ height is 12.7 mm, the resonant frequency occurs at 6.2824 GHz, while when the height of the pins is 18.7 mm, the resonant frequency rests at 5.9872 GHz.

The proposed sensor is fabricated using an aluminium alloy with an electric conductivity of σ = 3.56 × 107. The simulation conditions of the modelled sensor are listed in Table 2. The measured outcomes reveal that the advanced sensor produces an unloaded Q-factor of 2650 at 6.0950 GHz. To validate the proposed sensor, the measurement is carried out for ethanol and methanol at room temperature of 23.4°C. Figure 8 presents the predicted and measured transmission coefficient (S21) of the designed resonator with/without LUT. From Figure 8, it is clear that there is a slight difference between simulation and measurement where the unloaded resonator for the measured one appears at 6.0950 GHz. This shifting in S21 is because of the fabrication tolerance, variation of substrate’s dielectric constant, and the difference of the surrounding environment between the simulation and measurement. The noise that appears in the measured curve is due to external effects such as microwave signals that can affect the performance of the VNA [10].

In comparison, the simulated one seems at 6.1000 GHz, which is due to the manufacturing tolerance.

Based on the measurement results, the polynomial diagram is schemed based on the permittivity and frequency, as exhibited in Figure 9. The permittivity equation is obtained from the curve as follows:

On the other hand, the permittivity equation can be extracted using matrix. By expressing the relationship between variables in the form of an equation:

Thus, the three data points in Figure 9 can be matched with the x and y coordinates, where the frequency matches with the x coordinate and the permittivity matches with the y coordinate; thus:

This gives a system of three equations and three unknowns as suggested:

Using Cramer’s rule to solve the linear system of equations to find the regression coefficients using the determinants of the square matrix, each of the coefficients may be determined using the following equation:where k = N-1, and N = 3 (data points). Mi is the matrix with the ith column replaced with the column of the permittivity. For example:

Therefore, . Thus, the polynomial coefficient is as follows:

The same goes with a1 and a2, and then the permittivity equation can be expressed as follows:

The number of the simulated data points in Figure 9 has been increased and then the figure has been constructed as shown in Figure 10.

Figures 9 and 10 show the linearity of the frequency of resonance is decreasing by increasing the permittivity.

The measured polynomial fitting curve is plotted according to the loss tangent and the loaded Q-factor to obtain the loss tangent formula, as presented in Figure 11.

From Figure 11, the measured polynomial fitting curve of the loss tangent extracted equation is given by

By using matrix, the three data points in Figure 10 can be matched with the x and y coordinates, where the QL-factor matches with the x coordinate and the loss tangent matches with the y coordinate; thus:

This gives a system of three equations and three unknowns as follows:

As discussed previously, the loss tangent equation can be expressed as follows:

The number of the simulated data points in Figure 11 has been increased and then the figure has been constructed as shown in Figure 12.

Figures 11 and 12 show the linearity of the loaded quality factor is decreasing sharply by increasing the loss tangent.

The number of the simulated data points in Figures 10 and 12 has been increased to 7 data points, with the permittivity of the MUT at range from 1 to 13 and the loss tangent of the MUT at range from 0.043 to 0.901. The linearity helps achieve higher effective sensitivity at various ranges of MUT. Therefore, the measured data points have been chosen to fit the parabolic functions for both permittivity (15) and tan delta (24) since the measured data points are reliable for usage in real life considering the real environmental conditions; thus the proposed sensor is suitable for liquid characterisation and a good candidate for industrial applications. However, increasing the data points in simulation produces a different polynomial fitting curve and equation from the polynomial fitting curve and equation of the measurement. This is due to the dielectric properties of the additional MUT in the simulation. The differences between simulations and measurements are also affected by the SMA connector losses in practice, as an ideal connector is modeled in simulations.

A comparison of the suggested sensor with the other detailed studies is depicted in Table 3.

From the previous results, it can be noted that, by placing the quartz capillary in the proposed sensor, the resonant frequency shifts to 6.0475 GHz due to the dielectric properties of the capillary. The empty capillary is filled with air, which means that the resonant frequency at 6.0475 GHz is considered for the dielectric properties of the air where the measured results show a permittivity of 0.9600 and loss tangent of 0.0430. The permittivity of the LUT causes a frequency shift into the lower frequencies, where the measured permittivity of the ethanol sample is 4.4600. The measured resonant frequency resonates at 6.0100 GHz, while the measured permittivity is 11.0600 for the methanol sample, and the measured resonant frequency rises at 5.9675 GHz. Also, it can be noted that the Q-factor decreases by increasing the loss tangent where the measured loss tangent for ethanol is 0.8730, the measured QL-factor is 66.8150, and when the loss tangent for methanol is 0.9010, the measured QL-factor is 14.4100. From Figure 8, it is clear that the transmission coefficient is decreasing by increasing the loss tangent where the transmission coefficient for ethanol is better than the transmission coefficient of the methanol. The relationship between the loss tangent, QL-factor, and transmission coefficient is illustrated in Figure 13.

From Table 3, there is a good matching among the proposed sensor and lately published works of literature for ethanol and methanol at 6.1000 GHz. Nevertheless, the slight change of the permittivity values at 6.1000 GHz is principally due to temperature variation. The relative permittivity of ethanol was 8.9100 and 24.0200 for methanol at 2.45 GHz and 30°C as proposed in [43]. Moreover, at 5 GHz and 25°C, the ethanol and methanol relative permittivity were 5.2200 and 13.2300, respectively, as detailed in [44].

As aforementioned above, the results confirmed that the stated sensor has a high Q-factor and miniature size, which means that it has high performance and low cost. In addition, the small size of the microwave sensor required only a small volume of the LUT to detect the dielectric characteristics.

To recognise the current work in context, the performance of the proposed GWCR sensor is compared with the other planar resonator sensors in the literature. Table 4 shows the performance attributes of several typical microwave sensors with different structures in terms of sensing material, sensor technique, resonance frequencies (GHz), measured permittivity, and mean sensitivity. The experimental outcomes in the existing literature are offered for different ranges of permittivity; therefore, as a result of a fair comparison, the mean sensitivity values are calculated by the following expression:where is the incremental frequency shift, denotes to the resonant frequency (bare), and indicates the permittivity shift. From the mentioned earlier formula we can conclude that the proposed GWCR sensor provides a good mean sensitivity when exposed to known high-permittivity materials compared with the conventional highly sensitive sensor [45].

4. Conclusions

A novel design of a GWCR resonator sensor with a highly accurate Q-factor for liquid characterisation is introduced in this study. The main advantages of the proposed GWCR resonator are (1) more economical to build than the traditional metal waveguides and having (2) a lower loss with a higher Q-factor than microstrip structures. The modelled GWCR sensor used two liquid samples to validate the proposed sensor in the measurement setup: ethanol and methanol. The evaluated permittivity and loss tangent for ethanol and methanol at 6.1 GHz and temperature of 23.4°C have comparable results with the reference marks. The proposed GWCR shows the highest Q-factor compared to the recently published sensors. Besides, the GWCR sensor is shown to exhibit very high sensitivity and dynamic range versus the conventional MS sensors. Therefore, the finalised sensor is reliable, which makes it preferable for pharmaceutical applications.

Data Availability

All the data are included within the manuscript.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors would also like to thank Universiti Teknikal Malaysia Melaka (UTeM) for supporting this work under research Grants MTUNC/MC0011 and PJP/2019/FTKEE(2B)/S01654.