Fractal Interpolation Function based Thin Wire Antennas

Gizem Kalender1, Emine Yeşim Zoral2, and Serkan Günel2

1The Graduate School of Natural and Applied Sciences
2Engineering Faculty, Dept. of Electrical & Electronics Engineering
Dokuz Eylül University, Izmir, 35800, Turkey
gizem.kalender@deu.edu.tr, yesim.zoral@deu.edu.tr, serkan.gunel@deu.edu.tr

Submitted On: October 17, 2023; Accepted On: January 9, 2024

ABSTRACT

This paper presents an approach for the design of wire antennas based on fractal interpolation functions (FIFs). The interpolation points and the contraction factors of the FIFs are chosen as free parameters to modify the antenna geometry. The proposed structures’ gain and radiation pattern can be optimized using FIF parameters. Producible prefractal antennas obtained in the intermediate iterations of fractal generation have compact sizes compared to classical counterparts. The error in prefractal geometry and the original fractal is bounded, and can be determined in terms of the finest producible detail’s dimensions. The emerging structures have multiband behavior due to their self-similar and symmetric nature. To illustrate the approach, we have provided finite element based simulations for several prefractal antennas. |S11|, the gain, the radiation efficiency, the radiation patterns, and feed point impedances for the demonstrated antennas are calculated numerically. The results indicate that produced antennas can be used in applications that require limited mechanical size, multiple operating bands, and controlled radiation patterns.

Keywords: Fractal antenna, fractal interpolation functions, iterated function systems.

I. INTRODUCTION

Recent developments in wireless communications systems require more compact, wider bandwidth, multiband, and low-cost antennas. Fractal antennas can fulfill these requirements due to scale invariance, self-similarity, and space-filling properties of the fractals [1]. These properties enable the miniaturization of antenna structures [2, 3, 4, 5]. The fractal structures can be designed to increase the effective physical length of the antennas to achieve multiband behavior in a limited space [6, 7, 8]. Basic fractal geometries such as Koch and Hilbert curves, and Sierpinski carpet have been studied for their radiation characteristics in the literature widely [9, 10, 11, 12, 13, 14, 15]. Comprehensive and up-to-date reviews can be found in [16, 17].

However, the studies rarely relate the mathematical properties of fractals to the antenna radiation characteristics. One approach is to optimize the antenna geometry over the fractal dimension using genetic algorithms [18, 19] directly. The authors present the relation between the resonant frequencies and the fractal dimension of the parameterized Koch curves in [20]. In general, the studies in the literature focus on predefined well-known fractal templates such as variants of the infamous Koch curve or the Sierpinski carpet. On the other hand, restricting the geometry a priori limits the practical applications.

As a novel approach, we present fractal wire antenna geometries based on the FIFs. In contrast to the literature, we don’t assume a predefined topology in this study, and the designer is in full control of the antenna’s shape by setting a few interpolation points and contraction factors. The interpolation points and the contraction factors of the FIFs can be used to optimize the antennas for a specific purpose. Then, we investigate the effects of fractal parameters on antenna radiation properties, namely the resonant frequencies, the bandwidth, the radiation patterns, gain, and input impedance.

Fractal interpolation is a technique used to construct continuous functions whose graphs are fractals based on iterated function systems (IFS) [21, 22]. Following the pioneering research, FIFs have been applied in geometric design, signal processing, and wavelet theory in the context of engineering, physics, and chemistry [23, 24]. FIFs provide non-smooth alternatives to traditional smooth interpolation techniques and are more suitable for irregular curves that display self-similarity.

Fractal interpolation is an iterative procedure, and each iteration can be considered a prefractal. The various antennas can be constructed associated with each of the prefractals. The skeleton of the antenna geometry can be determined by the given interpolation points. Additionally, the FIFs have free parameters that can be used to manipulate the geometry to alter the fractal dimension and the symmetry of the structure. The antenna can be optimized by changing the interpolation points and the free parameters. Therefore, it can be constructed without a predefined fractal template in order to optimize the antenna performance. Note that several structures such as the Koch curve can also be obtained by specific choice of FIF parameters.

