Next Article in Journal
Flexible Self-Powered Low-Decibel Voice Recognition Mask
Previous Article in Journal
Human Activity Recognition Algorithm with Physiological and Inertial Signals Fusion: Photoplethysmography, Electrodermal Activity, and Accelerometry
Previous Article in Special Issue
Data-Monitoring Solution for Desalination Processes: Cooling Tower and Mechanical Vapor Compression Hybrid System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electric Susceptibility at Partial Coverage of a Circular One-Side Access Capacitive Sensor with Rigid Polyurethane Foams

Institute for Mechanics of Materials, University of Latvia, LV-1004 Riga, Latvia
Submission received: 18 March 2024 / Revised: 22 April 2024 / Accepted: 6 May 2024 / Published: 9 May 2024
(This article belongs to the Special Issue Sensors in 2024)

Abstract

:
The capability of dielectric measurements was significantly increased with the development of capacitive one-side access physical sensors. Complete samples give no opportunity to study electric susceptibility at a partial coverage of the one-side access sensor’s active area; therefore, partial samples are proposed. The electric susceptibility at the partial coverage of a circular one-side access sensor with cylinders and shells is investigated for polyurethane materials. The implementation of the relative partial susceptibility permitted us to transform the calculated susceptibility data to a common scale of 0.0–1.0 and to outline the main trends for PU materials. The partial susceptibility, relative partial susceptibility, and change rate of relative partial susceptibility exhibited dependence on the coverage coefficient of the sensor’s active area. The overall character of the curves for the change rate of the relative partial susceptibility, characterised by slopes of lines and the ratio of the change rate in the centre and near the gap, corresponds with the character of the surface charge density distribution curves, calculated from mathematical models. The elaborated methods may be useful in the design and optimization of capacitive OSA sensors of other configurations of electrodes, independent of the particular technical solution.

Graphical Abstract

1. Introduction

In the non-destructive evaluation of dielectric materials, in the frequency band up to 10 MHz, the capacitance method is one of the main testing methods [1]. The capability of dielectric measurements was significantly increased with the development of capacitive one-side access (OSA) physical sensors, which permitted non-destructive testing and the characterization of materials without preparation of specimens [2,3,4,5], by merely pushing the sensor against one side of the test object. Several problems exist when working with capacitive OSA sensors: they excite test objects with a non-homogeneous electric field, have low sensing (geometric) capacitances, and high stray capacitances; bulk reference specimens are needed, whose minimum thickness (height) is limited by the effective depth of the penetration of the electrical field. Contrary to the parallel-plate capacitive sensors, working with a homogeneous electric field, there are no simple formulas for calculating the sensing capacities of OSA sensors at a complete or partial coverage with a dielectric [3,6,7].
The intensity maximum of a low-frequency excitation field is situated in the direct vicinity of the active area of a capacitive OSA transducer’s electrodes [2]. The penetration depth of the excitation field into the test object determines the lateral dimensions (height/thickness) appropriate for permittivity measurements [2,3,7,8]. In [7], the penetration depth of a concentric, coplanar, and capacitive OSA sensor was estimated by identifying the test object’s height, for which the capacitance is 3% or 10% smaller than its value when in contact with a similar, but infinitely high, test-piece. For a given sensor configuration, the sensor penetration depth increases as the test object’s permittivity increases. The measured permittivity of polyurethane foams and PTFE was investigated in [8] with a dependence on the OSA sensor’s working volume, when the latter is filled partly along the lateral dimension (height) of the samples. The height of the samples was gradually decreased, and a 3% criterion to detect the penetration depth was applied.
Mathematical modelling and numerical calculations in [7] outlined an optimal two-electrode OSA sensor configuration for a given width of the gap, which ensured the maximum transcapacitance. Utilizing the spatial domain Green’s functions and the method of moments [9,10], the sensor surface charge density was calculated, and then the sensor’s transcapacitance was determined from its surface charge. In [11], the charge density on circular electrodes of a two-electrode OSA sensor was assumed to be constant, and the average potential of the electrodes was calculated. In [12,13], a circular three-electrode OSA sensor is considered. A jump in the normal component of the electric field across an interface is defined, which would be equal to the surface charge density in the absence of ohmic conduction in the medium. The jump is zero at every interface except in the plane of the electrodes. For every Fourier–Bessel mode [14], the jump is related to the difference in the surface capacitance density.
The regions of an OSA sensor that contribute most to the sensor’s capacitance are the outer edge of the inner electrode and the inner edge of the outer electrode, where the surface charge density is the highest [6,7,15]. This results in a low capacitance of the sensor. The most sensitive zone of the sensor closely corresponds to the location of the gap between its two electrodes, and for sensors with a relatively large radius of the inner electrode, there will be an insensitive zone at the centre [6]. Empirical knowledge of a certain physical property of dielectric material, dependent on the surface charge distribution on the active area, would facilitate a performance evaluation, modelling, and the design of optimal proportions of electrodes’ dimensions [16,17,18,19,20,21,22].
PU foam is a “Polymer-gas” composite, which has the same matrix—monolithic polyurethane, in the whole range of densities, i.e., 30–1280 kg/m3 [23,24,25]. For light-weight PU foams (density of 30 kg/m3), permittivity ε ≈ 1.065, for medium weight (550 kg/m3), ε ≈ 2.10, and for monolithic polyurethane (1280 kg/m3), ε ≈ 3.35 (all at a frequency of 1 kHz) [8,26,27]. Rigid polyurethane (PU) foams [23,24] exhibit low dielectric interference and nearly non-dispersive permittivity; therefore, excellent dielectric performance can be ensured in a wide frequency range. The foams are comparatively easy to process mechanically, which makes them appropriate for a study of permittivity and electric susceptibility at a partial coverage of an OSA sensor.
A complete dielectric sample covers the total active area of the OSA sensor, and the permittivity, measured for it, is the true permittivity, which characterises the overall effect of the electric charge on the sample. Complete samples provide no opportunity to study permittivity and electric susceptibility at a partial coverage of the active area and to estimate their correspondence to the charge density distribution. Therefore, partial samples are proposed for measurements of permittivity and the study of electric susceptibility. This study aims to investigate electric susceptibility at a partial coverage of a circular capacitive OSA physical sensor with samples—cylinders and shells, made from PU foams and monolithic polyurethanes. A transition to the relative partial electric susceptibility permitted us to transform the calculated data of the partial susceptibility to a common scale of 0.0–1.0 and to outline the main trends for PU materials of significantly differing true permittivity. The change rate of the relative partial susceptibility with a dependence on the coverage coefficients of the active area and the subsample’s radius is evaluated. A comparison with the results of other authors is made.

2. Materials and Methods

2.1. Polyurethane Materials

Rigid closed-cell petrochemical PU foams were made at apparent core densities of 30 kg/m3 < ρ < 80 kg/m3, in blocks, according to the formulation in [8,26,27]. The apparent core density (density) was determined in correspondence with ISO 845:2006 [28]. The difference in densities was achieved by varying the amount of physical or chemical blowing agents. To make monolithic polyurethane, a liquid mixture of the same formulation as for PU foams was poured into polyethylene ampoules with an inner diameter of 25.4 mm and a length of 115 mm; no foaming agent was added. The ampoules were centrifuged for ~20 min at 5000 rpm in a centrifuge Sigma 3-30KS (Sigma Laborzentrifugen GmbH, Osterode am Harz, Germany) [26]. The technology ensured that ≈2/3 of the length of the cylindrical part of the PU rods was free from gaseous inclusions.
Industrially manufactured rigid petrochemical PU foam SikaBlock, with densities of 85 kg/m3 ≤ ρ ≤ 415 kg/m3, monolithic polyurethane SikaBlock-M945, and a comparative material—monolithic epoxy LAB 975 NEW—were acquired from Sika JSC (Baar, Switzerland).
Depending on the production technology, PU foams can exhibit anisotropy of structure and properties [29,30], when the permittivity and electric susceptibility εij and χij are the 2nd-rank tensors. Only isotropic PU foams are considered further.

2.2. Capacitive One-Side Access Sensor

Permittivity was measured with an experimental dielectric spectrometer equipped with a capacitive sensor of the one-side access (OSA) type [3,8,26,27]. The dielectric sample was placed on the active area of the sensor, Figure 1, and excited via electrodes, by an electrical field generated by sinusoidal voltage signals. The amplitude value of the sinusoidal excitation signals was U0 = 20 V. The signals were generated at discrete frequencies, increasing in a geometric progression as follows:
fn = f1, 2f1, …, 2(n − 1)f1 Hz, where f1 = 10 Hz, n = 1, 2, …, 16,
f = 10, 20, …, 327,680 Hz,
where n is the ordeal number of a frequency. In the capacitive sensor system (CSS), the driven electrode №1 and sensing electrode №2 are connected to the multi-frequency excitation generator; the driven electrode is connected directly to the sensing electrode through a reference capacitor Cref [3].
To carry out stray-immune capacitance measurements, the sensing electrode, the reference capacitor, and the unity gain buffer amplifier are covered with a screen, which forms the guard electrode №3 on the active area of CSS around the sensing electrode, as shown in Figure 1. The guard electrode is fed by the same voltage as the sensing electrode (active guarding) through a voltage follower, thus suppressing the electric field between the driven electrode and the sensing electrode outside the sensor’s active surface.
The accuracy of the spectrometer in permittivity measurements, in conditions of repeatability, was estimated with the expanded uncertainty US95.45% = ±0.01 [8,31,32]. The calibration of the spectrometer was made before each measurement, with respect to permittivity, delivered by the sensor in the air. Three consecutive measurements of permittivity spectra were made for each data point.

