1. Introduction and Background
Domestic electric water heaters (DEWHs) are considered to be one of the largest energy-consuming devices in a typical domestic household. In first-world regions, such as Australia, the European Union and USA, water heating makes up 23%, 14% and 18% of the total residential load, respectively [
1]. In South Africa, it can contribute to anything in a range of 30% to 35% of a household’s total energy consumption, and accounts for approximately 7% of South Africa’s total grid load [
2,
3,
4].
It is no secret that electrical energy is a scarce commodity, especially in many developing countries that struggle to meet the ever-increasing energy demands. In Sub-Saharan Africa, approximately 600 million people still live without access to electricity, which is more than in any other part of the world [
5].
Electric water heaters (EWHs) can be used for grid stabilisation purposes using demand side management (DSM) for peak shaving and energy saving [
6]. There are at least 5.4 million of these heaters in the country [
3], making them ideal candidates for demand response (DR) [
7]. This provides the opportunity for struggling electricity generators, such as South Africa’s parastatal utility Eskom, to implement large-scale DSM techniques, such as ripple control, to control and limit the load on the electrical grid.
Xu et al. stated on DSM and DR strategies that “any reliable strategies aimed at controlling the demand from aggregated power usage of multiple water heaters should be based on accurate and representative energy models, which require a thorough and comprehensive understanding of the thermal behaviour for each individual water heater” [
7]. An accurate model of the EWH is needed to simulate its thermal behaviour and power usage. These are often termed “thermal” models because they primarily simulate the electrical energy supplied to the heater and the amount of energy extracted from the heater (exergy). The objective of using such models is firstly to simulate grid demand based on individual (usually statistical) hot water draws, and secondly, to allow centralised demand management to reduce the load on the grid without (a) sacrificing individual EWH energy savings, (b) increasing cold events experienced by the EWH user, and (c) prevent bacterial growth due to low temperatures in the tank [
3,
8]. From this, the load on the grid can be simulated and optimised when a network of these devices is connected. The scale of these simulations necessitates the implementation of thermal models that are computationally inexpensive.
One significant phenomenon that contributes to variations in energy consumption in an EWH is thermal stratification. It can be described as the vertical separation of water regions due to density differences. Since the density of water is relatively sensitive to temperature changes, lower density regions would tend to rise above regions that have higher densities. This means that water regions with higher temperatures (lower density) would rise to the top, while lower temperatures (higher density) would tend to rest at the bottom of the tank. This phenomenon is also known as the buoyancy effect and also occurs naturally in most large bodies of water, such as dams and rivers. The thermal stratification of the EWH device is an important factor to be considered when characterising the energy consumption of these devices [
9,
10,
11,
12].
There are many models and experimental data in the literature that aim to characterise EWH devices. Most are designed for vertically oriented tanks, and not many for horizontal tanks. EWHs are mounted horizontally in many developing countries, which makes it important to also focus on models representing horizontal tank orientation.
3. Experimental Platform Setup and Design
The experimental platform, shown in
Figure 1, is designed to control the variables that have an influence on the stratification in the tank, and consequently the energy characteristics of the EWH. Nel et al. [
28] observed that a 5% increase in the inlet water temperature can produce energy savings up to 13%, and 5% energy savings for the same change in ambient temperature. The results from Fernandez-Seara et al. [
18] and Castell et al. [
10] show that lower water flow rates improve the degree of stratification and thus improving exergy efficiency. Booysen et al. [
22] found that 29% energy savings can be achieved by using intelligent scheduled control and by lowering the thermostat set temperature. This is corroborated by Kizilors et al. [
24] who observed that discharge efficiencies increase when thermostat set points are lowered. Clearly, these environmental conditions, water usage characteristics and switching patterns of the heating element have a significant influence on EWH thermal behaviour and energy characteristics.
Therefore, to characterise EWH thermal behaviour, the platform is designed to emulate and control the ambient temperature using a climatic chamber, emulate the inlet water temperature using a controllable, in-line water heat exchanger, emulate the behaviour of the user by controlling water usage patterns and control the electrical switching frequencies and conditions for the heating element. Another significant objective of the experimental platform is data acquisition of the thermal stratification inside the tank, outlet water temperature for exergy analysis, environmental temperatures; ambient and inlet water, energy usage and water usage.
3.1. Environmental Emulation
The ambient temperature and inlet water temperature are two dynamic factors that make up the environmental conditions that influence the thermal behaviour and energy usage patterns of the EWH.
The ambient temperature is controlled using a custom climatic regulation chamber equipped with two 500-W air heaters, circulation fans and venting fans. These actuators work together to regulate a user-defined chamber temperature. The platform has the ability to thermally emulate ambient temperatures up to 50 °C.
The temperature of the inlet water is regulated using a servo-controlled, three-port mixer valve. One of the inlet ports to the mixer is connected to a separate 100-L EWH device that is placed inside a standard chest freezer. The insulation of this tank is completely removed to decrease the thermal resistance between the chest freezer cavity and the water inside the tank. Therefore, if cold water is required for the experiment, the element of the tank will remain off and the chest freezer will switch on and cool the water inside the tank. Conversely, if warm water is required, the heating element of the tank switches on and heats up the water while the chest freezer is off. The second inlet port of the mixer valve is connected to the municipal water supply. A PI (proportional/integral) control system is used to regulate the temperature at the outlet port of the mixer valve. The temperature set point used in this case is the desired inlet water temperature to the main 150-L EWH device under test.
3.2. User Emulation
The water usage patterns of the household user(s) have a substantial influence on the thermal behaviour of the EWH. Many studies focus on the development of stochastic water usage prediction algorithms, such as Heidari et al. [
29] and Ritchie et al. [
30]. A way to improve the accuracy of these models is to generate more water usage profiles and the corresponding EWH energy response from the platform developed in this study.
This is accomplished by regulating the volumes, frequencies, flow rates and times of water usage events with an actuator (an electric ball valve) and two sensors: a digital flow meter that produces a digital pulse after sensing a specific volume of passing water and an outlet temperature sensor. The pulses of the flow meter are accumulated using an interrupt service routine (ISR) on the controller. This accumulation is sampled and processed every second. The resolution of the volumetric flow sensor used is ±2.5 mL, 2.5 mL/pulse and the associated error is 2.5 mL/s.
The position of the electric ball valve is at the outlet of the tank to ensure that the EWH tank is pressurised to the regulated water line pressure when it is closed. In this setup, the water pressure is regulated to 100 kPa using a standard pressure regulator valve.
3.3. Electrical Utility and Thermostat Emulation
Most electric water heaters utilise a standard thermostat that controls the electrical input to the heating element based on a sensed temperature and a user-defined set temperature. If the sensed temperature within the tank is below the set temperature, the thermostat will deliver electrical energy to the heating element, causing the water to heat up. The heating element will remain on and only switch off when the sensed temperature of the water has reached the desired set temperature with a small hysteresis band.
The literature has shown that the thermostat set point has a significant influence on the EWH energy characteristics [
22] and the thermal stratification [
10]. The experimental platform therefore incorporates the use of an electronically controlled thermostat. A custom digital controller receives the sensed temperature from the thermostat and determines the state of the heating element based on the desired set temperature defined in the software. This allows for dynamic set point changes during an experiment and can be utilised by different heating control strategies.
The emulation of the electrical utility is also important. In South Africa, there are rolling blackouts, a severe DM strategy, that cause many EWH units to be without electricity for up to 4 h. A possible objective could be to observe the effects of these blackouts on the energy and stratification characteristics of the EWH. This is done by introducing an emulated power availability schedule in the controller.
Therefore, the platform’s controller does two software checks for changing the state of the heating element. The power availability schedule is checked to see if power is available from the emulated electrical utility and the sensed thermostat temperature is compared to the desired set temperature. This means that the heating element of the EWH could be off even when the sensed temperature is below the set temperature. This happens when there is no power available based on the emulated utility availability schedule.
In the event of a possible controller fault and for safety reasons, the thermostat has a hardware cut-off temperature of 90 °C. This ensures that the temperature within the sensing region of the thermostat never exceeds this temperature.
In addition to this, a safety protocol is established to prevent the mid-region of the tank to exceed a temperature limit of 5% above the EWHs current set point temperature in the case of a faulty thermistor reading. This is done by implementing the temperature reading of a few of the mid-region DS18B20 sensors in the EWHs element actuating logic.
3.4. Data Acquisition and Sensor Selection
An important objective of the platform is to sample experimental data that can aid in thermal and energy characterisation of the horizontal EWH. For this, the following measurements need to be taken: thermal stratification inside the tank using multiple internal temperature sensors, outlet water temperature for exergy analysis, environmental temperatures: ambient and inlet water, energy usage and volumetric flow rate of water.
In most studies, a uni-directional thermal stratification measurement strategy is used for vertical EWH experiments. However, the positional arrangement of the temperature sensors in this system allows for the measurement of vertical temperature variation and horizontal-radial temperature variation. In addition, the temperature variation along the length of the tank is also measurable. The three-dimensional sensor arrangement can be seen in
Figure 2. This arrangement allows for the measurement of stratification in multiple regions of the horizontal tank. Thus, the thermal variation inside the tank can be measured in a three-dimensional space.
The stratification measurement system is designed to be placed inside of the 150-L tank. The tank is initially sectioned near the one end, closest to the heating element. This is followed by a flange assembly installation that allows the tank to be bolted shut after the stratification measurement system is positioned and secured inside. The metal flange of the sectioned tank can be seen in
Figure 3.
The thermal stratification is measured using nine temperature sensors that are mounted on modified and waterproofed rectangular aluminium tubes. They are equally-spaced and vertically arranged through the centre of the tank. This number is chosen based on similar work done in the literature: Farooq et al. [
12] placed eight equally spaced temperature sensors along the height of a vertical tank wall, Fernandez-Seara et al. [
11,
18] mounted 11 equally spaced thermocouples into their vertical tank; the vertical tank used by Kepplinger et al. [
31] used 11 non-uniformly spaced thermocouples and Castell et al. [
10] used 6 equally spaced sensors positioned through the centre of their vertical tank. The resolution of the stratification data was suitable in all above-mentioned studies. The application of 9 sensors is therefore considered to be sufficient, especially since it is for a horizontally oriented tank.
The DS18B20 temperature sensor is specifically chosen based on its convenient digital 1-Wire communication protocol, sensing accuracy, measurement range and resolution. These sensors have a measuring range of −55 °C to 125 °C and have a rated measurement error of ±0.5 °C. However, for the range in which we measure, from 10 °C to 70 °C, the mean error in the datasheet is only ±0.2 °C. An added advantage of this sensor is that the 1-Wire communication protocol can support multiple devices on the same data line. This is ideal since there are fewer wires going into the pressurised tank. The design of the stratification measurement system implements eight separate data lines for the sake of modularity and fault detection convenience.
The sensor chip is encapsulated within a silicon-filled, stainless steel tube. The tube has a diameter of 6 mm with a length of 50 mm. The material properties of the stainless steel tube prevents the development of rust to a large extent. This is essential since these sensors would be exposed to heated and pressurised water for long periods of time.
3.5. Structural Design Considerations and Suitability
The thermal stratification measurement system is designed to withstand the thermal fluctuations and pressure of the internal environment of the tank. Normal tank operating pressures can range from 100 kPa to 600 kPa in South Africa. The water temperature can fluctuate dramatically with temperatures typically ranging from 20 °C to 70 °C in different regions of the tank. This makes it difficult to place electronic sensors in this type of environment.
The support frame material is selected to be aluminium since it is light-weight, corrosion-resistant in water and has suitable thermal properties, such as high thermal conductivity and lower specific heat capacity than that of water. The thermal conductivity and specific heat capacity of water are typically 0.598 W/mK and 4200 J/kgK, respectively. In contrast, the same properties of aluminium are typically 239 W/mK and 900 J/kgK, respectively. The thermal conductivity property is a parameter used to quantify how well a material conducts heat and affects the rate at which heat is transferred. This parameter, along with the exposed surface area, has a direct influence on the transient thermal response time of the material. This relationship is described by Fourier
’s law of heat conduction shown in Equation (
1),
where
is the rate of heat transfer through the material,
k is the thermal conductivity of the material,
is the exposed surface area subjected to the heat transfer and
is the temperature differential through the material in the x-dimension.
The aim of the frame is to provide rigid support for the sensors without influencing the true thermal response of the water inside the tank. Energy transfer between the water and the support frame material is inevitable. However, the aluminium frame does not act as a thermal reservoir which would cause slow thermal transients. The ability for water to store thermal energy is far greater than that of aluminium and conversely means that the aluminium tubes gain and lose thermal energy much faster than that of water. The presence of the aluminium thus has no significant impact on the slower thermal response of the body of water.
3.6. Digital System Control
A custom controller for the system is designed to incorporate the use of an “over-the-counter” microcontroller such as an Arduino Due. This device is responsible for the control of the experimental platform. The controller is responsible for the control of environmental emulation, user behaviour emulation, electrical utility emulation and EWH thermostat control. It is also responsible for the data acquisition of the platform.
The data sampling parameters are configured and sent to the controller prior to the experiment. These parameters include sampling frequency, duration of experiment, set ambient air temperature, set geyser thermostat temperature and the set water inlet temperature. In addition, there is also the option of providing two types of schedules for EWH heating element power availability and water usage patterns. After the controller receives the experiment parameters, it then starts to set up the EWH environment, such as the set ambient air temperature if required. Once the environment is ready, the experiment starts and the time-stamped data are recorded to an SD card and streamed to a computer via a serial port.
The physical and finalised version of the stratification measurement system is shown in
Figure 3. Each temperature bus module is connected to a three-core silicone cable that runs to the outside of the tank. For this to be possible, the existing anode rod had to be removed for the wires to have an entry point.
Figure 3 shows how the cables are positioned in the tank. It was important to minimise the disturbance of typical flow into the tank; therefore, the cables and the sensor busses were positioned appropriately to provide enough space from the inlet water diffuser.
3.7. Sensor Referencing and Geometry Considerations
A referencing convention is established to make sense of the recorded data from the test station. There are 67 temperature sensors in total that are positioned at specific locations inside of the tank.
Figure 2 shows the naming convention and sensor referencing notation for the sensors used for temperature variation measurement within the tank.
Most documented experimental stratification data from the literature are captured and analysed for vertically oriented tanks [
9,
10,
11,
12,
18,
19,
20,
31]. In these cases, the cross-sectional area as a function of height of the tank remains constant. This means that, assuming that the installed sensors are equally spaced, the volume of the water measured for each sensor node will be the same. This is not the case for a horizontally oriented tank since the cross-sectional area as a function of vertical height varies in a sinusoidal manner. The volume of water measured by the central node on the vertical plane (sensor position 4 on all sensor busses) is the largest as compared to the top and bottom layer volumes, which are equal and are the smallest.
This influences the way the average temperature of the tank along the vertical plane is calculated. Since the nodal volumes for each vertically positioned sensor are different, a weighted average calculation is more appropriate when the average temperature is required for analyses. Each nodal temperature that is measured in the vertical plane of the tank is multiplied by the volume of the measured region and divided by the total volume to act as a contribution factor. These 9 weighted values are added together to obtain the overall weighted average of the temperature in the tank at a specific point in time. The weighted average expression used is shown below in Equation (
2).
where
and
are the nodal volume and measured temperature at sensor
j, respectively.
For the physical design shown in
Figure 2, the nodal volumes for each sensed region in the vertical plane is tabulated below in
Table 1. It is important to note that when these volumes are used to calculate the weighted average of the vertical temperature variation in the tank, the assumption is made that the longitudinal and transverse temperature variation along the length, and along the horizontal radius of the tank is negligible and can be ignored. The volumetric-weighted temperature average is determined and visualised in
Section 4 for each longitudinal position in the tank.
5. Conclusions
Demand management strategies that rely on the direct control of electric water heaters, with their inherent potential for storing thermal energy, aim to reduce electrical load, reduce thermal losses, and to ensure user satisfaction (hot water) without fostering bacterial growth that occurs at lower temperatures. Their efficacy therefore relies on an accurate representation of the internal temperature distributions of these heaters. Although the internal temperature distribution of electric water heaters have been characterised for the vertical orientation, they have not been characterised for the horizontal orientation, which is more common in developing countries.
We presented a fully controllable platform with which the thermal behaviour of a horizontal water heater can be characterised under a wide range of environmental conditions (e.g., ambient temperature and inlet water temperature), control parameters (e.g., heating schedule, thermostat set temperature), and user behaviour (flow rate and draw volume). We presented the internal temperature measurements under static conditions (no water drawn) when the heater was heating from cold and kept at a set temperature. These were repeated under dynamic conditions (water drawn).
The static results show the dominant impact of stratification on the temperatures inside the tank, with the higher layers heating substantially faster than the lower ones. The thermal response of the three bottom layers were considerably slower than the layers immediately above. This is in contrast with what has been documented in literature for vertically oriented tanks [
10,
11,
12,
34]. It is believed that the location of the heating element and tank geometry have an influence on this recorded difference.
This thermal stratification also had an impact on the element switching frequency, which is higher shortly after the first heating cycle, and then slower as the bottom layers heat up. The longitudinal differences in stratification were observed to be minor, with only the middle layer varying by more than 2 °C from the inlet to the outlet. The differences between the two sides (left and right) of the water heater were mostly symmetrical for the middle and top layers. However, the inlet pipe, which is placed at a slight offset to the right, had a marked cooling thermal exchange impact on the right side of the bottom layer in static conditions, which extends even to the middle region along the longitudinal axis.
The dynamic results showed the introduction of the thermocline and its subsequent movement upwards in the tank to the top layer. Interestingly, the inlet region’s upper layer took the longest to cool down. An interesting phenomenon was observed in which the temperature at the outlet rose for a few minutes after valve was closed. We believe this may be a thermal exchange from the trapped heat in the fluid-driven vortex at the top layer, near the inlet side of the tank.
We have shown that there is more going on inside a horizontal water heater than what could be explained by considering vertical water heaters. We have shown that there exist interesting phenomena in the longitudinal and traversal axes, for both static and dynamic conditions. Finally, we have shown results that highlight the slight differences in thermocline dynamics between horizontal- and vertical tanks, specifically at the top and bottom layers, owing to different layer volumes and tank geometry. However, explaining these fully will require further work. Two examples of open questions are: the extent and exact mechanism by which a potential vortex affects the temperatures at the outlet, and the energy transfer characteristics between the inlet pipe coupling and the warm water in its vicinity.