See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/333021016

Design, Test and Optimization of Inductive Coupled Coils for Implantable

Biomedical Devices

Article in Journal of Low Power Electronics · March 2019

DOI: 10.1166/jolpe.2019.1597

CITATIONS READS

15 775

6 authors, including:

Mengyao Yuan Himmy Tam

University of Glasgow Imperial College London
12 PUBLICATIONS 115 CITATIONS 1 PUBLICATION 15 CITATIONS

SEE PROFILE SEE PROFILE

Muhammad Ali Imran Hadi Heidari

University of Glasgow University of Glasgow
1,030 PUBLICATIONS 21,307 CITATIONS 284 PUBLICATIONS 3,476 CITATIONS

SEE PROFILE SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Healthcare informatics View project

Towards Miniaturised Magnetomyography View project

All content following this page was uploaded by Himmy Tam on 27 May 2019.

The user has requested enhancement of the downloaded file.

 Copyright © 2019 American Scientific Publishers Journal of
All rights reserved Low Power Electronics
Printed in the United States of America Vol. 15, 76–86, 2019

Design, Test and Optimization of Inductive Coupled

Coils for Implantable Biomedical Devices
Jinwei Zhao, Rami Ghannam, Mengyao Yuan, Himmy Tam, Muhammad Imran, and Hadi Heidari∗
School of Engineering, University of Glasgow, G12 8QQ, Glasgow, United Kingdom

(Received: 20 December 2018; Accepted: 10 February 2019)

Biomedical implant devices are fast becoming a growing part of the healthcare industry. Providing
power to these devices in such a confined area is a critical challenge. Consequently, resonance-
based wireless power delivery provides a harmless yet effective way for powering these implantable
biomedical devices. This technique relies on transferring power via the inductive coupling technique.
In this regard, optimizing the quality factor and matched resonant frequency is required to achieve
high efficiency. However, the efficiency depends on the space available for the coil and the sepa-
ration distance between the two coils. In our case, the minimum separation distance between the
two coils needs to be at least 2 cm. Therefore, we demonstrate the design, simulation and exper-
imental procedure of an optimized wireless power delivery system for bio-implantable applications
with various considerations for size limitations. Our design delivers 68 mW output power to a 50-
load with an efficiency of 67% in vitro test and 74.8% in the FEM simulation.
Keywords: Implantable Electronics, Wireless Power Transfer (WPT), FEM Device Modelling,
Coil Design, High-Efficiency Inductive Power Coupling.
IP: 155.198.9.87 On: Mon, 27 May 2019 20:09:35
Copyright: American Scientific Publishers
Delivered by Ingenta
1. INTRODUCTION harvesting mechanism, heat is considered as “wasted”
The application of implantable and wearable electronic energy. However, Thermal Electric Generator (TEG) can
devices is rising rapidly in modern medicine due largely be used to make use of this heat energy. They can be com-
to their capability of both real-time monitoring and stim- bined with other energy harvesting devices to detect differ-
ulating organs locally.1–4 The power supply in most ence in human body temperature.22 Photovoltaic (PV) cell
conventional biomedical implant devices relies on batter- can deliver high output power from a small area, which
ies, which will require occasional surgical replacement or reduces system complexity in comparison to WPT and
recharging via a transdermal or percutaneous approach. PUEU. Compared with other power harvesting methods,
These procedures are either time-consuming or prone WPT is the only U.S. FDA (Food and Drug Administra-
to infection. Emerging power harvesting and generation tion) approved method.23 The high conversion efficiency
devices such as photovoltaic cells, piezoelectric generators and power delivered makes it more and more applied
or wireless power transfer (WPT) systems have previously in implantable devices.23 However, short power transfer,
5–20 tissue limitation, and complexity of tuning and circuitry
been investigated in the literature. There are several
factors affecting the device performance and our ability are still huge challenges.24 Near-Field Inductive Coupling
to surgically implant these device under the skin, which method is the viable candidate for implantable power har-
mainly are: (1) size, (2) harvested (or delivered) power vesting because it is easy to tune, advantageous of power
and (3) conversion efficiency.21 Table I shows a com- delivered and higher efficiency.23
parison between our work and other implantable power Magnetically coupled resonant wireless power delivery
harvesting topologies in terms of the key factors men- systems consist of two coupled air-core coils, which act
tioned above. The Piezoelectric Ultrasound Energy Har- as a transmitter (primary) coil and a receiver (secondary)
vester (PUEH) can deliver power at a long separation coil. When connecting an AC power source to the primary
distance and enable small geometry of the implants.9 How- coil (outside of tissue), a time-varying magnetic field will
ever, the amount of power is low (less 1 mW). In all energy induce a voltage in the secondary coil. To achieve maxi-
mum power delivery frequency, a strong magnetic field is
∗
Author to whom correspondence should be addressed. required, which leads to considerable energy loss in addi-
Email: [email protected] tion to various safety problems.25 There are a variety of

76 J. Low Power Electron. 2019, Vol. 15, No. 1 1546-1998/2019/15/076/011 doi:10.1166/jolpe.2019.1597

Zhao et al. Design, Test and Optimization of Inductive Coupled Coils for Implantable Biomedical Devices

Table I. The implantable power harvesting methods comparison state

of arts (primary coil, secondary coil).

Powering Power
scheme Size (mm) Specifications (mW) Efficiency Ref.

PUEG L:0.7 W:0.7 T:1.4 1 MHZ 0.18 50% [17]

