This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).

The paper presents a numerical energy harvesting model for sensor nodes, SIVEH (Simulator I–V for EH), based on I–V hardware tracking. I–V tracking is demonstrated to be more accurate than traditional energy modeling techniques when some of the components present different power dissipation at either different operating voltages or drawn currents. SIVEH numerical computing allows fast simulation of long periods of time—days, weeks, months or years—using real solar radiation curves. Moreover, SIVEH modeling has been enhanced with sleep time rate dynamic adjustment, while seeking energy-neutral operation. This paper presents the model description, a functional verification and a critical comparison with the classic energy approach.

The use of wireless sensor networks (WSNs) is widely extended in monitoring applications. Sensor nodes are usually battery-powered devices with finite energy. However, in many applications, nodes are required to operate for very long periods of time. Energy harvesting (EH) can extend node lifetime indefinitely, although there is a limit on daily power consumption rates. For example, if solar energy is used to recharge batteries or supercapacitors, energy harvested varies along the day and will probably be unavailable at night or in bad weather conditions. Thus, for nodes that operate without penalty on performance during sunny periods, the same power consumption at night could imply a quick loss of energy without being able to recharge.

Energy-aware management policies—which lead to energy-saving operation—are essential to avoid power-off failures on devices or unforeseen disconnection gap time periods with no energy to communicate. Both hardware component consumption and the sensor duty cycle are important factors to determine super- or ultra-capacitor available energy. Although an orthogonal design of the physical layer and power management policies is possible, a cross-layer design has proven to be more efficient, especially on energy constrained systems, such as WSNs [1–4]. Therefore, an accurate energy model is necessary in a very early design stage, before prototyping, to help the efficient design of the physical layer. Moreover, this model should also consider dynamically adjusted packet reception and transmission rates, depending on available energy, to simulate energy-neutral operation (ENO).

Energy harvesting circuits include hardware components for which modeling requires I–V (Current and Voltage) tracking to be reliable. However, as far as we know, the most recent implementations of the classic energy approach in the field lack numerical models of components with an energy level-dependent behavior. As an example, a resistor-capacitor (RC) network presents a different behavior at different energy levels. To address this problem, this paper presents a mathematical I–V model for the design of EH-enabled WSN. The model is called SIVEH (Simulator I–V for EH) and brings about a flexible solution for simulating most EH current circuits.

Moreover, as solar radiation values depend on the final WSN location around the world and along the season, SIVEH implementation allows the use of long solar radiation vectors (
Jm2) that can be directly extracted from public databases. Although the prediction of future solar radiation values is out of the scope of this work, the use of registered values provides very realistic information about the variation of the buffered energy during hours, weeks, months or years, according to sunlight progression in a specific world location.

The SIVEH model increases precision compared with the energy approaches used in the literature. Moreover, as SPICE is assumed to be a reference in electronics simulation, SIVEH computational implementation has also been compared with it. SIVEH has shown a maximum derivate error of less than 1% with a significant improvement on simulation speed. The aim is to propose a new model of energy harvesting nodes for wireless sensor network simulation with the following features: high precision, high simulation speed and allowing for complex networks simulation. As a case study, the implementation of an algorithm to reach energy neutral operation (ENO) conditions is proposed.

This paper is organized as follows. Section 2 gives a general related literature overview. Section 3 presents the mathematical approach of SIVEH. Section 4 presents the SIVEH framework with a generic description of energy harvesting circuits and an example describing the implementation of a prototype. Section 5 evaluates SIVEH error and simulation time vs. both an energy approach and SPICE. Section 6 describes the implementation of the SIVEH model. Section 7 carries out a verification by comparing SIVEH outputs with a prototype. Section 8 presents a case study: The ENO condition in WSN nodes. Finally, Section 9 concludes this paper.

2.Motivation

Techniques to harvest energy are used to convert energy from ambient sources into electrical energy stored in buffers [5,6], which are usually supercapacitors or rechargeable batteries [7]. Some of the literature has focused on sunlight energy-harvesting circuits, although this research does not usually address networking aspects. Some authors include rechargeable batteries as the main buffer [8,9]. However, batteries have a finite number of charge-discharge cycles, meaning that the use of another solution based on supercapacitors is advisable [9–12]. Moreover, an EH circuit can include both solutions—supercapacitors and backup battery [13,14]. Solar EH circuits currently already incorporate Maximum Power Point Tracking (MPPT) [15,16] and power management hardware features [11].

