Homogenization of the vertically stacked medium frequency magnetic metamaterials with multi

Scientific Reports volume 12, Article number: 20333 (2022) Cite this article

721 Accesses

1 Citations

1 Altmetric

Metrics details

The paper presents a homogenization method of the magnetic metamaterials, made of perpendicularly oriented resonators consisting of multi-turn planar coils. A resulting composite, in the form of parallel stripes with metamaterial cells, exhibits extraordinary properties in the medium frequency magnetic field, such as zero permeability. To identify an effective permeability of this metamaterial, two models were presented, i.e., a three-dimensional numerical model with current sheet approximation as well as Lorentz oscillator model, where individual coefficients are based on the lumped circuit parameters and directly related with a geometry of the unit cell. The accuracy of the second approach is improved by taking into account mutual inductances in a metamaterial grid. Then, a comparison is made with numerical model results to show adequacy of the adopted analytical attempt, and properties of this type of metamaterial are discussed. It is shown that discussed metamaterial structure can achieve negative permeability as well as its values, at identical resonant frequency, are dependent on number of turns of the planar coil.

Metamaterials were introduced at the end of the twentieth century and gained great attention in the scientific and industrial environment. Due to an ability to exhibit unordinary properties, for instance a zero or negative permittivity and permeability1,2,3,4, a potential interest of creating new functional structures had risen. The synthesis of superlenses5,6, cloaks7,8,9, new types of antennas2,10 and absorbers11,12,13 and devices for medical diagnosis14,15 are some of many considered applications. Recent advances in terahertz metamaterials also show an accelerating progress of tunable metastructures16,17, where a resonance wavelength can be tuned using appropriate external fields. A dynamic development of metamaterials in a high frequency electromagnetic field also encouraged scientists to analyze stationary as well as low frequency fields, and resulted in Laplacian metamaterials18,19. The scope of usage had extended to extraordinary devices, such as invisibility shells for sensors7, artificial wormholes20, magnetic hoses21, lenses increasing the resolution of magnetic resonance imaging22,23,24 and effective shielding of a low frequency magnetic field25,26. Nevertheless, some of these were intended to use in the stationary magnetic field, e.g., magnetic hose, through the utilization of anisotropic structures: layers of ferromagnetic material and superconductor with nearly zero permeability19. Still, it may be possible to mimic these abilities in alternating magnetic field at a specific frequency, if metamaterials characterized by identical properties would be synthesized. To effectively design such composite materials, a fast and accurate homogenization method to estimate their effective properties would be strongly desired.

Despite the low and high frequency structures, the medium frequency magnetic field metamaterials were discussed as well. Their topology is similar to the high frequency metamaterials, for example the split-ring resonators27,28, hence unit cells consist of an inductive element (in this case a planar coil) made of a conductive material and connected with capacitance (parasitic, or in a form of the lumped capacitor), which are then arranged periodically to form the structure with effective properties. Aside from the classical horizontally stacked unit cells, some examples of vertically stacked resonators were shown1,11,29,30,31, to create a composite possessing identical outer dimensions but able to interact with other component of an external field. In most cases, applications in the relatively high frequencies were considered, e.g., in MHz range32,33,34,35. Including widely developed stationary field meta-structures a peculiar gap, covering frequencies higher than several kHz but less than MHz band, in the used frequency spectrum of magnetic metamaterials can be observed.

It is important, since promising application of the medium frequency metamaterials was recently found in the wireless power transfer (WPT) systems36,37,38. The metamaterial can act as the device enhancing a magnetic coupling between transmitting and receiving coil38,39,40, which led to an increase of an efficiency of the WPT system34,37,38. There are several different types of metamaterials used for this purpose, where many of them were based on the printed circuit board planar coils with on-chip capacitors32,33,34,36. They are also varying in terms of positioning with relation to transmitting and receiving coil, e.g., parallel or perpendicularly oriented in-between them, behind these coils or in front/behind one of them35,39,40. This shows that both resonators and their position in space has a crucial impact on WPT system performance, yet properties of the medium frequency magnetic field metamaterials are not well known, and multi-turn resonators were not fully characterized either.

The extraction of medium frequency metamaterials effective properties is the most important part of their analysis. The homogenization of the unit cell can be performed using several methods, such as field averaging41, Bloch wave42,43,44 or scattering parameters (S-parameters)13,45. These methods, however, are based on a spatial field distribution, which is obtained predominantly from three-dimensional (3D) numerical model solved in a frequency domain. Since the number of degrees of freedom (NDOF) quickly becomes large in 3D models, the demand for computing memory rises sharply, putting pressure on computer memory, and the calculation time is also prolonged. To avoid disadvantages of numerical approach the analytical methods were used28,46 to estimate effective properties. Unfortunately, due to a complex derivation of the final formulae they are limited to the simplest or the best-known geometries of meta-cells. Techniques combining general solution of the analytical formulation with numerically (or experimentally) identified coefficients were also introduced. For example, the Lorentz oscillator model was used in47 and the quality factor and resonant frequency were found based on a modified 3D numerical model. To the best of the author's knowledge, no purely analytical solution was presented for the medium frequency metamaterials made of planar coils, especially of a multi-turn type, with resonators stacked vertically in space.