PUEG L:1 W:1 T:1.4 1 MHZ 0.36 50% [17]
PV L:1.3 W:1.3 T:1 Visible-NIR 25.4 20% [9]
PV L:2.7 W:3.3 NIR 3.44 27.7% [10]
Fig. 1. The schematic of two-coil near-field inductively coupled power
TEG L:2.4 W:3.3 T:1.1 27  C 0.15 N/A [19]
TEG D:9.3 T:1.4 20–25  C 0.03 N/A [19]
transfer system.
NFIC D:(60, 25) T:(10, 5) 13.56 MHZ 10 58.2% [14]
NFIC D:(30, 20) 8 MHZ 68.7 63.9% and SAR safety limits, while enabling a high efficiency
[16]
NFIC D:(60, 20) T:(10, 2) 3.25 MHZ 68 67% and a more compact device size as shown in Table I.
This work
The two-coil near-field inductive coupled system is
shown in Figure 1. The transmitting and receiving coils
implantable applications that use WPT as the main source can be regarded as two inductances coupled by a mutual
of power, such as cochlear implants, neurostimulators and inductance, M. Since the coils are resistive and due to stray
cardiac pacemakers.26 capacitance, the losses attributed by these components can-
In fact, WPT can be achieved using three main methods: not be ignored. The effect of these unwanted parameters
(1) Inductively Coupled Method (ICM); (2) Capacitively can be predicted using COMSOL, which is a macroscopic
Coupled Method (CPM); and (3) Ultrasonic Power Trans- finite element simulations tool. Thus, the power transfer
fer Method (UPM). Compared with other methods, the can be optimized by changing the relative parameters of
inductively coupled method enables the transfer of power inductive coils. Impedance matching is one of the tech-
outputs ranging from microwatts to several watts at rel- niques for optimizing the output efficiency. Here, reactive
atively high efficiency in the near-field range.26 The aim power is always lagging in power transfer systems based
of this paper is twofold: (a) to demonstrate the size con- on inductive coils. Thus, capacitors are used to compen-
straints of implantable near-field inductively coupled coils, sate and reduce the reactive power in impedance matching,
and (b) report the optimized power delivery efficiency which will increase the transfer efficiency and stability of
IP: 155.198.9.87 On: Mon, 27 May 2019 20:09:35
through human tissue, which is achieved by improving the the system.
Copyright: American Scientific Publishers
quality factor and resonant matching frequency inDelivered
the coil. by Ingenta
2.2. Numerical Analysis of kQ Product
Three major requirements must be fulfilled for achieving
2. PRINCIPLES OF INDUCTIVE COIL
high power transfer efficiency. Firstly, the product of the
2.1. Near-Field Inductively Coupled Coils
Coupling Coefficient (k) and Quality Factor (Q) must be
Near-field inductive power transfer can be achieved if the optimized. Secondly, the system must be operated at the
coils are aligned in close proximity to one another. There appropriate resonant frequency. Finally, impedance match-
are different coil technologies for WPT. For instance, the ing is necessary to improve the power transmission effi-
2-coil, 3-coil and 4-coil systems are examples of multi-coil ciency when the quality is already optimal.
systems. Multi-coil systems enable greater power trans- The transmission efficiency  can be determined using
fer efficiency and a larger separation distance between Eq. (1):26 30
the coils. However, this is predominantly useful for
small power delivery in the mid- and far-field frequency k2 Q Q
26–29 link =  1 2 (1)
ranges. In our case, we adopt the 2-coil system due to 1 + 1 + k2 Q1 Q2 2
its low cost, simple design, easy tuning and large power
delivery.27 However, to prevent tissue loss, we must ensure Consequently, the product of k and Q affect the magni-
that our system transfers less than 100 mW through the tude of the power transmittance efficiency. The Coupling
body and that the operating frequency is much greater than coefficient can be determined using:26
23
100 kHz. Consequently, due to specific absorption rate M
(SAR) safety limits, our work focuses on the design and k= √ (2)
L1 L2
optimization of a 2-coil system that is operated in the MHz
range, which will transfer 68 mW of power with 67% high Where M is the mutual inductance between the two coils
efficiency.23 As shown in Table I, increasing the operating with inductance L1 and L2 , respectively. The k is normal-
frequency enables us to induce the coil geometry. How- ized and can be compared in different coil geometries.
ever, the amount of declined power and the efficiency will Nevertheless, it is very hard to improve the value of k
also be reduced. The opposite is true if we reduce the oper- when the geometry is limited in some scenario, such as
ating frequency. Our design satisfies the required power implantable application. In this case, increasing the value

J. Low Power Electron. 15, 76–86, 2019 77

Design, Test and Optimization of Inductive Coupled Coils for Implantable Biomedical Devices Zhao et al.