The energy lifetime of self-powered wireless communication devices depends not only on energy-harvesting capabilities, but also on energy demands, which are rarely static and change as a function of radio operating modes.

The main motivation to work on a new EH model is the need to accurately describe some energy-dependent effects frequently observed in these circuits. I–V tracking is demonstrated to be accurate when components present different power dissipation at either different operating voltages or drawn currents. Previous energy-based approaches calculate capacitor available energy by adding and subtracting harvested and consumed energy, respectively [11,17–19]. However, there are some objections to these works. In [17], energy management components, such as DC-DC converters or low drop-out regulators (LDO) are not considered, so the model is not representative in ultra low power modes. In [19], a constant energy buffer charging efficiency and a constant leakage power are considered with some boundaries. However, as [18] demonstrates, these parameters are not constant, but depend on instant voltage and current conditions. In [18], DC-DC circuits are included, and I–V operating conditions are concluded to be needed. However, the mathematical formulation is very complex and will require many computational resources, which can be efficiently reduced, as this paper shows. Empirical methods and prototyping are used in [11]. However, empirical models cannot always be generalized.

In conclusion, a reliable numerical energy model implies I–V tracking. Nevertheless, in this paper, an alternative is presented with the aim of reducing computational load. The model is called SIVEH and is I–V tracking-based. Computational complexity is reduced, and simulation time is significantly improved using both variables.

Energy management policies have already been studied for networks with energy-harvesting sensor nodes [19–22], etc. Power management systems contribute to adjust sensor node consumption to available energy. Thus, SIVEH has also been designed to enable the simulation of energy-neutral operation (ENO) conditions.

However, numerical computing is limited in sensor network simulation, and event-based computer network simulation tools are powerful [22]. This motivated an early implementation of SIVEH in ns3 [23]. This paper addresses the full description of SIVEH for a more general model, which can be implemented in both network simulation tools and numerical computation tools.

3.Mathematical Approach

I–V tracking is strongly necessary for an accurate circuit simulation. The SPICE simulator [24], which is based on this principle, is widely used and has been used as a reference in previous EH systems design [8] and energy estimation tools [25]. The objective of SIVEH is to reduce computational load, increasing simulation speed and decreasing memory requirements, but still with high accuracy. To understand the model, basics on EH circuits based on supercapacitors are described in this section.

3.1.Supercapacitor Modeling

Supercapacitors are non-ideal rechargeable buffers. Energy stored in them is given by expression (1), where Vscap is the voltage in the capacitor and C is the rated capacitance. Therefore, a supercapacitor voltage decreases, as long as its remaining energy is depleted.

(1)Escap(t)=12C(Vscap(t))2

An energy model estimates remaining energy Escap at time (t + 1) using expression (2). The energy either consumed or harvested at time (t) depends on every component in the system. ΔE calculation increases complexity, which can be reduced with a I–V formulation without accuracy loss.

(2)Escap(t+Δt)=Escap(t)−ΔEconsumed(t)+ΔEharvested(t)

Alternatively, capacitor voltage at time (t) is calculated in expression (3). Vscap is directly dependent on current flowing through capacitor (Iscap).

(3)Vscap(t)=1C∫Iscap(t)dt

Iscap is the supercapacitor current at a given time instant (t) and it is calculated as the difference between demands (Iconsumed) and harvested current (Iharvested) in expression (4).

(4)Iscap(t)=Iconsumed(t)−Iharvested(t)

Circuit currents are estimated according to the circuit description and expressed in Equation (5), where circuit components have been divided into four groups discussed in Section 3.2.

(5)Iconsumed(t)=∑IR(t)+∑Isink(t)+∑IDC−DC(t)+∑Ilinear(t)

The presented approach assumes a constant current during each time slot (Δt). Thus, a key issue is the definition of the time slot. This assumption introduces an error to be bounded below an acceptable threshold. As this threshold can be set, as will be later demonstrated, SIVEH uses expression (6) instead of Equation (3) to calculate supercapacitor remaining voltage at the next time slot (Δt).

(6)Vscap(t+Δt)≃Vscap(t)−1CIscap(t)Δt

3.2.Circuit Components

Current demand, Iconsumed, is the addition of current drained in every component in the node. Reviewing the most recent EH circuits and the overall sensor node architectures, hardware usually consists of components of the following nature (Figure 1): (1a) resistors; (1b) current sinks; (1c) DC-DC energy management blocks; and (1d) Linear regulators.

3.2.1.Resistors