2.3. Dielectric Losses

The dielectric loss part ε″(f) of the complex permittivity ε ~ (jf) = ε′(f) − jε″(f) was evaluated for PU foams and monolithic polyurethane, using data from direct measurements and scientific information sources [8]. The dielectric losses of PU materials were measured with a broadband dielectric spectrometer BDS-50 (Novocontrol Technologies GmbH & Co. KG, Montabaur, Germany), comprising a parallel plate capacitor, for samples with a diameter of 30 mm and a thickness of 2 mm. Five measurements were made for each data point. The real and the imaginary parts ε′(f) and ε″(f) were compared.

2.4. Lateral Dimensions of a Complete Sample

In order to study the special cases of the partial coverage of the OSA sensor’s active area, several benchmark experiments were carried out (Points 2.4., 2.5., 2.6., and 2.8.), and their results were analysed (Points 4.2., 4.3., 4.4., and 4.6.).
Let us consider a PU foam sample with a diameter d = D0 = 45 mm and thickness t = 65 mm. The cross-sectional surface of the sample covers the whole active area of the circular OSA sensor, and the thickness considerably exceeds the penetration depth of the low-frequency electric field into the PU foams [1 d c, PTFE]. The impact of the lateral dimension (the thickness) of such a sample on the permittivity was studied. The intensity maximum of the dielectric spectrometer’s low-frequency excitation field is situated in the direct vicinity of the active area of the OSA sensor’s electrodes [2,8]. Therefore, first, the permittivity of the samples was determined. Then, a thin layer, thickness Δt = 0.5–2.0 mm, was cut from the top of the sample, and the permittivity of the remaining part was measured. Then, the next layer with a thickness Δt = 0.5–2.0 mm was cut, etc. The experimental data were plotted against the thickness values of the sample. Then the penetration depth of the electric field of a certain frequency f was defined as such a thickness t3% of a sample, at which the electric susceptibility is 3% less than the true susceptibility χt = εt − 1.0 of an infinitely thick sample, χ = 97% χt [7,8].
A sample is considered complete if (1) its cross-sectional surface covers the whole active area of the sensor and (2) its thickness equals or exceeds the penetration depth of the electric field of a certain frequency into the given dielectric. A measurement on a complete sample provides the true permittivity εt of the dielectric material.
The reading of the dielectric spectrometer for any sample, which either (1) covers the active area of the sensor partially and/or (2) has a thickness, smaller than the penetration depth of the electric field of a certain frequency, is denoted further as “The measured permittivity εm”.

2.5. Transversal Dimensions of Samples

The impact of the transversal dimensions of a sample on the permittivity was studied for parallelepiped-shaped samples of dimensions (1) 20 × 65 × 65 mm and (2) 20 × 45 × 45 mm and for (3) cylindrical samples with a diameter d = D0 = 45 mm and a thickness of 20 mm, Figure 2a. The “(2)” parallelepipeds were made from the “(1)” ones, and then the cylinders were made from the “(2)” parallelepipeds. All three kinds of samples are considered complete, since covering the whole active area is ensured. The permittivity of samples “(1)”, “(2)”, and “(3)” was determined for lab-made PU foam (50 kg/m3), PU foam SikaBlock-M150 (144 kg/m3), and SikaBlock-M450 (415 kg/m3).
Another experiment was performed for three other parallelepipeds of the kind “(1)” (20 × 65 × 65 mm), made from the same PU foams. First, the permittivity was measured for each parallelepiped, and then a concentric cylindrical hole, diameter d ≈ D0 = 45 mm, was made in each sample. The sample was placed on the isolated circular case of the sensor, Figure 2b, and a reading of the spectrometer was made.

2.6. Inner Vertical Surfaces

The impact of inner vertical surfaces in a complete, circular cylindrical sample on the measured permittivity was estimated experimentally. The task arises due to the limited transversal and lateral dimensions of the lab-made monolithic PU rods, when only one semi-cylinder (D0 = 45 mm and h = 12 mm) can be made from a PU rod. Side by side, the two semi-cylinders form a combined sample with two inner vertical surfaces along the diagonal. The semi-cylinders were pushed together, fixed firmly with an elastic band along the outer circular perimeter, Figure 3a, and the reading of the spectrometer was made.
The experiment was carried out also for a monolithic PU SikaBlock-M960. The dielectric SikaBlock-M960 was available in blocks of sufficiently large dimensions to make complete cylindrical samples of thickness t = 12 mm and diameter D0 = 45 mm. First, the permittivity of a complete, continuous sample was measured, and then it was cut into two semi-cylindrical parts. The parts were pushed together and fixed firmly with the elastic band, and the reading of the spectrometer was made. The experiment was made for the PU foam SikaBlock-M450 as well; the diameter of a sample was dc = D0 = 45 mm, and the thickness was t = 25 mm, as shown in Figure 3b.

2.7. Coverage Coefficients and the Corresponding Cylinders and Shells

Cylinders and shells of circular symmetry were used to study the partial susceptibility of rigid PU foams and monolithic polyurethanes, depending on the coverage coefficients of the active area of the OSA capacitive sensor. When a cylinder or a shell is placed concentrically on the circular active area, its coverage can be characterised by coverage coefficients kc and ks as follows:
k c = S c S 0 = d c D 0 2 = r c R 0 2 , k s = S s S 0 = π D 0 2 d s i n 2 π D 0 2 = 1.0 d s i n D 0 2 = 1.0 r s i n R 0 2 ,   d s i n = d sout t , 0.0 r c ,   r sin R 0  and  0.0 k c ,   k s 1.0 ,
where dc—the diameter of the cylinder (dc ≤ D0), t—the thickness of the shell wall, dsin and dsout = D0 = 45 mm—the inner and the outer diameters of the shell (dsin ≤ dsout), Sc—the cross-sectional area of a cylinder, Ss—the cross-sectional area of a shell and the total active area of the sensor, S0 = π/4D02 = 1590.4 mm2.
A shell of outer diameter dsout = D0 = 45 mm and inner diameter dsin is denoted as corresponding to a cylinder of diameter dc, if dsin = dc. The summary cross-sectional area of a cylinder and a corresponding shell equals the sensor’s active area, providing 100% coverage, Figure 4, as follows:
Sc + Ss = S0 and kc + ks = 1.0.
The radius of a cylinder and the inner radius of a shell are related to the coverage coefficients kc and ks as follows:
r c = 1 / 2 D 0 k c   and   r sin = 1 / 2 D 0 1 k s i n .
The dependence of kc and ks and their change rates kc′ = d k c d r c = 2 r c R 0 2 and ks′ = d k s d r s i n = 2 r s i n R 0 2 on the radial dimension of the subsamples was evaluated.

2.8. The Measured Permittivity of Concentric Circular Subsamples

A complete cylindrical sample of PU foam SikaBlock-M240 (density of 230 kg/m3, thickness t = 20 mm, diameter d = D0 = 45 mm, and cross-sectional area S = S0 = 1590.4 mm2) was made, and its true permittivity was measured. The sample was cut into four concentric circular subsamples—a cylinder and three cylindrical shells, Figure 5. The dimensions of the subsamples were calculated to ensure an equal cross-sectional area S = S0/4 = 397.6 mm2, neglecting the finite width of the cutting lines. Then, the actual dimensions of the subsamples were measured, and the coverage coefficients were calculated. The measured permittivity εm was determined for each subsample. To estimate the impact of PU foams, lost as sawdust, the four subsamples were put one into another to form a combined, quasi-complete sample, and its measured permittivity εm was determined.
The true and the partial electric susceptibility χt and χp were calculated from experimental data of the true and the measured permittivity εt and εm.

2.9. The Measured Permittivity at a Partial Coverage

To study the partial electric susceptibility of rigid PU foams and monolithic polyurethanes, with a dependence on the coverage coefficients of the sensor, two identical right, circular, and cylindrical complete samples (i = 0 and j = 0, thickness t = 20 mm, and diameter D0 = 45 mm) were made from each dielectric material, as shown in Figure 6.
The true permittivity εt of both complete samples was measured. Then a concentric, circular, and cylindrical subsample (i = 1) was made by removing a shell of wall thickness t from outside of the first complete sample, as in Figure 6a, and the measured permittivity of the subsample was determined. The next subsample (i = 2) was made, and its measured permittivity was determined, etc. A concentric, circular, and cylindrical shell was made by removing a cylinder of radius t from the centre of the second complete sample, as in Figure 6b, and the measured permittivity of the subsample j = 1 was determined, etc. Then the next shell (j = 2) was made by removing a shell of wall thickness t from the centre of the shell j = 1, and the permittivity of the shell, j = 2, was measured, etc. Altogether, 5 … 6 subsamples were made from each complete sample.
The measured permittivity of lab-made PU foams (density 32 and 78 kg/m3) and industrial PU foam SikaBlock-M80 (85 kg/m3), SikaBlock-M150 (144 kg/m3), SikaBlock-M240 (226 kg/m3), and SikaBlock-M450 (415 kg/m3) was determined for cylinders. The measured permittivity of the shells was determined for light- to medium-weight PU foam SikaBlock-M80 (85 kg/m3), SikaBlock-M150 (144 kg/m3), and SikaBlock-M450 (415 kg/m3).
The measured permittivity of monolithic dielectrics was determined for the cylinders, made of the lab-made monolithic polyurethane (1280 kg/m3), the industrial monolithic polyurethane SikaBlock-M945 (1352 kg/m3), and epoxy Lab-975New (708 kg/m3). The shell-shaped subsamples of these high-density materials were not made due to technical challenges. The thickness of lab-made monolithic PU complete samples was t = 12 mm.
The spectrometer and the test materials were situated in the same premises for the entire study to attain thermodynamic equilibrium; no conditioning was made for the samples. The permittivity spectra were approximated according to methodology [8]. The true and the partial electric susceptibility, χt and χp, were calculated from experimental data of the true and the measured permittivity values, εt and εm, of the listed dielectric materials.
Data, corresponding to a frequency f = 1/2(f7 + f8) = 960 Hz ≈ 1 kHz, are given in this paper. For two density cases, (1) light-weight PU foam SikaBlock-M80 (85 kg/m3) and (2) Monolithic PU SikaBlock-M945 (1352 kg/m3) data, corresponding to f13 = 40,960 Hz, are given as well to evaluate the dielectric dispersion [26].

