Methodology for measuring dielectric and magnetic material properties

Abstract

A methodology for measuring the dielectric and magnetic properties of materials is provided. The method comprises a multi-section transmission line (TL) and multi-section fluidic channels. Two of the TL sections are incorporated with channel sections. The length-ratio of the two line-channel sections (LCS) is L:L. A network analyzer coupled to the TL can be configured to measure scattering parameters when the LCS have no material-under-measurement (MUM), one LCS has MUM, both LCSs have MUM. Then, the S-parameters are transformed to ABCD matrices. When L:L is 2:1, the propagation constant of the LCSs can be calculated with analytical formulas over the measured frequency range. When L:L is not 2:1, the propagation constant can be obtained with numerical root-searching algorithms. The properties of MUM can be extracted by use of models which relate the propagation constants with the given geometry and substrate of LCSs as well as the unknown MUM properties.

Claims

1. A method for measuring the dielectric properties of liquids using a transmission line, wherein the transmission line comprises: a signal line; and two fluidic channel sections, wherein said fluidic channel sections are perpendicular to said signal line, wherein each of said fluidic channel sections have rectangular cross sections, wherein the widths of said fluidic channel sections are parallel to said signal line, wherein the width ratio of said two fluidic channel sections is 2:1, wherein said fluidic channel sections intercept said signal line, and wherein said transmission line has an L section on the side that the fluidic channel has the larger width and a 0.5 L section on the side with the smaller width; wherein said transmission line is connected to a device for measuring an amplitude and phase of reflected and transmitted electromagnetic signals and for measuring two-port scattering parameters of said transmission line, the method comprising: (a) measuring a plurality of scattering parameters of said transmission line at selected frequencies while said two fluidic channel sections are empty; (b) converting the scattering parameters to ABCD matrices M1; (c) measuring a plurality of scattering parameters of said transmission line at the selected frequencies while one of said fluidic channel sections is filed with said liquid and the other said fluidic channel section is empty; (d) converting the scattering parameters to ABCD matrices M2; (e) measuring a plurality of scattering parameters of said transmission line at the selected frequencies while said two fluidic channel sections are filed with said liquid; (f) converting the scattering parameters to ABCD matrices M3; (g) determining the propagation constant .sub.liquid of said transmission line sections while said fluidic channel sections are filled with said liquid by using the expression: ${}_{liquid}^{1,2}=\frac{2}{L}.Math.\cos h^{-1}\left(\frac{B.Math.\cos h\left({}_{air}.Math.\frac{L}{2}\right)\sqrt{\left(B^2-4\right)\cos h^2\left({}_{air}.Math.\frac{L}{2}\right)-A+2}}{2}\right)$ A=Tr(M1M2) and B=Tr(M2M3), while .sub.air, is the propagation constant of said empty transmission line L and 0.5 L sections; (h) determining the relative complex permittivity, 1, for said liquid based on .sub.liquid.