An examples of resistor abstract-modeled components are resistor network voltage dividers (e.g., feedback networks and comparator threshold). Current leakage in the supercapacitors can be also modeled as an ideal resistor [26].

Integrated circuits directly connected to supercapacitors dissipate some current, even inactive (quiescent) current. Moreover, components, such as the control unit (MCU), attached sensors or a wireless RF circuit are also modeled as a current basic type described in Figure (1b).

As the SIVEH model is formulated in terms of I–V, dissipated currents in expression (10) are directly obtained from vendor data-sheets.

(10)Isink(t)=f(Datasheet)

3.2.3.Voltage Regulators

Voltage regulators are frequently observed in harvesting circuits. These electronic components can be divided into two main separate groups: DC-DC converters (Figure 1c) and linear devices (Figure 1d).

DC-DC converters

Energy consumed by a DC-DC converter depends on the converter efficiency; expression (11), where Eout is the energy consumed by every circuit element attached to the DC-DC output.

(11)ΔEDC−DC(t)=Eout(t)effDC−DC

As a I–V formulation is used in SIVEH, input current to DC-DC converters is described using expression (12):

Vin denotes the voltage measured in the converter input. This case corresponds to the energy buffer voltage.

Vout denotes the back module supply voltage. This value is fixed beforehand, according to node demands (MCU, sensors, etc.).

Iout denotes the current dissipated by the modules connected to the converter output.

On the other hand, linear regulators and LDO voltage regulators can be considered linear devices (Figure 1d). Exceptionally, diodes can also be considered a part of this group under especial circumstances (e.g., supercapacitor charge using a diode, as in [14]).

Energy consumed using such devices is calculated in expression (13).

(13)ΔElinear(t)≃Vin(t)Iout(t)Δt

In the SIVEH approach, these elements are ideally modeled by setting input current (Iin) equal to the devices connected to its output current consumption (Iout), as expressed in Equation (14)(14)Ilinear(t)=Iout(t)

3.3.Energy Harvesters

In case of solar EH, supercapacitor recharge depends on solar radiation conditions, solar cells and the hardware between supercapacitors and solar cells. The SIVEH model contemplates these three items to calculate supercapacitor supplied current.

Firstly, solar energy values can be directly imported from solar radiation databases. This off-line process obtains more realistic scenarios with respect to world location conditions and meteorological variability. For example, some national or state local weather services record hourly-integrated solar energy
Jm2 levels that can be downloaded into comma separated values (.csv) format.
Jm2 vectors contain enough data to simulate from some hours to several years.

Secondly, solar power average and cell harvested currents are calculated by expressions (16), the ratio between hourly-integrated energy and average power, and (15), cell supplied power mean, respectively. Pcell limit in expression (15) is Pcellmax.

Icell is calculated following expression (17) and is limited to 300 mA.

(17)Icell(t)=P¯cell(t)Vscap(t)1effBlockA

4.Generic EH Circuits Description

Figure 2 shows a generic EH circuit block diagram. The figure tries to summarize all the recent advances in energy harvesters in a general architecture, embedding the most significant blocks. Additionally, each block is linked with the corresponding model described in Section 3.

The circuit includes two energy buffers. Supercapacitors are preferred as a primary energy buffer. They can be recharged several times very fast and without significant degradation. Since supercapacitor energy is limited and self-discharge is quite high, a secondary buffer is usually included. Both rechargeable and non-rechargeable batteries have been found in the literature.

Since the two types of energy buffers present different specifications, independent A, B and comparer (COMP) blocks have been considered in the design.

Block A handles buffer recharge using the energy harvested. Although a solar cell is shown in Figure 2, wind and even vibrations can be used as the primary energy source to recharge energy buffers. This block is usually implemented by DC-DC converters or linear devices. Resistors can also be considered in closed-loop feedback networks. The current sink model describes leakage currents.

Block B allows matching energy buffer output voltage and sensor node operating voltage. This block can also be implemented by means of DC-DC converters and/or linear devices. As in the previous block, some quiescent current is drawn that can be modeled using the current sink model.

An energy threshold comparer (COMP) is useful for estimating buffered remaining energy and handling flags in the sensor node. This block dissipates some current, modeled as the current sink, and a resistor network, modeled as a resistor.

Block C is a mechanism designed for switching the input supply line to the sensor node between the primary and the secondary energy buffers without MCU piggybacking. It guarantees sensor node stable supply voltage and optimizes the use of energy, optimizing energy consumed from the secondary buffer as long as there is remaining energy in the primary buffer. This block is usually modeled as a linear device with an associated quiescent current, modeled as a current sink.