To demonstrate the proposed approach, we generated a simple curved wire dipole antenna using FIFs based on affine transformations. The parameters of affine transformations consist of contraction factors on the horizontal axis and scaling factors on the vertical axis, which simply controls the antenna geometry.

images

Figure 1: FIF’s corresponding to prefractals for the parameters in Table 1. The red dots indicate interpolation points. The fractal dimensions for γi=γ0.25, γ=0.4, and γ=0.8 are D=1, D1.34, and D1.84 for the graph of h*(x) respectively.

Table 1: The parameters of FIF wire prefractal antennas for γ=0.2 and 0.4 in Fig. 1

Parameter Value Description
N 5 Number of interpolation points
[ziyi], i=1,2,,5 [0.50]+[0 5 1015200 10 0-100] Interpolation points (mm), yi=h[k](zi)
αi, i=1,2,3,4 0.25 z-scaling factor
βi, i=1,2,3,4 yi-yi-120 y-scaling factor
ui, i=1,2,3,4 zi-1 z-translation
vi, i=1,2,3,4 yi-1 y-translation
r0 0.005 Radius of the antenna wire (mm)

The scattering parameter |S11|, the input impedance, the gain, and the bandwidth of the constructed structure are calculated via extensive numerical simulations.

We have observed that even the simple structure can show multiband behavior for prefractals obtained at each iteration. Using the proposed procedure, the designer has flexibility in the determination of the skeletal structure of the antenna first. Afterwards, the vertical scaling parameters that are particularly significant on fractal properties can be used to optimize the antenna for a specific application. By means of this flexibility, the technique can be extended to design effective antennas confined in a limited spaceespecially.

II. FRACTAL INTERPOLATION

Let the set of interpolation points, {[xi,yi]T2:i= 0,1,2,,N} where x0<x1<<xN be given, and the continuous function h:[x0,xN] that satisfies h(xi)=yi, be the interpolation function.

We can construct an IFS from a set of contractive shear transformations wi:22, i=1,2,,N, of the form

wi([xy])=[αi 0βiγi][xy]+[uivi], (1)

such that its attractor is the graph of continuous function h. Clearly, 0|αi|,|γi|<1,i. The contraction factor of wi is than σi=max{|αi|,|γi|}, and the contraction factor of the IFS is σ=maxiσi.

Following the steps in [23] and choosing γi’s as free parameters, one can construct wi’s in such a way that the line segment between [x0,y0]T and [xN,yN]T is mapped to the line segment between [xi-1,yi-1]T and [xi,yi]T. Therefore, the parameters must be chosen to satisfy

αi =xi-xi-1xN-x0,ui=xNxi-1-x0xixN-x0,
βi =yi-yi-1xN-x0-γiyN-y0xN-x0, (2)
vi =xNyi-1-x0yixN-x0-γixNy0-x0yNxN-x0.

Denoting as the space of continuous functions h:[x0,xN] such that h(x0)=y0 and h(xN)=yN with a metric d(h,g)=max{|h(x)-g(x)|,h,g}, lets us define a transformation T: that satisfies

(Th)(x) =βili-1(x)+γih(li-1(x))+vi, (3)
li(x) =αix+ui  i=1,2,,N,

for x[xi-1,xi]. T is a contraction in the metric space with contraction factor σT=max{|γi|} and has a unique fixed point h*, i.e. (Th*)(x)=h*(x),x[x0,xN] [22]. For any h[0], the sequence of functions for k=1,2,

h[k](x)=(Th[k-1])(x)x[x0,xN], (4)

converges to h*, i.e.,

limkh[k](x)=h*(x),x[x0,xN]. (5)

Furthermore, the points on the attractor of the IFS is determined by the function h* since