2. A method for measuring the dielectric properties of liquids using a transmission line, wherein the transmission line comprises: a signal line; and two fluidic channel sections, wherein said fluidic channel sections are perpendicular to said signal line, wherein each of said fluidic channel sections have rectangular cross sections, wherein the widths of said fluidic channel sections are parallel to said signal line, wherein the width ratio of said two fluidic channel sections is 2:1, wherein said fluidic channel sections intercept said signal line, and wherein said transmission line has an L section on the side that the fluidic channel has the larger width and a 0.5 L section on the side with the smaller width; wherein said transmission line is connected to a device for measuring an amplitude and phase of reflected and transmitted electromagnetic signals and for measuring two-port scattering parameters of said transmission line, the method comprising: (a) obtaining a material, wherein said material is a liquid or gas, wherein said material has known complex permittivities; (b) measuring a plurality of scattering parameters of said transmission line at selected frequencies while said two fluidic channel sections are filled with said material; (c) converting the scattering parameters to ABCD matrices M1; (d) measuring a plurality of scattering parameters of said transmission line at the selected frequencies while one of said fluidic channel sections is filed with said liquid and the other said fluidic channel section is filled with said material; (e) converting the scattering parameters to ABCD matrices M2; (f) measuring a plurality of scattering parameters of said transmission line at the selected frequencies while said two fluidic channel sections are filed with said liquid; (g) converting the scattering parameters to ABCD matrices M3; (h) determining the propagation constant .sub.liquid of said transmission line sections while said fluidic channel sections are filled with said liquid by using the expression: ${}_{liquid}^{1,2}=\frac{2}{L}.Math.\cos h^{-1}\left(\frac{B.Math.\cos h\left({}_{material}.Math.\frac{L}{2}\right)\sqrt{\left(B^2-4\right)\cos h^2\left({}_{\mathrm{material}}.Math.\frac{L}{2}\right)-A+2}}{2}\right)$ A=Tr(M1M2) and B=Tr(M2M3), while .sub.material, is the propagation constant of said transmission line L and 0.5 L sections while said fluidic channel sections are filled with said material; (i) determining the relative complex permittivity, 1, for said liquid based on .sub.liquid.

Description

(3) BRIEF DESCRIPTION OF THE DRAWINGS

(1) A full and enabling disclosure of the present invention, including the best mode thereof, directed to one of ordinary skill in the art, is set forth in the specification, which makes reference to the appended figures, in which:

(2) FIG. 1(a) illustrates a schematic block diagram of a coplanar-waveguide transmission line with a single microfluidic channel as an exemplary embodiment of the methodology;

(3) FIG. 1(b) illustrates a schematic block diagram of a coplanar-waveguide transmission line with two disconnected microfluidic channels as an exemplary embodiment of the methodology;

(4) FIG. 1(c) illustrates a schematic block diagram of a meandered coplanar-waveguide transmission line with a single microfluidic channel as an exemplary embodiment of the methodology;

(5) FIG. 1(d) illustrates a cross-section view of a coplanar waveguide with microfluidic channel according to an exemplary embodiment of the presented method;

(6) FIG. 2 illustrates a diagram representing a cascaded network of the transmission line comprising five dual-port networks according to an exemplary embodiment of the presented method;

(7) FIG. 3(a) illustrates a top view of a coplanar waveguide based measurement device for liquid property measurement as an exemplary embodiment of the methodology;

(8) FIG. 3(b) illustrates a cross-section view of a coplanar waveguide (CPW) with microfluidic channel according to an exemplary embodiment of the presented methodology;

(9) FIG. 4(a)-(d) illustrates the graphs of experimentally obtained complex permittivity of (a) IPA, (b) ethanol, (c) methanol, and (d) DI water with comparison to the values from Cole-Cole model according to an exemplary embodiment of the presented method;

(10) FIG. 5(a) illustrates a graph depicting the real-time measurement results from 100 MHz to 9 GHz when DI water and IPA flows in succession in an exemplary embodiment of the methodology;

(11) FIG. 5(b) illustrates a real-time measurement result of DI water and IPA flowing in succession at 2 GHz along the arrow in FIG. 5(a);

(12) FIG. 6(a)-(b) illustrates the experimentally obtained complex permittivity of (a) DI water and (b) IPA with comparison to the values from Cole-Cole model according to an exemplary embodiment of the presented methodology;

(13) FIG. 7 illustrates a flow diagram of a method according to an exemplary embodiment of the presented invention.

(14) Repeat use of reference characters throughout the present specification and appended drawings is intended to represent same or analogous features or elements of the invention.

(4) DETAILED DESCRIPTION OF THE INVENTION

(15) Reference now will be made in detail to embodiments of the invention, one or more examples of which are illustrated in the drawings. Each example is provided by way of explanation of the invention, not limitation of the invention. In fact, it will be apparent to those skilled in the art that various modifications and variations can be made in the present invention without departing from the scope or spirit of the invention. For instance, features illustrated or described as part of one embodiment can be used with another embodiment to yield a still further embodiment. Thus, it is intended that the present invention covers such modifications and variations as come within the scope of the appended claims and their equivalents.

