1. Introduction
Pyranometers are devices used to measure solar radiation that are composed of direct and diffuse incident radiation [
1]. It is used, for example, to measure the energy transfer rate per unit area to produce the Global Horizontal Irradiation (GHI), which, in turn, is used to determine the energy potential for photovoltaic systems in a region [
2].
According to the ISO Norm 9060:2018 [
3], there are two pyranometer types: (i) the silicon-cell pyranometer, which provides a fast time response but with limited spectral range; and (ii) the thermopile pyranometer, which provides a high and flat spectral range but with a slow time response. Alternatively, there is a class of pyranometers based on the Electrical Equivalence Principle (EEP) or Electrical Substitution Principle (ESP). The EEP uses a feedback architecture to self-balance the sensor heat dissipation with the sensor incoming power, for which the latter is composed of the electrical power and the incident radiation power. The incident radiation can, then, be estimated through the thermal balance equation [
4,
5,
6,
7].
The majority of pyranometers based on the EEP operate at a constant temperature, which has the main benefit of reducing the response time [
8,
9]. The classical analog architecture is composed of a two-analog Wheatstone Bridge in a differential voltage configuration, employing two sensors [
6]. One sensor is painted in black and is used mainly for sensing the incident radiation, while the other one is painted in white and serves to compensate the influence of the ambient temperature. Similarly, other constant temperature architectures, like the PWM [
10] and the Sigma-Delta [
7,
11], also employ two sensors.
The ambient temperature directly affects the sensor thermal balance, which, for the architecture with the sensor operating at a constant temperature, has an important effect on the sensitivity of the system, among other parameters [
1]. In this paper, we propose an analog pyranometer architecture based on the EEP that operates with a controlled temperature difference, which has the main advantage of mitigating the influence of ambient temperature variation of the system performance. The analysis performed show five times sensibility improvement, reduction of the ambient temperature influence, and
increase in the output voltage useful span.
The remainder of this paper is organized as follows:
Section 2 presents the principle of electrical equivalence applied to the classical architecture of the Wheatstone bridge. Limitations of the bridge in relation to the variation of the ambient temperature are demonstrated. A method for changing the sensor operation and the necessary modifications of the classical architecture to transform them into the proposed architecture is also proposed. In
Section 3, an analytical analysis of the proposed architecture is carried out, followed by a validation by computer simulation. In
Section 4, general considerations about the work developed are presented, as well as the main contributions of the work.
Section 5 contains the final conclusions and, in
Section 6, patents deposited from this research are presented.
2. Material and Methods
The classical pyranometers based on EEP control the temperature of a thermoresistive sensor through the Joule effect to keep it constant. As consequence, the incident radiation can be estimated by the thermal balance, while the feedback control reduces the system response time.
2.1. State of the Art
The electrical equivalence principle is historically used for measure incident radiation and fluids speeds. In recent years, there has been a growing trend to replace analog architecture with digital for the greater ease of correction of ambient temperature variation. In addition, many works have contributed to the study of uncertainties and the improvement of sensitivity.
Table 1 summarizes the State of the Art in the last 6 years.
2.2. Thermal Balance in Electrical Equivalence Pyranometers
A thermoresistive sensor varies its resistance with the temperature, and can present either a negative or a positive temperature coefficient, called NTC and PTC, respectively. NTC sensors have a nonlinear resistance–temperature relationship but generally a much higher sensitivity than PTC, and, for this reason, they are chosen. The NTC sensor resistance can be approximated by the
-model equation as [
17,
18]:
where
[
] represents the electrical resistance at 25
C,
[K] is the intrinsic thermistor temperature, and
[
C] is the sensor temperature.
The sensor thermal balance equation relates its temperature with the sensor incoming and lost power, and can be given by [
6,
7,
14,
19,
20,
21]:
where
[W] represents the electrical power,
[W] represents the absorbed radiation power,
[W
C
] is the thermal dissipation factor,
[J
C
] the thermal capacity,
[
C] the sensor temperature, and
[
C] the ambient temperature.
The thermal dissipation factor
is nonlinearly dependent on the fluid dynamics, and to avoid the influence of the latter, a transparent dome is used. Hence, the thermal dissipation factor can be approximated by a constant [
22,
23]. The sensor module, composed of a transparent dome, the sensor, and a supporting base, is illustrated in
Figure 1.
The power absorbed by sensor due to the incident radiation,
, expressed as [
24,
25]:
where
H is the incident radiation power per unit area [W/m
], and
can be viewed as a calibration constant [m
], which is composed of the dome transmittance, sensor absorption coefficient, and exposed sensor area to the incident radiation.
In steady-state, (
2) can be written as:
where
represents the voltage across the sensor and
the sensor resistance.
2.3. Constant Temperature Wheatstone Bridge Architecture
The Constant Temperature (CT) Wheatstone bridge architecture, shown in
Figure 2, is responsible for controlling the electrical power delivered to the sensor to keep it at a constant temperature [
6,
15]. In steady-state, the difference between the sensor heat dissipation and electrical power is equivalent to its absorbed power, and
H can be estimated from (
4).
As shown in
Figure 2, the Wheatstone bridge branches are classified as: (i) sensing branch and (ii) the reference branch. As the name suggests, the sensing branch contains a sensor that heats up with electrical energy and incident radiation. The reference branch is used as a reference for the sensor temperature. An operational amplifier (OpAmp) is used to provide a high negative feedback gain making the uncertainties associated with measurements parameters negligible [
28]. It may also include a bipolar transistor to provide sufficient current gain to drive the bridge [
8,
9,
29].
Any variation in the incident radiation power tends to change the temperature and resistance of the sensor, and consequently to unbalance the bridge. However, this variation is sensed by the OpAmp, which, due to its high open-loop gain, increases or decreases the electrical output power delivered to the bridge to bring it back to equilibrium, keeping the sensor temperature and resistance constant [
14,
27].
From
Figure 2, neglecting the effect of the OpAmp input voltage offset, in the thermal balance steady-state, the following can be written:
where
and
Considering that, in steady state, the voltages
and
are approximately equal, it becomes possible to determine the sensor resistance
that is greater than zero, as a ratio of the resistances of the Wheatstone bridge:
and this relationship is used to implement our proposed architecture.
Recombining (
5) and (
6), the OpAmp output voltage can be found as:
in which
is used to calculate the output voltages needed to keep the sensor warm in any condition of ambient temperature and incident radiation.
2.4. Constant Temperature Difference Wheatstone Bridge Architecture
The CT architecture features a high feedback gain in the analog domain, which results in a rapid compensation of the electrical power delivered to the sensor. This compensation is directly influenced by the ambient temperature variation that drives the architecture to increase or decrease the output power influencing directly the incident radiation estimation. To preserve the architecture advantage and mitigate the effect caused by ambient temperature variation, it is proposed to use a constant temperature difference instead of constant temperature in Wheatstone bridge architecture. The proposed architecture is called the Constant Temperature Difference (CTD) Wheatstone bridge.
The proposed CTD Wheatstone bridge architecture minimizes the problem of dependence of the output voltage range on the ambient temperature presented by the classical architecture, controlling the sensor temperature so that it is different from the ambient temperature, i.e.,
, is constant. Rewriting (
9) using
gives:
with
The best value of
in (
10) is the one that ensures no output signal saturation for the maximum incident radiation and maximizes the architecture sensitivity. For maximizing the sensitivity, it is desirable for the output voltage (
) to be zero for the maximum incident radiation. Hence,
can be defined as:
Keeping the sensor at above the ambient temperature requires some modifications in the Wheatstone bridge reference branch. The main change is the addition of a second sensor used to sense the ambient temperature.
The sensor responsible for sensing the ambient temperature is called cold sensor (
), while the sensor responsible for receiving incident radiation is called hot sensor (
). The constant
represents the temperature difference between the hot and cold sensors. The branch that contains the hot and cold sensors are referred to as the hot branch and cold branch, respectively.
Figure 3 shows the proposed Wheatstone CTD bridge architecture using NTCs.
Like CT architecture, the OpAmp also provides a positive voltage to the transistor base, which is used to increase the current capacity. However, the transistor may be reverse-biased if the heated sensor temperature saturates and/or the OpAmp offset voltage is negative. Both cases can cause the circuit to stop functioning correctly. It is also important that the transistor threshold be greater than the differential power supply to avoid damage.
The problem of saturation can be avoided by keeping the OpAmp offset positive and adding a reference resistor between nodes and that provides a minimum current for the OpAmp inputs to function properly. The guarantees a minimum voltage reference so that the OpAmp returns to provide a positive voltage output when the sensor returns to the normal temperature condition.
The hot sensor in the CTD Wheatstone bridge architecture is chosen in a similar manner to the one in the CT Wheatstone bridge architecture. The value of
is optimized to provide the maximum voltage variation in the Wheatstone-Bridge over the change of
. This guarantees the maximum possible linearization over the ambient temperature range of the nonlinear sensor-based transfer function. The CTD Wheatstone bridge architecture value for (
) is chosen to be the geometric mean of the extreme sensor values. The CTD Wheatstone bridge architecture optimal value for (
) is:
where
is the hot sensor static response.
To keep the hot sensor (
) heated up
over the ambient temperature is necessary to satisfy the bridge balance that, neglecting the
effect, is:
Resistors and serve to compensate for the intrinsic temperature of the cold sensor. A good choice of presents a high value of and similar to . The high value of prevents the cold sensor from self-heating.
One methodology to find
,
and
values employs a nonlinear system approximated by three or more points of ambient temperature. A good practice is to use the maximum, average and minimum points of ambient temperature. Other constraints of the nonlinear system are the values of
,
and
which cannot be negative and the cold branch must have enough resistance to prevent
from self-heating. This way, the nonlinear system is given by:
where
[W/
C] are the thermal dissipation factor for hot and cold sensors, respectively, and
[
C] the maximum allowable self-heat of the cold sensor. The solution of (
15) can be achieved using an optimization algorithm like the trust-region or Levenberg–Marquardt.
2.5. Experimental Setup
We prototyped the proposed CTD Wheatstone bridge using a glass dome, a metal case, and a printed circuit board for the electronic components. The printed circuit board contains the connectors to input/output and the conditioning circuit that is composed of: 2PC4081R as a bipolar NPN transistor, the TL081 as operational amplifier with offset adjustment, and two NTCs thermistors are employed as hot and cold sensors.
Afterwards, to validate our proposed architecture, we compared the experimental results with the Hukseflux SR-05 thermopile pyranometer.
Figure 4 shows the commercial SR-05 pyranometer and proposed CTD prototype from a perspective view.
The data of each pyranometer were sampled using a Data Acquisition Device (DAQ), and processed in a computer to compute the incident radiation, as illustrated in
Figure 5. The DAQ used is a NI USB-6361 and the computer presented the following configuration: AMD A8-5500 @ 3.2 GHz Quad-core 32 nm, 4 GB of RAM, 512 GB HD and a 3.0 USB port. Two 12 volt @ 7 ampere-hour batteries were used as a pyranometer power supply.
The data acquisition device was configured to read the differential mode with a resolution of 16 bits, −10 to 10 V input range and a sampling rate of 100 Sa/s. The SR-05 pyranometer signal was acquired through the 4–20 mA analog output, and the CTD prototype signal acquired by nodes and .
A Fluke 1502A thermometer and a Keysight DMM 34470 were used to calibrate the CTD prototype sensors’ static response. The dynamics response calibration was performed applying a Pseudo-Random Multilevel Sequence (PRMS) in a current source with Daq [
30,
31,
32]. The dynamics response obtained was used in a system identification toolbox to obtain the dynamics parameters of the system [
24]. A search method of Levenberg–Marquardt was applied with the current range between 300
A to 40 mA. Experimental results and CTD parameters obtained are shown in
Section 3.
3. Results and Analysis
3.1. Analysis of CT and CTD Wheatstone Bridge Architecture
To analyze the performance of the CT and CTD Wheatstone bridge architecture, consider that the ambient temperature varies from 0 to 60
C, and the incident radiation varies between 0 and 1600 kW/m
and
V. The sensor and parameters for the CT and CTD architectures are summarized in
Table 2 and
Table 3, respectively. The parameters presented were adapted from [
1,
7,
11] with computed values from Equation (
12), which avoid any temperature saturation, and (
13), which optimize the value of
for CTD architecture.
Applying the parameters from
Table 3 to (
15), it is possible to obtain the relation of resistances that allow a constant temperature difference between hot and cold sensors. The resistance values obtained for the parameters used were:
k
,
and
M
. These resistances values allow a fit of 99.6% in the comparison of the transfer functions of the hot branch and the cold branch.
Figure 6 shows the hot and cold branch transfer function of CTD Wheatstone bridge architecture. In order for
be constant, it is necessary that hot and cold branches’ transfer functions, represented respectively by
and
, are the same for the entire ambient temperature range, as described in (). If the hot branch and cold branch transfer functions are not same, it will lead to a
different from the one projected for the heated sensor and, consequently, to an error in the incident radiation estimation.
The values from
Table 2 and
Table 3 were applied in the simulation of CT and CTD Wheatstone bridge architectures. From (
9) and (
11), the output voltage span was obtained as a function of ambient temperature and incident radiation, as shown in
Figure 7.
From
Figure 7, it can be observed that the output voltage varies significantly between the architectures. For comparing the performance of both architectures, the metrics of average sensitivity
, sensitivity variation
, and useful voltage range
, defined as:
Ideally, the sensitivity should be as high and constant as possible, while the useful voltage range should be 100%. The following metrics were obtained from results (
Table 4):
Compared to the CT architecture, the proposed architecture has a 5-fold increase in average sensitivity, approximately 3-fold decrease in sensitivity variation, and a 76.3% increase in the useful output voltage range.
Considering that the electrical power is mostly dissipated in the branch that contains the hot-sensor, one can estimate this branch electrical power consumption considering only the heated sensor and the resistor
[
13].
Figure 8 illustrates the hot branch highest electrical power consumption in steady-state by ambient temperature.
Figure 8 shows that CT architecture power consumption is inversely proportional to the ambient temperature while the CTD architecture increases nonlinearly with ambient temperature. In this case, the CT is more efficiently in ambient temperature above 50
C due to the low value of
.
3.2. SPICE Simulation
A simulation was also performed using a SPICE environment with the classic architecture (CT) and the proposed architecture (CTD), considering varying both the ambient temperature and the incident radiation. The objective of this simulation was to verify if the temperature of the heated sensor would follow the intended reference and to explore the effects of using the temperature difference in the estimation of H.
The Wheatstone CT and CTD bridge architecture, shown in
Figure 2 and
Figure 3, were simulated using a thermistor model provided in [
1]. The simulation parameters are summarized as follows:
Table 5 contains the CT Wheatstone bridge architecture sensor parameters;
Table 6 contains the CT Wheatstone bridge architecture fixed resistances that allow
to be
C;
Table 7 contains the CTD Wheatstone bridge architecture sensor parameters;
Table 8 contains the CTD Wheatstone bridge architecture fixed resistances that allow
to be
C;
Table 9 contains simulation parameters shared between CT and CTD Wheatstone bridge architecture. The parameters
represent the offset voltage and
the open loop gain of the operational amplifier.
The ambient temperature and incident radiation are defined, respectively, in simulation as and , where is the step function.
Figure 9 shows the ambient temperature
, the sensor temperature
in CT architecture, the hot
and cold
sensor temperature in CTD architecture, the incident radiation H, and the difference of temperature
.
As evidenced by
Figure 9, the sensor temperature
of the classical CT Wheatstone bridge architecture remains practically constant regardless of the ambient temperature, while the temperature of the hot sensor
in the CTD Wheatstone bridge architecture remains
degrees warmer than the ambient temperature.
To estimate the incident radiation, (
5) can be used as:
where
[W/m
] is the estimated incident radiation. The estimated incident radiation for CT and CTD Wheastone bridge architectures is shown in
Figure 10. The
sensing in CT architecture simulation followed the cold sensor time constant for fair comparison.
Small fluctuations in sensor temperature (CT) or temperature difference (CTD) lead to a small error that becomes noticeable due to amplification caused by . The mean absolute error of simulated architectures was W/m for CT and W/m for CTD.
3.3. Experimental Results
The developed prototype was submitted to a field experiment with the Hukseflux SR-05 pyranometer as a reference. Before the measurement, the sensor calibration was performed as described in
Section 2.5.
Table 10 presents the static and dynamics sensor parameters used in CTD prototypes.
The pyranometer was exposed to the sun and observed the maximum temperature difference between the hot and cold sensors, without control, in order to find the maximum temperature variation and
parameter, according to (
12). The maximum difference observed was
C, so the project parameter was set to
C. From
and the sensors’ static response, a new calculus of the bridge resistances was performed. The values of resistances are summarized in
Table 11, and other project parameters are summarized in
Table 12.
All measurements were performed for 20 min in three different weather conditions, at the Federal University of Rio Grande do Norte, RN, Brazil, with the following geospatial coordinates: −55032.6034 (lat.), −351150.9274 (long.).
The SR-05 pyranometer and the CTD prototype results were acquired by DAQ and digitally processed to estimate incident radiation, the measurement error, and the sensors’ temperature. A low-pass filter, with a cutoff frequency of 10 hertz, was applied to the prototype CTD signal to remove high-frequency components. The filter was not applied to the SR-05 as this device already has a low frequency response.
Figure 11 shows the “A” measurement taken on 16 February 2022 at 11:10 a.m. with clear sunny conditions. During the experiments, few cumulus clouds in the sky were noticed.
Figure 12 shows the measurement difference between the CTD prototype and the pyranometer in “A” measurement. The absolute error measured was
W/m
.
Figure 13 shows the temperatures of the hot sensor (
), cold sensor (
), and the temperature difference (
) during “A” measurement.
Figure 14 shows the “B” measurement taken on 14 February 2022 at 11:10 a.m. with partially cloudy conditions. During the experiments, cumulus and stratocumulus clouds were present.
Figure 15 shows the measurement difference between the CTD prototype and the pyranometer in “B” measurement. The absolute error measured was
W/m
.
Figure 16 shows the temperatures of the hot sensor (
), cold sensor (
), and the temperature difference (
) during “B” measurement.
Figure 17 shows the “C” measurement taken on 18 February 2022 at 11:30 a.m. with cloudy weather conditions. Nimbostratus clouds were present but without rain.
Figure 18 shows the measurement difference between the CTD prototype and the pyranometer in “C” measurement. The absolute error measured was
W/m
.
Figure 19 shows the temperatures of the hot sensor (
), cold sensor (
), and the temperature difference (
).
3.4. Enhancement and Post-Processing Information
The enhancements can be performed in two ways: (i) in calibration process; (ii) in post-processing. In the calibration process, several measurements can be made simultaneously between calibrated pyranometers and the CTD prototype. The optimization algorithms can be employed to determine the best constant values, , and that reduce the measurement difference between the calibrated device and the prototype. This method must be used carefully, and it is recommended only to refine the sensor parameters.
Post-processing methods include smoothing the response using digital filters and applying correction factors. In general, the correction factor is only used if the incident radiation is constructed using only the signal. Using signal node with signal node provides greater reliability in signal acquisition. The voltage signal at nodes and can be used to determine if the sensor has already reached its steady state, avoiding measurement errors.
From
,
, and
, it is possible to select the best values to reconstruct
H that happens when
is equal to
. This prevents any dynamic interference caused by the hot and cold sensor constants. For this, a comparator can be added to nodes
and
to generate a trigger signal for the sampling device. The comparator output signal must be one if
and zero if
.
Table 13 shows the absolute error using only
, using
and
, and using all nodes to reconstruct H.
4. Discussion
This paper presented the constant temperature difference Wheatstone bridge architecture and shows its analytical model, a SPICE simulation, and a set of field experiments carried out with a developed prototype. Information and methodologies to find the ideal temperature difference for proposed architecture are also provided.
The proposed architecture is intended to mitigate the effects of ambient temperature variation on incident radiation measurement, while maintaining the advantage of fast response of the Wheatstone bridge based architecture. The constant temperature difference use abruptly changes the output voltage range of the circuit, promoting a greater range of useful voltage, with greater constancy without DC level as shown in
Figure 7. It is a great advantage especially if it is planned to add an analog-to-digital converter.
The increase of useful output range reflects directly on the device sensitivity that presents more constancy in all ambient temperature range, and in power consumption. The power consumption of our proposed architecture is significantly at lower temperatures when compared with classical architecture, as shown in
Figure 8. The energy consumption of the proposed architecture is only higher at high temperatures due to the
that is generally great in CTD than CT architecture. The power consumption is an important parameter for this type of device as it can be powered by batteries and deployed in remote locations.
A simulation was performed and showed that CTD Wheatstone bridge architecture can increase or decrease the delivered power to the sensor in order to keep it higher than ambient temperature. An equation was presented to find out the optimized value. From the parameters, it is possible to estimate the incident radiation simulated. Simulated values were compared and showed that our proposal provides lower uncertainty than classical architecture. This happens because there is a greater margin for temperature variation compensation by OpAmp when using the difference instead constant temperature.
The simulation results corroborate the theory by keeping the difference between the sensors almost constant. There is also a small fluctuation that occurs in the sensor temperature due to the sensors thermal capabilities (
). Usually this variation is minimal but eventually becomes noticeable due to amplification by the calibration constant (
). This effect can be reduced by considering the actual temperature instead of the constant in the reconstruction equation (
19).
A set of field experiments were performed with the proposed CTD pyranometer that is compared to SR-05 as a reference pyranometer. The main difference between the SR-05 pyranometer and proposed CTD pyranometer is the response time. The SR-05 has a response time of ≈18 s, while our proposed architecture shows a response time in milliseconds due to the high proportional feedback gain.
The CTD pyranometer variation is directly associated with its dynamics constants. While the temperature presents high precision, due to the high gain in closed-loop, the incident radiation is amplified in an open loop by the low value of constant. In order to smooth the response and reduce the uncertainty, it is desirable to increase the constant and reduce the dissipation factor . This may allow for increasing the temperature difference and, consequently, reducing the variation.
One way to reduce the variation in the proposed prototype response is adding a low-pass filter to the OpAmp output voltage or by increasing the calibration constant. The calibration constant can be improved by using a thin layer of special inks with high light absorption capacity.
Figure 11,
Figure 14 and
Figure 17 showed a good similarity between the CTD prototype and the SR-05 pyranometer in different weather conditions. In some specific points, there is a noticeable variation only in our prototype. This contrast is explained by the different spectral device responses, which can be adjusted by covering the sensor with the same special black ink used in the commercial pyranometer. Despite this, all valleys and peaks are perceived to a similar degree and intensity in both devices.
The difference between the CTD prototype and the SR-05 pyranometer shown in
Figure 12,
Figure 15, and
Figure 18 is similar to a random variable with a non-zero mean. The main challenge is to keep its dispersion as low as possible with an average close to zero. To reach this objective, this error can be employed in minimization algorithms to improve calibration of CTD pyranometers.
Regarding the temperature differences shown in
Figure 13,
Figure 16, and
Figure 19, it can be concluded that, as the incident radiation increases, the ambient temperature also increases. In all the cases presented, the controller managed to keep the temperature difference stabilized at 12
C and thus drastically reduce the influence of the ambient temperature variation on the measurement.
The main contributions of this work were to analyze classical architecture from a different perspective, proposing to use the temperature difference instead of the constant temperature. A method was developed to find parameters that allow the Wheatstone bridge to keep the sensor warm to degrees above ambient temperature and to carry out simulations and field experiments to validate the proposed architecture.
In this work, the fundamentals were applied to Wheatstone bridge architecture, but can be adapted to other architectures, analog or digital. This method can also be adapted to a Wheatstone bridge anemometer.