(Th)(αix+ui)=βix+γih(x)+vi,x[xi-1,xi]. (6)

We consider each set {[x,h[k](x)]T2,x[x0,xN]} associated with h[k] as a prefractal and a candidate antenna. Given the transformations wi, and h[0](x)0, we can construct the geometry of the antenna using (3) and the random iteration algorithm for IFS [22]. The convergence rate to final attractor depends on the contraction factor σT. Given 0<ϵ1, the convergence can be assumed if

d(h[k],h[k-1])σTk-1d(Th[0],h[0])
=σTk-1maxx[x0,xN]{|βix-uiαi+vi|}i=1Nϵ, (7)

is satisfied. Clearly, σT depends on the number of interpolation points and chosen γi; hence the designer has two means of controlling how fast the convergence to h* is. ϵ can be chosen according to the finest detail that can be manufactured in practice, and the necessary number of iterations, k, in (3) determined accordingly. We also have

d(h*,h[k])σT1-σTd(h[k-1],h[k]), (8)

in order to measure how close the prefractal associated with h[k] is to the fractal associated with h*.

Note that the selection γi’s has a significant impact on the overall topology of the FIF, as depicted in Fig. 1.

The fractal dimension, D, of the final attractor of the associated IFS satisfies

D=1+{log(i=1N|γi|)log(N),i=1N|γi|>10,otherwise. (9)

Hence 1D<2 if the interpolation points are spaced equally. Clearly, we have absolute control of the fractal’s dimension and the complexity [22].

If γi=γ=0,i, h* corresponds to the linear interpolator. Besides, the small contraction factor (σT0) yields 1D fractals without much detail, and associated prefractals are not of much interest. The choice of i=1N|γi|>1 results in fractals with self similar and symmetric structures with fine details. On the other hand, |γi|1 implies i=1N|γi|N and as a result D2. The attractor associated with the underlying IFS has finer details in this case. However, the details of the associated prefractals may not be suitable for manufacturing for large k.

The parameter ϵ in (7) can be chosen with respect to the finest producible detail. Hence, k can be determined automatically. Then the error estimate between the prefractal at kth iteration and the attractor can be estimated by (8). Although the fractal structures are more interesting in terms of radiation properties as D2, the practical realization can be cumbersome for a given accuracy due to large number of required iterations, yielding extremely fine details.

The main advantage of using FIF as the basis for the structure of the antenna is its flexibility. The designer can optimize the structure by a few number of points on the structure ({xi,h(xi)}i=1N) and altering the free parameters ({γi}i=1N) for desired radiation properties.

images

Figure 2: HFSS 3D model for γ=0.4, k=10. (b) Corresponding mesh. yz-plane is set to be a perfect electric symmetry boundary to for a smaller FEM model. The other outer surfaces of the mesh region is set to be radiation boundaries.

Table 2: The properties of the simulated FIF wire prefractal antennas

γ k Frequency[GHz] Bandwidth [MHz] Gain Impedance [Ω]
- 1 1.807 126 0.47 67.2-j2.0
5.277 146 1.36 53.3-j0.0
0.2 3 1.653 108 0.40 64.4-j1.4
4.838 115 1.27 45.5-j0.1
10 1.602 90 0.38 64.2-j0.5
4.636 94 1.01 46.6-j0.1
0.4 3 1.088 56 0.2 57.4-j0.8
3.043 37 0.42 34.4-j0.1
4.832 49 2.09 78.9-j0.4
6.484 47 2.09 83.5+j0.1
10 0.866 42 0.13 55.2-j0.9
2.468 28 0.29 35.4+j0.3
3.824 39 1.40 62.3-j2.8
5.058 45 1.90 46.8-j1.2

images

Figure 3: |S11| versus the frequency f for varying k and γ. and indicates the prefractal and the reference bend wire antenna corresponding to k=1 case, respectively.

images

Figure 4: The H-plane (0,0) – plot[mark=*] coordinates (0.25cm,0) – (0.5cm,0);: ϕ=0), and E-plane (0,0) – plot[mark=*] coordinates (0.25cm,0) – (0.5cm,0);:  ϕ=90) normalized radiation patterns for the first four modes in frequency range 0.8GHzf8GHz with γ=0.4.