(16) Generally, the present application relates to automatic and self-calibration of devices for quantitative material measurement. The device described herein simultaneously addresses the need for calibration standards or reference material and the uncertainty of multiple connections. In actual embodiment, the transmission line does not need to be straight. Any type of transmission lines, including waveguides, can be used. The materials to be measured can be solid, such as thin films, gas, such as volatile organic chemicals (VOC, which can be absorbed with polymers for stronger signals), and liquids, including liquids with suspended particle.

(17) In operation, a vector network analyzer (VNA) can be configured to measure the four scattering parameters or S-parameters (S.sub.11, S.sub.12, S.sub.21, S.sub.22), under three different states: there is no material-under-measurement (MUM, e.g. liquid), there is MUM on one transmission line section, and there is MUM on both transmission line sections. Then, the S-parameters are transformed to ABCD matrices. After processing the ABCD matrices, the propagation constant of the MUM line section can be calculated with analytical formulas when the length ratio of the two line-channel sections is 2:1 (or 1:2). When the ratio is not 2:1, numerical root-searching algorithms can be used to obtain the propagation constant of the MUM line section. The use of two line-channel sections, instead of one, can facilitate the searching process. Finally, the permittivity or permeability of MUM can be extracted by use of geometry and material specific models.

(18) Referring to the drawings, FIG. 1(a) illustrates a schematic block diagram of a coplanar-waveguide (CPW) transmission line (TL) with a single microfluidic channel as an exemplary embodiment of the methodology. The CPW comprises sections 100, 102, 104, 106 and 108. The microfluidic channel comprises sections 110, 112 and 114. Section 102 intercepts with section 110, and section 106 intercepts with section 114. The length ratio between 102 and 104 is L:L. During MUM measurements, the CPW experiences three different states. Before MUM reaches the section 102, air fills both 102 and 106. The second state is when MUM fills 102 but not 106. The third state is when MUM fills both 102 and 106. The order to fill TL line sections with MUM can be reversed.

(19) Full S-parameters for the three states can be recorded in a single-connection measurement while the vector network analyzer performs measurement operations continuously.

(20) Referring to the drawings, FIG. 1(b) illustrates a schematic block diagram of a coplanar-waveguide transmission line with two disconnected microfluidic channels as an exemplary embodiment of the methodology. The CPW comprises sections 120, 122, 124, 126 and 128. The two microfluidic channel sections 130 and 132 are connected together with external tubing or connected to separate material delivery devices. Section 122 intercepts with section 130, and section 126 intercepts with section 132. The length ratio between 122 and 126 is L:L. During MUM measurements, the CPW experiences three different states. Before MUM reaches the section 122, air fills both 122 and 126. The second state is when MUM fills 122 but not 126. The third state is when MUM fills both 122 and 126. Full S-parameters for the three states can be recorded in a single-connection measurement while the vector network analyzer performs measurement operations continuously.

(21) FIG. 1(c) illustrates a schematic block diagram of a meandered coplanar-waveguide transmission line with a single microfluidic channel as an exemplary embodiment of the methodology; The CPW comprises sections 140, 142, 144, 146 and 148. The microfluidic channel has two sections 150 and 152. Section 142 intercepts with section 150, and section 146 intercepts with section 152. The width ratio between 142 and 146 is L:L. During MUM measurements, the CPW experiences three different states. Before MUM reaches the section 142, air fills both 142 and 146. The second state is when MUM fills 142 but not 146. The third state is when MUM fills both 142 and 146. Full S-parameters for the three states can be recorded in a single-connection measurement while the vector network analyzer performs measurement operations continuously.