In this article such technique of homogenization was presented. The main focus was put on numerical modeling and calculations, but mostly on determining formulae and procedure for analyzing discussed composite material by analytical approach. Accordingly, the vertically stacked stripes of resonators, forming metamaterial intended for the medium frequency bandwidth, were analyzed numerically (3D model) and by the Lorentz oscillator model. The novelty included in the second approach can be found in taking into account magnetic couplings in the periodic grid and derived equations, where all coefficients of the Lorentz model were identified and directly related with the geometry of the unit cell and a type of planar coil. This allows estimating permeability of metamaterial without numerical model or experimentally identified coefficients, reducing computational requirements. While the numerical approaches were not able to find properties of each individual unit cell in a finite grid, the presented analytical model opens this possibility. An analysis of the effective permeability was performed for different resonant frequencies, distances between layers and filling factor of the resonator. Both models were compared and the results as well as accuracy of the proposed analytical model were discussed.

The medium frequency metamaterial is intended to operate at frequencies higher than 10 kHz but less than 10 MHz. In this bandwidth it is possible to create a metamaterial unit cell consisting of the periodic grid of resonators. These resonators are constructed of planar coils, interacting with external magnetic field, and a capacitance to achieve a resonance at some specific resonant frequency f0. Typical approach of realizing this metamaterial is to form layers of a rectangular 2D periodic grid of horizontally stacked resonators. Yet, another possibility is to orient these unit cells perpendicularly and utilize only one column of the periodic grid. As the result, instead of 2D array, a pack of many stripes with linearly organized resonators is used to build the composite, as shown in Fig. 1a. An elementary unit cell, presented in Fig. 1b, has centrally located planar coil connected in series with a lumped capacitor. The cell has width dw, length dl and the distance between stripes is h. Hence, the thickness of the metamaterial is equal to the width of unit cell.

(a) Considered metamaterial composite consisting of an array of vertically stacked planar coils connected with lumped capacitors. Resonators are perpendicularly oriented with respect to magnetic field H of the incident plane wave S. (b) Unit cell of the metamaterial structure with a length dl, width dw and height h.

The resonator is made in a form of multi-turn coil with attached lumped capacitor, as shown in Fig. 2a. Although in Figs. 1 and 2 the square planar coils were presented, the shape of coil can be changed to, e.g., circular or octagonal and even deformed. For example, the width of cells and, simultaneously, the thickness of the metamaterial may be reduced utilizing planar coil with rectangular shape, where dw ≠ dl. The main idea is to wound the turns of any type of the planar coil using ultra-thin wire, where a diameter of conductor is dc, separation between wounds is ds and the number of all turns is n. When the outer size of the planar coil is ro, then the inner resection will be expressed as

(a) A view on the resonator wounded using thin wires with diameter dc and separation ds. (b) The 3D numerical model of the metamaterial (unit cell) with boundary conditions and homogenized planar coil.

An application of ultra-thin wire has two main advantages: firstly, many turns can be wounded in an area restricted by outer size ro which will result in a higher self-inductance and secondly, if the diameter of the conducting part is less than skin depth, than the skin effect can be omitted in a mathematical model. On the other hand, the lesser diameter will result in a smaller cross-section area of the wire and hence a greater equivalent resistance of the coil.

The numerical model of the presented metamaterial structure can be reduced to a single cell. This approach simplifies the modeling task and decreases NDOF with respect to the model of the entire composite material. The considered three-dimensional model, presented in Fig. 2b, is a cuboid consisting (in the central part) of the unit cell with homogenized multi-turn resonator as well as two parts of an air (ambience) and boundary conditions. To mimic the infinite periodic arrangement of vertically stacked stripes, the periodic boundary condition is assigned to all outer surfaces, except for bases of the cuboid where lumped ports are imposed. Port 1 is a source of an external magnetic field with the only vector component directed perpendicularly to the surface of the resonator, so in this case H = [0 Hy 0]. Port 2 does not excite any external field since it acts as a perfectly absorbing surface which imitates infinite space behind the metamaterial.

The distribution of the magnetic field in a frequency domain can be solved using magnetic vector potential A = [Ax Ay Az], defined as μrμ0H = ∇ × A, and formulation of the phenomena through the Helmholtz equation48

where: μr is the relative permeability; μ0 is the permeability of air (H/m); ω is the angular frequency (rad/s); σ is the electrical conductivity (S/m). The solution of the Eq. (2) in 3D space can be obtained using, e.g., the finite element method (FEM)4.

A crucial part of the model is a homogenized multi-turn coil. Without any homogenization the complex structure, shown in Fig. 2a, would have to be recreated in the model, along with the thin wires (dc) and even thinner separations (ds). From the modeling point of view, it is troublesome task leading to the greater NDOF, since each thin part must be reconstructed by extremely fine mesh. As the result, the model may not be solved, or it will require long computation time. A solution to this problem is a current sheet approximation49,50 or an equivalent foil approach51. In the software used for the computations of the introduced metamaterial it is known as the Homogenized multi-turn boundary condition52. This approach implements a homogenized model of a coil consisting of many conducting wires or printed paths, where the number of wires and their cross-section area are taken into account. Due to this homogenization, coils containing many wires can be modeled as the negligibly thin conducting sheets.

The magnetic field distribution will be determined based on the magnetic vector potential. The distribution of magnetic field on the surfaces of Port 1 and Port 2 is the most important since magnetic field values are used to find the scattering parameters S11 and S2145. These parameters and equations allow identifying the effective permeability of the metamaterial13.