4.1.EH Prototype for Evaluation Purpose Description

As final implementation of the different blocks depends on the actual hardware, an example is proposed: an improved version of [23]. The main innovations are MPPT and a backup battery. Energy storage is mainly carried out by supercapacitors. Two-point-three-volt Panasonic supercapacitors and a Li-Ion non-rechargeable CGR18650DA battery have been included as the main and secondary buffers, respectively. In future work, rechargeable batteries will also be considered. However, for verification purposes, the Li-Ion battery model is available [27] and has been added to SIVEH.

The energy-harvesting circuit is designed with solar cells rated at 4 V open-circuit and 3.5 peak voltage, 48.5 mA short-circuit and 45 mA peak current. The circuit connects either one cell or several, depending on energy demands. Block A1 in Figure 2 collects energy from solar cells and stores it in two 2.3 V supercapacitors connected in series. This block incorporates a Texas Instruments (TI) BQ25504 integrated circuit [28] that is a PWM converter to implement the MPPT algorithm for harvesting energy from low input sources (80 mV) with an efficiency above 80%.

The behavior of BQ25504 is as follows:

If supercapacitor voltage is below 1.8 V, a cold start procedure runs, and supercapacitor recharge current (Icell) is calculated using Equation (18) according to [23].

(18)Icell(t)=IcellmaxP¯solar(t)(Wm2)1000(Wm2)