(22) Referring now to FIG. 1(d), a cross-section view of the TL section with the microfluidic channel in FIG. 1(a)-(c) is illustrated. The conductor backed coplanar waveguide (CPW) is demonstrated as an example, which comprises metallic signal line 160, ground 162, and substrate 164. Other types of transmission line, such as microstrip line, slot line or rectangular waveguide, are also applicable. The microfluidic channel cover 166 covering the TL can be filled with air or material-under-measurement (MUM) 168. A dielectric layer 170 can be optionally coated on the metallic line 160, 162 for different MUM measurements, such as polymer for gas measurement.

(23) Referring now to FIG. 2, a diagram representing a cascaded network of the devices in FIG. 1 with ABCD matrices is presented. The network between ports 1 and 2 is comprised of left network 200, section 202, center network 204, Section 206 and right network 208. Therein, left and right networks 200, 208 include all the connections between ports 1 and 2 and sections 202, 206. The center network 204 isolates sections 202, 206 and provides or enables a time delay needed for measurements and data recording. For sections 202 and 206 transmission line sections, their ABCD matrices are:

(24) $\begin{array}{cc}{M}^{I\left(\mathrm{II}\right)}=\left[\begin{array}{cc}\cosh \left({}^{I\left(\mathrm{II}\right)}{L}^{I\left(\mathrm{II}\right)}\right)& Z^{I\left(\mathrm{II}\right)}\sinh \left({}^{I\left(\mathrm{II}\right)}{L}^{I\left(\mathrm{II}\right)}\right)\\ \frac{1}{Z^{I\left(\mathrm{II}\right)}}\sinh \left({}^{I\left(\mathrm{II}\right)}{L}^{I\left(\mathrm{II}\right)}\right)& \cosh \left({}^{I\left(\mathrm{II}\right)}{L}^{I\left(\mathrm{II}\right)}\right)\end{array}\right]& \left(1\right)\end{array}$

(25) The two TL sections have the same cross section geometry, thus identical propagation constant, .sup.I=.sup.II==+j, and characteristic impedance Z.sup.I=Z.sup.II=Z when air or MUM filled both sections and when MUM is at the same temperature as the system in FIG. 1. The overall ABCD matrix of the 2-port network is
M.sub.LeftM.sub.i.sup.IM.sub.CenterM.sub.i.sup.IIM.sub.Right=M.sub.i(2)

(26) where subscript i denotes one of the three states of sections 202, 206.

(27) At (700), the measurement procedure for an exemplary embodiment in the present disclosure is demonstrated. The measured scattering parameters of the device in three states are described as the following:

(28) 1) At (702), both section 202 and section 206 are full of air.
M.sub.LeftM.sub.air.sup.IM.sub.CenterM.sub.air.sup.IIM.sub.Right=M.sub.1(3)

(29) 2) At (704), section 202 is filled with MUM while section 206 with air.
M.sub.LeftM.sub.liquid.sup.IM.sub.CenterM.sub.air.sup.IIM.sub.Right=M.sub.2(4)

(30) 3) At (710), both sections 202, 206 are full of MUM.
M.sub.LeftM.sub.liquid.sup.IM.sub.CenterM.sub.liquid.sup.IIM.sub.Right=M.sub.3(5)

(31) Multiply matrix M.sub.1 (M.sub.2) by the inverse matrix of M.sub.2 (M.sub.3), we obtain:
M.sub.LeftM.sub.air.sup.I(M.sub.liquid).sup.1(M.sub.Left).sup.1=M.sub.1M.sub.2.sup.1(6)
(M.sub.LeftM.sub.liquid.sup.IM.sub.Center)M.sub.air.sup.II(M.sub.liquid.sup.II).sup.1(M.sub.LeftM.sub.liquid.sup.IM.sub.Center).sup.1=M.sub.2M.sub.3.sup.1(7)

(32) Notice that M.sub.1M.sub.2.sup.1(M.sub.2M.sub.3.sup.1) and M.sub.air.sup.I(M.sub.liquid.sup.I).sup.1(M.sub.air.sup.II(M.sub.liquid.sup.II).sup.1) are similar matrices. According to trace matrix theory, we can obtain:

(33) $\begin{array}{cc}\begin{array}{c}A=\mathrm{Tr}\left(M_1{M}_2^{-1}\right)\\ =\mathrm{Tr}\left({M_{\mathrm{air}}^I\left(M_{\mathrm{liquid}}^I\right)}^{-1}\right)\\ =2\cosh \left({}_{\mathrm{air}}{L}^I\right)\cosh \left({}_{\mathrm{liquid}}{L}^I\right)-\\ \left(\frac{Z_{\mathrm{air}}}{Z_{\mathrm{liquid}}}+\frac{Z_{\mathrm{liquid}}}{Z_{\mathrm{air}}}\right)\sinh \left({}_{\mathrm{air}}{L}^I\right)\sinh \left({}_{\mathrm{liquid}}{L}^I\right)\end{array}& \left(8\right)\end{array}$$\begin{array}{cc}\begin{array}{c}B=\mathrm{Tr}\left(M_2{M}_3^{-1}\right)\\ =\mathrm{Tr}\left({M_{\mathrm{air}}^{\mathrm{II}}\left(M_{\mathrm{liquid}}^{\mathrm{II}}\right)}^{-1}\right)\\ =2\cosh \left({}_{\mathrm{air}}{L}^{\mathrm{II}}\right)\cosh \left({}_{\mathrm{liquid}}{L}^{\mathrm{II}}\right)-\\ \left(\frac{Z_{\mathrm{air}}}{Z_{\mathrm{liquid}}}+\frac{Z_{\mathrm{liquid}}}{Z_{\mathrm{air}}}\right)\sinh \left({}_{\mathrm{air}}{L}^{\mathrm{II}}\right)\sinh \left({}_{\mathrm{liquid}}{L}^{\mathrm{II}}\right)\end{array}& \left(9\right)\end{array}$

(34) For generality, we define the length ratio of section 202 and section 206 is always larger or equal to 1, i.e.

(35) $\frac{L^I}{L^{\mathrm{II}}}1.$
After collecting the S-parameters of the first two states, we need to determine if the length ratio is equal to 2 at (706). For the case

(36) $\frac{L^I}{L^{\mathrm{II}}}2$
at (708), root-searching algorithm techniques can be used to solve eq. (8) or (9) for .sub.liquid. For the case

(37) $\frac{L^I}{L^{\mathrm{II}}}=2$
at (710), the S-parameters of the third state allows the removal of term

(38) $\left(\frac{Z_{\mathrm{air}}}{Z_{\mathrm{liquid}}}+\frac{Z_{\mathrm{liquid}}}{Z_{\mathrm{air}}}\right)$
in eqs. (8) and (9) by using the sums of arguments rule of hyperbolic function:

(39) $\begin{array}{cc}4{\cosh }^2\left({}_{\mathrm{liquid}}\frac{L}{2}\right)-4B\cosh \left({}_{\mathrm{air}}\frac{L}{2}\right)\cosh \left({}_{\mathrm{liquid}}\frac{L}{2}\right)+A+4{\cosh }^2\left({}_{\mathrm{air}}\frac{L}{2}\right)-2=0& \left(10\right)\end{array}$

(40) Since

(41) $\cosh \left({}_{\mathrm{air}}\frac{L}{2}\right)$
can be calculated with eq. (13) below when device dimensions and air permittivity are known, eq. (10) can be considered as a quadratic equation of

(42) $\cosh \left({}_{\mathrm{liquid}}\frac{L}{2}\right).$
Then, the transmission-line propagation constants can be obtained from the measured S-parameters in three states with the following formulas at (712):