III. ILLUSTRATIVE EXAMPLE

To illustrate the approach, we present FIF perfectly conducting thin wire prefractal antennas oriented along z-axis embedded in yz-plane as a proof of concept.

The simulations were carried out using Ansoft High Frequency Structural Simulator (HFSS)™, on an Intel Xeon based workstation with 32 physical cores and 256 Gb of memory.

The antennas have been fed through a gap of 1 mm located at the origin with a 50 Ω lumped port. The common parameters for the design are listed in Fig. 1. The prefractal curves that form the antennas are obtained by running iterations in (4) with Julia programming language [25]. The generated curves are imported to HFSS for further processing. A circle of radius r0= 5 μm has been extruded along the imported path to create the 3D model. The 3D model is simplified to exclude irrelevant details with respect to operating wave length. The 3D models are simulated using finite element method (FEM). The mesh used in FEM has been fine-tuned with adaptive meshing. Only half z0 plane is considered with an electric symmetry boundary at xy-plane (Fig. 2 (a)). The largest FEM model had 852,603 mixed order tetrahedral elements for the case with γ=0.4 and k=10 (Fig. 2 (b)).

The resonant frequencies, 10 dB bandwidths, peak gains, and feed point impedances are listed in Table 2. The radiation efficiencies have been confirmed to be unity in all cases listed, as the antennas have been assumed to be perfect electric conductors. The frequency sweep analysis in the range of 0.8GHzf8GHz is shown in Fig. 3. Note that the case with k=1 corresponds to linear interpolation over the set {xi,h(xi)}i=1N, and its shape is independent of the contraction factors {γi}i=1N. It is a simple bend wire dipole antenna. The decrease in the first mode’s frequency and the emergence of several other modes is apparent with respect to the reference bend wire dipole.

When the contraction factor is close to 0, the prefractals in each iteration converge to a simple wire antenna with slight decrease in the resonant frequencies for increasing k (Table 2). This is expected since the physical length of the antenna increases with successive iterations of FIF as well. More interesting results are observed when γ=0.4. Several new bands of operation with excellent matching emerge as k increases. Besides, the deviation in resonant frequencies is more pronounced compared to γ=0.2. For larger contraction factors, the antenna is still confined to the same space compactly, although it is electrically longer.

The normalized radiation patterns for γ=0.4 are presented in Fig. 4. The multi-directional radiations patterns emerge for various operating frequencies. Note that although the antenna’s orientation was kept fixed along z-axis, the radiation patterns are almost perpendicular in mode 1 (f=0.866 GHz) and mode 4 (f=5.058 GHz) for k=10 and γ=0.4. The fractal structure allows such possibilities, which would not be available in classical wire antennas.

IV. CONCLUSION

In this article, we propose an approach based on FIFs to design fractal wire antennas. In this approach the geometry need not be predetermined, but can be altered flexibly, in contrast to the fractal antenna studies in the literature. The geometry of the antenna can be controlled by the free parameters, i.e., the interpolation points and the contraction factors of the FIF. Therefore, the antenna’s radiation properties can be controlled directly. The FIF parameters can be adjusted to optimize the performance for the desired antenna properties in terms of gain, radiation pattern, and matching. Furthermore, the optimization can be carried out for multiple bands of operation under spatial constraints.

One of the advantages of the proposed approach is the possibility to bound errors between the prefactals obtained in the intermediate iterations of fractal generation and ideal fractals. This is particularly important because it is impossible to manufacture the infinite details of the ideal fractal. We can determine the required number of iterations a priori for a given manufacturing tolerance based on the constructiveness of the underlying transformations leading to FIF for given antenna performance measures.

