1 s2.0 S2211285523006985 Main
1 s2.0 S2211285523006985 Main
1 s2.0 S2211285523006985 Main
Nano Energy
journal homepage: www.elsevier.com/locate/nanoen
A R T I C L E I N F O A B S T R A C T
Keywords: Off-grid connected objects are widely used for environmental monitoring applications such as air and water
Triboelectric nanogenerators quality, soil salinity, weather forecasting and more. In recent studies, mechanical energy harvesting using
Free-standing balls energy harvesting triboelectric nanogenerators (TENGs) have gained in efficiency, even though the design of these TENGs was kept
Norton equivalent model
simple. In this context, the durability of the mechanical contact between triboelectric materials is challenged.
Output current multiplication
Brushless energy conversion
This is why, free-standing balls have been proposed to reduce the friction in TENGs instead of foils, which are
commonly used. Thus, a non-contact gravitational TENG at state-of-the-art performances is proposed in this
work. This two-stage design is based on free-standing balls in a rotor and a contactless capacitive coupler. The
power generation has been investigated for these two stages according to analytical models. A new methodology
to extract the Norton model parameters from a moving TENG is suggested. Through an optimized interdigitated
design of the capacitive coupler, a 27-fold increase of the primary stage power is demonstrated and a state-of-the-
art peak output power density of 24.45 W m− 2 is reached. This study is a leading-edge example of how TENGs
could be adopted to power small remote apparatus, in addition of battery.
* Corresponding author.
E-mail address: blayac@emse.fr (S. Blayac).
https://doi.org/10.1016/j.nanoen.2023.108861
Received 22 May 2023; Received in revised form 30 July 2023; Accepted 2 September 2023
Available online 9 September 2023
2211-2855/© 2023 Published by Elsevier Ltd.
E. Kharbouche et al. Nano Energy 117 (2023) 108861
geometrical shapes, where airflows can be efficiently harvested. In cylinder along the surface or falling top-down between electrodes. Next,
2015, TENGs were used to harvest low-frequency water wave energy the electrical performances of the primary stage are studied as a function
thanks to a rolling ball [14]. Since then, performances of free-standing of ball and triboelectric materials, ball diameter, number of balls and
triboelectric balls in sphere unit were investigated for blue energy har rotating speed. Afterwards, we study and model the secondary stage
vesting [15–18]. This technique allows to reduce frictions between dedicated to the electrical power amplification and transfer to a stator.
triboelectric layers in order to maximize the TENGs durability [19]. In The principle of operation is similar to that of rotary-sliding disk TENGs
return, the effective contact surface is less important, which weakens the [29–31], except that there is no contact between the disks and thus no
triboelectric effect. Using soft balls to achieve enhanced power density contact electrification at this stage, but also no friction loss. Indeed, the
may be a good option [20]. It is also possible to multiply the units for the main drawback of a contact coupling is that sparks occur between the
same TENG in order to achieve high output power density [21–25]. electrodes, which reduces the output performance. Therefore, a capac
Another solution to increase TENGs power density is to use a two- itive coupling has been chosen to avoid brushes in the secondary stage as
stage architecture. In this configuration, a pump TENG injects charges implemented in Fig. S1. Moreover, the capacitive coupler is designed
into a main TENG through a floating layer [26] or a buffer capacitor with interdigitated electrodes [32] to amplify the output current. The
[27]. Thus, the pump TENG initiates the charges by triboelectric effect. output electrical performance is then investigated at the secondary
These charges are then accumulated in the main TENG that conse stage, where we demonstrate and explain a 27-fold increase of electrical
quently converts mechanical to electrical power through electrostatic power. A peak output power of 20.54 mW (24.45 W m− 2) at 500 RPM is
induction in output electrodes. Since output power has a quadratic obtained, which ranks our newly developed device the best in class of its
dependence on accumulated charges [28], a strong multiplication of type (Tables S1 and S2). Using these performance results, an application
output power is obtained. scenario has been suggested for air quality monitoring.
In this study, we report the performance of a two-stage gravitational
energy harvester based on the triboelectric effect and a capacitive 2. Materials and methods
coupling between interdigitated electrodes. The non-contact gravita
tional triboelectric nanogenerator (NCG-TENG) takes advantage of free 2.1. Fabrication of the primary stage
balls in motion into a rotating pipe (primary stage) and a circular con
tactless capacitive coupler (secondary stage), as illustrated in Fig. 1a. We A poly(methyl methacrylate) (PMMA) pipe (internal diameter =
first investigate how ball movements allow conversion of mechanical 26 mm, external diameter = 30 mm, length = 185 mm, Polydis) was
into electrical energy thanks to a complete theoretical analysis. Then, we chosen as a core cylinder to support the rotor (Fig. S2a). The flexible
expose our methodology to extract Norton equivalent model parameters electrodes were etched on an 81.5 mm × 75 mm copper (Cu)-polyimide
and relate them to physical parameters of the TENG. Two new operating (PI) film (thickness = 100 µm, Goodfellow Ltd.) to obtain 1-complemen
modes are presented for the primary stage, with balls rolling inside the tary pair of electrodes (surface of 1 electrode = 0.84E-3 m2). The
Fig. 1. Structure of the non-contact gravitational triboelectric nanogenerator (NCG-TENG). a) Exploded view of the TENG showing the two-stage architecture. b)
Photograph of the TENG in the test bench showing the DC motor and the micro-positioner. c) Focused photo revealing the two stages: the main TENG with the
rectifier and the primary capacitor (primary stage in green), and the capacitive coupler (secondary stage in red).
2
E. Kharbouche et al. Nano Energy 117 (2023) 108861
triboelectric layer is an 81.5 mm × 40 mm film directly stuck on the PI also measured using the digital oscilloscope through a low-noise current
surface with a spray adhesive (3M Company). The tested film materials preamplifier (Model SR570). The secondary open-circuit voltage was
are fluorinated ethylene propylene (FEP) (thickness = 250 µm, ADDEV measured using the digital oscilloscope through a 5 GΩ resistive load.
Materials) and polycarbonate (PC) (thickness = 250 µm, Goodfellow For all experiments, the real-time acquisition was done by Python lan
Ltd.). At one end, the film was wrapped around a PMMA end-cap guage scripts. All finite element method (FEM) simulations were
(external diameter = 25.8 mm, length = 10 mm) illustrated in computed by COMSOL Multiphysics® software. The slow motion movies
Fig. S2a and it was dragged into the pipe in order to create a container in supplementary material were filmed using a compact high-speed
for balls. The balls shown in Fig. S2b were purchased from the CIMAP camera (Sony RX0 II) at 1000 frames per second with a 1920 × 1080
company and selected in different materials (zirconium oxide (ZrO2), pixels resolution.
glass (SiO2) and polyacetal (POM)). These ball materials were chosen for
their ability to form efficient triboelectric pair with FEP [33–37] and to 3. Working principle of the NCG-TENG
be commercially available in spherical form in order to propose an
easy-to-make TENG. The material density and triboelectric charge The NCG-TENG is a two-stage rotational energy harvester as depic
density (TECD) [35] of various practical materials are compared in ted in Fig. 1c. The primary stage is the main TENG where small balls are
Table S3, in terms of power density normalized to frequency. At the free to move into a closed rotating pipe (Fig. S2a). The secondary stage is
opposite side, the container was sealed with a resin part printed with a a circular capacitive coupler based on non-contact coupling between
stereolithography (SLA) 3D printer (Formlabs Form 3 +). The electrodes interdigitated pairs of electrodes (Fig. S3).
were connected to a rectifier circuit made with low-leakage diodes
(Nexperia BAV199) and a high-voltage capacitor (Kyocera 3.1. Operation of the primary stage
2225WC103KAT1A) for charge storage. The circuit is flexible and made
on a Cu-PI 100-µm thick film. The electrical connections between the The generation of electrical power in the primary stage is based on
TENG electrodes and the rectifier were done with wires noticeable in the successive ball movement inside the rotating pipe, tribo
Fig. S2a. electrification process and induction of charges at the electrical termi
nals. In order to properly explain the conversion of mechanical into
2.2. Fabrication and integration of the capacitive coupler as the electrical power, we present in the following section the physical
secondary stage modeling of the primary stage, and then derive the Norton equivalent
model of the triboelectric generator.
The PMMA pipe was cut after the rectifier in order to integrate a
rotor disk (diameter = 108 mm) sliced in a 5-mm thick PMMA sheet by a 3.2. Triboelectrification phase
laser cutting machine (Universal PLS6.75). A complementary pair of
electrodes (EP1 and EP2) was laser-cut (LPKF ProtoLaser S4) in a copper The movement of the balls inside the pipe induces multiple contacts
tape (thickness = 66 µm, 3M Company) to be bounded on the PMMA between balls surface and polymer film. This induces a tribo
surface (Fig. S3). Each electrode consists of (nd) digits patterned in pe electrification process, implying the increase of positive charges at the
riodic angular sectors (αT), the angle between two digits of the same ball surface, and negative charges at the polymer film surface. After a
electrode. Three capacitive couplers have been manufactured with 8, 16 certain time, the amount of generated surface charges reaches a satu
and 32 digits as shown in Fig. S4. A void sector was added between each ration value. The total amount of positive charges generated at balls
pair of electrodes to avoid air breakdown sparks. In order to isolate the surface (QB), is equal to the total amount of negative charges at polymer
electrodes, a 90-mm diameter PI film (thickness = 13 µm, Goodfellow film surface (QF). One can then deduce the relationship between space
Ltd.) was stuck on it with spray adhesive. The stator disk (diameter = charge densities at balls and polymer film surfaces (σB) and (σF),
124 mm) was also laser-cut in 5-mm thick PMMA sheet. The same respectively. Once a stable triboelectrification state has been reached,
interdigitated electrodes were bounded on the PMMA surface and a 90- stable periodic waveforms are generated at the contacts, which is
mm diameter PI film (thickness = 13 µm) was also stuck on the ES1 and studied in the next section.
ES2 electrodes. The stator disk can be adjusted toward the rotor disk
thanks to a micro-positioner (Fig. 1b) to control the gap between them 3.3. 2D electrostatic modeling of the primary stage
and a specific design was done to ensure a very precise gap spacing
between the two disks. To guarantee the concentricity, the stator disk is In order to explain and quantify the mechanical to electrical power
free to move along 4 screws thanks to 4 polyethylene (PE) foams (in transduction, the physical modeling of charge transfer induced at the
ternal diameter = 3 mm, external diameter = 10 mm, height = 10 mm, electrodes by the moving ball is investigated. First, the physical 2D
RS Pro). For the experiments, the NCG-TENG is connected to 3D printed rotating model developed within COMSOL Multiphysics®, comprises a
supports (Raise3D E2) by means of bearings (SKF 6001-Z) as explained layer of 12 balls that rotate around the pipe as depicted in Fig. 2a, where
in Fig. S5. The rotation was driven by a DC motor (Fig. 1b) through a 3D the rotation direction is represented by an arrow. The number of balls is
printed coupler, and a contact tachymeter (RS Pro) was used to monitor selected to sufficiently cover the surface of an electrode. As presented
the TENG rotating speed (N) noticeable in Movie S1. before, an equal amount of fixed charges are present at ball and polymer
Supplementary material related to this article can be found online at film surfaces and the system is globally neutral. The displacement of
doi:10.1016/j.nanoen.2023.108861. balls inside the pipe will thus result in the modification of local charge
balance and electrostatic forces acting on the device terminals. The
2.3. Measurement setup initial model considers electrostatic laws based on 2D Poisson equation
solution, dielectric properties of materials, surface charge state of ma
The methodology to extract Norton equivalent model parameters terials, and boundary conditions at the terminals. The amount on charge
from the primary stage is presented in Fig. S6, where the primary (Q1) and (Q2) is induced at the electrodes by moving charge QB. In this
capacitor voltage (VCP) across the primary capacitor (CP) is acquired short-circuit condition, local neutrality can be satisfied near the balls
using a source measure unit (SMU) (Keithley 2636) through a resistance through electron transfer from one electrode to the other. As a result, the
ladder configuration. The coupling capacitors (CC) are monitored by a ball layer induces an equivalent negative charge in the nearest electrode.
capacitance meter (Agilent 4288 A). The root mean square (RMS) The rotating movement of balls inside the pipe thus induces an alter
voltage across resistive loads was acquired using a digital oscilloscope native movement of charge between the electrodes, thus resulting in a
(LeCroy waveRunner 6100 A). The secondary short-circuit current was periodic and complementary time variation of Q1 and Q2 (Fig. 2b). In a
3
E. Kharbouche et al. Nano Energy 117 (2023) 108861
Fig. 2. Models and FEM simulation results of the primary stage. a) Physical model of electrostatics with a layer of 12 rotating balls within the pipe (2D rotating
model). b) Simulated charge in short-circuit for the three models: 2D rotating, 2D vertical and 1D analytic. c) Simulated voltage in open-circuit for the three models.
d) Physical model using a single ball with concentrated charge in vertical movement (2D vertical model). This model is equivalent to the 2D rotating model when a
homogeneous balls layer moving around the pipe.
second simulation, the open-circuit operation is investigated. In this equations system, allowing calculation of output current and voltage
case, no current transfer can occur between the electrodes and the ball in any load condition. By simulating the short and open conditions
layer induces a voltage variation in the nearest contact. Connecting across the contacts of the 1D cross-section using COMSOL Multi
electrode 1 to a zero reference potential, an alternating sinewave in physics®, the charge and voltage values shown in Fig. 2b-c are obtained.
electrode 2 is observed in Fig. 2c. The parameters of the simulation were The ball material carries the full positive charge at balls QB. Each side of
calibrated on actual measurements, in order to obtain a maximum the polymer film has a surface charge QF in much lower concentration
transferred charge of 35 nC, and an open-circuit voltage of +/− 4500 V. since it is distributed over a large area around the inner pipe. Calculation
The 2D vertical model considers a unique ball, as shown in Fig. 2d, of the electrostatic equilibrium of this system allows quantification
concentrating the full charge of the balls layer simulated in the 2D induced charge balance at the contacts Q1 and Q2. The electrical
rotating model. In this model, the ball moves along the vertical Y-axis as boundary condition V21 at the contacts is set by the external circuit,
depicted in Fig. 2d. The vertical position of the ball represents exactly either open, short or load. We consider here a fringing capacitance (Cfg),
the actual vertical position of the rotating ball layer. One can observe accounting for 2D spreading of the field lines, and the influence of pipe
that this simplification of the model allows a fair approximation of the plastic material dielectric properties. The study in load condition is
actual system (Fig. 2b-c). We further confirm this by calculating the made with the capacitance CP fed by a full rectifier.
impact of a horizontal displacement of the ball along the X-axis, which Connecting contact 1 to reference voltage, we can express V21
has a very limited impact on both transferred charge and voltage fluc voltage across the two metal contacts with all the parameters described
tuation. Thus, deriving a simplified 1D analytic electrostatic model of in Fig. S7a.
the primary stage generator in order to develop a simplified Norton [ ]
1 YB
model seems to be a relevant approach. V21 (YB ) = ⋅ QB ⋅ + (Q2 + QF + QB )⋅aT − QB ⋅aF (1)
A.ε0 εair
3.4. Derivation of Norton equivalent model of the primary stage generator In expression (1), all parameters are fixed, except voltage V21, ter
minal charge Q2 and the average vertical ball layer position (YB). The
From the physical electrostatic model, we want to derive the Norton rotation movement will induce a sinusoidal time variation of YB, and
equivalent model in order to relate the potential output power of the subsequent time variation of (Q2, V21) couple. Calculating V21 and
device to the physical parameters. Fig. S7a shows a cross-section of the subtracting it at extreme values of YB when R is the internal radius of the
device, obtained from the 2D vertical model, which is used to elaborate pipe, we can derive the amplitude of intrinsic open-circuit voltage
the governing equations of the 1D analytic model. Solving Poisson (VOCi):
equation in this system results in Charge-Electric and Field-Voltage
4
E. Kharbouche et al. Nano Energy 117 (2023) 108861
Fig. 3. Output current multiplication mechanism through the capacitive coupler of the secondary stage. a) 2D schematic describing the characteristic steps of the
capacitive coupling mechanism between EP and ES electrodes for a periodic angular sector (αT). b) Evolution of the simulated charge on the ES1 electrode (QES1) as a
function of the angle when the secondary electrodes are short-circuited and the gap between the rotor and stator is 100 µm. (I), (II), (III) and (IV) correspond to the
characteristic steps during an αT rotation. c) Output current model resulting in QES1 fit derivation over time when the air gap is 100 µm and the rotating speed (N) is
125 RPM. If the number of digits (nd) is doubled, the peak value and the frequency are also doubled. d) Evolution of the electrical potential on the ES1 electrode (VES1)
over time when the secondary electrodes are in open-circuit condition and ES2 electrode is grounded. While the frequency doubles with the number of digits, the peak
value remains constant.
5
E. Kharbouche et al. Nano Energy 117 (2023) 108861
stator. This is shown in Fig. 3b, where we consider the maximum function of time. It appears that a sinusoidal model can fit the simulated
transferred charge of 35 nC in the primary stage to be QEP1 and an 100- QES1 (Fig. S10a) as follows,
µm air gap between the rotor and the stator. This result is extracted from
N.nd
the 3D FEM simulation of the simplified structure illustrated in Fig. S9a. QES1 (t) = Qmax ⋅sin(2.π.fd .(t − t0 ) ) where fd = (7)
60
The smaller the air gap, the closer the QES1 maximum (Qmax) tends to
QEP1, due to better coupling between the primary and the secondary where Qmax is the maximum charge induced in ES1 electrode when EP
electrodes as observed in Fig. S9b. First, EP1 is aligned to ES1 and QES1 is and ES electrodes are aligned at the time t0 and nd is the number of digits
minimum (I). When EP1 moves to ES2 (II), electrons are injected into ES2 of an electrode. We deduce the current in ES1 electrode (IES1) from (7) by
to counterbalance the positive charge QEP1. Then, the charge QES1 is derivation over time,
maximum when EP2 is aligned to ES1 electrode (III). After the step (III),
the negative charge on EP2 pushes back electrons from ES2 to ES1 elec IES1 (t) =
dQES1 (t)
= Imax ⋅cos(2.π .fd .(t − t0 ) ) where Imax = 2.π.fd .Qmax
trode (IV). dt
At a constant rotating speed N, the angle modification is proportional (8)
to the time. Therefore, we can consider the simulated charge QES1 as a For one period, the proposed model for IES1 is consistent with the
Fig. 4. Influence of the mechanical parameters on the electrical performances of the primary stage, where the dimensions of an electrode are 20.4 mm × 41 mm. a-b)
For each material, 266 balls of 2-mm diameter are investigated at 125 RPM. When a FEP film is used, the ZrO2 balls achieve the highest power of 91.5 µW. c) Effect of
ZrO2 balls diameter on the electrical characteristic at 125 RPM. The more the ball diameter is small, the more the output power is high. A maximum power of 152 µW
is reached for 1.587-mm diameter balls. d) The optimal number of 1.587-mm diameter balls in ZrO2 is searched at 125 RPM. The maximum power is attained when
the number of balls covers the surface of an electrode. e-f) Several rotating speeds are tested on 266 2-mm diameter ZrO2 balls. In the same operating mode, the
output power grows linearly over the rotating speed.
6
E. Kharbouche et al. Nano Energy 117 (2023) 108861
short-circuit current measured on the generator as demonstrated in 4.1. Influence of ball and triboelectric layer materials
Fig. S10b. Finally, the sinusoidal model can also fit the electrical po
tential measured on the ES1 electrode (VES1) when the ES2 electrode is First, Fig. 4a illustrates the investigation of three different materials
grounded (Fig. S10c) and it can be expressed as, of 2-mm diameter balls at 125 RPM when a FEP film is used. ZrO2 balls
demonstrate the highest output power of 91.5 µW while POM balls show
VCP
VES1 (t) = Vmax ⋅sin(2.π .fd .(t − t0 ) ) where Vmax = (9) the lowest output performance of 12.3 µW. According to the triboelec
2
tric series [35–37], glass outlines the biggest TECD of 38.63 µC m− 2
In summary, both current and voltage frequencies are proportional whereas for ZrO2 is 0.09 µC m− 2. However, ZrO2 has the highest density
to the rotating speed N and the number of electrode digits nd. The (6000 kg m− 3) among the three materials tested. It can be seen in
maximum current (Imax) can be multiplied by nd as demonstrated in Movies S6–7, for similar rotation speed, glass and POM balls stay stuck
Fig. 3c. In contrast, the maximum voltage (Vmax) only depends on the on the FEP film due to the electrostatic force (FE) and the centrifugal
voltage between the primary electrodes (EP1 and EP2). Thus, Vmax is force (FC) influence, but this is not observed with ZrO2 balls. For them,
constant when VCP tends to the primary open-circuit voltage VOCP and the weight along the Y-axis (FGy) is sufficient to overcome FE and FC
this is illustrated in Fig. 3d. forces. Consequently, all ZrO2 balls follow at least one of the two modes
of operation (sliding or falling mode presented in Fig. S8), resulting in
4. Results and discussions better performance. Another explanation for the performance difference
between ZrO2 balls and others could be the operating modes. Indeed, the
In order to characterize the performance of the NCG-TENG, the falling balls mode is not observed in Movies S6–7 for glass and POM
primary stage was studied before being connected to the secondary stage materials. In contrast, a mix of the two operating modes is noticeable in
(Fig. 4). Different mechanical parameters (ball and triboelectric mate Movie S3 when ZrO2 balls are used. This visual observation is confirmed
rials, ball diameter, number of balls and rotating speed) are investigated by I-V curves in Fig. 4a as for glass and POM balls follow a linear trend,
thanks to I-V curves. For each experiment, the maximum power is whereas the I-V curve for ZrO2 balls is clearly non-linear.
extracted from these curves (Fig. S6) and the number of balls is chosen to Supplementary material related to this article can be found online at
cover only one electrode surface. This means that more the balls diam doi:10.1016/j.nanoen.2023.108861.
eter is small, more the number of balls to maintain the same covered The selection of FEP material for the triboelectric layer was based on
surface is high. After connection to the secondary stage, N and the its high TECD of -249.42 µC m− 2 [37]. In opposition, a film material of
number of electrode digits nd are considered to highlight the strong PC with a low TECD of -24.4 µC m− 2 [37] was tested and the results
frequency dependence of the TENG output power (Fig. 5). were compared to the FEP film I-V curve in Fig. 4b. An output power of
Fig. 5. Output performance of the NCG-TENG after the secondary stage when 1.587-mm diameter ZrO2 balls are used in the primary stage. a-b) Effect of the
electrode digits number of the capacitive coupler on the output performance at 125 RPM. The electrical circuit (I) was used to investigate the load impedance
matching. At the same time, the short-circuit current (VOUT = 0 V) is directly proportional to the number of electrode digits (nd) and the open-circuit voltage (IOUT =
0 µA) is constant. c) Effect of the rotating speed on the output performance for an 8-digit capacitive coupler. d) Investigation of the TENG performance with 32-elec
trode digits for different capacitor sizes. At 500 RPM, a power of 20.54 mW is produced by the TENG when a 1-µF capacitor (CS) is connected to the circuit (II).
7
E. Kharbouche et al. Nano Energy 117 (2023) 108861
79.5 µW is observed, which closely resembles that obtained with an FEP 4.5. Average power density of the primary stage
film. One possible explanation is a positive charge saturation on the balls
surface prior to that of the film. Finally, the material of the triboelectric The literature review in Table S1 compares performances of free-
layer becomes less crucial when ZrO2 balls are used, as its TECD is standing balls TENGs to those of the primary stage when the surface
relatively weak. In all the following experiments, ZrO2 and FEP are reference is 1-electrode area (equal to 20.4 mm × 41 mm). The rotating
respectively adopted as the ball and polymer film materials. ball-mill-inspired TENG (BMI-TENG) [19] appears to be as efficient as
the primary stage, but its operating frequency is more than 3 times
4.2. Influence of the ball diameter higher. Moreover, the literature tends to focus on peak power density,
which is always far superior to average power density. The average
Three different diameters are examined at 125 RPM and results are power density of the primary stage is 181.95 mW m− 2 from an average
presented in Fig. 4c. Best output power performance is recorded with the power of 152 µW at 125 RPM. When the reference is changed to a vol
smallest balls (1.587-mm diameter). Under this condition, the TENG ume reference as presented in Table S2, the average power density of the
produces 152 µW whereas only 48.2 µW is produced with 3-mm diam primary stage is 6.98 W m− 3. Here, the volume reference is the inner
eter balls. It has been observed in Movies S3, S8, and S9 that the smaller cylinder (diameter = 26 mm, height = 41 mm) where balls are placed.
the ball diameter is, the further the ball climbs up the pipe surface, due In our work, the output power density of the primary stage exceeds the
to their lightness. The 1.587-mm diameter balls can climb until the top one reported for the BMI-TENG by at least a factor of two.
of the pipe and then fall along the entire pipe diameter, which means
that balls velocity increases and consequently enhancing the current 4.6. Output performance of the NCG-TENG
flowing between the electrodes. The primary short-circuit current (ISCP)
is directly equal to the current ICP when VCP is zero. In Fig. 4c, ISCP is The primary stage is now coupled to the capacitive coupler to form
equal to 164 nA for 1.587-mm diameter balls. The primary open-circuit the NCG-TENG. In the primary stage, 422 ZrO2 balls of 1.587-mm
voltage (VOCP) is limited to 2000 V by current leakages in the rectifier diameter are placed as they gave the best performances. The capaci
and by the measurement range of the SMU. As illustrated in Fig. S6, VOCP tive coupler gap is tuned in order to maximize and set the values CC1 and
can be extracted from the intersection of the linear fit of the data and the CC2 (Fig. 3a) at 100 pF during the experiments. First, a resistive load (RL)
X-axis. Finally, VOCP is extrapolated to 4509 V. is connected between the ES1 and ES2 electrodes to investigate the load
Supplementary material related to this article can be found online at impedance matching. The results presented in Fig. 5a, confirm that the
doi:10.1016/j.nanoen.2023.108861. output power (POUT) is proportional to the number of electrode digits
(nd). Using an 8-digit capacitive coupler at 125 RPM, POUT is 1.24 mW
compared to 152 µW produced by the primary stage. When N is con
4.3. Influence of the number of balls stant, nd is a multiplying factor for the output current (IOUT) as illustrated
in Fig. 5b, where doubling nd doubles IOUT. This experimental observa
In all experiments, the number of balls was determined to ensure full tion confirms the theoretical analysis carried out in Fig. 3. Furthermore,
coverage of the electrode surface, thereby maximizing the amount of the output voltage (VOUT) remains constant even if nd increases, since
charges QB rotating in the primary stage. However, an excessive number VOUT only depends on the voltage VCP. Similarly, POUT raises linearly
of balls may cause them to extend beyond the electrode surface, leading with N due to the relationship between the current and the frequency
to undesirable charge transfers. In addition, for optimum efficiency, the demonstrated in (8). This is confirmed by the experiment in Fig. 5c. At
balls must form a unique layer, enabling free balls movements in either 600 RPM, the NCG-TENG can deliver an average power of 6.39 mW,
the sliding or falling mode. A significant number of balls can result in which is 5 times higher than the power recorded at 125 RPM. Finally, RL
multiple layers stacking, leading to compromise ball performance. To is substituted by a full-bridge rectifier to investigate the power across a
validate this hypothesis, experimental tests are conducted using 1.587- storage capacitor (CS) as highlighted in Fig. 5d. In this configuration, the
mm diameter balls, which have demonstrated the highest output power is the derivative of the energy stored in CS with respect to time, as
power. The outcomes, shown in Fig. 4d, confirm that the optimal explained in Fig. S11. Three different sizes of CS were tested and a peak
number of balls is directly dependent of the electrode surface. In our power of 20.54 mW at 500 RPM was observed for 1 µF, which is the
specific case, employing 422 balls of 1.587-mm diameter covered an matching capacitance. In summary, a 32-digit capacitive coupler is able
electrode and resulted in the most favorable performances. to boost the NCG-TENG output power density to 24.45 W m− 2 when the
reference is 1-electrode area of the primary stage (Table S1), which is to
4.4. Influence of the rotating speed (N) our knowledge the best-published performance. A stability test was
conducted to verify the sustainable energy collection capability of the
The linearity between N and the output current has been highlighted NCG-TENG in Fig. S12. After running for 36 h at 300 RPM (equivalent to
in the analytic model. Therefore, several N are investigated with 2-mm 648,000 cycles), more than 90% of the peak current is sustained across a
diameter ZrO2 balls in Fig. 4e-f. This ball diameter has been chosen to 10 MΩ resistive load. This result is satisfactory, especially considering
offer more freedom in terms of rotating speed. Several tests have been the reliability tests conducted by Qin et al. [19], where 80.13% of the
carried out with 1.587-mm diameter balls, but they tend to stick directly transferred charge was preserved after 5000 cycles.
to the pipe surface as the centrifugal force FC increases under the in
fluence of the rotating speed. As a result, the proportionality between N 4.7. Application scenario
and ISCP (VCP = 0 V) is observed in Fig. 4e, as well as a constant VOCP (ICP
= 0 nA). As explained in the analytic model, the voltage VOCP depends This study emphasizes the great potential of the NCG-TENG in
on the charge QB, which saturates at a certain time. Finally, the output achieving high power density. Hence, using wind energy to power an air
power (PCP) shown in Fig. 4f grows linearly according to N until 300 quality monitoring application is proposed. A custom-made plenum
RPM. The observed operating mode is a mix between sliding and falling chamber substitutes the DC motor to generate a controlled airflow as
balls mode. Above 300 RPM, a change in operating mode to a stationary shown in Movie S10. In this scenario, the propeller is driven by an
mode is observed (Movie S5). In this situation, balls accumulated in the airflow of 12 m s− 1, leading to a rotating speed of the pipe of 280 RPM
bottom of the pipe and electrodes turn around them. The stationary and charging a 1-µF capacitor to 415 V in 60 s. Following the method
mode goes on up to 650 RPM when the centrifugal force sticks the balls ology outlined in Fig. S11, a power of 6 mW is generated.
on the pipe surface (Movie S4). In this situation, PCP drops to zero since Supplementary material related to this article can be found online at
the balls do not induce charge movement in the electrodes. doi:10.1016/j.nanoen.2023.108861.
8
E. Kharbouche et al. Nano Energy 117 (2023) 108861
Subsequently, the feasibility of powering the air quality monitoring administration, Funding acquisition, Writing – review & editing.
application (as illustrated in Fig. S13a) is explored. Based on power
measurement (Fig. S13b) conducted on the object, the average power
required for a measurement and data transmission operation is Declaration of Competing Interest
13.85 mW. This power (PON) is consumed during the time (TON) lasting
for 6.5 s. After an operation, the system enters in sleep mode and con The authors declare that they have no known competing financial
sumes an average power (POFF) of 0.0036 mW during the time (TOFF). interests or personal relationships that could have appeared to influence
This time can be adjusted to reduce the average power consumption the work reported in this paper.
(PMEAN) expressed as,
Data availability
PON • TON + POFF • TOFF
PMEAN = (10)
TON + TOFF Data will be made available on request.
To achieve power balance, PMEAN must be set equal to 6 mW. In this
case, TOFF is calculated as follows, Acknowledgments
PON − PMEAN
TOFF = TON • (11) This work was supported by IPCEI Flexeo project funded by the
PMEAN − POFF
French Ministry of Industry. Part of this work was done with the support
Finally, the connected object should remain in sleep mode for 8.5 s to of ID-Fab (Prototyping platform: project funded by the European
attain PMEAN for a subsequent operation. However, this outcome can Regional Development Fund, the French state and local authorities).
only be considered when using an impedance matching circuit to extract
the maximum power available at 50–70% of the open-circuit voltage as Appendix A. Supporting information
illustrated in Fig. 5d. Zhang et al. [38] proposed a high voltage condi
tioning circuit based on a switching strategy, achieving an efficiency of Supplementary data associated with this article can be found in the
35%. Consequently, TOFF is increased from 8.5 s to 36.4 s. With an online version at doi:10.1016/j.nanoen.2023.108861.
operation conducted every minute, it would be possible to accomplish
1440 operations per day, which is suitable for an air quality monitoring References
application.
[1] A.J. Cohen, M. Brauer, R. Burnett, H.R. Anderson, J. Frostad, K. Estep,
5. Conclusion K. Balakrishnan, B. Brunekreef, L. Dandona, R. Dandona, V. Feigin, G. Freedman,
B. Hubbell, A. Jobling, H. Kan, L. Knibbs, Y. Liu, R. Martin, L. Morawska, C.
A. Pope, H. Shin, K. Straif, G. Shaddick, M. Thomas, R. van Dingenen, A. van
Energy resources and energy harvesting play a vital role in powering Donkelaar, T. Vos, C.J.L. Murray, M.H. Forouzanfar, Estimates and 25-year trends
IoT devices, extending device lifespan, enabling connectivity and data of the global burden of disease attributable to ambient air pollution: an analysis of
data from the Global Burden of Diseases Study 2015, Lancet 389 (2017)
transmission, facilitating real-time monitoring and control, promoting 1907–1918, https://doi.org/10.1016/S0140-6736(17)30505-6.
energy efficiency and sustainability in IoT applications. For instance, in [2] World Health Organization, Ambient Air Pollution: A Global Assessment of
order to operate energy-intensive air quality sensor nodes, the devel Exposure and Burden of Disease, World Health Organization, 2016. 〈https://apps.
who.int/iris/handle/10665/250141〉 (Accessed April 3, 2023).
opment of high-performance energy harvesters is requested. In this
[3] D. Banister, Unsustainable Transport: City Transport in the New Century,
study, the capabilities of a non-contact triboelectric nanogenerator to Routledge, London, 2005, https://doi.org/10.4324/9780203003886.
harvest the gravitational energy of sliding and falling balls are pre [4] D. Kandris, C. Nakas, D. Vomvas, G. Koulouras, Applications of wireless sensor
networks: an up-to-date survey, Appl. Syst. Innov. 3 (2020) 14, https://doi.org/
sented. We propose a new methodology to extract Norton model pa
10.3390/asi3010014.
rameters from a rotating TENG. Moreover, a simplified analytic model [5] IoT Devices Installed Base Worldwide 2015–2025, Statista, n.d. 〈https://www.stat
has been suggested to relate these electrical parameters to the physical ista.com/statistics/471264/iot-number-of-connected-devices-worldwide/〉
parameters of the TENG. Finally, we have explained and demonstrated (Accessed April 3, 2023).
[6] P. Das, S. Ghosh, S. Chatterjee, S. De, Energy harvesting-enabled 5G advanced air
from a 500 RPM rotating speed that the peak power density of our two- pollution monitoring device, in: Proceedings of the IEEE 3rd 5G World Forum
stage TENG is 24.45 W m− 2. Based on these performances, an applica 5GWF, 2020, pp. 218–223. 〈https://doi.org/10.1109/5GWF49715.2020.9221
tion scenario has been proposed for air quality monitoring. This state-of- 330〉.
[7] F.-R. Fan, Z.-Q. Tian, Z. Lin Wang, Flexible triboelectric generator, Nano Energy 1
the-art power density value takes benefit from the capacitive coupling (2012) 328–334, https://doi.org/10.1016/j.nanoen.2012.01.004.
between a rotor and a stator, which increases the output power by means [8] Z.L. Wang, Triboelectric nanogenerators as new energy technology for self-
of thin interdigitated electrodes. By using advanced manufacturing powered systems and as active mechanical and chemical sensors, ACS Nano 7
(2013) 9533–9557, https://doi.org/10.1021/nn404614z.
techniques to increase the number of digits and the rotating speed range, [9] E. Kharbouche, W.L. Ferreira, D. Garcia, F. Bernier, S. Blayac, Highly integrated
the reported performance could certainly be improved. With optimized planar airflow energy harvester for self-powered air quality monitoring, in:
mechanical parameters, triboelectric free-standing balls in rotation are a Proceedings of the 2022 IEEE International Smart Cities Conference ISC2, 2022,
pp. 1–5. 〈https://doi.org/10.1109/ISC255366.2022.9921813〉.
promising concept for increasing the durability of TENGs and main
[10] D.A. Howey, A. Bansal, A.S. Holmes, Design and performance of a centimetre-scale
taining high performance. Further investigations are ongoing to provide shrouded wind turbine for energy harvesting, Smart Mater. Struct. 20 (2011),
a power management circuit to adjust the load impedance for maximum 085021, https://doi.org/10.1088/0964-1726/20/8/085021.
[11] X. Wu, D.-W. Lee, An electromagnetic energy harvesting device based on high
power extraction. A real-world test is also being considered to validate
efficiency windmill structure for wireless forest fire monitoring application, Sens.
the concept for an air quality application. Actuators Phys. 219 (2014) 73–79, https://doi.org/10.1016/j.sna.2014.09.002.
[12] F.-R. Liu, W.-M. Zhang, Z.-K. Peng, G. Meng, Fork-shaped bluff body for enhancing
the performance of galloping-based wind energy harvester, Energy 183 (2019)
CRediT authorship contribution statement 92–105, https://doi.org/10.1016/j.energy.2019.06.044.
[13] J. Zhang, Z. Fang, C. Shu, J. Zhang, Q. Zhang, C. Li, A rotational piezoelectric
Elias Kharbouche: Conceptualization, Methodology, Validation, energy harvester for efficient wind energy harvesting, Sens. Actuators Phys. 262
(2017) 123–129, https://doi.org/10.1016/j.sna.2017.05.027.
Formal analysis, Investigation, Data curation, Writing – original draft,
[14] X. Wang, S. Niu, Y. Yin, F. Yi, Z. You, Z.L. Wang, Triboelectric nanogenerator based
Visualization, Writing – review & editing. Marc Ramuz: Validation, on fully enclosed rolling spherical structure for harvesting low-frequency water
Investigation, Supervision, Resources, Writing – review & editing. wave energy, Adv. Energy Mater. 5 (2015) 1501467, https://doi.org/10.1002/
François Bernier: Conceptualization, Methodology. Cyril Calmes: aenm.201501467.
[15] H. Wang, L. Xu, Z. Wang, Advances of high-performance triboelectric
Conceptualization, Methodology, Visualization. Sylvain Blayac: nanogenerators for blue energy harvesting, Nanoenergy Adv. 1 (2021) 32–57,
Conceptualization, Formal analysis, Supervision, Resources, Project https://doi.org/10.3390/nanoenergyadv1010003.
9
E. Kharbouche et al. Nano Energy 117 (2023) 108861
[16] L. Xu, T. Jiang, P. Lin, J.J. Shao, C. He, W. Zhong, X.Y. Chen, Z.L. Wang, Coupled Marc Ramuz received his PhD in 2010 from the Centre Suisse
triboelectric nanogenerator networks for efficient water wave energy harvesting, d′Electronique et de Microtechnique (CSEM) in Zurich. Then he
ACS Nano 12 (2018) 1849–1858, https://doi.org/10.1021/acsnano.7b08674. joined Zhenan Bao group in Stanford for two years post
[17] Y. Pang, S. Chen, Y. Chu, Z.L. Wang, C. Cao, Matryoshka-inspired hierarchically doctoral. In 2012, he undertook a postdoctoral fellowship in
structured triboelectric nanogenerators for wave energy harvesting, Nano Energy the Bioelectronics department of the School of Mines Saint
66 (2019), 104131, https://doi.org/10.1016/j.nanoen.2019.104131. Etienne under the supervision of Róisín Owens. He is currently
[18] Y. Wang, A. Matin Nazar, J. Wang, K. Xia, D. Wang, X. Ji, P. Jiao, Rolling spherical an associate professor in Flexible Electronics Department at the
triboelectric nanogenerators (RS-TENG) under low-frequency ocean wave action, Campus Georges Charpak Provence of Mines Saint Etienne. He
J. Mar. Sci. Eng. 10 (2022) 5, https://doi.org/10.3390/jmse10010005. is co-author of more than 30 peer-reviewed scientific publica
[19] Q. Qin, X. Cao, N. Wang, Ball-mill-inspired durable triboelectric nanogenerator for tions, including 10 as first author. His research interests are
wind energy collecting and speed monitoring, Nanomaterials 13 (2023) 939, related to multi-stack sensors for artificial skin, soft/stretch
https://doi.org/10.3390/nano13050939. able electronic and investigation of nanomaterials for energy
[20] P. Cheng, H. Guo, Z. Wen, C. Zhang, X. Yin, X. Li, D. Liu, W. Song, X. Sun, J. Wang, conversion and storage, including nanogenerators.
Z.L. Wang, Largely enhanced triboelectric nanogenerator for efficient harvesting of
water wave energy by soft contacted structure, Nano Energy 57 (2019) 432–439,
https://doi.org/10.1016/j.nanoen.2018.12.054.
François Bernier is currently responsible of the prototyping
[21] X. Yang, L. Xu, P. Lin, W. Zhong, Y. Bai, J. Luo, J. Chen, Z.L. Wang, Macroscopic
platform ID-Fab in the Campus Georges Charpak Provence of
self-assembly network of encapsulated high-performance triboelectric
Mines Saint Etienne in Gardanne (France). He provides tech
nanogenerators for water wave energy harvesting, Nano Energy 60 (2019)
nical support for different academic and research projects. His
404–412, https://doi.org/10.1016/j.nanoen.2019.03.054.
areas of interests are rapid prototyping process, embedded
[22] M. Xu, T. Zhao, C. Wang, S.L. Zhang, Z. Li, X. Pan, Z.L. Wang, High power density
systems development and electronic platform design. In the
tower-like triboelectric nanogenerator for harvesting arbitrary directional water
past, he has developed various systems related to medical and
wave energy, ACS Nano 13 (2019) 1932–1939, https://doi.org/10.1021/
astronomical imaging systems.
acsnano.8b08274.
[23] L. Liu, X. Yang, L. Zhao, H. Hong, H. Cui, J. Duan, Q. Yang, Q. Tang, Nodding duck
structure multi-track directional freestanding triboelectric nanogenerator toward
low-frequency ocean wave energy harvesting, ACS Nano 15 (2021) 9412–9421,
https://doi.org/10.1021/acsnano.1c00345.
[24] J. Lee, H. Yong, Y.I. Choi, J. Ryu, S. Lee, Stackable disk-shaped triboelectric
nanogenerator to generate energy from omnidirectional wind, Int. J. Precis. Eng.
Manuf. Green Technol. 9 (2022) 557–565, https://doi.org/10.1007/s40684-021-
00340-9. Cyril Calmes completed his studies at University Montpellier II
[25] H. Yong, J. Chung, D. Choi, D. Jung, M. Cho, S. Lee, Highly reliable wind-rolling in material sciences and received his PhD in 2003 at the
triboelectric nanogenerator operating in a wide wind speed range, Sci. Rep. 6 Institute of Fundamental Electronics at the University of Paris
(2016) 33977, https://doi.org/10.1038/srep33977. XI – ORSAY (France). He joined the Provence Microelectronics
[26] L. Xu, T.Z. Bu, X.D. Yang, C. Zhang, Z.L. Wang, Ultrahigh charge density realized Center, Mines Saint-Etienne, Gardanne, France, in 2008. He is
by charge pumping at ambient conditions for triboelectric nanogenerators, Nano currently working as a research engineer at the Flexible Elec
Energy 49 (2018) 625–633, https://doi.org/10.1016/j.nanoen.2018.05.011. tronics Department at the Campus Georges Charpak Provence
[27] H. Wang, L. Xu, Y. Bai, Z.L. Wang, Pumping up the charge density of a triboelectric of Mines Saint Etienne. He is also responsible for end-of-study
nanogenerator by charge-shuttling, Nat. Commun. 11 (2020) 4203, https://doi. projects in electronics and energy for engineering students on
org/10.1038/s41467-020-17891-1. the ISMIN course at Mines Saint Etienne.
[28] A.N. Ravichandran, C. Calmes, J.R. Serres, M. Ramuz, S. Blayac, Compact and high
performance wind actuated venturi triboelectric energy harvester, Nano Energy 62
(2019) 449–457, https://doi.org/10.1016/j.nanoen.2019.05.053.
[29] T. Jiang, X. Chen, K. Yang, C. Han, W. Tang, Z.L. Wang, Theoretical study on
rotary-sliding disk triboelectric nanogenerators in contact and non-contact modes,
Nano Res. 9 (2016) 1057–1070, https://doi.org/10.1007/s12274-016-0997-x.
[30] Y. Bai, L. Xu, S. Lin, J. Luo, H. Qin, K. Han, Z.L. Wang, Charge pumping strategy for Sylvain Blayac was born in 1972. He received the M.S. degree
rotation and sliding type triboelectric nanogenerators, Adv. Energy Mater. 10 in electronics from the University of Montpellier, Montpellier,
(2020) 2000605, https://doi.org/10.1002/aenm.202000605. France in 1997, and the Ph.D. degree in electronics in 2001. He
[31] A. Ghaderiaram, A. Bazrafshan, K. Firouzi, M. Kolahdouz, A multi-mode R-TENG then joined Alcatel, Paris, France, as a Device Physical and
for self-powered anemometer under IoT network, Nano Energy 87 (2021), 106170, Electrical Modeling Engineer for high-speed fiber-optic
https://doi.org/10.1016/j.nanoen.2021.106170. communication circuits. He joined the Provence Microelec
[32] A.N. Ravichandran, F. Depoutot, E. Kharbouche, M. Hamand, M. Ramuz, S. Blayac, tronics Center, Ecole de Mines de Saint-Etienne, Gardanne,
Demonstration of friction-based triboelectric nanogenerator and integration in a France, in 2003, where he became a Teacher of semiconductor
power-balanced fully autonomous system, Nano Energy 83 (2021), 105796, devices physics and conducted research on embedded stress
https://doi.org/10.1016/j.nanoen.2021.105796. sensors. His current research interests include flexible sensors
[33] J. Henniker, Triboelectricity in polymers, 474–474, Nature 196 (1962), https:// and energy harvesting devices. Professor Blayac authored and
doi.org/10.1038/196474a0. co-authored more than 50 communications in refereed inter
[34] A.F. Diaz, R.M. Felix-Navarro, A semi-quantitative tribo-electric series for national journals and conferences and holds 7 patents.
polymeric materials: the influence of chemical structure and properties,
J. Electrost. 62 (2004) 277–290, https://doi.org/10.1016/j.elstat.2004.05.005.
[35] H. Zou, Y. Zhang, L. Guo, P. Wang, X. He, G. Dai, H. Zheng, C. Chen, A.C. Wang,
C. Xu, Z.L. Wang, Quantifying the triboelectric series, Nat. Commun. 10 (2019) Elias Kharbouche is presently a PhD student in the Flexible
1427, https://doi.org/10.1038/s41467-019-09461-x. Electronics Department at the Campus Georges Charpak Pro
[36] H. Zou, L. Guo, H. Xue, Y. Zhang, X. Shen, X. Liu, P. Wang, X. He, G. Dai, P. Jiang, vence of Mines Saint Etienne in Gardanne, France. His research
H. Zheng, B. Zhang, C. Xu, Z.L. Wang, Quantifying and understanding the focuses on self-sufficient connected objects using triboelectric
triboelectric series of inorganic non-metallic materials, Nat. Commun. 11 (2020) nanogenerators. He holds an engineering degree in Micro
2093, https://doi.org/10.1038/s41467-020-15926-1. electronics and Telecommunications from Polytech Marseille,
[37] D. Liu, L. Zhou, S. Cui, Y. Gao, S. Li, Z. Zhao, Z. Yi, H. Zou, Y. Fan, J. Wang, Z. France, which he attained in 2016. Subsequently, he joined the
L. Wang, Standardized measurement of dielectric materials’ intrinsic triboelectric Secure Architectures and Systems Department in Gardanne as a
charge density through the suppression of air breakdown, Nat. Commun. 13 (2022) research engineer, where he specialized in developing test
6019, https://doi.org/10.1038/s41467-022-33766-z. benches to characterise the integrity of electronic components.
[38] H. Zhang, F. Marty, X. Xia, Y. Zi, T. Bourouina, D. Galayko, P. Basset, Employing a His research interests encompass large-scale environmental
MEMS plasma switch for conditioning high-voltage kinetic energy harvesters, Nat. monitoring, energy transduction and harvesting, power con
Commun. 11 (2020), https://doi.org/10.1038/s41467-020-17019-5. ditioning circuits and self-powered sensors.
10