In another case, the MPPT algorithm is carried out, and Icell is calculated using expression (19, a variation of the proposed current equation for DC-DC converters (12).

(19)Icell(t)=P¯cell(t)Vscap(t)effDC−DC

In any case, current is limited to 300 mA.

Since supercapacitor voltage decreases when stored energy is drained, block B1 provides a stable 3.3 V MCU supply voltage, while supercapacitor voltage is between 1 V and 4.6 V. This block includes, for testing every basic circuit described in Figure 1, a DC-DC step-up converter (L6290 ST Microelectronics [29]) together with a TI TPS782333 LDO [28]. Block B behavior is modeled as follows:

If supercapacitor voltage is below 1 V, block B does not work properly, and only quiescent current (the current sink model) is dissipated according to datasheet values.

If supercapacitor voltage is between 1 V and 3.3 V, block B acts as a 3.3 V step-up converter. The DC-DC converter model is used (expression (12)), adding the correspondent quiescent current (the current sink model).

If supercapacitor voltage is above 3.3 V, block B acts as a linear regulator. The linear device model (expression (14)) is used in this case.

A flag is set when supercapacitor voltage is below 1 V. The COM block consists of a voltage comparator that drains 400 nA (the current sink model) and a 1 MΩ voltage divider resistor network (resistor model).

Since a 3.6 V non-rechargeable battery has been included into the prototype, blocks A2 and B2 are disabled. Iinbat and Ioutbat are equal in Figure 2. Both are renamed as Ibat.

As already mentioned in Section 3, time slot approximation introduces a numerical error. Errors can be quantified using Maclaurin series and the Taylor error theorem, as shown in the Appendices. Errors are summarized in Table 1. SIVEH error is bounded by defining a simulation time slot (Δt) according to the capacitor. For example, a DC-DC converter, C = 10 mF, eff = 92% and Vout = 2.3 V, attached to a resistor of 100 Ω, shows an error below 1% for any Δt below 3.9 ms.

5.1.Simulation Analysis

The use of analog electronic circuit simulators, such as SPICE, is widely accepted. The basic circuits in Figure 1 have been simulated using a SPICE-based simulator (LT-SPICE [30]), the SIVEH model and the energy approach. The objective of this section is to evaluate both the accuracy and simulation performance.

DC-DC comparison is carried out using L6920 of ST [29] and outputs plotted in Figure 3 that shows a model deviation from the SPICE reference less than 0.7%. The mean error and mean square error [Me, Mqe] in Isink and resistors are [0.04%, 0.04%] and [0.008%, 0.012%], respectively. Energy tracking and SIVEH present barely the same accuracy for all cases: 0.1% maximum deviation measured in the resistor case in which both solutions present the some error.

However, SPICE simulation time is very lengthy with DC-DC converters. For example, LT-SPICE takes around 20 h to simulate the complete discharge of a 10 F capacitor. However, numerical computing is almost constant, independent of component complexity. Table 2 shows a practical comparative example. The higher performance speed of the SIVEH approach is shown in both SPICE and the numerical energy approach. For example, simulation speed increases, on average, 240% compared with the energy models.

6.SIVEH Model Implementation

The SIVEH model has been implemented using a numerical computing language for fast evaluation. It has also successfully been implemented in discrete-event network simulators, such as ns-3 [31]. Figure 4 presents the overall abstraction flow diagram of SIVEH.

Energy data depends on the sensor node operating mode at the current time slot (t) and includes three tasks:
(a)

Solar energy is read from the solar database, and cell power is calculated using expression (15).

(b)

Supercapacitor self-discharge, as well as resistor networks, IC quiescent current, etc., directly attached to the supercapacitors, are calculated using expressions (8) and (10).

(c)

Sensor node current consumption is set according to the current operating mode.

Supervisor: Block C in Figure 2 checks for a minimum threshold condition in Vscap and determines whether the sensor node will be powered by the supercapacitor (Iscapout) or by the battery (Ibat). Some quiescent current is drawn by the non-connected energy source.

Energy dissipated by DC-DC components themselves is also calculated with expression (10).

(c)

Block A1 calculates Icell(t).

Energy buffers models: Buffering voltages are updated.

(a)

Supercapacitor voltage is updated with expressions (6) and (4), as the addition of Ir, Isink, Icell and Iscapin, shown in Figure 4.

(b)

Battery voltage (Vbat) is updated using the model proposed in [27].

New simulation time: At the end of the flow diagram, time is also updated to the next period (t + 1).

7.Model Verification

SIVEH model estimations have been compared with the prototype. For these experiments, the sensor node was programmed with the state diagram in Figure 7: sleep time: 7 s and 10 μW; duty cycle 7.6%: Rx reception 512 ms and 24 mW; Tx transmission: 64 ms and 108 mW. Consumption values have been collected from [23,32].

Two experiments have been carried out. Firstly, indoor experiments to force fast recharge and discharge (860
Wm2 were measured in recharge and 0 W/m2 in discharge). Figure 5 shows the indoor experimental results. Both SIVEH and prototype outputs match with a slide deviation of 2%.

Secondly, the prototype was located outdoors from March 1 to 4. The experiment was carried out in a river 7.8 km away from a AEMet weather station placed at 39° 29′ 7” N–0° 28′ 28”. The AEMet solar record is shown in Figure 6a.

The experiment started on February 29, 2012, at 8 pm. The initial Vscap was 2.8 V. Discharge reached 1 V before sunrise (see Figure 6b), leaving the node without enough energy to communicate. Both simulated and real curves show the same effect. However, the simulation curve starts the supercapacitor recharge before the prototype. This effect is given by the one-hour time slot grain of the database. While the value observed by the prototype varies continuously, the simulation uses a sample hourly-integrated solar energy (
Jm2) level, which means an error experimentally observed at sunrise and sunset of 0.1 V in the worst case; Figure 6b. From sunrise to sunset, both lines match. Figure 6a shows solar conditions during the experiment.

The low consumption of the modem architecture and the energy harvesting circuit, as well as the quick recharge of supercapacitors at sunrise make it possible to maintain high voltage levels in the supercapacitor during the whole day, even in cloudy weather (see March 2 in Figure 6a). It should be noted that all simulation average values, obtained from manufacturer datasheets, have been used for all components. Statistically distributed results can also be obtained, introducing randomness in the nominal values of components, which can also be obtained in the same datasheets. Experiments should then be repeated as needed.

8.Case Study: Energy Neutral Operation

Energy neutral operation (ENO) is the use of harvested energy at an appropriate rate, such that the system continues to operate perennially [19]. The definition of energy neutral operation needs to be evaluated from a very early stage of the design. This implies simultaneous cross layer work on physical design decisions and management power policies.

This section shows an example of the ENO condition in a system with the basic state diagram shown in Figure 7, which will be dynamically re-adjusted (Tsleep, TRX or TTX).

Following the state diagram in Figure 7, node average current consumption is calculated with Equation (20), whose reformulation supports Figure 7 extensions.

(20)I¯node=ITXTTX+IRXTRX+IsleepTsleepTTX+TRX+Tsleep

ENO is based on the law of conservation of energy: harvested energy must be equal or greater (excess would be stored) than dissipated energy; Equation (21).

(21)Eharvested≥EdissipatedPharvested≥Pdissipated

Assuming Vscap is higher than a minimum threshold, Pharvested is given in expression (22). This threshold is 1.8 V and assures that the BQ25504 energy harvester operates with a MPPT algorithm.

Pdissipated depends on sensor node power consumption requirements (Pnodeef f, since Block B1 efficiency has to be considered), quiescent and leakage currents (Isinkeq) and the circuit resistor equivalent value (Req). Pnodeeff is usually the most energy-demanding element. Therefore, Pdissipated can be approximated with expression (20). Thus, the initial ENO premise (21) is expanded in (23).

Node sleep time is calculated from Equation (23) and is dynamically adjusted using the algorithm, defined below as the ENO Algorithm. Sleep time variation increases or decreases the packet transmission rate in order to consume energy proportionally to the harvested energy. Additionally, the ENO Algorithm sets the energy safe maximum boundary (Vscap = 4.5 V) and minimum and maximum sleep times.

Listing 1. Node sleep time basic decision algorithm using ENO

int ENO_algorithm (int Vscap, int Pcell,

int Tsleep max, int Tsleep min)

{

if (Vscap>4.5 V & Pcell >=0.2*Pcell_max)

Tsleep = (eq 22);

else

Tsleep = Tsleep / rho;

if (Tsleep > Tsleep_max)

Tsleep=Tsleep max;

else

if (Tsleep < Tsleep_min)

Tsleep=Tsleep min;

return (Tsleep);

}

The following simulations analyze the use of ENO condition. Simulations were set on December 8. The ENO settings were: Tsleepmin = 0 s and Tsleepmax = 40 s.

The evolution of the energy in the supercapacitors is shown in Figure 8.

Simulation without the ENO condition: 2,614 cycles (1.36 cycles per minute). The output showed that the node needed to be powered by the backup battery for 3 h.

Simulation with ENO condition: 67,942 cycles (35.3 cycles per minute). The output showed that the node needed to be powered by the backup battery for 3.5 h.

In the SIVEH numerical model, it is possible to introduce power management policies based on probabilities. For instance, sleep time can be conditioned to an external variable, such as Icell. When Icell decreases below 20% of Icell maximum value, the packet transmission rate decreases into (1 — ρ), for example, with a ρ = 10%:

Simulation with ENO condition + ρ 10%: 68,571 received messages and 67,419 transmitted messages (35.7 cycles per minute in rechargeable conditions and 35.1 cycles per minute at night). With this adjustment, no battery energy was required.

These results show that a static duty cycle is not really efficient. The ENO condition enables a dynamic adjustment of transmission and reception packet rates for harvested energy which becomes more profitable when taking into account solar radiation differences observed between summer and winter in many world locations.

9.Conclusions

This paper presents a numerical model to simulate capacitor energy storage with solar energy-harvesting applied to WSN nodes (SIVEH). The mathematical basis and error boundaries have been accurately described. The proposed model is computationally faster than both an energy-based formulation and a SPICE model. For example, SIVEH reaches a simulation speed up to 54,678 times faster than a SPICE simulation of DC-DC components. Moreover, SIVEH presents an accuracy of 98% of a SPICE reference model.

The model has been designed to support any state-of-the-art energy harvesting platforms, after an extensive analysis, and verified by comparing simulations with a real prototype with less than 2% maximum experimental error.

The model versatility allows for introducing additional features, such as the ENO condition. An example demonstrates that an efficient use of harvested energy could increase communication cycle rates up to 25.9% on average.

Beyond that, the model is enhanced with the ability to work with (
Jm2) solar radiation vectors. Thus, the model is a tool that enables an efficient design of the hardware architecture, according to the power management policies specified for the system and the specific world location of the final system.

As a summary, the model is generic and wide enough to be able to simulate any energy harvesting circuit composed of the typical elements of these circuits.

Thus, with the accuracy observed in SIVEH, many architectural design decisions and predictions can be analyzed without requiring a prototype.

The authors gratefully acknowledge financial support from CICYT. ANDREA: Automated Inspection and Remote Performance of Marine Fish Farms (CTM2011-29691-C02-01); and RIDeWAM: Research on Improvement of the Dependability of WSN-based Applications by Developing a Hybrid Monitoring Platform. (TIN2011-28435-C03-01).

Differential Equation (24) describes the behavior of the basic RC circuit shown in Figure (1a). The exact solution for a certain time slot (t) is shown in Equation (25).

(24)0=1C∫i(t)dt+Ri(t)(25)vc(t+Δt)=vc(t)e−ΔtRC

This expression can be approximated by a finite number of terms using Taylor series. In fact, since Δ t is very small (the Δ t≪RC condition must be fulfilled), the Maclaurin series is centered on zero (expression (26)).