3. Theoretical

3.1. Model Functions

The true and the partial susceptibility of a dielectric material, χt and χp, are related to the true and the measured permittivity, εt and εm, as follows:
χt = εt − 1.0 and χp = εm − 1.0.
The experimental data curves εm = εm(r) of cylinders and shells were recalculated to χp = χp(r), and the relative partial susceptibility RΧ was implemented by normalising the partial susceptibility χp(r) with χt:
R X r c = χ p c r c χ t = ε m c r c 1.0 ε t 1.0   and   R X r s i n = χ p s r s i n χ t = ε m s r s i n 1.0 ε t 1.0 ,
where rc or rsin are the radiuses of cylinders and shells. With such a definition, RX characterises the fraction, which a subsample’s partial susceptibility forms from the true susceptibility of a complete sample. The values of the relative partial susceptibility lie in a common scale of 0.0 ≤ RΧ(rc) and RΧ(rsin) ≤ 1.0. Analysis of calculated RΧ(rc) and RΧ(rsin) data curves identified an inflection point at re2 < rinfl < re1in, where re1in is the inner radius of the electrode №1 and re2 is the radius of the electrode №2.
To model the calculated data curves RX(rc) and RΧ(rsin), different functions were tested, but none of them described the relationships properly. Therefore, the model function is proposed as a combination of two normalized functions: (1) a cumulative normal distribution and (2) a cumulative lognormal distribution. For cylinders, the model function Φ(rc) ≈ RX(rc) equals the following:
f 1 ( r c ) = 1 2 π σ 1 0 r c e ( r μ 1 ) 2 2 σ 1 2 d r  at  r c < r infl  and f 2 ( r c ) = 1 2 π σ 2 r c r c R 0 e [ l n ( r ) μ 2 ] 2 2 σ 2 2 d r  at  r c r infl ,
where parameters μ1, μ2 are mean values and σ1, σ2 are standard deviations. The functions f1(rc) and f2(rc) are joined at the intersection point, which coincides with the inflection point. To join f1(rc) and f2(rc) smoothly, their average value can be assigned to RX at the joining point, when necessary.

3.2. Change Rate of the Relative Partial Susceptibility

In order to characterise the change rate of the relative partial susceptibility RX, with a dependence on the radius of dielectric cylinders, ΔRX(rc)/Δrc can be estimated from the model function Φ(rc) as follows:
ΔRX(rc)/Δrc ≈ ΔΦ(rc)/Δrc.
But, with an increase in the cylinder’s radius rc, the coverage coefficient kc increases nonlinearly, as in Equation (2), and the increments of the coverage coefficient, which correspond to constant increments of radius Δrc = const., are radius-dependent:
k c Δ k c Δ r c = 2 r c R 0 2   and   Δ k c = 2 r c r c R 0 2 .
Therefore, estimation (8) is not fit for the purpose, and the change rate of RX has to be determined at constant increments of the OSA sensor’s coverage. Substituting r c = R 0 k c into Equation (7), we obtain the following:
ΔRX(kc)/Δkc ≈ ΔΦ(kc)/Δkc.
This corresponds to a constant increase in the cross-sectional area of the cylinders: when Δkc = const., then ΔSc = S0Δkc = const., as shown in Equation (2).
Analysis of the calculated data curves RX(rc) and RΧ(rsin) of low- to medium-density PU foams (85 kg/m3, 144 kg/m3, and 415 kg/m3) identified a relationship between the relative partial susceptibilities of corresponding cylinders and shells:
RX(rc) + RX(rsin) = 1.0.
Then, taking into account rsin = rc for the corresponding cylinders and shells, the model function for shells Ψ(rsin) ≈ RX(rsin) can be constructed as follows:
f 3 ( r sin ) = 1.0 f 1 ( r sin ) = 1.0 1 2 π σ 1 0 r s i n e ( r μ 1 ) 2 2 σ 1 2 d r  at  r sin < r infl  and f 4 ( r sin ) = 1.0 f 2 ( r sin ) = 1.0 1 2 π σ 2 r s i n r s i n R 0 e [ l n ( r ) μ 2 ] 2 2 σ 2 2 d r  at  r sin r infl .
For shells, constant increments of radius Δrsin = const. provide radius-dependent increments of coverage coefficient ks as follows:
k s Δ k s Δ r s i n = 2 r s i n R 0 2 ,   then   Δ k s = 2 r s i n r s i n R 0 2 .
Therefore, the change rate of RX has to be estimated at constant increments of the OSA sensor’s coverage with dielectric shells:
ΔRX(ks)/Δks ≈ ΔΨ(ks)/Δks.
A constant change in the cross-sectional area of a shell is ensured: when Δks = const., then ΔSs = S0Δks = const., as shown in Equation (2). The ΔRX(ks)/Δks curves for shells, cut from low- to medium-density PU foams (85 kg/m3, 144 kg/m3, and 415 kg/m3), were calculated from model functions according to Equation (14).
Due to a lack of experimental data on the measured permittivity for shells of the high-density PU materials, the corresponding data curves RΧ(rsin) could not be calculated. Consequently, Equation (14) could not be used to calculate the change rate of the relative partial susceptibility for shells, i.e., ΔRX(ks)/Δks. As an alternative, with ΔΦ(kc)/Δkc known for cylinders, its relation to ΔΨ(ks)/Δks for shells was found. Taking into account RX(rc) ≈ Φ(rc) and RX(rsin) ≈ Ψ(rsin), we obtain the following from Equation (11):
Φ(rc) + Ψ(rsin) = 1.0.
Substituting r c = R 0 k c and r s i n = R 0 1.0 k s into Equation (15) leads to the following:
Φ(kc) = 1.0 − Ψ(ks).
Let us give both sides of Equation (16) a finite increment, ΔΦ(kc) = −ΔΨ(ks), and then divide by Δkc, Φ ( k c ) k c = Ψ ( k s ) k c . If the coverage coefficient of a cylinder increases, the coverage coefficient of the corresponding shell decreases for the same amount: Δkc = −Δks. Then Φ ( k c ) k c = Ψ ( k s ) k s , which permits us to calculate ΔΨ(ks)/Δks when ΔΦ(kc)/Δkc is known.
In numerical calculations, the functions f1(rc) and f2(rc) were determined as normalised, cumulative standard functions NORMDIST(rc, μ1, σ1, TRUE) and LOGNORM.DIST(rc, μ2, σ2) of the MS EXCEL software v.12 (Microsoft Corporation; Redmond, WA, USA). A transformation of the radial coordinate rc′ = rc + rT permitted us to translate the function f2(rc) for a distance rT to the best-fitting position. The parameters of f1(rc) and f2(rc), ensuring the best fitting of the calculated RX data, were determined for each dielectric material. Similar calculations were made for shells.

3.3. Increments of Radius at Constant Increments of Coverage Coefficient

To illustrate the distribution of the change rate of the relative partial susceptibility over the radius of the sensor’s zone, covered with a dielectric cylinder, the curves “ΔRX(kc)/Δkc − kc” were recalculated to such rc values, which correspond to constant increments of the coverage coefficient Δkc = kc(i+1) − kci = const. Here “kci” and “kc(i+1)” are points on the kc − axis, i = 1, …, I + 1 and I = kc/Δkc = 1.0/Δkc. Taking into account that kci = (rci/R0)2, the corresponding increments of the radius of a cylinder were calculated as follows:
Δ r ci ( r c ) = r c i ± r c i 2 + R 0 2 k c ,  where i = 1 ,   2 ,   ,   I .
Since rc ≥ 0.0 mm for all 0.0 mm ≤ rc ≤ R0 mm, the positive square root is chosen. Then, the values of Δrci depend on rc, but in a way that ensures Δkc = const., as in Table 1. Similar calculations were made for the curves “ΔRX(ks)/Δks − ks” in the case of shells.
Values of the change rate ΔRX(kc)/Δkc were calculated numerically at Δkc = 0.05 (5%) for cylinders of PU materials with significantly differing permittivity.

4. Results and Discussion

4.1. Dielectric Losses