An alternative for numerical modeling of metamaterials is an analytical approach, based on the general formulation of the phenomena occurring in the meta-cell. In various works2,28,46,47,53 the formula for complex effective permeability μeff = μRe + jμIm was studied

where: f0 is the resonant frequency (Hz); F is the fractional volume factor; Q is the quality factor. These three main parameters of the function expressed by Eq. (3) must be found to describe the properties of the unit cell. The fractional volume factor is47

where: Leff is the effective inductance of the unit cell (H); Si is the area of i-th loop of the coil (m2); V is the volume of the cell (m3), which is

The analyzed unit cell is in a specific position within the metamaterial structure and surrounded by the other cells, as shown in Fig. 3. The magnetic coupling between this unit cell and each resonator appears and affects the effective magnetic properties. Similar situation was previously considered in the periodic arrays of wireless power transfer systems54 and it was found that inductive magnetic couplings are directly affecting the inductance of a coil by adding to its self-inductance all mutual inductances between considered inductor and all the others. Hence, at this point it will be assumed that effective inductance of the unit cell can be expressed as

where: Lself is the self-inductance of the planar coil (H); M is the sum of mutual inductances (H); a, b is the row and column number in the periodic grid, a ≠ 0 ˄ b ≠ 0; Ma,b is the mutual inductance between an arbitrary unit cell and the cell in a-th row and b-th column (H). The mutual inductance Ma,b is55

where: g = (dc + ds)/(2π); Φi = ri/g, Φo = ro/g; φ1, φ2 are the elementary angles (rad); ρ is the shape factor related with the type of coil, e.g., for circular ρ = 1 while for a square ρ varies from 1.1 to 1.2756,57. After applying the rectangle rule Eq. (7) takes the following form54

where: ΦK = (Φo – Φi)/K is an integration step; K is assumed number of integration subintervals, K ≥ ro/g and K ∈ N.

Front view on the periodic structure of resonators with indicated mutual inductances between considered unit cell and other cells in the array.

The self-inductance of a coil can be estimated using formula58

where davg is a mean diameter, e.g., for a circular or square coil it will be

and v is a fill factor

while c1, c2, c3, c4 are shape coefficients which differs for a square, circular, hexagonal, octagonal or any other shape of a planar coil58. The fractional volume in Eq. (4) also requires the sum of areas of each loop of the coil. The area of i-th loop is

Assuming that square coils are used to build a metamaterial, the sum of all areas from Eq. (12) may be calculated as

Finally, using Eqs. (5), (6) and (13) the fractional factor can be expressed in a form

From Eqs. (1), (8)–(11) it can be seen that for any type of a planar coil the coefficients ρ, c1–c4 can be taken from the literature and then, for a given outer size (ro), number of turns (n) and used conductor (with diameter dc and separation ds), the effective inductance from Eq. (6) can be found, and then substituted into Eq. (14) to calculate the fractional factor.

The second parameter is a quality factor (Q). A resonator in the unit cell acts as a RLC circuit, hence the quality factor is

where R is a resistance of the coil defined as

where: ltotal is the total wire length (m); σc is the electrical conductivity of a wire (S/m), ac is the wire cross-section area (m2). The length of the wire can be simply obtained as a product of the number of turns and an average perimeter (half length of the outer edges and edges of the inner resection). Hence, for a square shape planar spiral, wounded by a wire with a circular cross-section, the resistance can be obtained using formula

For a printed spiral coil with rectangular cross-section and copper path thickness t, simply ac = dct must be substituted in Eq. (16). Finally, the resonant frequency (f0) of the resonator will be found based on a series RLC circuit resonance condition

If in Eq. (18) the resonant frequency will be imposed, then a capacitance (C) can be derived. In further investigations this second approach was applied to show that if mutual inductances M would be omitted in Eq. (6), then the resonant frequency will not be properly estimated in Lorentz model, presented in Eq. (3). Nevertheless, the formulae shown above indicates that all parameters in this model are directly related with the unit cell geometry, such as distance between layers (h), pitch size (dl), size of the planar coil (ro) and inner resection (ri). As a result, the effective permeability can be estimated without numerical model and much less computational effort (i.e., computer memory and calculation time) must be used. All equations are also valid for coils with any number of turns; both printed and wounded using a wire. Additionally, the extremes of a real part of the Eq. (3) can be now easily found as

which are determined by the fill and quality factor. Both these factors are dependent on effective inductance and geometry of a coil as well as the quality factor, as shown in Eq. (15), makes them also dependent on a resonant frequency.

An exemplary metamaterial was studied to show a validity of the analytical Lorentz oscillator model. The metamaterial with vertically stacked stripes of meta-cells, shown in the previous chapter, was analyzed from a point of view of a single cell located in the center of the infinitely expanding periodic grid of resonators. The square planar coil with outer size ro = 20 mm was used. Two cases of the distance between layers were examined: h = dw = 44 mm and h = dw/2 = 22 mm. Electrical conductivity of a wire (copper) was σc = 5.7·107 S/m. In analytical models, the coefficients used for estimating a self-inductance were c1 = 1.26, c2 = 2.08, c3 = 0.14, c4 = 0.115. For a mutual inductance a shape coefficient was ρ = 1.27 and the number of rows and columns was limited to a, b ∈ < –10; 10 > .