(43) 0$\begin{array}{cc}{}_{\mathrm{liquid},1}^{1,2}={\cosh }^{-1}\left[\frac{\begin{array}{c}B\cosh \left({}_{\mathrm{air}}\frac{L}{2}\right)-\\ \sqrt{\left(B^2-4\right){\cosh }^2\left({}_{\mathrm{air}}\frac{L}{2}\right)-A+2}\end{array}}{2}\right].Math.\frac{2}{L}& \left(11\right)\end{array}$$\begin{array}{cc}{}_{\mathrm{liquid},2}^{1,2}={\cosh }^{-1}\left[\frac{\begin{array}{c}B\cosh \left({}_{\mathrm{air}}\frac{L}{2}\right)+\\ \sqrt{\left(B^2-4\right){\cosh }^2\left({}_{\mathrm{air}}\frac{L}{2}\right)-A+2}\end{array}}{2}\right].Math.\frac{2}{L}& \left(12\right)\end{array}$

(44) where .sup.1 and .sup.2 are two eigenvalue solutions from

(45) $\cosh \left({}_{\mathrm{liquid}}\frac{L}{2}\right).$
They describe and incident and a reflected wave with .sup.1=.sup.2. Furthermore, the physically meaningful propagation constant needs to be identified from .sub.liquid,1.sup.1,2 and .sub.liquid,2.sup.1,2. Some educated, but often obvious, rules can be exploited from the obtained and || and ||. For a reasonable solution, its .sub.liquid is expected to be larger than .sub.air and have the same sign as .sub.liquid for each frequency.

(46) With the obtained values, MUM permittivity can be obtained by use of geometry and material specific models at (714). Many of which are available in literature, such as the following one:

(47) $\begin{array}{cc}={}_c+j\frac{2f}{c}\sqrt{{}_{eff}}={}_c+j\frac{2f}{c}\sqrt{q{}_1+C_0}& \left(13\right)\end{array}$

(48) where c is the speed of light in vacuum, .sub.eff is the effective permittivity, .sub.1 is the relative permittivity of MUM, and q and C.sub.0 are constants that can be derived from conformal mapping (CM) or single layer reduction techniques. Conductor loss (.sub.c) can also be considered for better accuracy.

(49) The equations to calculate from measured S-parameters are applicable to any type of transmission lines even though different transmission lines have different models to extract MUM permittivity from .

(50) Referring to FIG. 3(a) a top view of the experimental device constructed according to an exemplary embodiment of the presented methodology is illustrated. In accordance with one embodiment, the two microfluidic channels have a width ratio of 2:1 and are connected together with a capillary tube to form a one-channel device for single-channel and single-connection measurement demonstration without the need for liquid reference materials or calibration standards. FIG. 3(b) shows the cross section of the microfluidic channels (sections 202 and 206) on the device. In accordance with one embodiment, the upper .sub.3 layer may be constructed by polydimethylsiloxane (PDMS). The lower .sub.2 layer can be constructed from various materials, such as glass, quartz and/or silicon. The coplanar waveguide is illustrated with signal line width w and gap g.

(51) Referring to FIG. 4(a)-(d), plots illustrating the comparison between measured .sub.liquid* and Cole-Cole model are presented for deionized (DI) water, methanol, ethanol and IPA, respectively. The samples are measured separately, using air as reference in accordance with an embodiment of the present method. The plots reveal that, for all the MUMs, the obtained .sub.liquid* agree with Cole-Cole model reasonably well. The calculated correlation coefficients R.sup.2 further confirm the agreement since R.sup.2 are higher than 0.937 for both real and imaginary permittivity components. Besides, the average relative errors (AREs) of the MUM permittivity are lower than 7.28%. ARE is defined as

(52) $\begin{array}{cc}\mathrm{ARE}=\frac{1}{N}\underset{i=1}{\overset{N}{.Math.}}.Math.\frac{{}_i^{\mathrm{Measured}}-{}_i^{\mathrm{reference}}}{{}_i^{\mathrm{reference}}}.Math.100\%& \left(14\right)\end{array}$

(53) to describe spectrum accuracy, where N is the number of measured frequency points over the spectrum.