For PU foams with densities of 95–222 kg/m3, dielectric losses were measured as ε″(f) = 0.0022–0.0063 (1 kHz) and ε″(f) = 0.0032–0.0084 (0.1 MHz). For monolithic lab-made polyurethane, ε″(f) = 0.042 at 1 kHz and ε″(f) = 0.088 at 0.1 MHz. The acquired data are in good correspondence with the experimental data in [26]. For the monolithic PU, ε′(f) = 3.42 (1 kHz) and the dielectric loss tangent tgδ = ε″(f)/ε′(f) = 0.012 << 1.0. At higher frequencies of the considered range, f = 10 Hz–0.33 MHz, as in Equation (1), the dielectric losses are even smaller. The loss part of the PU materials, ε″(f), is small compared to the real one, ε′(f). Then ε ~ (jf) ≈ ε′(f), and the real part ε′(f) = ε(f) is referred to as permittivity.

4.2. Lateral Dimensions of Complete Samples

For PU foams of densities 50–228 kg/m3 and εt = 1.14–1.42 (1 kHz), the penetration depth was determined as 5.72 mm ≤ t3% ≤ 5.87 mm ± 0.02 mm. That corresponds to the conclusions in [7]:, according to the definition of penetration depth, which is proposed in [7], the penetration depth of the OSA sensor increases as the permittivity of the sample increases. In order to achieve the same 3% difference, samples with high values of the true permittivity εt need to have a bigger penetration depth, t3%, whereas samples with low εt values can have smaller t3% values to achieve the same percentage of difference [8]. The samples have to be thick enough to provide the true permittivity; therefore, 3–4 times the thickness of the penetration depth t3% was taken as appropriate for PUR foams’ samples with densities of 50–1280 kg/m3 and t ≈ 20–25 mm [8]. Due to the limited dimensions of the lab-made monolithic PU rods, the semi-cylindrical samples were made with a thickness of 12 mm.

4.3. Transversal Dimensions of Complete Samples

The measured permittivity of PU foams and the expanded uncertainty U95% of the experimental data point are practically equal for all three considered kinds of samples: “(1)”, “(2)”, and “(3)”. The differences lie in the limits of the uncertainty of the spectrometer US95.45% = ±0.01, as shown in Table 2. No correlation was identified between the size/shape of the samples and the measured permittivity.
It is concluded that the measured permittivity εm, determined for parallelepiped-shaped samples, which ensure a 100% coverage of the sensor’s active area, differs insignificantly from the true permittivity value εt, measured for complete cylindrical samples, exactly matching the sensor’s active area. All three kinds of samples can be used to measure the true permittivity of PU foams.
For the parallelepiped-shaped PU foam samples with a concentric cylindrical hole (D0 = 45 mm), the readings of the spectrometer were equal to those acquired in measurements with no sample on the sensor (“Sensor in air”): ε ≈ 1.00 ± 0.01. PU foams, located outside the cylindrical zone above the active area of the sensor, have practically no effect on the measured permittivity value.

4.4. Inner Vertical Surfaces

For both PU materials, monolithic and cellular, the results at each frequency fn (data at f1 = 10 Hz, f8 = 1280 Hz, and f16 = 327,680 Hz are displayed) showed a small relative difference between the permittivity of a complete cylindrical sample εt and the measured permittivity of two semi-cylindrical samples εmc: R = (εt − εmc)/εt; |R| < 0.5%, Table 3.
It was concluded that a similar result could be expected for the lab-made monolithic PU as well. The inner vertical surfaces between the semi-cylinders have a negligible impact on the measured permittivity of a combined sample.

4.5. Coverage Coefficients

At an increase in the radius of the cylinders rc and the inner radius of the shells rsin from 0.0 mm to R0 = 22.5 mm, the coverage coefficient of the cylinders kc increases nonlinearly from 0% to 100%, while that of shells ks decreases in the same way, as in Figure 7. The change rate kc′ increases linearly, and ks′ decreases by the same amount.

4.6. The Measured Permittivity of Concentric Circular Subsamples

The theoretical parameters of the PU foam SikaBlock-M240 (230 kg/m3) sample and subsamples are given in Table 4. The actual dimensions, in Table 5, can differ slightly. Around 5% of the cross-sectional area of a complete sample is changed into sawdust while cutting the three circular lines, each of width w ≈ 0.3 mm. As a result, the actual coverage coefficients kc and ks differ from those calculated with the assumption w = 0.0 mm.
Analysis of experimental data showed that the true permittivity of a complete cylindrical sample can be calculated from the measured permittivity of the four concentric circular subsamples, each of them covering the sensor’s active area partially, as follows:
ε t ε = ε mc + j = 1 3 ε m s j 3.0 = 1.38 .
where εmc and εms are the measured permittivities of the subsamples. The relative difference between the true permittivity of the complete cylinder and (1) the calculated permittivity ε′ is estimated as R1 = (εt − ε′)/εt = (1.40 − 1.38)/1.14 ≈ 1.4% and (2) the measured permittivity of the combined sample εm as R2 = (εt − εm)/εt < 1.0%. The true susceptibility χt of the complete cylindrical sample was determined from the calculated partial susceptibilities of the four concentric circular subsamples as follows:
χ pc = ε mc 1.0 ;   χ ps = ε ms 1.0 ,  then χ t χ = χ pc + j = 1 3 χ p s j + 1.0 = 1.38 .
It can be seen that, although the subsamples have similar coverage coefficients kc ≈ ks, the partial susceptibilities differ considerably, as shown in Table 5. The cylindrical subsample provides ~90% of the true susceptibility, while the three shells only ~10%, indicating an uneven distribution of the electric charge over the active area of the sensor. In a general case, for J subsamples (J − 1 shells and a cylinder), the following is calculated:
ε t ε = ε mc + j = 1 J 1 ε m s j J   and   χ t χ = χ pc + 1 J 1 χ p s + 1.0 .

4.7. The Measured Permittivity and the Partial Susceptibility at a Partial Coverage

4.7.1. PU Foams with Similar True Permittivity

The experimental data of the measured permittivity of low- to medium-density PU foams (32–415 kg/m3), with similar true permittivity 1.00 < εt < 1.80, are depicted in Figure 8 with a dependence on the radius of the cylinders.
At an increase of rc from 0.0 mm to R0 = 22.5 mm, the measured permittivity εmc increases from 0.0 to the value of the true permittivity εt. The concave–onvex curves “εmc − rc” have inflection points, located at the limits of re2 ≤ rc ≤ re1in, where re1in is the inner radius of electrode №1, and re2 is the radius of electrode №2. The dashed verticals in Figure 8 and other Figures mark the location of the radiuses re1in and re2 of the sensor’s electrodes №1 and №2 as well as the radius, re3, of the thickness centre of the guard electrode, as in Table 6. Coverage coefficients are given for cylinders in Table 6, when the radius of a cylinder rc = re1in, rc = re1out, rc = re2, and rc = re3, as well as for shells, when the inner radius of a shell rsin = re1in, rsout, rsin = re2, and rsin = re3.
The cross-sectional area of a cylinder increases, as its radius increases, from the centre of the cylinder to its outer perimeter. The cross-sectional area of a shell increases as the inner radius of the shell decreases, in the direction from the outer perimeter of the shell to its centre. The ring-shaped area between the electrodes №1 and №2 corresponds to ≈7.9% of the active area (for a cylinder 21.8% − 13.9% = 7.9% as well as for a shell, 86.1% − 78.2% = 7.9%).
Figure 9 gives experimental data of the measured permittivity εmc with a dependence on radiuses rc and rsin of the cylinders and shells.
The curves “εmc − rc” and “εms − rsin are symmetrical about a straight line, drawn parallel to the Orc (Orsin) axis and passing through the intersection point of the curves. Then, for each radius of a cylindrical subsample, 0.0 mm ≤ rc ≤ 22.5 mm, the following relation holds:
εc(rc) + εs(rsin) = εt.
It can be seen that the true permittivity of PU foams can be determined as a sum of the measured permittivity of the corresponding circular subsamples: a cylinder and a shell. The intersection points of the curves practically coincide with the inflection points.
The transition to the calculated relative partial susceptibility RX, as in Equations (6), brings the experimental data from Figure 8 and Figure 9 to a common scale of 0.0 ≤ RΧ(rc) and RΧ(rsin) ≤ 1.0, as in Figure 10 and Figure 11.
The parameters of the best-fitting model functions Φ(rc) and Ψ(rsin) of the calculated relationships “RX(rc) − rc” and “RX(rsin) − rsin” are given in Table 7. Theoretically, the parameters of Ψ(rsin) have to be the same as those of Φ(rc), as in Equations (7) and (10), but practically, they can differ due to irregularities in the subsamples’ shape, density, etc. For PU foams of a similar true permittivity, 1.00 < εt < 1.80, the curves “RX(rc) − rc” and “RX(rsin) − rsin” have a similar shape, are close together, and partially overlap, as in Figure 10 and Figure 11, which hinders the determination of the parameters and the identification of the main trends.
It can be seen that the model curves “Φ(rc) − rc” and “Ψ(rsin) − rsin” for cylinders and shells of both PU foams are symmetrical about a straight line, which is drawn parallel to the Orc (Orsin) axis and passes through the intersection of the curves at Φ(rc) = Ψ(rsin) = 50%, as in Figure 12. Then, for each radius of a cylinder, 0.0 mm ≤ rc ≤ 22.5 mm, and the radius of the cylinder’s corresponding shell, rsin = rc, the following hold:
Φ(rc) + Ψ(rsin) = 1.0 and RX(rc) + RX(rsin) = 1.0.
For each cylinder of a coverage coefficient kc 0.0% ≤ kc ≤ 100% and the cylinder’s corresponding shell of a coverage coefficient ks, kc + ks = 100% is valid, as shown in Equation (3). Then, from Figure 13 we obtain the following:
Φ(kc) + Ψ(ks) = Φ(kc) + Ψ(100% − kc) = 100%.