The frequency bandwidth was limited to f ∈ < 100 kHz; 1 MHz > . Within this limit five desired resonant frequencies were imposed fz = {200, 375, 550, 725, 900} kHz. To achieve a resonance in meta-cell a proper capacitor has to be used. Based on Eq. (18) it can be concluded that effective inductance has an impact on a resonant frequency. To show the shift of this frequency, resulting from various Leff in different approaches (numerical and analytical), the same capacitor should be used. The simplest approach to find its capacitance is to assume identical "basic" inductance for each model, e.g., one of the possibilities is a self-inductance, effortless to find using Eq. (9). Therefore, the lumped capacitance adopted for the models was estimated as

The resulting geometrical parameters and capacitances connected to planar resonators in the unit cells are presented in Table 1. The numerical model of a meta-cell, presented in Fig. 2b, was solved by FEM. COMSOL Multiphysics 4.3b software was used, where in AC/DC module (Magnetic Fields physics) the frequency domain calculations were performed. The model was complemented by an additional electrical circuit, where an ideal lumped capacitor was connected in series to a homogenized multi-turn inductor.

The main parameter varying during tests was the number of turns (n) to study the characteristics and properties of metamaterial. More wounds help to improve a self-inductance (Lself), but also lead to a greater resistance (R), which limits possible values of the real part (μRe) and increases imaginary part (μIm). To make the analysis more general, the influence of n on the effective permeability (μeff) was presented in function of fill factor (v) from Eq. (11), dependent directly on n, and thereby representing the extent to which the available area of a planar coil is occupied by the turns. Abovementioned cases were calculated using three models:

3D numerical model with homogenized planar coil (FEM),

proposed Lorentz oscillator model with mutual inductances M (LMM),

Lorentz model without M (LM).

The frequency characteristics for cells at a distance h = dw were shown in Fig. 4. In all cases, a typical shape of μRe and μIm function was obtained, leading to a dynamic change from positive to negative values of the real part and a maximum value of μIm at resonant frequency. In Fig. 4a numerical model (FEM, red lines) has indicated that lower actual resonant frequency (f0 = 195,7 kHz) was obtained than the one imposed (fz = 200 kHz), if the capacitance from Eq. (21) was used in a resonator. Nearly identical characteristics were achieved using Lorentz model with a sum of mutual inductances (LMM, blue lines) where actual resonant frequency was f0 = 196,2 kHz, and the extreme values of the real and imaginary part are consistent with the one from numerical model. However, the analytical model which was not considering magnetic couplings (LM, green lines) preserved unchanged resonant frequency (f0 = fz). In case of LM the effective inductance is equal to self-inductance, hence frequencies from Eqs. (18) and (21) are the same.

Real and imaginary parts of the effective permeability obtained using three models at h = dw, for the specified resonant frequencies: (a) fz = 200 kHz. (b) fz = 900 kHz.