of Q is the best way to increase the power transmittance magnetic field. The coil’s self-resonant frequency (fself 
efficiency. can be calculated using:29
The coil quality factor can be determined using
1
Eq. (3):26 29 fself =  (9)
2 Lself Cself
Leff 2Leff × 1 − f 2 /fself
2

Q= = (3) The coil’s self-resonant frequency should be 3 or 4 times
ESR Rdc × 1 + f /fh 
2 2
greater than the operating frequency to ensure the coil is
where Leff is effective inductance of the coil, f is the still dominated by inductance. Since the coil design is both
operating frequency, fself is the coil’s self-resonant fre- multi-turn and multi-layer solenoid coil, the power loss
quency, fh is the frequency at which active power is half in the coil is predominantly caused by proximity effects
of reactive power in the coil. Lself is the self-inductance rather than skin effect. Thus, fh can be determined using:
and RDC is the DC resistance of the coil. From Eq. (1), √
2 2
increasing the operating frequency improves the coil qual- fh =  (10)
ity factor, which maximizes the power-transfer efficiency. rs2 0 Nt Ns 
However, when the operating frequency is increased to the
Where is the conductivity of the material, Nt is the
coil’s self-resonant frequency, the transmission efficiency
number of turns, Ns is the number of strands,  and are
drops to nearly zero. Since the coil turns are very close to
ratios related to coil geometry, rs is the radius of strands.
each other, they suffer from parasitic capacitance.25 29 31–33
Similarly, the mutual inductance, M is significantly
There are two types of parasitic capacitance in solenoid
dependent on k, as previously demonstrated in Eq. (2). The
coil: (a) capacitance between turns (Cm  and (b) capaci-
mutual inductance can also be determined using:34
tance between turns (Cb . The parasitic capacitances in the
coils can be determined using:29 V2
M12 = (11)
 /4 Dc Rw
jI1
Cb = 0 r × d (4)
0 T + r Rw 1 − cos  Where V2 is the voltage between the secondary coil and
 /4 Dc Rw I 1 is the current across the primary coil and  is the angu-
Cm = 0 r × IP: 155.198.9.87 dOn: (5) larMay
Mon, 27 frequency
2019 of the system.
20:09:35
0 T + r Rw 1 − cos  + 0 5r h
Copyright: American Scientific Publishers
In conclusion, the coil quality factor can be maximized
Delivered by Ingenta
by increasing the operating frequency to several megahertz
Where 0 , r are the permittivity of material, T is the
thickness of coil, Dc is the diameter of the coil, Rw is the and by using a greater number of turns, which increases
radius of wire and h is the separation between two layers. the effective inductance. Increasing the inductance will
The total parasitic capacitance can be calculated using: 29 33 decrease k, which results in a poorer link efficiency.
In design procedure, it is always required to consider the
 
1 Nt trade-off point between k and Q. However, implantable
Cself = 2
× Cb N t −1N l +C m 2i −12
N l −1 biomedical devices suffer from coil size limitations, fre-
Nt Nl  i=1
(6) quency limitation due to safety concerns. Consequently,
Where Nt is the number of turns and Nl is number of the aim of this paper is to demonstrate the design process
layers. The analytical expression of self-inductance is: 29 that is required in order to determine a coil geometry that
yields optimal transmission efficiency.
 
N Dc
L= D ln (7)
2 0 c Dw 3. DESIGN METHODOLOGY
Where the 0 is permeability of vacuum, Dw is the The design procedure of the whole system is shown in
diameter of a single turn and Leff is the effective induc- Figure 2. The first step of the design process involves list-
tance related to frequency. The effective inductance is ing the design constraints based on the limitation of coil
size. Next, we aim to achieve as many coil turns as pos-
higher than the self-inductance at a low frequency, which
sible for better quality factor. Next, we determine whether
is computed as:29
the coil’s self-resonant frequency is higher than the oper-
L ating frequency. This is the most important step, which
Leff = (8) ensures that the coil will still behave like an inductor when
1 − 4 2f 2 LCself
operating at that frequency. The coil’s self-resonant fre-
In essence, the coils can be approximated as a simple quency can be increased by using fewer turns and increas-
RLC circuit, as shown in Figure 1. When the operating fre- ing the separation between the wires. When designing the
quency is higher than the self-resonant frequency, the coil primary coil, there are generally no limitations of the trans-
behaves as a capacitor, which can only set up a very weak mitter coil being placed in vitro.

78 J. Low Power Electron. 15, 76–86, 2019

Zhao et al. Design, Test and Optimization of Inductive Coupled Coils for Implantable Biomedical Devices

No Yes
Inductor?

Reduce Turns

Extract Coil
Geometry: Circuit
Constraints: Properties: Link Performance:
Turns Performance:
Coil size L and R vs f Power
Layers Tuning C
Frequency Stray C kQ product
Thickness Coil L
Coil Separation M Link efficiency
Separation Efficiency
Frequencies

Change Separation
and wire radius
No Yes
Q Optimized?

Fig. 2. The design procedure of implantable inductively coupled coils.

Based on the electromagnetic theory of a circular coil, is the sum of the coil’s parasitic capacitance and the tun-
the magnetic field strength between two coils can be deter- ing capacitance. If the value of tuning capacitance is much
mined using:31 higher than the parasitic capacitance, the effect of para-
sitic capacitance can be ignored. However, if the inductor
Ir2 has inductance within one hundred micros Henrys and is
Hx r =  (12)
2 r 2 + x2  being tuned to resonate at a few megahertz, the parasitic
capacitance will usually cause a resonant frequency shift
where r is the coil radius, x is the separation between two and thus reduce the transmission efficiency at the desired
coils and Ir is the magnitude of current. When the distance operating frequency. Accurately determining the coil para-
x is increased a weaker magneticIP: 155.198.9.87
field On: Mon, 27
will be experienced siticMay 2019 20:09:35
capacitance is important for boosting the transmission
by another coil. Based on the designCopyright:
constraintsAmerican
summa- Scientific Publishers
efficiency.
rized in Table II, the separation between the two Delivered
coils will by Ingenta
The total impedance of the primary coil circuit can be
be 2 cm. To maximize the magnetic field strength that is determined using:
experienced by the implantable coil, the radius of the pri-
  
mary coil should be the same as the distance separating L 1
Z = R + j −
the two coils multiplied by root 2. Therefore, the outer total
1 − 2 LC C1
diameter of the coil is approximately 6 cm. After both  
2 LC2 C1
coils have an acceptable quality factor a tuning capacitor 1− + + jwRC2 (14)
should be connected to both coils, which will make both 1 − 2 LC C2
the primary and implantable coils resonate at the same fre- At resonant frequency, the imaginary part of the
quency. The resonant frequency of the LC tank circuit can impedance is equal to zero. Thus:
be determined using:25
  
1 C
1 −1−w 2 LC2 w 2 R2 C2 + 1+ 2
f=  (13) C1 C1
2 LCtotal  
C
+1−w 2 LC w 2 L+2w 2 L 2 −w 4 L2 C2 = 0 (15)
where L is the inductance of the coil, C is the total capaci- C1
tance of the circuit. The inductance of the coil can be accu-
rately determined by the LCR meter. Ctotal of the circuit where L, C, C1 , C2 is the same component which shown
as Figure 1. C1 is the value of tune capacitor and C2
is the BNC wire capacitance = 7.60 pF, which can be
Table II. Design constrains.
measured precisely. The only unknown parameter is the
Parameters Design value parasitic capacitance of the coil C. By placing differ-
Receiver coil outer diameter 2 cm ent tuning capacitances, different resonant frequencies can
Receiver coil thickness 0.20 cm be observed. In Eq. (13) is the resonant frequency × 2.
Transmitter coil outer diameter 6 cm By solving this equation with a different tuning capacitor
Operating frequency 3–5 MHz and corresponding resonant frequency, the same root can
Power delivered Few mW
be obtained. This is the parasitic capacitance. The resonant

J. Low Power Electron. 15, 76–86, 2019 79

Design, Test and Optimization of Inductive Coupled Coils for Implantable Biomedical Devices Zhao et al.

frequency of the coil circuit when including the parasitic Table IV. The coil properties at 3.25 MHz compared between simula-
tion (S) and measurement (M).
capacitance can be calculated as Eq. (16):
Parameters Primary coil Secondary coil
1
freal =  (16)
2 C + C1 Inductance (M) 141.36 H 58.42 H
Inductance (S) 137.06 H 52.65 H
And the comparison of measured resonant frequency to Resistance (M) 5.60  6.68 
predicted one from Eq. (16) can be plotted as Figure 4. In Resistance (S) 3.50  5.33 
Quality factor (M) 750 230
addition, the coils are tested by powering through a piece Quality factor (S) 761.10 167.77
of ham in vitro analysis, which is shown in Figure 3.

is advantageous in both weakly coupling and strong cou-

4. EXPERIMENTAL RESULTS pling. Compared to the series connection method, the par-
AND COMPARISON allel connection topology suffers from strong coupling
Our experimental results are compared with COMSOL regime. The capacitors applied in our circuit are shown
simulations for verification purposes. The coil dimen- in Figure 4. The analytical solution is simulated by MAT-
sions were shown in Table III. The simulation process is LAB using Eq. (13). The magnitude of the tuning capac-
explained in greater detail in Section 5. Firstly, the coil is itor cannot be an arbitrary value, since it depends on its
wounded by the wounding machine precisely. The primary market value. As previously mentioned, with the excep-
coil currently contains 35 turns only and the coil size is tion of the capacitance, the inductance and resistance of
slightly smaller. The real coils are shown in Figure 3. The the coil can be clearly measured by an LCR meter. Three
coil parameters such as inductance, resistance and para- capacitors are involved in the system, which potentially
sitic capacitance can be measured by LCR Meter, while influences the value of self-resonant frequency. The capac-
the capacitance is difficult to test because the meter yield itance of the BNC wire is from 7.60 pF to 130 pF in
nonsensical reading and the parameter capacitance is too the instrumentation, which did not cause any noticeable
small in picofarad scale. The coil properties are shown in resonant frequency effect, which means that it can be
Table IV with a frequency of 3.25 MHz. neglected. Equation (16) is included to solve self-resonant
Secondly, the impedance matching is necessary to com- frequency mismatch. The analytical solution is perfectly
pensate the reactive power in the IP:circuit. The series
155.198.9.87 On:con-
Mon, 27 May 2019
matched with 20:09:35
real-circuit measurement. Thus, it is obvi-
nection method is utilized because the Copyright: American Scientific
series connection Publishers
ous the Eq. (16) is improved the capacitator approxima-
Delivered by Ingenta
tion compared with Eq. (13) in real circuit design. As all
Table III. The coil geometry that yield optimal quality factor.
the steps are completed, the power transmission efficiency
Parameters Primary coil Secondary coil can be achieved, which is shown in Figure 5. The data
Coil diameter 6 cm 2 cm
from the measurement is triggered by applying the “curve
Wire radius 0.10 mm 0.05 mm fitting” tool in Matlab. The efficiency from the simula-
Total number of turns 35 40 tion is slightly higher than the measurement result because
Turns per layer 35 20 the heating loss from the instrument. The main reason
Number of layers 1 2

5.0 Eq.15
Eq.12
4.5 Measurement

4.0
Frequency (MHz)

3.5

3.0

2.5

2.0

1.5

1.0
0 10 20 30 40 50
Capacitance (pF)

Fig. 4. The resonant frequency as a function of tuning capacitor. The

Fig. 3. The in-vitro test of near-field power transfer system. predicted (Eq. (15)) and experimental results show the good agreement.

80 J. Low Power Electron. 15, 76–86, 2019

Zhao et al. Design, Test and Optimization of Inductive Coupled Coils for Implantable Biomedical Devices

80 selected to be lower than 4 MHz, the capacitance of tun-

ing capacitor is much higher than the parasitic capacitance
contributed by the tissue, thus the tissue effect becomes
60
minuscule.
Efficiency (%)

5. FEM SIMULATION AND DEVICE

40
OPTIMIZATION
COMSOL was used for simulating the magnetic field dis-
tribution and magnetic resonance between the coils. The
coil model was developed in the 2D axisymmetric formu-
20 Simulation lation. 3D simulations were also developed to show the
Measurement
magnetic field distribution. Since the processing time of
Curve fit of measurement
2D simulations is considerably less than 3D, 2D simula-
0 tions were used to process important calculations such as
0 1 2 3 4 5 6 7 8 9 parameter sweep. The axisymmetric principle in 2D mod-
Frequency (MHz) elling is shown in Figure 6. A circle (sphere in 3D) with
Fig. 5. The comparison between COMSOL simulation and experimen- a diameter of 15 cm is aligned as the FEM magnetic field
tal result for the power transmission efficiency. boundary, which means that only the objects inside the cir-
cle will be analyzed otherwise will be neglect. The mate-
rial of the circle is aligned with air, and the coils’ material
is the self-resonant capacitance because the value is not is set up with copper. Now turning point to the geometry
evaluable in measurement and difficult to simulate the of the coils, the geometry parameters are set up as same
frequency-dependent capacitance accurately. as experimental works. According to the design constraint,
The transmission efficiency is determined by the power the separation distance between two coils is less than 2 cm,
delivered to the 50-Ohm secondary resistor divided by where the distance is regarded as the length between the
the input power in the primary side. A 0.5 cm thick circle’s centers of two coils. Next, the Magnetic Field
chunk of ham was used to imitate human tissue while
IP: 155.198.9.87 On: Mon, 27 Physics is added
May 2019 in the model to set up the magnetic to
20:09:35
the receiver coil was covered. WithCopyright: the high coil qual- the electrical interface.
American Scientific Publishers It is based on the Ampere’s Law,
ity factor and accurate resonant frequency prediction, the by and
Delivered Ingenta
the magnetic boundary is set up by the sub-section
two coil systems with the structure mentioned in Table IV called Magnetic Insulation. The initial condition is set by a
can achieve around 67% of transmission efficiency when standstill magnetic field, where the magnetic vector poten-
the relative distance is 2 cm at 3.28 MHz. A total of tial is zero at x-, y- and z-axis respectively. The opera-
68 mW can be delivered to the load resistor which meets tional temperature is set to room temperature 293.15 K,
the power requirement of delivering about few tens’ mW and some material properties such as conductivity permit-
for the implantable device.11 In this experiment, the tis- tivity and permeability are aligned to the built-in mate-
sue is nearly transparent to the system. As the permeabil- rial: air (Ciddor) and copper (AC–DC). The sub-module
ity of tissue is close to the free space, the inductance of called Coil ports are involved in exciting the transmitter
the coil remains unchanged. As the resonant frequency is coil and reviver coil. The relative parameters are utilized

Axis
Coil 1

2D Cross section (Tx)

2D FEM
2D Cross section (Rx)
Coil 2

Fig. 6. The 2D axisymmetric model of the 2-coil near-field inductively coupled power transfer system, which was simulated using COMSOL.

J. Low Power Electron. 15, 76–86, 2019 81

Design, Test and Optimization of Inductive Coupled Coils for Implantable Biomedical Devices Zhao et al.

in this setting, such as turns number, conductivity, cross- because the inductance is increased. The stray capacitance
section area of the coil. Moreover, the conductor model is decreased with higher frequency. With the self-resonant
is set to Homogenize multi-turn coil and the coil can be degrading, the quality will meet a trade-off point where
excited by voltage, current and power. The self-inductance the maximum quality factor is achieved. The maximum
can be determined by exciting a coil with a current, which quality factor for the primary coil is 1480 with 160 turns,
can be instrumented by the software automatically. The while the one for the secondary coil is 380 with 35 turns.
relative results can be achieved when the sending coil is According to the turn numbers, it is easy to find out the
supplied with 30 mA AC current, the receiving is induced relative self-resonant frequency of both coils. Where the
and supply a 50 ohms resister by involving electrical cir- parameter values of the optimized coil can be obtained,
cuit physics in AC–DC module. The quality factor of which is L1 = 2 30 H, C1 = 0 10 pF and L2 = 164 10 H,
the primary coil and the secondary coil are shown in C2 = 6 50 pF. As it is an implantable device, the dimen-
Figures 7(a) and (b), which vary with the turns per layer. sions of the coils are required to be considered, where
With the turn numbers rise, the inductance of coil rises the thickness of the primary coil is 1.64 cm, the thick-
according to the Eq. (7) because of the product of Cself , ness of the secondary coil is 0.17 cm. The width of the
Leff and  are much smaller than 1. With the increase coils can be only achieved if the number of layers can
of turn numbers, the self-resonant frequency is decreased be determined, which will be shown in the next step. The
fh also affect the magnitude of quality factor in Eq. (3)

(a)
(a)
70 Q peak 1500 f h(TX)
fself(TX) 5 2xM
Q TX
Q TX
60 0

1000 4 Q TX>0 Q TX<0

Frequency (MHz)

50 –2
Frequency (MHz)
Quality factor

Quality factor
40 –4
3
500
–6
30 IP: 155.198.9.87 On: Mon, 27 May 2019 20:09:35
Copyright: American Scientific Publishers
20 Delivered by Ingenta –8
2
0
–10
10

50 100 150 200 1 –12

1 2 3 4 5
Number of turns
Number of layers
(b) (b)
30 400 8
fself(RX) Q peak fh(RX) 5000
Q RX Q RX Q RX>0 Q RX<0
25 7
0
Frequency (MHz)
Frequency (MHz)

20 200 6
Quality factor

Quality factor

–5000
15 5

10 0 4 –10000

5 3
–15000

0 –200 2
10 20 30 40 50 0 2 4 6
Number of turns Number of layers

Fig. 7. Quality factor optimization by sweeping number of turns. Fig. 8. Quality factor optimization by sweeping number of layers.
(a) Quality factor and self-resonance frequency of primary coil accord- (a) The quality factor and half-power frequency of primary coil according
ing to number of turns (b) quality factor and self-resonance frequency to number of layers. (b) The quality factor and half-power frequency of
of secondary coil according to number of turns. The number of turns is secondary coil according to number of layers. The dashed line indicates
defined by the maximum quality factor point, which are 160 (primary the number of layers for primary coil and secondary coil, which are 1
coil) and 35 (secondary coil) respectively. and 2 respectively.

82 J. Low Power Electron. 15, 76–86, 2019

Zhao et al. Design, Test and Optimization of Inductive Coupled Coils for Implantable Biomedical Devices

and it can be optimized by changing layer numbers shown For instance, the Resistance for 1-layer primary coil at
in Figures 8(a) and (b). For a great inductive coupling 3.25 MHz is 6.72 Ohms, and 1-layer secondary coil at the
coil, a large value of fh is always required, which means same frequency is 5.34 Ohms. The bandwidth of induc-
the change in coil AC resistance according to the proxim- tive coil is also essential for the quality of designs, which
ity effect and the skin effect is not changed rapidly with regulates the range of operational frequency. To obtain
the frequency increased. The dashed lines in Figures 8(a) the bandwidth and compare between the different number
and (b) provide the maximum layers for the coils because of turns, it is necessary to sweep the relative parameters,
the quality factor cannot less than 0, which are 1 layer for where the results are shown in Figures 9(a) and (b). It
the primary coil and 2 layers for the secondary coil. How- is clear to see that the bandwidth starts rapidly narrowed
ever, the layer should be as small as possible, since the when the turn numbers are greater than a threshold value.
quality factor is inversely proportional to the number of This phenomenon appears when the number of turns of
layers. The self-resonant frequency is slightly influenced
primary and secondary coils are higher than 130 and 10
by layer numbers, but the change is only a 1% difference
respectively. It is always necessary to guarantee the band-
between one layer and two layers. As the half power fre-
width greater than the operational frequency. Due to the
quency is achieved, it is easy to obtain the AC resistance.
constraints (3–5 MHz), the turns number of secondary
coils are limited below 30.
(a)
N1=10
1600
N1=30
N1=50
1400 (a)
N1=70 0.01995 k (frequency)
N1=90
1200
N1=110 0.01990
Quality Factor

N1=130
1000
N1=150 0.01985
N1=170
800 0.01980
Bandwidth

600 0.01975
k

400 IP: 155.198.9.87 On: Mon, 27 0.01970

May 2019 20:09:35
Copyright: American Scientific Publishers
200 Delivered by Ingenta
0.01965
Bandwidth 2 Bandwidth 1
0 0.01960
0 1 2 3 4 5 6 7 8 9 10
0.01955
Frequency (MHz)
(b) 0 2 4 6 8 10
600 N2=10 Frequency (MHz)
N2=20
(b) Number of turns
N2=30
500 20 40 60 80 100 120 140 160 180
N2=40 0.022 0.035
N2=50
Quality Factor

400 Bandwidth 0.021 kTX 0.030

kRX
0.020
Coupling coefficient

Coupling coefficient
0.025
300
0.019
0.020
200
0.018
0.015

100 0.017
0.010

0.016
0 0.005
0 1 2 3 4 5 6 7 8 9 10 0.015
Frequency (MHz) 0 10 20 30 40 50
Number of turns
Fig. 9. The quality factor and bandwidth with different turn numbers.
(a) The quality factor and half-power frequency of primary coil versus Fig. 10. The coupling coefficient optimization. (a) The coupling coef-
frequency for different turn numbers. (b) The quality factor and half- ficient according to the operation frequency. (b) Coupling coefficient
power frequency of secondary coil versus frequency changes and swept according to primary coil number of turns (kTX ) and secondary coil num-
by different turn numbers. The band gap becomes narrowed as the num- ber of turns (kRX ). Increase the number of turns and operating frequency,
ber of turns increases. degrade the coupling coefficient.

J. Low Power Electron. 15, 76–86, 2019 83

Design, Test and Optimization of Inductive Coupled Coils for Implantable Biomedical Devices Zhao et al.

Table V. Ideal design based on FEM simulation. indicate how good the design, and a kQ product needs to
Parameters Primary coil Secondary coil be involved. The value can be achieved by multiplying
the data in Figures 7 and 8 properly. The final optimized
kQ product 9.84 6.47
design is shown in Table V. With this design, efficiency
Turn numbers 20 30
Layer numbers 1 2 of 74.8% can be finally achieved, which is improved by
Diameter (cm) 6.04 2.02 7.80% compared with the initial design. However, this is
Thickness (cm) 0.4 0.2 an ideal design from the simulation, the practical results
Inductance ( H) 60.93 126.65 also depend on some technical conditions, such as the coil
Capacitance (pF) 0.78 5.94
separations in the winding, the gap between the strands.
Resistance ( 2.00 10.67
The efficiency variation due to separation distance is
shown in Figure 11. This shows the experimented simula-
Now turning point to the coupling coefficient, it varies tion results of an unoptimized WPT system. Similarly, we
with the operational frequency, number of turns, num- have used these experiment results to develop a COMSOL
ber of layers and distance between two coils. First, the simulation program that enabled us to optimize a WPT
influence of layers numbers can be neglected, because the system. The optimized system obtained a 74.8% maximum
primary coil is fixed at 1, and 3% difference between conversion efficiency at 2 cm as mentioned above. With
2-layer and 1-layer secondary coil (0.0196 and 0.0190 the distance increased, the flux linkage between two coils
respectively). Secondly, the effect from the various dis- was degraded. This makes the conversion efficiency and
tance can be neglected as well, because the distances for power transfer capability diminished.23 The separation dis-
the implantable devices are limited. The k over operational tance is always constrained by the implantable applications
frequency is shown in Figure 10(a), and the coupling coef- in practice.23
ficients of the primary coil and secondary coil over the
number of turns are shown in Figure 10(b). Based on the 6. CONCLUSION
bandwidth in the previous analysis, the coefficient varies This paper demonstrated a highly efficient WPT system
only with 1.50% difference. Thus, the effect due to oper- design and determined its constraints. The iterative design
ation frequency can be also neglected. There are signifi- procedure discussed in Figure 2 enabled to improve the
cant changes of coefficient when the number of turns is performance
IP: 155.198.9.87
changed. To be specific, the coupling coefficient On: Mon, 27
of pri- May 2019 of a WPT system in implantable applica-
20:09:35
tions. The experimentally
Copyright: American Scientific Publishers achieved transmission efficiency
mary coil drops from 0.0315 to 0.005 when the Delivered
turn num- by Ingenta
of 67% at 3.28 MHz. Moreover, 68 mW of power were
ber is increased up to 180, while the one for secondary delivered to the secondary coil when the relative distance
coil decreases from 0.0215 to 0.0158 when the turn num- between the coils is 2 cm. 68 mW power delivered can
ber rise to 50. In this case, the Quality factor itself cannot easily drive implants with massive consumption, such as
cochlear hearing aids. Device optimization was achieved
100
Experiments
using COMSOL, which was used to obtain an improved
Curve Fit efficiency reaching 74.8%.
unoptimized Simulation An advanced WPT system requires additional compo-
80 optimized Simulation nent to convert and manage the harvested energy in order
to power the load. For example, we need DC–AC, AC–DC
and DC–DC converter. Moreover, implantable sensors pro-
Efficiency (%)

60 duce reactive power. Consequently, in order to optimize the

design of a WPT system, we will need to consider these
components in addition to the reactive power produced by
40 implantable sensors. These issues will be investigated and
reported in our next publication.

20 Acknowledgment: This work was supported by the

Royal Society under grant RSG\R1\180269 and UK Engi-
2 3 4 5 6 7
neering and Physical Sciences Research Council (EPSRC)
Distance (cm) under grant EP/R511705/1.

Fig. 11. The variation of efficiency according to different separation

distances. The experimental data was evaluated from the in-vitro test in References
Section 4, and the data was refined using the built-in curve fit tool in 1. H. Li, A. Shrestha, H. Heidari, J. Le Kernec, and F. Fioranelli,
MATLAB. The unoptimized simulation is developed using COMSOL A multisensory approach for remote health monitoring of older
with the same experimental setup in Section 4. The optimized simulation people. IEEE Journal of Electromagnetics, RF and Microwaves in
is built up with the coil properties in Section 5. Medicine and Biology 2, 102 (2018).

84 J. Low Power Electron. 15, 76–86, 2019

Zhao et al. Design, Test and Optimization of Inductive Coupled Coils for Implantable Biomedical Devices

2. N. Alduais, J. Abdullah, A. Jamil, and H. Heidari, Performance eval- and a hybrid bi-directional data link. IEEE Journal of Solid-State
uation of real-time multivariate data reduction models for adaptive- Circuits 50, 1741 (2015).
threshold in wireless sensor networks. IEEE Sensors Letters 1, 1 19. K. O. Htet, J. Zhao, R. Ghannam, and H. Heidari, Energy-efficient
(2017). start-up power management for batteryless biomedical implant
3. X. Liang, R. Ghannam, and H. Heidari, Wrist-worn gesture sensing devices, IEEE Int. Conference on Electronics, Circuits and Systems
with wearable intelligence. IEEE Sens. J. 19, 1082 (2018). (ICECS) (2018), pp. 597–600.
4. K. O. Htet, H. Fan, and H. Heidari, Switched capacitor DC–DC 20. J. Zhao, R. Ghannam, Q. H. Abbasi, M. Imran, and H. Heidari,
converter for miniaturised wearable systems, IEEE Int. Symposium Simulation of photovoltaic cells for implantable sensory applications.
on Circuits and Systems (ISCAS) (2018), pp. 1–5. IEEE Sensors (2018), pp. 1–4.
5. S. Wen, H. Heidari, A. Vilouras, and R. Dahiya, A wearable fabric- 21. B. Shi, Z. Li, and Y. Fan, Implantable energy-harvesting devices.
based RFID skin temperature monitoring patch. 2016 IEEE Sensor Adv. Mater. 30, 1801511 (2018).
(2016), pp. 1–3. 22. E.-J. Yoon, J.-T. Park, and C.-G. Yu, Thermal energy harvesting
6. K. Oo Htet, R. Ghannam, Q. H. Abbasi, and H. Heidari, Power circuit with maximum power point tracking control for self-powered
management using photovoltaic cells for implantable devices. IEEE sensor node applications. Frontiers of Information Technology and
Access 6, 42156 (2018). Electronic Engineering 19, 285 (2018).
7. K. O. Htet, J. Zhao, R. Ghannam, and H. Heidari, Energy-efficient 23. K. Agarwal, R. Jegadeesan, Y.-X. Guo, and N. V. Thakor, Wire-
start-up power management for batteryless biomedical implant less power transfer strategies for implantable bioelectronics. IEEE
devices. IEEE Int. Conference on Electronics Circuits and Systems Reviews in Biomedical Engineering 10, 136 (2017).
(ICECS) (2018), pp. 597–600. 24. J. Chen, R. Ghannam, M. Imran, and H. Heidari, Wireless power
8. M. Kiani and M. Ghovanloo, An RFID-based closed-loop wireless transfer for 3D printed unmanned aerial vehicle (UAV) systems,
power transmission system for biomedical applications. IEEE Trans- IEEE Asia Pacific Conf. on Postgraduate Research in Microelectron-
actions on Circuits and Systems II: Express Briefs 57, 260 (2010). ics and Electronics (PrimeAsia) (2018), pp. 72–76.
9. Q. Shi, T. Wang, and C. Lee, MEMS based broadband piezoelec- 25. M. Kiani and M. Ghovanloo, The circuit theory behind coupled-
tric ultrasonic energy harvester (PUEH) for enabling self-powered mode magnetic resonance-based wireless power transmission. IEEE
implantable biomedical devices. Scientific Reports 6, 24946 (2016). Transactions on Circuits and Systems I: Regular Papers 59, 2065
10. T. Tokuda, T. Ishizu, W. Nattakarn, M. Haruta, T. Noda, (2012).
K. Sasagawa, M. Sawan, and J. Ohta, 1 mm3 -sized optical neural 26. M. Schormans, V. Valente, and A. Demosthenous, Practical inductive
stimulator based on CMOS integrated photovoltaic power receiver. link design for biomedical wireless power transfer: A tutorial. IEEE
AIP Advances 8, 045018 (2018). Transactions on Biomedical Circuits and Systems 12, 1112 (2018).
11. Y. J. Hung, M. S. Cai, J. F. Chen, H. W. Su, P. C. Jen, P. Chen, and 27. M. Kiani, U.-M. Jow, and M. Ghovanloo, Design and optimization
T. C. Chang, High-voltage backside-illuminated CMOS photovoltaic of a 3-coil inductive link for efficient wireless power transmission.
module for powering implantable temperature sensors. IEEE Journal IEEE Transactions on Biomedical Circuits and Systems 99, 1 (2011).
of Photovoltaics 8, 342 (2018). IP: 155.198.9.87 On: Mon, 27 C. Xiao,
28. May D. Cheng,
2019 and K. Wei, An LCC-C compensated wireless
20:09:35
12. E. Moon, D. Blaauw, and J. D. Phillips, Small-area
Copyright: Si photovoltaics
American Scientificcharging system for implantable cardiac pacemakers: Theory, exper-
Publishers
for low-flux infrared energy harvesting. IEEE Transactions on Elec- by Ingenta
Delivered iment, and safety evaluation. IEEE Transactions on Power Electron-
tron Devices 64, 15 (2017). ics 33, 4894 (2018).
13. E. Moon, D. Blaauw, and J. D. Phillips, Subcutaneous photovoltaic 29. Z. Yang, W. Liu, and E. Basham, Inductor modeling in wireless
infrared energy harvesting for bio-implantable devices. IEEE Trans- links for implantable electronics. IEEE Transactions on Magnetics
actions on Electron Devices 64, 2432 (2017). 43, 3851 (2007).
14. Y. Hung Jr., T.-Y. Chuang, C.-L. Chun, M.-S. Cai, H.-W. Su, and 30. S. Li and C. C. Mi, Wireless power transfer for electric vehicle appli-
S.-L. Lee, CMOS-enabled interdigitated back-contact solar cells for cations. IEEE Journal of Emerging and Selected Topics in Power
biomedical applications. IEEE Transactions on Electron Devices Electronics 3, 4 (2015).
61, 4019 (2014). 31. U.-M. Jow and M. Ghovanloo, Design and optimization of printed
15. R.-F. Xue, K.-W. Cheng, and M. Je, High-efficiency wireless power spiral coils for efficient transcutaneous inductive power transmis-
transfer for biomedical implants by optimal resonant load trans- sion. IEEE Transactions on Biomedical Circuits and Systems 1, 193
formation. IEEE Transactions on Circuits and Systems I: Regular (2007).
Papers 60, 867 (2013). 32. M. Kesler, Highly Resonant Wireless Power Transfer: Safe, Efficient,
16. Y. Guo, D. Zhu, and R. Jegadeesan, Inductive wireless power trans- and Over Distance, Witricity Corporation (2013), pp. 1–32.
mission for implantable devices, 2011 International Workshop on 33. A. K. RamRakhyani, S. Mirabbasi, and M. Chiao, Design and
Antenna Technology (iWAT) (2011), pp. 445–448. optimization of resonance-based efficient wireless power delivery
17. D. Ahn and S. Hong, Wireless power transmission with self- systems for biomedical implants. IEEE Transactions on Biomedical
regulated output voltage for biomedical implant. IEEE Transactions Circuits and Systems 5, 48 (2011).
on Industrial Electronics 61, 2225 (2014). 34. F. Flack, E. James, and D. Schlapp, Mutual inductance of air-cored
18. J. Charthad, M. J. Weber, T. C. Chang, and A. Arbabian, A mm-sized coils: Effect on design of radio-frequency coupled implants. Medical
implantable medical device (IMD) with ultrasonic power transfer and Biological Engineering 9, 79 (1971).

J. Low Power Electron. 15, 76–86, 2019 85

 Design, Test and Optimization of Inductive Coupled Coils for Implantable Biomedical Devices Zhao et al.

Jinwei Zhao
Jinwei Zhao received a B.Eng. degree in Electrical and Electronics Engineering, University of Edinburgh, 2016 and received an M.Sc.
degree in Electric Power, Newcastle University, UK, 2017. He is currently a Ph.D. student in University of Glasgow, from January
2018. The research interest is Energy harvesting, Implantable Systems and Solar cells.

Rami Ghannam
Rami Ghannam (B.Eng., DIC, M.Sc., Ph.D., MIET, SMIEEE) is a lecturer in Electronic and Electrical Engineering with research
interests in the broad field of photonics. Following his Ph.D. from Cambridge University in 2007, Dr. Ghannam has spent the past
ten years in the field of photovoltaics. He has held previous industrial appointments at Nortel Networks, IBM Research GmbH and
Starkon S.A.E. He received his B.Eng. degree from King’s College (top First-Class Honours), as well as his DIC and M.Sc. degrees
from Imperial College London. Dr. Ghannam is also a member of the Microelectronics Lab, where he is investigating the use of PV
cells for energy harvesting applications in wearable and implantable electronic device.

Mengyao Yuan
Mengyao Yuan received a B.Eng. degree in Electrical and Electronics Engineering, University of Glasgow College UESTC, 2018.
She is currently studying the M.Sc. degree in Electrical and Electronics Engineering, University of Glasgow. The research interest is
Energy harvesting for Implantable and Wearable Systems.

Himmy Tam
Himmy Tam received a B.Eng. Degree. In Electrical and Electronics Engineering, University of Glasgow, 2017. Currently studying
the M.Sc. degree in Machine Learning in Imperial College London. The research interest is Robotics and deep learning.

Muhammad Imran
Muhammad Imran (B.Sc., M.Sc., DIC, Ph.D., CENG, FIET, SMIEEE) Fellow IET, Senior Member IEEE, Senior Fellow HEA is a
Professor of Wireless Communication Systems with research interests in self organised networks, wireless networked control systems
and the wireless sensor systems. He heads the Communications, Sensing and Imaging CSI research group at University of Glasgow.
He is an Affiliate Professor at the University of Oklahoma, USA and a visiting Professor at 5G Innovation Centre, University of Surrey,
UK. He has over 18 years of combined academic and industry experience with several leading roles in multi-million pounds funded
projects. He has been awarded 15 patents; has authored/co-authored over 400 journal and conference publications; was editor of
2 books and author of more than 15 book chapters; has successfully supervised over 40 postgraduate students at Doctoral level. He
has been a consultant to international
IP: projects and localOn:
155.198.9.87 companies
Mon, 27in the
May area of self-organised
2019 20:09:35 networks.
Copyright: American Scientific Publishers
Hadi Heidari Delivered by Ingenta
Hadi Heidari is a Lecturer (Assistant Professor) in the school of Engineering and head of the Microelectronics Lab (meLAB) at the
University of Glasgow, UK. He has authored over 90 articles in peer reviewed journals (e.g., IEEE Solid-State Circuits and Systems I
and IEEE Trans. Electron devices) and in international conferences. He has organised several conferences, workshops and special
sessions, e.g., he is founder chairs of UK-China Emerging Technology (UCET) workshop, UK Circuits and Systems (UKCAS) workshop,
and member of organising committee of SENSORS’17–’18, BioCAS’18, PRIME’19 and ICECS’20. He is an IEEE Senior Member, an
Editor for the Elsevier Microelectronics Journal and lead Guest Editor for four journal special issues. He is member of the IEEE
Circuits and Systems Society Board of Governors (BoG), IEEE Sensors Council and IEEE Solid-State Circuits Society Administrative
Committee (AdCom). He has received several best paper awards from IEEE international conferences including ISCAS’14, PRIME’14,
ISSCC’16, and travel scholarship from IEEE NGCAS’17. He has grant portfolio of £1 million funded by major research councils
and funding organizations including European Commission, UK’s EPSRC, Royal Society, Royal Academy of Engineering and Scottish
Funding Council.

86 J. Low Power Electron. 15, 76–86, 2019

View publication stats