4.7.2. PU Materials with Significantly Differing True Permittivity

The experimental data of the measured permittivity of dielectrics (PU foams, monolithic polyurethane, and an epoxy) with significantly differing true permittivity (1.0 < εt < 9.0) are given in Figure 14 with a dependence on the radius of the cylinders. With an increase of rc from 0.0 mm to R0 = 22.5 mm, the measured permittivity εmc increases from 0 to εt.
The εm = εm(rc) curves, registered for the light-weight PU foam SikaBlock-M80 (85 kg/m3) at frequencies 1000 Hz and 40,960 Hz, practically overlap, because the differences Δεm between the εm values at both frequencies are less than 0.001. For the monolithic polyurethane SikaBlock-M945 (1352 kg/m3), the εm = εm(rc) curve at 40,960 Hz is given in Figure 14f and Δεm < 0.20. When the εm = εm(rc) curves are recalculated to the relative partial susceptibility RX, as in Equations (5) and (6), the RX = RX(rc) curves at 1000 Hz and 40960 Hz coincide closely for both PU materials and provide no new information. Therefore, the curves RX = RX(rc) at f13 = 40,960 Hz are not depicted in Figure 15. It is reasonable to expect the same for PU foams of density values below the density of monolithic PU SikaBlock-M945. It is concluded that the studied PU materials exhibit low dielectric dispersion.
The calculated data curves of the relative partial susceptibility RX(rc) follow a similar, concave/convex pattern, with inflection points at the limits, re2 ≤ rc ≤ re1in, as in Figure 15. For dielectrics with higher true permittivity, the inflection points correspond to higher values of the radius of the cylinders, rc.
The εt values of the considered dielectrics are sufficiently different to prevent overlapping of the curves RX(rc) and to facilitate a determination of the main trends and parameters of model functions, as shown in Table 8.
The main features of the calculated data curves RX(rc), as in Figure 15, are valid for the model functions Φ(rc) of the relative partial susceptibility as well, as in Figure 16. Compared to the PU foam SikaBlock-M80 (εt = 1.15), the inflection point of the Φ = Φ(rc) curve for epoxy LAB975 (εt = 8.95) is located ≈2 mm further in the direction of the higher values of rc. This can be explained by a higher input of the outer circular zones (covered with the dielectric) of the electrode №1 into the total transcapacitance. At higher values of the true permittivity, the input becomes sufficient to be detected by a spectrometer of a certain sensitivity.
For PU foams with densities of 78–85 kg/m3, the true permittivity 1.13 ≤ εt ≤ 1.15 is only ≈13–15% higher than that of vacuum εt0 = 1.0. The outer circular zones of the electrode №1 (rc > 12 … 15 mm) have a lower charge density, compared to that of the inner zones [2]; as a result, their input is comparatively small.
The higher the true permittivity, the higher the coverage coefficient, at which a certain value of the relative partial susceptibility RX is reached: (a) 50%, (b) 75%, and (c) 90%, as in Figure 17.
The mentioned trend is displayed in Figure 18.
This may be explained by a higher input of the outer circular zones of electrode №1 into the total transcapacitance, when at a higher true permittivity, the input becomes sufficient to be detected by a spectrometer of a certain sensitivity.
The coverage of the sensor’s active area increases from 0% to 100% in the direction from (1) the centre to the perimeter at an increase of a cylinder’s radius from 0 to R0, and when Δrc > 0, then Δkc > 0, and in the direction from (2) the perimeter to the centre at a decrease of a shell’s inner radius from R0 to 0, and when Δrsin < 0, then Δks > 0. The values of the change rate ΔRX(kc)/Δkc were calculated numerically at Δkc = Δks = 0.05 (5%).
In the centre of the OSA sensor, in Figure 19, at kc = 0% for cylinders, the change rate is ΔRX/Δkc ≈ 4–10%/5%, because the change rate value ΔRX/Δkc, calculated for each coverage interval Δkc = kc(i+1) − kci = 5%, is assigned to its left point kci. ΔRX/Δkc has the highest values of 20–25%/5% at 14% < kc < 22%, which corresponds to the zone between electrodes №1 and №2. In the centre of the sensor, ΔRX/Δkc is 2–5 times smaller, and it decreases to < 5%/5% as kc increases above 35%. For shells, in the centre of the OSA sensor at ks = 100%, the corresponding value of the change rate ΔRX/Δks ≈ 4–10%/5% is assigned to the right point of each Δks interval. ΔRX/Δks has the highest values of 20–25%/5% at 78% < ks < 86% and decreases to less than 5%/5% at ks < 65%, as shown in Table 6.
Figure 20 shows the change rate of the relative partial susceptibility ΔRX/Δkc over the radius of the sensor’s circular concentric zone, which is covered with a dielectric cylinder of radius rc. The increments Δrci, corresponding to the markers on the curves, depend on rc in a way that ensures constant increments of the coverage coefficient Δkc = 5%, as in Equation (17). It can be seen that the change rate is the highest, ΔRX/Δkc = 20–25%/5%, in the ring-shaped zone of radius 8.4 mm < rc < 10.5 mm between electrodes №1 and №2. In the centre of the sensor, the change rate is 2–5 times smaller. As the radius rc increases above 12 mm, ΔRX/Δkc decreases to <5%/5%.
The numerical results of several mathematical models of other authors are available for circular OSA sensors with two electrodes in [2,7,13,15]. The overall character of the curves in Figure 20 like slopes of lines and the ratio of ΔRX/Δkc values in the centre and near the gap corresponds with the character of the surface charge density distribution curves, calculated from mathematical modelling. In [11], a simplification is made by assuming an even distribution of the surface charge on the sensor’s electrodes. In [12,13], a complex mathematical model for a circular three-electrode OSA sensor is elaborated, but numerical results for the surface charge are not displayed.
A shortage of experimental data for the surface charge density distribution over the active area of circular OSA sensors hinders a relevant comparison.

5. Conclusions

The electric susceptibility at the partial coverage of a circular OSA sensor with cylinders and shells is investigated for PU foams and monolithic polyurethane. It is shown experimentally that the true susceptibility can be determined from the partial susceptibilities of the corresponding subsamples. The implementation of the relative partial susceptibility permitted us to transform the calculated susceptibility data to a common scale of 0.0–1.0 and to outline the main trends for PU materials.
The partial susceptibility, the relative partial susceptibility, and the change rate of the relative partial susceptibility exhibit a marked dependence on the coverage coefficient of the OSA sensor’s active area as well as a correlation with the coordinates of the sensor’s electrodes and the gap. Numerical calculations showed that, for cylindrical subsamples, the change rate of the relative partial susceptibility is the highest in the zone of the active area between electrodes №1 and №2: 20–25%/5%. In the centre of the sensor, it is 2–5 times smaller and decreases below 5%/5% as the coverage coefficient increases above 35% in correspondence with the surface charge density distribution, reported in scientific information sources. The curves for shells are symmetric to those for the cylinders with respect to the location of the gap. The overall pattern of curves for the change rate of the relative partial susceptibility, characterised by slopes of lines and the ratio of ΔRX/Δkc values in the centre and values near the gap, corresponds with the character of the surface charge density distribution curves, calculated from mathematical models.
The uncertainty of the spectrometer, geometric irregularities of subsamples, and varying humidity of the ambient atmosphere can be named as the main sources of measurement uncertainties for PU foams. The elaborated methods can be applied in the design and optimization of capacitive OSA sensors of other configurations of electrodes, independent of the particular technical solution.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within this article.

Conflicts of Interest

The author declares no conflicts of interest.