The fractal nature of the designed geometries allows the apparent electrical length of the antenna to be larger than the equivalent dipole fitted to the same limited space. In other words, the antenna size can be miniaturized relative to classical structures while operating at low frequencies. Additionally, the self-similarity of the generated fractals results in multiband behavior. These properties render proposed antennas suitable for mobile and wearable wireless applications that require long-range communication especially.

ACKNOWLEDGMENT

This study has been fully supported by Dokuz Eylül University, Scientific Research Fund under project 2021.KB.FEN.004 and partially supported under project 2021.KB.FEN.007.

REFERENCES

[1] B. B. Mandelbrot, The Fractal Geometry of Nature, vol. 1, WH Freeman New York, 1982.

[2] D. H. Werner and S. Ganguly, “An overview of fractal antenna engineering research,” IEEE Antennas and Propagation Magazine, vol. 45, no. 1, pp. 38-57, 2003.

[3] W. J. Krzysztofik, “Fractal geometry in electromagnetics applications from antenna to metamaterials,” Microwave Review, vol. 19, no. 2, 2013.

[4] N. Sharma and V. Sharma, “A journey of antenna from dipole to fractal: A review,” J. Eng. Tech., vol. 6, no. 2, pp. 317-351, 2017.

[5] P. Prabhu and S. Malarvizhi, “Compact dual-band hybrid-fractal MIMO system for UMTS and LTE mobile applications,” Applied Computational Electromagnetics Society (ACES) Journal, pp. 135-140, 2019.

[6] D. H. Werner, R. L. Haupt, and P. L. Werner, “Fractal antenna engineering: The theory and design of fractal antenna arrays,” IEEE Antennas and Propagation Magazine, vol. 41, no. 5, pp. 37-58,1999.

[7] D. Werner and P. Werner, “Frequency-independent features of self-similar fractal antennas,” Radio Science, vol. 31, no. 6, pp. 1331-1343, 1996.

[8] J. Anguera, C. Puente, C. Borja, and J. Soler, “Fractal shaped antennas: A review,” Encyclopedia of RF and Microwave Engineering, 2005.

[9] H. Sagan, Space-filling Curves, Springer Science & Business Media, 2012.

[10] A. Ismahayati, P. J. Soh, R. Hadibah, and G. Vandenbosch, “Design and analysis of a multiband koch fractal monopole antenna,” in 2011 IEEE International RF & Microwave Conference, pp. 58-62, IEEE, 2011.

[11] C. P. Baliarda, J. Romeu, and A. Cardama, “The Koch monopole: A small fractal antenna,” IEEE Trans. on Antennas and Propagation, vol. 48, no. 11, pp. 1773-1781, 2000.

[12] C. Puente-Baliarda, J. Romeu, R. Pous, and A. Cardama, “On the behavior of the Sierpinski multiband fractal antenna,” IEEE Trans. on Antennas and Propagation, vol. 46, no. 4, pp. 517-524, 1998.

[13] A. Bhattacharya, B. Roy, S. K. Chowdhury, and A. K. Bhattacharjee, “Design and analysis of a Koch snowflake fractal monopole antenna for wideband communication,” Applied Computational Electromagnetics Society (ACES) Journal, pp. 548-554, 2017.

[14] P. Chowdary, A. M. Prasad, P. M. Rao, and J. Anguera, “Design and performance study of Sierpinski fractal based patch antennas for multiband and miniaturization characteristics,” Wireless Personal Communications, vol. 83, no. 3, pp. 1713-1730, 2015.

[15] K. Vinoy, K. Jose, V. Varadan, and V. Varadan, “Hilbert curve fractal antenna: A small resonant antenna for VHF/UHF applications,” Microwave and Optical Technology Letters, vol. 29, no. 4, pp. 215-219, 2001.

[16] J. Anguera, A. Andújar, J. Jayasinghe, V. V. S. S. S. Chakravarthy, P. S. R. Chowdary, J. L. Pijoan, T. Ali, and C. Cattani, “Fractal antennas: An historical perspective,” Fractal and Fractional, vol. 4, no. 1, 2020.

[17] A. Karmakar, “Fractal antennas and arrays: A review and recent developments,” Int. J. of Microwave and Wireless Tech., vol. 13, no. 2, pp. 173-197, Mar. 2021.

[18] H. Oraizi and S. Hedayati, “Miniaturized UWB monopole microstrip antenna design by the combination of Giusepe Peano and Sierpinski carpet fractals,” IEEE Antennas and Wireless Propagation Letters, vol. 10, pp. 67-70, 2011.

[19] D. H. Werner, P. L. Werner, and K. Church, “Genetically engineered multiband fractal antennas,” Electronics Letters, vol. 37, no. 19, pp. 1150-1151, 2001.

[20] K. Vinoy, K. Jose, V. Varadan, and V. Varadan, “Resonant frequency of Hilbert curve fractal antennas,” in IEEE Antennas and Propagation Society Int. Symp. 2001 Digest. Held in conjunction with: USNC/URSI National Radio Science Meeting (Cat. No. 01CH37229), vol. 3, pp. 648-651, IEEE, 2001.

[21] K. Vinoy, J. K. Abraham, and V. K. Varadan, “On the relationship between fractal dimension and the performance of multi-resonant dipole antennas using Koch curves,” IEEE Trans. on Antennas and Propagation, vol. 51, no. 9, pp. 2296-2303, 2003.

[22] M. F. Barnsley, Fractals Everywhere, Academic Press, 2014.

[23] M. F. Barnsley, “Fractal functions and interpolation,” Constructive Approximation, vol. 2, no. 1, pp. 303-329, 1986.

[24] M. A. Navascués, A. K. B. Chand, V. P. Veedu, and M. V. Sebastián, “Fractal interpolation functions: A short survey,” Applied Mathematics, vol. 2014, 2014.

[25] J. Bezanson, A. Edelman, S. Karpinski, and V. B. Shah, “Julia: A fresh approach to numerical computing,” SIAM Review, vol. 59, no. 1, pp. 65-98, 2017.

BIOGRAPHIES

images

Gizem Kalender received a degree in electrical and electronics eng. from Çukurova University, Turkiye, in 2013 and her M.Sc. degree from Dokuz Eylül University in 2015. She is a Ph.D. candidate and holds a research assistant position at Dokuz Eylül University. Her main research interests are antennas, fractal theory, and metamaterials.

images

Emine Yeşim Zoral graduated from Dokuz Eylül University, Izmir, Turkiye, in 1990. She received the Ph.D. from the Illinois Inst. of Tech., Chicago, USA, in 1999, in electrical and electronics eng. She has been working as a professor at the Department of Electrical and Electronics Eng. at Dokuz Eylül University since 2000. Her main research interests are perturbation techniques in electromagnetic theory, microwave circuits, dielectric resonators, and antennas.

images

Serkan Günel completed his Ph.D. at Dokuz Eylül University, Inst. of Natural Sciences, Dept. of electrical and electronics eng., in 2006. His main research areas include nonlinear system analysis, information theory, and electromagnetic theory. His current research topics cover analysis of nonlinear dynamics via entropy measures, design of fractal antennas based on iterated function systems, analysis of resonant electromagnetic structures via perturbation of matrices, and the applications of cluster synchronization of chaotic systems. Currently, he is employed as an associate professor in Dokuz Eylül University.

ABSTRACT

I. INTRODUCTION

images

II. FRACTAL INTERPOLATION

images

images

images

III. ILLUSTRATIVE EXAMPLE

IV. CONCLUSION

ACKNOWLEDGMENT

REFERENCES

BIOGRAPHIES