(54) Referring to FIG. 5(a), a plot illustrating measurement results of DI water and IPA in succession in an embodiment are shown. FIG. 5(b) shows |S.sub.21| vs. time at 2 GHz. The five phases of different channel conditions are clearly shown. By using the proposed method, the permittivity of DI water can be extracted from Phase I to Phase III with air as reference. Then the permittivity of IPA can be extracted from Phase III to Phase X by considering DI water as a reference material.

(55) Referring now to FIG. 6(a)-(b), plots illustrating the comparison between measured .sub.liquid* and Cole-Cole model for DI water and IPA in succession is presented. The extracted .sub.liquid* agrees with results from Cole-Cole model reasonably well. The correlation coefficients R.sup.2 are higher than 0.948 and the AREs are lower than 7.11% for both MUMs.

Experimental Measurement

(56) A coplanar waveguide (CPW) with two microfluidic channels in FIG. 3, is built to verify the proposed method. The microfluidic channels made by polydimethylsiloxane (PDMS) are bonded to the CPW to hold MUMs. The CPW and glass substrate surface serve as the bottom wall of the microfluidic channels. The device is fabricated with standard microfabrication techniques on a 700-m-thick quartz wafer. A 200-nm-thick gold film was deposited on a 20-nm-thick Ti adhesion layer to form a CPW structure. The signal line width of the CPW is 40 um and the gap is 5 um. The microfluidic channels are placed perpendicular to transmission line, with a length of 200 um and 100 um, respectively, along the CPW. The channel height is 100 um and the PDMS thickness is 3 mm. The device is settled on a brass holder for better resistance to mechanical vibration and parasitic modes that can dominate at high frequencies. A vector network analyzer (VNA) is used to measure the 2-port S parameters from 100 MHz to 9 GHz.

(57) In embodiment I, where the two channels are not connected together, de-ionized (DI) water, methanol, ethanol and 2-propanol (IPA) from Sigma-Aldrich (US), are measured separately. For each of the sample liquids, the S-parameters of the three states are obtained. Each measurement is repeated three times. The measured data are then processed by use of eqs. (11) and (12) to obtain .sub.liquid. The complex permittivity of the MUM is obtained from eq. (5). FIG. 4 shows the obtained .sub.liquid* for each MUM. The results agree with Cole-Cole model reasonably well. The calculated correlation coefficients R.sup.2 further confirm the agreement since R.sup.2 are higher than 0.937 for both real and imaginary permittivity components. Besides, the average relative errors (AREs) of the MUM permittivity are lower than 7.28%.

(58) In embodiment II, where the two microfluidic channels are connected together with a silicone tubing to form a single channel, DI water and IPA are measured in succession. The 100 mm long tubing has a 0.51 mm inner diameter and provides 60 s delay time. The obtained S parameters, such as S.sub.21 magnitude, are plotted in FIG. 5. FIG. 5(b) shows |S.sub.21| vs. time at 2 GHz, where the five phases of different channel conditions are clearly shown. By using the proposed method, the permittivity of DI water can be extracted from Phase I to Phase III with air as reference. Then the permittivity of IPA can be extracted from Phase III to Phase X by considering DI water as a reference material. FIG. 6 shows the extracted .sub.liquid*, which agrees with results from Cole-Cole model. The correlation coefficients R.sup.2 are higher than 0.948 and the AREs are lower than 7.11% for both MUMs. The measurement accuracy is comparable with the results in separate-channel measurements.

(59) In addition to measuring material properties, the present methodology has significant applicability in the development of portable, quantitative and fast RF/microwave probes that are easy to operate.

(60) Such probes can be used to characterize various materials in-situ and on-line for real-time monitoring.

(61) While the present subject matter has been described in detail with respect to specific embodiments thereof, it will be appreciated that those skilled in the art, upon attaining an understanding of the foregoing may readily produce alterations to, variations of, and equivalents to such embodiments. Accordingly, the scope of the present disclosure is by way of example rather than by way of limitation, and the subject disclosure does not preclude inclusion of such modifications, variations and/or additions to the present subject matter as would be readily apparent to one of ordinary skill in the art.