References

  1. Kremer, F.; Schonhals, A. (Eds.) Broadband Dielectric Spectroscopy; Springer: Berlin/Heidelberg, Germany, 2003; 730p. [Google Scholar] [CrossRef]
  2. Matiss, I. Capacitive Transducers for Non-Destructive Testing; Zinatne: Riga, Latvia, 1982; 302p. (In Russian) [Google Scholar]
  3. Kalpinsh, A.; Shtrauss, V.; Lomanovskis, U. Digital Emulation of Dielectric Relaxation Functions for Capacitive Sensors of Non-Destructive Dielectric Spectrometry. Comput. Methods Exp. Meas. 2019, 125, 111–119. [Google Scholar] [CrossRef]
  4. Baxter, L.K. Capacitive Sensors: Design and Applications; Wiley Institute IEEE Press, Inc.: New York, NY, USA, 1997; 277p, 302p with Appendices. [Google Scholar]
  5. Fraden, J. Handbook of Modern Sensors, 4th ed.; Springer: Berlin/Heidelberg, Germany, 2010; 663p. [Google Scholar] [CrossRef]
  6. Chen, T.; Bowler, N. Analysis of a concentric coplanar capacitive sensor for nondestructive evaluation of multi-layered dielectric structures. IEEE Trans. Dielectr. Electr. Insul. 2010, 17, 1307–1318. [Google Scholar] [CrossRef]
  7. Chen, T. Capacitive Sensors for Measuring Complex Permittivity of Planar and Cylindrical Structures. Ph.D. Thesis, Iowa State University, Ames, IA, USA, 2012. Available online: https://lib.dr.iastate.edu/etd/12294 (accessed on 11 March 2024).
  8. Beverte, I.; Cabulis, U.; Gaidukovs, S. Polytetrafluoroethylene Films in Rigid Polyurethane Foams’ Dielectric Permittivity Measurements with a One-Side Access Capacitive Sensor. Polymers 2021, 13, 1173. [Google Scholar] [CrossRef] [PubMed]
  9. Chen, T.; Bowler, N. Analysis of concentric coplanar capacitor for quantitative dielectrometry. Stud. Appl. Electromagn. Mech. 2010, 33, 61–68. [Google Scholar] [CrossRef]
  10. Chen, T.; Song, J.M.; Bowler, J.R.; Bowler, N. Analysis of a concentric coplanar capacitive sensor using a spectral domain approach. Rev. Prog. Quant. Nondestr. Eval. 2011, 30, 1647–1654. [Google Scholar]
  11. Cheng, H.; Zhang, Y.; Huang, X.; Rogers, J.A.; Huang, Y. Analysis of a concentric coplanar capacitor for epidermal hydration sensing. Sens. Actuator A Phys. 2013, 203, 149–153. [Google Scholar] [CrossRef]
  12. Sheiretov, Y. Deep Penetration Magnetoquasistatic Sensors. Ph.D. Thesis, Massachusetts Institute of Technology, Cambridge, MA, USA, 2001. Available online: https://web.mit.edu/yksheire/www/thesis.pdf (accessed on 11 March 2024).
  13. Shay, I.C.; Zahn, M. Cylindrical geometry electroquasistatic dielectrometry sensors. IEEE Trans. Dielectr. Electr. Insul. 2005, 12, 41–49. [Google Scholar] [CrossRef]
  14. Sheiretov, Y.; Zahn, M. Modeling of spatially periodic dielectric sensors in the presence of a top ground plane bounding the test dielectric. IEEE Trans. Dielect. Electr. Insul. 2005, 5, 993–1004. [Google Scholar] [CrossRef]
  15. Bowler, N.; Chen, T. Concentric Coplanar Capacitive Sensor System with Quantitative Model. Patent US 8791707 B2, 29 July 2014. [Google Scholar]
  16. Chu, Z.; Peng, C.; Ren, R.; Ling, B.; Zhang, Z.; Lei, H.; Xia, S. A High Sensitivity Electric Field Microsensor Based on Torsional Resonance. Sensors 2018, 18, 286. [Google Scholar] [CrossRef] [PubMed]
  17. Guo, J.; Hu, P.; Tan, J. Annular Coplanar Capacitive Tilt Sensor. Sensors 2016, 16, 133. [Google Scholar] [CrossRef] [PubMed]
  18. Hu, X.; Yang, W. Planar capacitive sensors–designs and applications. Sens. Rev. 2010, 30, 24–39. [Google Scholar] [CrossRef]
  19. Gong, Y.; Zhu, H.; Li, N. Research on Coplanar Capacitive Sensor Design for Film Thickness Measurement. Adv. Mater. Res. 2014, 945–949, 2030–2036. [Google Scholar] [CrossRef]
  20. Zaitsev, I.O.; Levytskyi, A.S. Determination of response characteristic of capacitive coplanar air gap sensor. In Proceedings of the 5th IEEE Microwaves, Radar and Remote Sensing Symposium, Kyiv, Ukraine, 29–31 August 2017. article number 8075034, code 131255. [Google Scholar] [CrossRef]
  21. Zaitsev, I.O.; Levytskyi, A.S.; Novik, A.I.; Bereznychenko, V.O.; Smyrnova, A.M. Research of a capacitive distance sensor to grounded surface. Telecommun. Radio Eng. Elektrosvyaz Radiotekhnika 2019, 78, 173–180. [Google Scholar] [CrossRef]
  22. Abdollahi-Mamoudan, F.; Savard, S.; Ibarra-Castanedo, C.; Filleter, F.; Maldague, X. Influence of different design parameters on a coplanar capacitive sensor performance. NDT Int. 2022, 126, 102588. [Google Scholar] [CrossRef]
  23. Hilyard, N.C. Mechanics of Cellular Plastics; MacMillan: New York, USA, 1982; 360p. [Google Scholar]
  24. Klempner, D.; Frisch, K.C. Handbook of Polymeric Foams and Foam Technology; Hanser Publishers: Munich, Germany, 1991; 413p. [Google Scholar]
  25. Kirpļuks, M. Development of Renewable Feedstock Based Rigid Polyurethane Foam and Nanoclay Composites. Ph.D. Thesis, Riga Technical University, Riga, Latvia, 2020. Available online: https://ortus.rtu.lv/science/en/publications/31088 (accessed on 11 March 2024).
  26. Beverte, I.; Shtrauss, V.; Kalpinsh, A.; Lomanovskis, U.; Cabulis, U.; Sevastyanova, I.; Gaidukovs, S. Dielectric Permittivity of Rigid Rapeseed Oil Polyol Polyurethane Biofoams and Petrochemical Foams at Low Frequencies. J. Renew. Mater. 2020, 8, 1151–1170. [Google Scholar] [CrossRef]
  27. Beverte, I.; Gaidukovs, S.; Andersons, J.; Skruls, V. The Impact of Atmospheric Parameters on the Dielectric Permittivity Values of SikaBlock Block®-M150 and Other Rigid Polyurethane Foams Measured with a Capacitive One-Side Access Sensor. Sensors 2022, 22, 7859. [Google Scholar] [CrossRef] [PubMed]
  28. ISO 845:2006; Cellular Plastics and Rubbers—Determination of Apparent Density. ISO: Geneve, Switzerland, 2006.
  29. Beverte, I. Determination of Highly Porous Plastic Foams’ Structural Characteristics by Processing LM Images Data. J. Appl. Polym. Sci. 2014, 131, N39477. [Google Scholar] [CrossRef]
  30. Beverte, I.; Cabulis, U.; Andersons, J.; Kirpluks, M.; Skruls, V.; Cabulis, P. Anisotropy and Mechanical Properties of Nanoclay Filled, Medium Density Rigid Polyurethane Foams Produced in a Sealed Mold, from Renewable Resources. Polymers 2023, 15, 2582. [Google Scholar] [CrossRef]
  31. Evaluation of Measurement Data—Guide to the Expression of Uncertainty in Measurement. JCGM 100:2008 GUM 1995 with Minor Corrections, 120p. Available online: https://www.bipm.org/utils/common/documents/jcgm/JCGM_100_2008_E.pdf (accessed on 18 March 2024).
  32. International Vocabulary of Metrology—Basic and General Concepts and Associated Terms (VIM), 3rd Edition 2008 Version with Minor Corrections. JCGM 200:2012. Available online: https://www.bipm.org/utils/common/documents/jcgm/JCGM_200_2012.pdf (accessed on 18 March 2024).