In Fig. 4b it can observed that between FEM and LMM models differences became greater than in Fig. 4a (f0FEM = 887 kHz and f0LMM = 883 kHz), yet still smaller than comparing FEM to LM (f0FEM = 887 kHz and f0LM = 900 kHz). In this specific instance the resonant frequency shift from LMM was more significant than the one from FEM (f0LMM < f0FEM), even though numerical model considers all mutual inductances in the grid, hence FEM computed Leff should be the highest. As a result, according to Eq. (18), f0FEM always ought to be the lowest. The discrepancy in Fig. 4b is the example of potential issue which may come from the usage of Eq. (7). Numerical solution of this formula produces additional errors, finite number of rows and columns is considered during the calculations and ρ may be different (it strongly depends on coil's structure and varies from 1.1 to 4/π). Thus, the calculated M is an approximation and, in some cases, can produce higher Leff in LMM than in FEM. Such situation is not correct and to avoid it ρ has to be chosen carefully.

After reducing distance between stripes to h = dw/2 the divergences in favor of LMM approach were more significant, which can be seen comparing Figs. 4a and 5a as well as Figs. 4b and 5d. In Fig. 5 it was shown that resonant frequency and values of effective permeability from LMM were closer to FEM results than LM. Differences between the maximum values of real and imaginary parts were also visible, especially comparing FEM with LM (Fig. 5c). For example, in Fig. 5b max(μIm) from LM was 43.45, while from FEM it was 37.43 and from LMM it was 38.83. This example shows a dependency of both real and imaginary parts on mutual inductances appearing in the metamaterial made of vertically stacked unit cells. Hence, it follows that magnetic couplings in the periodic grid must be taken into account during calculations of the capacitance, to more accurately predict both resonant frequency and values of effective permeability at any frequency.

Real and imaginary parts of the effective permeability obtained using three models at h = dw/2, for specified resonant frequencies: (a) fz = 200 kHz. (b) fz = 375 kHz. (c) fz = 725 kHz. (d) fz = 900 kHz.

In Fig. 6 extremes of the real part were shown and LMM approach was again compared with FEM results. One can conclude that for a fill factor (v) increasing from 0 to 0.45, which was caused by the larger number of turns (n), maximum and minimum values of the effective permeability were increasing. For the higher fill factors a decrease of these values occurred, which was mainly observed for the highest resonant frequency (900 kHz). Metamaterial exhibited paramagnetic and ferromagnetic behavior since max(μRe) exceeds 26 at h = dw, as shown in Fig. 6a, and 46 at h = dw/2 (Fig. 6c). Moreover, the negative permeability was obtained and reached even −25 at h = dw, as shown in Fig. 6b, and −45 at h = dw/2 (Fig. 6d). The most interesting observation is that optimal values of max(μRe) or min(μRe), at specific resonant frequency, were found for v between 0.35 and 0.55. Further addition of turns had negative impact, since extreme values started to decrease and—from an economical point of view—more material was used, enlarging potential cost and mass of the metamaterial. It would be reasonable to prior perform theoretical investigation of potential impact of filling factor on the effective permeability and identify optimum n to avoid an unnecessarily large number of turns.

Effective permeability for specified resonant frequencies and two models: (a) Maximum value of the real part at h = dw. (b) Minimum value of the real part at h = dw. (c) Maximum value of the real part at h = dw/2. (d) Minimum value of the real part at h = dw/2.

The Eqs. (19) and (20) indicated that uneven absolute values of max and min μRe should be observed and both models have shown that asymmetry—the absolute max(μRe) is higher comparing to min(μRe). Moreover, these extremes were strongly dependent on the imposed resonant frequency, rather than on filling factor, and this is also in line with theoretical predictions. For identical v the optimum min(μRe) varied from approximately −5 to −25, as shown in Fig. 6b. In other words, the same planar coil will help to achieve not only various resonant frequencies but also effective permeabilities, due to an application of different capacitors.

In the finite grid of resonators each unit cell possesses distinct set of inductive couplings and, resulting from it, sum of mutual inductances, thus resonant frequency and fractional factor will differ within a grid, and therefore each resonator will exhibit various effective permeability at any arbitrary frequency. This explains differences in effective permeabilities, expected and observed, in terms of their distributions on a surface of the physically realized metamaterials. Approach presented in LMM enables to calculate these effective permeabilities for every unit cell independently, by an appropriate application of Eq. (8), i.e., selection of appropriate numbers of rows and columns. Furthermore, it can be also concluded from Fig. 5 that for f > > f0 imaginary part tends to 0, however real part tends to some constant value greater than 0 but less than 1, since based on Eq. (3) we have \(\mathop {\lim }\limits_{f \to \infty } \mu_{eff} = 1 - F\). This attribute of metamaterial suggests that it is possible to realize the composite structure with diamagnetic properties, where imaginary parts is negligible and real part of permeability will be approximately constant in some frequency range.

To compare quantitatively numerical (FEM) and analytical model (LMM and LM) a coefficient was proposed, e.g., for LMM given as

The coefficient from Eq. (22) expresses the sum of relative differences among extreme values of effective permeability for all five specified resonant frequencies fz. Additionally, to compare errors in estimating the actual resonant frequency similar relative difference was calculated

Cumulative error for LM approach (Δμ ≈ 14 ÷ 17%) was more than four times larger than the one for LMM (Δμ ≈ 2 ÷ 4%) at h = dw/2, as it can be seen from Fig. 7a. Relative difference was generally increasing with higher values of v. This suggests that the usage of LMM and planar coils with optimum filling factor should bring results to be more accurate. Nevertheless, the differences for LMM and LM for the case with h = dw were smaller (for v = 0.48 LMM resulted in nearly 0% error) which can be explained by the decreasing values of M with increasing distances between resonators in the periodic grid. In Fig. 7b similar findings can be listed, although LMM and LM for v < 0.29 gave comparatively low relative differences in estimating actual resonant frequency. For v > 0.29 the differences for LM were increasing but for LMM started to decrease. It was the only discrepancy comparing to Fig. 7a and other results in Fig. 7b.

(a) Relative difference in extreme values of the effective permeability between analytical and numerical model. (b) Relative difference in values of the resonant frequencies between analytical and numerical model.

In this article the medium frequency metamaterials with vertically stacked multi-turn resonators were presented and analyzed in terms of their computational models and effective properties. Numerical model with homogenized planar resonator and the analytical model, based on Lorentz oscillator, were characterized. The Lorentz model was formulated using electrical lumped parameters of the resonator and a periodic grid, where the mutual inductances between planar coils were taken into account to improve an accuracy of the analytical approach. As a result, the effective permeability may be estimated for any resonators with different number of turns, wires, and sizes of the unit cell, since all lumped parameters have been associated with geometric parameters.

The introduced Lorentz model with mutual couplings eliminates the need of performing numerical computations with the use of 3D model solved in a frequency domain. Hence, the pressure put on memory of a computational unit and the overall solution time can be dramatically reduced. Proposed analytical approach also enables to find the effective permeability of each unit cell in a finite grid, and thus to take into account edge effects appearing in real metamaterials.

An exemplary metamaterial with square planar coils was solved both numerically and analytically. The results showed that by including magnetic couplings of individual resonators in Lorentz model, the differences between numerical and analytical model had been reduced from several to only few percent (or nearly to zero in some cases). While omitting mutual inductances in Lorentz model had led to significant differences in the actual resonant frequency value. Additionally, further analysis had indicated that, in considered metamaterial, it was possible to achieve negative permeability reaching even min(μRe) = −45, but also the dependency of effective permeability on planar coil filling factor was non-linear, with an optimum value for filling factor between 0.35 and 0.55. The permeability of metamaterial was higher with higher values of the desired resonant frequency; hence it was strongly depending on the capacity of a series connected capacitor. Further research will focus on experimental verification of the abovementioned results, as well as modification of the unit cell geometry, where irregular planar coils will be used.

The data generated and analysed during the current study is available from the corresponding author.

Zhang, J. et al. Experimental characterization and cell interactions of a two-dimensional isotropic left-handed metamaterial. Appl. Phys. Lett. 92, 084108. https://doi.org/10.1063/1.2884528 (2008).

Article ADS CAS Google Scholar

Banerjee, B. An Introduction to Metamaterials and Waves in Composites 4th edn. (CRC Press Taylor & Francis Group, 2011).

Book Google Scholar

Moore, R. Electromagnetic Composites Handbook 5th edn. (McGraw-Hill Education, 2016).

Google Scholar

Zhong, S.-L., Dang, Z.-M. & Zha, J.-W. Prediction on effective permittivity of 0–3 connectivity particle/polymer composites at low concentration with finite element method. IEEE Trans. Dielectric Electric Insul. 25(6), 2122–2128. https://doi.org/10.1109/TDEI.2018.007193 (2018).

Article CAS Google Scholar

Pendry, J. B. Negative refraction makes a perfect lens. Phys. Rev. Lett. 85(18), 3966–3969. https://doi.org/10.1103/PhysRevLett.85.3966 (2000).

Article ADS CAS PubMed Google Scholar

Hegde, R., Hor, Y.L., Szabo, Z., Li, E.P., & Hoefer, W.J.R. Veselago-Pendry superlens imaging modeled with a spectral waveguide approach. in 2011 XXXth URSI General Assembly and Scientific Symposium. 1–4. https://doi.org/10.1109/URSIGASS.2011.6050622 (2011).

Navau, C. et al. Enhancing the sensitivity of magnetic sensors by 3D metamaterial shells. Sci. Rep. 7, 44762. https://doi.org/10.1038/srep44762 (2017).

Article ADS CAS PubMed PubMed Central Google Scholar

Steckiewicz, A. High-frequency cylindrical magnetic cloaks with thin layer structure. J. Magn. Magn. Mater. 534, 168039. https://doi.org/10.1016/j.jmmm.2021.168039 (2021).

Article CAS Google Scholar

Schurig, D. et al. Metamaterial electromagnetic cloak at microwave frequencies. Science 314, 977. https://doi.org/10.1126/science.1133628 (2006).

Article ADS MathSciNet CAS PubMed Google Scholar

Shaw, T., Samanta, G., Mitra, D., Mandal, B. & Augustine, R. Design of metamaterial based efficient wireless power transfer system utilizing antenna topology for wearable devices. Sensors 21, 3448. https://doi.org/10.3390/s21103448 (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Sichao, Q., Yuxiao, H. & Sheng, P. Conceptual-based design of an ultrabroadband microwave metamaterial absorber. Proc. Natl. Acad. Sci. USA 118(3), e2110490118. https://doi.org/10.1073/pnas.211049011 (2021).

Article Google Scholar

Xu, G. et al. 0.1–20 THz ultra-broadband perfect absorber via a flat multi-layer structure. Opt. Exp. 24(20), 23177–23185. https://doi.org/10.1364/OE.24.023177 (2016).

Article ADS CAS Google Scholar

Hannan, S. et al. A filling-factor engineered, perfect metamaterial absorber for multiple applications at frequencies set by IEEE in C and X bands. J. Mater. Res. Technol. 19, 934–946. https://doi.org/10.1016/j.jmrt.2022.05.071 (2022).

Article ADS CAS Google Scholar

Silalahi, H. M. et al. Floating terahertz metamaterials with extremely large refractive index sensitivities. Photon. Res. 9(10), 1970–1978. https://doi.org/10.1364/PRJ.433335 (2021).

Article Google Scholar

Wang, Z. et al. Plasmonic metasurfaces for medical diagnosis applications: A review. Sensors 22(1), 133. https://doi.org/10.3390/s22010133 (2022).

Article ADS CAS Google Scholar

Silalahi, H. M., Tsai, Y.-H., Lee, K.-L., Wei, P.-K. & Huang, C.-Y. Large shift of resonance wavelengths of silver nanoslit arrays using electrowetting-on-dielectric cells. Opt. Lett. 46(4), 705–708. https://doi.org/10.1364/OL.415500 (2021).

Article ADS PubMed Google Scholar

Shih, Y.-H., Lin, X.-Y., Silalahi, H. M., Lee, C.-R. & Huang, C.-Y. Optically tunable terahertz metasurfaces using liquid crystal cells coated with photoalignment layers. Curr. Comput.-Aided Drug Des. 11(9), 1100. https://doi.org/10.3390/cryst11091100 (2021).

Article CAS Google Scholar

Han, T. & Qiu, C.-W. Transformation Laplacian metamaterials: Recent advances in manipulating thermal and dc fields. J. Opt. 18(4), 1–13 (2016).

Article MathSciNet CAS Google Scholar

Sanchez, A., Navau, C., Prat-Camps, J. & Chen, D.-X. Antimagnets: Controlling magnetic fields with superconductor–metamaterial hybrids. New J. Phys. 13, 093034. https://doi.org/10.1088/1367-2630/13/9/093034 (2011).

Article ADS CAS Google Scholar

Prat-Camps, J., Navau, C. & Sanchez, A. A magnetic wormhole. Sci. Rep. 5, 12488. https://doi.org/10.1038/srep12488 (2015).

Article ADS CAS PubMed PubMed Central Google Scholar

Navau, C., Prat-Camps, J., Romero-Isart, O., Cirac, J. I. & Sanchez, A. Long-distance transfer and routing of static magnetic fields. Phys. Rev. Lett. 112, 253901 (2014).

Article ADS CAS PubMed Google Scholar

Duan, G. et al. Boosting magnetic resonance imaging signal-to-noise ratio using magnetic metamaterials. Commun. Phys. 2, 35. https://doi.org/10.1038/s42005-019-0135-7 (2019).

Article PubMed PubMed Central Google Scholar

Wu, K., Zhao, X., Bifano, T. G., Anderson, S. W. & Zhang, X. Auxetics-inspired tunable metamaterials for magnetic resonance imaging. Adv. Mater. 34(6), 2109032. https://doi.org/10.1002/adma.202109032 (2021).

Article CAS Google Scholar

Hurshkainen, A. et al. A novel metamaterial-inspired RF-coil for preclinical dual-nuclei MRI. Sci. Rep. 8, 9190. https://doi.org/10.1038/s41598-018-27327-y (2018).

Article ADS CAS PubMed PubMed Central Google Scholar

Boyvat, M. & Hafner, C. Magnetic field shielding by metamaterials. PIER 136, 647–664. https://doi.org/10.2528/PIER12121805 (2013).

Article Google Scholar

Wei, B. et al. Matrix metamaterial shielding design for wireless power transfer to control the magnetic field. Materials 15, 2678. https://doi.org/10.3390/ma15072678 (2022).

Article ADS CAS PubMed PubMed Central Google Scholar

Vallecchi, A. et al. Superdirective dimers of coupled self-resonant split ring resonators: Analytical modelling and numerical and experimental validation. Sci. Rep. 10, 274. https://doi.org/10.1038/s41598-019-56988-6 (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Chen, H. et al. A meta-substrate to enhance the bandwidth of metamaterials. Sci. Rep. 4, 5264. https://doi.org/10.1038/srep05264 (2014).

Article CAS PubMed PubMed Central Google Scholar

Maurya, S., Nyman, M., Kaivola, M. & Shevchenko, A. Highly birefringent metamaterial structure as a tunable partial polarizer. Opt. Exp. 27, 27335–27344. https://doi.org/10.1364/OE.27.027335 (2019).

Article CAS Google Scholar

Burokur, S. N. & de Lustrac, A. Negative index from asymmetric metallic cut wire pairs metamaterials. Int. J. Microw. Wirel. Technol. 1(6), 52–527. https://doi.org/10.1017/S1759078709990808 (2010).

Article Google Scholar

Gong, Z. & Yang, S. One-dimensional stacking miniaturized low-frequency metamaterial bulk for near-field applications. J. Appl. Phys. 127, 114901. https://doi.org/10.1063/1.5135974 (2020).

Article ADS CAS Google Scholar

Zhang, Y., Tang, H., Yao, C., Li, Y. & Xiao, S. Experiments on adjustable magnetic metamaterials applied in megahertz wireless power transmission. AIP Adv. 5, 017142. https://doi.org/10.1063/1.4907043 (2015).

Article ADS CAS Google Scholar

Rong, C. et al. A critical review of metamaterial in wireless power transfer system. IET Power Electron. 14, 1541–1559. https://doi.org/10.1049/pel2.12099 (2021).

Article Google Scholar

Shan, D., Wang, H., Cao, K. & Zhang, J. Wireless power transfer system with enhanced efficiency by using frequency reconfigurable metamaterial. Sci. Rep. 12, 331. https://doi.org/10.1038/s41598-021-03570-8 (2022).

Article ADS CAS PubMed PubMed Central Google Scholar

Lee, W. & Yoon, Y.-K. Wireless power transfer systems using metamaterials: A review. IEEE Access 8, 147930–147947. https://doi.org/10.1109/ACCESS.2020.3015176 (2020).

Article Google Scholar

Zhou, J., Zhang, P., Han, J., Li, L. & Huang, Y. Metamaterials and metasurfaces for wireless power transfer and energy harvesting. Proc. IEEE 110(1), 31–55. https://doi.org/10.1109/JPROC.2021.3127493 (2022).

Article Google Scholar

Wang, B., Yerazunis, W. & Teo, K. H. Wireless power transfer: Metamaterials and array of coupled resonators. Proc. IEEE 101(6), 1359–1368. https://doi.org/10.1109/JPROC.2013.2245611 (2013).

Article Google Scholar

Li, W., Wang, P., Yao, C., Zhang, Y. & Tang, H. Experimental investigation of 1D, 2D, and 3D metamaterials for efficiency enhancement in a 6.78 MHz wireless power transfer system. in 2016 IEEE Wireless Power Transfer Conference (WPTC). 1–4. https://doi.org/10.1109/WPT.2016.7498809 (2016).

Huang, X. et al. Development and prospects of metamaterial in wireless power transfer. IET Power Electron. 14, 2423–2440. https://doi.org/10.1049/pel2.12189 (2021).

Article Google Scholar

Zhang, Z., Zhang, B., Deng, B., Wei, X. & Wang, J. Opportunities and challenges of metamaterial-based wireless power transfer for electric vehicles. Wirel. Power Transfer 5(1), 9–19. https://doi.org/10.1017/wpt.2017.12 (2018).

Article Google Scholar

Smith, D. R. & Pendry, J. B. Homogenization of metamaterials by field averaging. J. Opt. Soc. Am. B. 23(3), 391–403. https://doi.org/10.1364/JOSAB.23.000391 (2006).

Article ADS CAS Google Scholar

Sista, S. G. & Tewary, V. Bloch wave approach to almost periodic homogenization and approximations of effective coefficients. Discrete Contin. Dyn. Syst. B 27(4), 1989–2024. https://doi.org/10.3934/dcdsb.2021119 (2022).

Article MathSciNet MATH Google Scholar

Rockstuhl, C., Menzel, C., Paul, T. & Lederer, F. Homogenization of metamaterials from a Bloch mode perspective. in 2012 International Conference on Electromagnetics in Advanced Applications. 178–181. https://doi.org/10.1109/ICEAA.2012.6328614 (2012).

Hossain, A. N. M. S. & Tsukerman, I. Non-asymptotic homogenization of 3-D periodic structures. Phys. Lett. A 398, 127278. https://doi.org/10.1016/j.physleta.2021.127278 (2021).

Article MathSciNet CAS Google Scholar

Arslanagić, S. et al. A review of the scattering-parameter extraction method with clarification of ambiguity issues in relation to metamaterial homogenization. IEEE Antennas Propag. Mag. 55(2), 91–106. https://doi.org/10.1109/MAP.2013.6529320 (2013).

Article ADS Google Scholar

Pendry, J. B., Holden, A., Robbins, D. J. & Stewart, W. J. Magnetism from conductors, and enhanced non-linear phenomena. IEEE Trans. Microw. Theory Tech. 47(11), 2075–2084. https://doi.org/10.1109/22.798002 (1999).

Article ADS Google Scholar

Liu, J., Gong, Z., Yang, S., Sun, H. & Zhou, J. Practical model for metamaterials in wireless power transfer systems. Appl. Sci. 10, 8506. https://doi.org/10.3390/app10238506 (2020).

Article CAS Google Scholar

Stankiewicz, J. M. & Choroszucho, A. Comparison of the efficiency and load power in periodic wireless power transfer systems with circular and square planar coils. Energies 14, 4975. https://doi.org/10.3390/en14164975 (2021).

Article CAS Google Scholar

Gyselinck, J. & Dular, P. Frequency-domain homogenization of bundles of wires in 2-D magnetodynamic FE calculations. IEEE Trans. Magn. 41(5), 1416–1419. https://doi.org/10.1109/TMAG.2005.844534 (2005).

Article ADS Google Scholar

Meeker, D. C. An improved continuum skin and proximity effect model for hexagonally packed wires. J. Comput. Appl. Math. 236(18), 4635–4644. https://doi.org/10.1016/j.cam.2012.04.009 (2012).

Article MATH Google Scholar

Meeker, D.C. Continuum Representation of Wound Coils Via an Equivalent Foil Approach. Finite element method magnetics website. http://www.femm.info/examples/prox/notes.pdf (2010).

COMSOL Inc. AC/DC Module User's Guide: Chapter 3. COMSOL multiphysics official website. https://doc.comsol.com/5.4/doc/com.comsol.help.acdc/ACDCModuleUsersGuide.pdf (2018).

Cummer, S. A., Popa, B.-I. & Hand, T. H. Q-based design equations and loss limits for resonant metamaterials and experimental validation. IEEE Trans. Antennas Propag. 56(1), 127–132. https://doi.org/10.1109/TAP.2007.912959 (2008).

Article ADS Google Scholar

Steckiewicz, A., Stankiewicz, J. M. & Choroszucho, A. Numerical and circuit modeling of the low-power periodic WPT systems. Energies 13(10), 2651. https://doi.org/10.3390/en13102651 (2020).

Article Google Scholar

Liu, S., Su, J. & Lai, J. Accurate expressions of mutual inductance and their calculation of Archimedean spiral coils. Energies 12(10), 1–14. https://doi.org/10.3390/en12102017 (2019).

Article Google Scholar

Jow, U.-M. & Ghovanloo, M. Design and optimization of printed spiral coils for efficient transcutaneous inductive power transmission. IEEE Trans. Biomed. Circuits Syst. 1(3), 193–202. https://doi.org/10.1109/TBCAS.2007.913130 (2007).

Article Google Scholar

Raju, S., Wu, R., Chan, M. & Yue, C. P. Modeling of mutual coupling between planar inductors in wireless power applications. IEEE Trans. Power Electron. 29(1), 481–490. https://doi.org/10.1109/TPEL.2013.2253334 (2014).

Article ADS Google Scholar

Mohan, S. S., del Mar Hershenson, M., Boyd, S. P. & Lee, T. H. Simple accurate expressions for planar spiral inductances. IEEE J. Solid-State Circuits 34(10), 1419–1424. https://doi.org/10.1109/4.792620 (1999).

Article ADS Google Scholar

Download references

This work was supported by the Ministry of Science and Higher Education in Poland at the Białystok University of Technology under research subsidy No. WZ/WE-IA/2/2020.

Bialystok University of Technology, Faculty of Electrical Engineering, 15-351, Białystok, Poland

Adam Steckiewicz

You can also search for this author in PubMed Google Scholar

The work reported in this paper is formulated by the corresponding and single author.

Correspondence to Adam Steckiewicz.

The author declares no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Steckiewicz, A. Homogenization of the vertically stacked medium frequency magnetic metamaterials with multi-turn resonators. Sci Rep 12, 20333 (2022). https://doi.org/10.1038/s41598-022-24809-y

Download citation

Received: 04 August 2022

Accepted: 21 November 2022

Published: 25 November 2022

DOI: https://doi.org/10.1038/s41598-022-24809-y

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.