Figure 1. The electrodes: (1) Driven (D1in = 21.0 mm, D1out = 45.0 mm), (2) Sensing (diameter D2 = 16.8 mm), and (3) Guard (width of 1.00 mm) as well as (4) Insulator (width of 0.55 mm) and (5) Grounded screen of the OSA sensor.
Figure 1. The electrodes: (1) Driven (D1in = 21.0 mm, D1out = 45.0 mm), (2) Sensing (diameter D2 = 16.8 mm), and (3) Guard (width of 1.00 mm) as well as (4) Insulator (width of 0.55 mm) and (5) Grounded screen of the OSA sensor.
Sensors 24 03003 g001
Figure 2. PU foam samples on the active area of the sensor: (a) A cylindrical complete sample with a diameter d = D0 = 45 mm and a thickness of 20 mm and (b) a parallelepiped with a concentric cylindrical hole.
Figure 2. PU foam samples on the active area of the sensor: (a) A cylindrical complete sample with a diameter d = D0 = 45 mm and a thickness of 20 mm and (b) a parallelepiped with a concentric cylindrical hole.
Sensors 24 03003 g002
Figure 3. Complete samples with inner vertical surfaces on the active area of the sensor: (a) Lab-made monolithic PU and (b) PU foam SikaBlock-M450.
Figure 3. Complete samples with inner vertical surfaces on the active area of the sensor: (a) Lab-made monolithic PU and (b) PU foam SikaBlock-M450.
Sensors 24 03003 g003
Figure 4. Cross-section of corresponding circular figures: cylinder 1 and shell 2.
Figure 4. Cross-section of corresponding circular figures: cylinder 1 and shell 2.
Sensors 24 03003 g004
Figure 5. Cross-section of the concentric circular subsamples.
Figure 5. Cross-section of the concentric circular subsamples.
Sensors 24 03003 g005
Figure 6. Complete samples (i = 0 and j = 0) and circular subsamples: (a) cylinders (i = 1 and 2) and (b) cylindrical shells (j = 1 and 2).
Figure 6. Complete samples (i = 0 and j = 0) and circular subsamples: (a) cylinders (i = 1 and 2) and (b) cylindrical shells (j = 1 and 2).
Sensors 24 03003 g006
Figure 7. Coverage coefficients kc and ks (blue) and their change rates k′c and k′s (red), depending on (a) the radiuses of the cylinders and (b) the inner radiuses of the shells.
Figure 7. Coverage coefficients kc and ks (blue) and their change rates k′c and k′s (red), depending on (a) the radiuses of the cylinders and (b) the inner radiuses of the shells.
Sensors 24 03003 g007
Figure 8. The measured permittivity εmc with a dependence on the cylinder’s radius rc: lab-made PU foams of density (a) 32 kg/m3 and (b) 78 kg/m3 and industrial PU foams (c) SikaBlock-M80 (85 kg/m3), (d) SikaBlock-M150 (144 kg/m3), (e) SikaBlock-M240 (226 kg/m3), and (f) SikaBlock-M450 (415 kg/m3) at f = 1 kHz.
Figure 8. The measured permittivity εmc with a dependence on the cylinder’s radius rc: lab-made PU foams of density (a) 32 kg/m3 and (b) 78 kg/m3 and industrial PU foams (c) SikaBlock-M80 (85 kg/m3), (d) SikaBlock-M150 (144 kg/m3), (e) SikaBlock-M240 (226 kg/m3), and (f) SikaBlock-M450 (415 kg/m3) at f = 1 kHz.
Sensors 24 03003 g008
Figure 9. The measured permittivity εmc of PU foams’ cylinders: (a,c,e) as well as εms of shells: (b,d,f), with a dependence on the subsamples’ radiuses rc (circles) and rsin (triangles). Densities: (a,b) SikaBlock-M80 (85 kg/m3), (c,d) SikaBlock-M150 (144 kg/m3), and (e,f) SikaBlock-M450 (415 kg/m3) at f = 1 kHz.
Figure 9. The measured permittivity εmc of PU foams’ cylinders: (a,c,e) as well as εms of shells: (b,d,f), with a dependence on the subsamples’ radiuses rc (circles) and rsin (triangles). Densities: (a,b) SikaBlock-M80 (85 kg/m3), (c,d) SikaBlock-M150 (144 kg/m3), and (e,f) SikaBlock-M450 (415 kg/m3) at f = 1 kHz.
Sensors 24 03003 g009
Figure 10. The relative partial susceptibility RX with a dependence on the cylinder’s radius rc: lab-made PU foams with densities of (a) 32 kg/m3 and (b) 78 kg/m3 and industrial PU foams (c) SikaBlock-M80 (85 kg/m3), (d) SikaBlock-M150 (144 kg/m3), (e) SikaBlock-M240 (226 kg/m3), and (f) SikaBlock-M450 (415 kg/m3) at f = 1 kHz.
Figure 10. The relative partial susceptibility RX with a dependence on the cylinder’s radius rc: lab-made PU foams with densities of (a) 32 kg/m3 and (b) 78 kg/m3 and industrial PU foams (c) SikaBlock-M80 (85 kg/m3), (d) SikaBlock-M150 (144 kg/m3), (e) SikaBlock-M240 (226 kg/m3), and (f) SikaBlock-M450 (415 kg/m3) at f = 1 kHz.
Sensors 24 03003 g010
Figure 11. The relative partial susceptibility RX of PU foams’ cylinders (a,c,e) as well as of the corresponding shells (b,d,f) with a dependence on the subsamples’ radiuses rc (circles) and rsin (triangles). Densities: (a,d) SikaBlock-M80 (85 kg/m3), (b,e) SikaBlock-M150 (144 kg/m3), and (c,f) SikaBlock-M450 (415 kg/m3) at f = 1 kHz.
Figure 11. The relative partial susceptibility RX of PU foams’ cylinders (a,c,e) as well as of the corresponding shells (b,d,f) with a dependence on the subsamples’ radiuses rc (circles) and rsin (triangles). Densities: (a,d) SikaBlock-M80 (85 kg/m3), (b,e) SikaBlock-M150 (144 kg/m3), and (c,f) SikaBlock-M450 (415 kg/m3) at f = 1 kHz.
Sensors 24 03003 g011
Figure 12. Model functions Φ and Ψ of the relative partial susceptibility RX of PU foams’ cylinders (ac) as well as of the corresponding shells (df) with a dependence on the radiuses of subsamples rc and rsin. Densities: (a,d) SikaBlock-M80 (85 kg/m3), (b,e) SikaBlock-M150 (144 kg/m3), and (c,f) SikaBlock-M450 (415 kg/m3) at f = 1 kHz.
Figure 12. Model functions Φ and Ψ of the relative partial susceptibility RX of PU foams’ cylinders (ac) as well as of the corresponding shells (df) with a dependence on the radiuses of subsamples rc and rsin. Densities: (a,d) SikaBlock-M80 (85 kg/m3), (b,e) SikaBlock-M150 (144 kg/m3), and (c,f) SikaBlock-M450 (415 kg/m3) at f = 1 kHz.
Sensors 24 03003 g012
Figure 13. Model functions Φ and Ψ of the relative partial susceptibility RX of PU foams’ cylinders (ac) as well as of the corresponding shells (df) with a dependence on the coverage coefficients kc and ks. Densities: (a,d) SikaBlock-M80 (85 kg/m3), (b,e) SikaBlock-M150 (144 kg/m3), and (c,f) SikaBlock-M450 (415 kg/m3) at f = 1 kHz.
Figure 13. Model functions Φ and Ψ of the relative partial susceptibility RX of PU foams’ cylinders (ac) as well as of the corresponding shells (df) with a dependence on the coverage coefficients kc and ks. Densities: (a,d) SikaBlock-M80 (85 kg/m3), (b,e) SikaBlock-M150 (144 kg/m3), and (c,f) SikaBlock-M450 (415 kg/m3) at f = 1 kHz.
Sensors 24 03003 g013
Figure 14. The measured permittivity εmc with a dependence on the cylinder’s radius rc: PU foams (a) SikaBlock-M80 (85 kg/m3), (b) SikaBlock-M450 (415 kg/m3), monolithic materials, (c) lab-made PU (1280 kg/m3), (d) industrial PU SikaBlock-M945 (1352 kg/m3), (e) industrial epoxy Lab-975 New (708 kg/m3) at f = 1 kHz, and (f) industrial PU SikaBlock-M945 (1352 kg/m3) at f = 40960 Hz.
Figure 14. The measured permittivity εmc with a dependence on the cylinder’s radius rc: PU foams (a) SikaBlock-M80 (85 kg/m3), (b) SikaBlock-M450 (415 kg/m3), monolithic materials, (c) lab-made PU (1280 kg/m3), (d) industrial PU SikaBlock-M945 (1352 kg/m3), (e) industrial epoxy Lab-975 New (708 kg/m3) at f = 1 kHz, and (f) industrial PU SikaBlock-M945 (1352 kg/m3) at f = 40960 Hz.
Sensors 24 03003 g014
Figure 15. The relative partial susceptibility RX with a dependence on the cylinder’s radius rc: (a) SikaBlock-M80 (85 kg/m3), (b) PU foam SikaBlock-M450 (415 kg/m3), monolithic materials, (c) lab-made PU (1280 kg/m3), (d) industrial PU SikaBlock-M945 (1352 kg/m3), and (e) industrial epoxy Lab-975 New (708 kg/m3) at f = 1 kHz.
Figure 15. The relative partial susceptibility RX with a dependence on the cylinder’s radius rc: (a) SikaBlock-M80 (85 kg/m3), (b) PU foam SikaBlock-M450 (415 kg/m3), monolithic materials, (c) lab-made PU (1280 kg/m3), (d) industrial PU SikaBlock-M945 (1352 kg/m3), and (e) industrial epoxy Lab-975 New (708 kg/m3) at f = 1 kHz.
Sensors 24 03003 g015
Figure 16. Model functions Φ of the relative partial susceptibility RX with a dependence on the cylinder’s radius rc: (a) SikaBlock-M80 (85 kg/m3), (b) PU foam SikaBlock-M450 (415 kg/m3), monolithic materials, (c) lab-made PU (1280 kg/m3), (d) industrial PU SikaBlock-M945 (1352 kg/m3), and (e) an epoxy Lab-975 New (708 kg/m3) at f = 1 kHz.
Figure 16. Model functions Φ of the relative partial susceptibility RX with a dependence on the cylinder’s radius rc: (a) SikaBlock-M80 (85 kg/m3), (b) PU foam SikaBlock-M450 (415 kg/m3), monolithic materials, (c) lab-made PU (1280 kg/m3), (d) industrial PU SikaBlock-M945 (1352 kg/m3), and (e) an epoxy Lab-975 New (708 kg/m3) at f = 1 kHz.
Sensors 24 03003 g016
Figure 17. Model functions Φ of the relative partial susceptibility RX with a dependence on the coverage coefficient kc of the cylinders: (a) SikaBlock-M80 (85 kg/m3), (b) PU foam SikaBlock-M450 (415 kg/m3), monolithic materials, (c) lab-made PU (1280 kg/m3), (d) industrial PU SikaBlock-M945 (1352 kg/m3), and (e) an epoxy Lab-975 New (708 kg/m3) at f = 1 kHz.
Figure 17. Model functions Φ of the relative partial susceptibility RX with a dependence on the coverage coefficient kc of the cylinders: (a) SikaBlock-M80 (85 kg/m3), (b) PU foam SikaBlock-M450 (415 kg/m3), monolithic materials, (c) lab-made PU (1280 kg/m3), (d) industrial PU SikaBlock-M945 (1352 kg/m3), and (e) an epoxy Lab-975 New (708 kg/m3) at f = 1 kHz.
Sensors 24 03003 g017
Figure 18. The dependence of the coverage coefficient kc on the true permittivity εt at the relative partial susceptibility RX (a) 50%, (b) 75%, and (c) 90%.
Figure 18. The dependence of the coverage coefficient kc on the true permittivity εt at the relative partial susceptibility RX (a) 50%, (b) 75%, and (c) 90%.
Sensors 24 03003 g018
Figure 19. The change rate of the relative partial susceptibility, with a dependence on the coverage coefficient of the cylinders (continuous lines) and of shells (dashed lines): (a,f) SikaBlock-M80 (85 kg/m3), (b,g) PU foam SikaBlock-M450 (415 kg/m3), monolithic materials (c,h) lab-made PU (1280 kg/m3), (d,i) industrial PU SikaBlock-M945 (1352 kg/m3) as well as (e,j) an epoxy Lab-975 New (708 kg/m3) at f = 1 kHz.
Figure 19. The change rate of the relative partial susceptibility, with a dependence on the coverage coefficient of the cylinders (continuous lines) and of shells (dashed lines): (a,f) SikaBlock-M80 (85 kg/m3), (b,g) PU foam SikaBlock-M450 (415 kg/m3), monolithic materials (c,h) lab-made PU (1280 kg/m3), (d,i) industrial PU SikaBlock-M945 (1352 kg/m3) as well as (e,j) an epoxy Lab-975 New (708 kg/m3) at f = 1 kHz.
Sensors 24 03003 g019
Figure 20. The change rate of the relative partial susceptibility ΔRX/Δkc with a dependence on the radius rc of a cylinder: (a) SikaBlock-M80 (85 kg/m3), (b) PU foam SikaBlock-M450 (415 kg/m3), monolithic materials, (c) lab-made PU (1280 kg/m3), (d) industrial PU SikaBlock-M945 (1352 kg/m3), as well as (e) an epoxy Lab-975 New (708 kg/m3) at f = 1 kHz.
Figure 20. The change rate of the relative partial susceptibility ΔRX/Δkc with a dependence on the radius rc of a cylinder: (a) SikaBlock-M80 (85 kg/m3), (b) PU foam SikaBlock-M450 (415 kg/m3), monolithic materials, (c) lab-made PU (1280 kg/m3), (d) industrial PU SikaBlock-M945 (1352 kg/m3), as well as (e) an epoxy Lab-975 New (708 kg/m3) at f = 1 kHz.
Sensors 24 03003 g020
Table 1. Increments of cylinder’s radius Δrc at a constant increment of the coverage coefficient Δkc = 5%.
Table 1. Increments of cylinder’s radius Δrc at a constant increment of the coverage coefficient Δkc = 5%.
iCoverage
Coefficient kc
Radius rc; mmIncrement Δrc; mmiCoverage
Coefficient kc
Radius rc; mmIncrement Δrc; mm
10.000.05.03120.5516.70.74
20.055.02.08130.6017.40.71
30.107.11.60140.6518.10.68
40.158.71.35150.7018.80.66
50.2010.11.19160.7519.50.64
60.2511.31.07170.8020.10.62
70.3012.30.99180.8520.70.60
80.3513.30.92190.9021.30.58
90.4014.20.86200.9521.90.57
100.4515.10.82211.0022.50.00
110.5015.90.78
Table 2. The measured permittivity and measurement uncertainties of PU foams.
Table 2. The measured permittivity and measurement uncertainties of PU foams.
NDensity ρ;
kg/m3
Kind of a SampleMeasured Permittivity εmExpanded Uncertainty U95%
f1f7f16f1f7f16
150(1)1.091.091.080.0100.0100.010
(2)1.101.091.080.0100.0100.010
(3)1.091.091.080.0100.0100.010
2144(1)1.281.261.240.0130.0100.100
(2)1.271.261.250.0120.0100.100
(3)1.271.261.240.0120.0100.100
3427(1)1.861.841.750.0110.0110.011
(2)1.871.841.760.0150.0110.010
(3)1.861.841.760.0110.0110.010
Table 3. The permittivity of cylindrical and semi-cylindrical samples.
Table 3. The permittivity of cylindrical and semi-cylindrical samples.
NDielectric MaterialDensity ρ; kg/m3SampleThickness t; mmPermittivity εt and εmc
f1f8f16
1PU foam SikaBlock-M450437A cylinder121.8551.8021.756
2Two semi-cylinders1.8521.8041.757
|R|; %0.160.110.06
3Monolithic PU SikaBlock- M9601180A cylinder123.7483.6723.507
4Two semi-cylinders3.7643.6813.520
|R|; %0.430.250.37
5Monolithic
PU, lab-made
1280Two semi-cylinders123.4743.4083.310
Table 4. Theoretical parameters of the sample and subsamples (din and dout—the inner and outer diameters).
Table 4. Theoretical parameters of the sample and subsamples (din and dout—the inner and outer diameters).
Sample/SubsampleDiameterThickness of a Shell Wall t; mmCross-Sectional Area
S; mm2
Coverage
Coefficients
kc and ks; %
d;
mm
din;
mm
dout;
mm
Cylinder, complete45.0---1590.4100
Subsamples
Cylinder22.50.022.5-397.625
Shell 1-22.531.84.7397.625
Shell 2-31.839.03.6397.625
Shell 3-39.045.03.0397.625
Sum:1590.4100
The combined sample
Cylinder + shell 1 + shell 2 + shell 345.0-- 1590.4100
Table 5. The parameters of the sample and subsamples, determined experimentally.
Table 5. The parameters of the sample and subsamples, determined experimentally.
Sample/SubsampleDiameterThickness of a Shell Wall t; mmArea S;
mm2
kc and
ks; %
εt and εm(1 kHz)Susceptibility χt and χp
(1 kHz)
χp/Sumχp;
%
d;
mm
din;
mm
dout;
mm
Cylinder, complete44.9---1583.499.61.400.40
Subsamples
Cylinder22.2---387.124.31.340.3489.3
Shell 1-22.831.54.35371.023.31.030.038.1
Shell 2-32.138.73.3367.023.11.010.012.1
Shell 3-39.144.92.9382.624.11.000.000.5
Sum:1507.794.84.380.38100.0
The combined sample
Cylinder + shell 1 + shell 2 + shell 344.9---1507.794.81.400.40100.0
Table 6. Radiuses of the sensor’s electrodes and corresponding coverage coefficients for cylinders and shells.
Table 6. Radiuses of the sensor’s electrodes and corresponding coverage coefficients for cylinders and shells.
Radius
re1in; mm
Coverage Coefficients kc and ks; %Radius
re1out; mm
Coverage Coefficients kc and ks; %Radius
re2; mm
Coverage Coefficients
kc and ks; %
Radius
re3; mm
Coverage Coefficients kc and ks; %
10.5021.822.50100.08.4013.99.4517.6
78.20.086.182.4
Table 7. Parameters of the model functions Φ(rc) and Ψ(rsin) for cylinders and shells.
Table 7. Parameters of the model functions Φ(rc) and Ψ(rsin) for cylinders and shells.
NPU FoamsDensity ρ;
kg/m3
True Permitt. εt
(1 kHz)
Model Functions Φ(rc) and Ψ(rsin)
Normal Cumulative
f1c(rc) or f1s(rsin)
Lognormal Cumulative
f2c(rc) or f2s(rsin)
μ1σ1Logicalμ2σ2rT
Cylinders
1SikaBlock-M80851.158.82.9TRUE1.320.56−4.8
2SikaBlock-M1501441.249.02.8---“---1.290.59−5.1
3SikaBlock-M4504151.789.13.1---“---1.320.57−5.3
Shells
1SikaBlock-M80851.158.82.9TRUE1.320.56−4.8
2SikaBlock-M1501441.249.02.7---“---1.100.67−5.3
3SikaBlock-M4504151.789.13.1---“---1.320.57−5.3
Table 8. The parameters of model function Φ(rc) (cylinders).
Table 8. The parameters of model function Φ(rc) (cylinders).
NDielectric MaterialDensity; kg/m3Permitt.
εt
(1 kHz)
Model Function Φ(rc)
Normal Cumulative f1(rc)Lognormal Cumulative f2(rc)
μ1σ1Logicalμ2σ2rT
1Monolithic PU; lab-made12803.409.83.5TRUE1.300.58−6.0
2Monolithic PU Sika M94513524.3410.83.6---“---1.150.63−7.0
3Epoxy LAB9757088.9512.04.0---“---1.100.73−7.6
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Beverte, I. Electric Susceptibility at Partial Coverage of a Circular One-Side Access Capacitive Sensor with Rigid Polyurethane Foams. Sensors 2024, 24, 3003. https://doi.org/10.3390/s24103003

AMA Style

Beverte I. Electric Susceptibility at Partial Coverage of a Circular One-Side Access Capacitive Sensor with Rigid Polyurethane Foams. Sensors. 2024; 24(10):3003. https://doi.org/10.3390/s24103003

Chicago/Turabian Style

Beverte, Ilze. 2024. "Electric Susceptibility at Partial Coverage of a Circular One-Side Access Capacitive Sensor with Rigid Polyurethane Foams" Sensors 24, no. 10: 3003. https://doi.org/10.3390/s24103003

APA Style

Beverte, I. (2024). Electric Susceptibility at Partial Coverage of a Circular One-Side Access Capacitive Sensor with Rigid Polyurethane Foams. Sensors, 24(10), 3003. https://doi.org/10.3390/s24103003

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop