# Simulating Peak Ice

This year ice in the tank was finally melted between March 5 to March 10 – as ‘visual inspection’ showed. Level sensor Mr. Bubble was confused during the melting phase; thus it was an interesting exercise to compare simulations to measurements.

Simulations use the measured ambient temperature and solar radiation as an input, data points are taken every minute. Air temperature determines the heating energy needed by the house: Simulated heat load is increasing linearly until a maximum ‘cut off’ temperature.

The control logic of the real controller (UVR1611 / UVR16x2) is mirrored in the simulation: The controller’s heating curve determines the set temperature for the heating water, and it switches the virtual 3-way valves: Diverting heating water either to the hygienic storage or the buffer tank for space heating, and including the collector in the brine circuit if air temperature is high enough compared to brine temperature. In the brine circuit, three heat exchangers are connected in series: Three temperatures at different points are determined self-consistently from three equations that use underground tank temperature, air temperature, and the heat pump evaporator’s power as input parameters.

The hydraulic schematic for reference, as displayed in the controller’s visualization (See this article for details on operations.)

The Coefficient of Performance of the heat pump, its heating power, and its electrical input power are determined by heating water temperature and brine temperature – from polynomial fit curves to vendors’ data sheet.

So for every minute, the temperatures of tanks – hot and cold – and the volume of ice can be calculated from energy balances. The heating circuits and tap water consume energy, the heat pump delivers energy. The heat exchanger in the tank releases energy or harvests energy, and the collector exchanges energy with the environment. The heat flow between tank and ground is calculated by numerically solving the Heat Equation, using the nearly constant temperature in about 10 meters depth as a boundary condition.

For validating the simulation and for fine-tuning input parameters – like the thermal properties of ground or the building – I cross-check calculated versus measured daily / monthly energies and average temperatures.

Measurements for this winter show the artificial oscillations during the melting phase because Mr. Bubble faces the cliff of ice:

Simulations show growing of ice and the evolution of the tank temperature in agreement with measurements. The melting of ice is in line with observations. The ‘plateau’ shows the oscillations that Mr. Bubble notices, but the true amplitude is smaller:

Simulated peak ice is about 0,7m3 greater than the measured value. This can be explained by my neglecting temperature gradients within water or ice in the tank:

When there is only a bit of ice yet (small peak in December), tank temperature is underestimated: In reality, the density anomaly of water causes a zone of 4°C at the bottom, below the ice.

When the ice block is more massive (end of January), I overestimate brine temperature as ice has less than 0°C, at least intermittently when the heat pump is turned on. Thus the temperature difference between ambient air and brine is underestimated, and so is the simulated energy harvested from the collector – and more energy needs to be provided by freezing water.

However, a difference in volume of less than 10% is uncritical for system’s sizing, especially if you err on the size of caution. Temperature gradients in ice and convection in water should be less critical if heat exchanger tubes traverse the volume of tank evenly – our prime design principle.

I have got questions about the efficiency of immersed heat exchangers in the tank – will heat transfer deteriorate if the layer of ice becomes too thick? No, according also to this very detailed research report on simulations of ‘ice storage heat pump systems’ (p.5). We grow so-called ‘ice on coil’ which is compared to flat-plate heat exchangers:

… for the coil, the total heat transfer (UA), accounting for the growing ice surface, shows only a small decrease with growing ice thickness. The heat transfer resistance of the growing ice layer is partially compensated by the increased heat transfer area around the coil. In the case of the flat plate, on the contrary, also the UA-value decreases rapidly with growing ice thickness.

__________________________________

For system’s configuration data see the last chapter of this documentation.

# Mr. Bubble Was Confused. A Cliffhanger.

This year we experienced a record-breaking January in Austria – the coldest since 30 years. Our heat pump system produced 14m3 of ice in the underground tank.

The volume of ice is measured by Mr. Bubble, the winner of The Ultimate Level Sensor Casting Show run by the Chief Engineer last year:

The classic, analog level sensor was very robust and simple, but required continuous human intervention:

So a multitude of prototypes had been evaluated …

The challenge was to measure small changes in level as 1 mm corresponds to about 0,15 m3 of ice.

Mr. Bubble uses a flow of bubbling air in a tube; the measured pressure increases linearly with the distance of the liquid level from the nozzle:

Mr. Bubble is fine and sane, as long as ice is growing monotonously: Ice grows from the heat exchanger tubes into the water, and the heat exchanger does not float due to buoyancy, as it is attached to the supporting construction. The design makes sure that not-yet-frozen water can always ‘escape’ to higher levels to make room for growing ice. Finally Mr. Bubble lives inside a hollow cylinder of water inside a block of ice. As long as all the ice is covered by water, Mr. Bubble’s calculation is correct.

But when ambient temperature rises and the collector harvests more energy then needed by the heat pump, melting starts at the heat exchanger tubes. The density of ice is smaller than that of water, so the water level in Mr. Bubble’s hollow cylinder is below the surface level of ice:

Mr. Bubble is utterly confused and literally driven over the edge – having to deal with this cliff of ice:

When ice is melted, the surface level inside the hollow cylinder drops quickly as the diameter of the cylinder is much smaller than the width of the tank. So the alleged volume of ice perceived by Mr. Bubble seems to drop extremely fast and out of proportion: 1m3 of ice is equivalent to 93kWh of energy – the energy our heat pump would need on an extremely cold day. On an ice melting day, the heat pump needs much less, so a drop of more than 1m3 per day is an artefact.

As long as there are ice castles on the surface, Mr. Bubble keeps underestimating the volume of ice. When it gets colder, ice grows again, and its growth is then overestimated via the same effect. Mr. Bubble amplifies the oscillations in growing and shrinking of ice.

In the final stages of melting a slab-with-a-hole-like structure ‘mounted’ above the water surface remains. The actual level of water is lower than it was before the ice period. This is reflected in the raw data – the distance measured. The volume of ice output is calibrated not to show negative values, but the underlying measurement data do:

Only when finally all ice has been melted – slowly and via thermal contact with air – then the water level is back to normal.

In the final stages of melting parts of the suspended slab of ice may break off and then floating small icebergs can confuse Mr. Bubble, too:

So how can we picture the true evolution of ice during melting? I am simulating the volume of ice, based on our measurements of air temperature. To be detailed in a future post – this is my cliffhanger!

>> Next episode.

# Ice Storage Hierarchy of Needs

Data Kraken – the tentacled tangled pieces of software for data analysis – has a secret theoretical sibling, an older one: Before we built our heat source from a cellar, I developed numerical simulations of the future heat pump system. Today this simulation tool comprises e.g. a model of our control system, real-live weather data, energy balances of all storage tanks, and a solution to the heat equation for the ground surrounding the water/ice tank.

I can model the change of the tank temperature and  ‘peak ice’ in a heating season. But the point of these simulations is rather to find out to which parameters the system’s performance reacts particularly sensitive: In a worst case scenario will the storage tank be large enough?

A seemingly fascinating aspect was how peak ice ‘reacts’ to input parameters: It is quite sensitive to the properties of ground and the solar/air collector. If you made either the ground or the collector just ‘a bit worse’, ice seems to grow out of proportion. Taking a step back I realized that I could have come to that conclusion using simple energy accounting instead of differential equations – once I had long-term data for the average energy harvesting power of the collector and ground. Caveat: The simple calculation only works if these estimates are reliable for a chosen system – and this depends e.g. on hydraulic design, control logic, the shape of the tank, and the heat transfer properties of ground and collector.

For the operations of the combined tank+collector source the critical months are the ice months Dec/Jan/Feb when air temperature does not allow harvesting all energy from air. Before and after that period, the solar/air collector is nearly the only source anyway. As I emphasized on this blog again and again, even during the ice months, the collector is still the main source and delivers most of the ambient energy the heat pump needs (if properly sized) in a typical winter. The rest has to come from energy stored in the ground surrounding the tank or from freezing water.

I am finally succumbing to trends of edutainment and storytelling in science communications – here is an infographic:

Using some typical numbers, I am illustrating 4 scenarios in the figure below, for a  system with these parameters:

• A cuboid tank of about 23 m3
• Required ambient energy for the three ice months is ~7000kWh
(about 9330kWh of heating energy at a performance factor of 4)
• ‘Standard’ scenario: The collector delivers 75% of the ambient energy, ground delivers about 18%.
• Worse’ scenarios: Either collector or/and ground energy is reduced by 25% compared to the standard.

Contributions of the three sources add up to the total ambient energy needed – this is yet another way of combining different energies in one balance.

Ambient energy needed by the heat pump in  Dec+Jan+Feb,  as delivered by the three different sources. Latent ‘ice’ energy is also translated to the percentage of water in the tank that would be frozen.

Neither collector nor ground energy change much in relation to the base line. But latent energy has to fill in the gap: As the total collector energy is much higher than the total latent energy content of the tank, an increase in the gap is large in relation to the base ice energy.

If collector and ground would both ‘underdeliver’ by 25% the tank in this scenario would be frozen completely instead of only 23%.

The ice energy is just the peak of the total ambient energy iceberg.

You could call this system an air-geothermal-ice heat pump then!

____________________________

Continued: Here are some details on simulations.

# Earth, Air, Water, and Ice.

In my attempts at Ice Storage Heat Source popularization I have been facing one big challenge: How can you – succinctly, using pictures – answer questions like:

How much energy does the collector harvest?

or

What’s the contribution of ground?

or

Why do you need a collector if the monthly performance factor just drops a bit when you turned it off during the Ice Storage Challenge?

The short answer is that the collector (if properly sized in relation to tank and heat pump) provides for about 75% of the ambient energy needed by the heat pump in an average year. Before the ‘Challenge’ in 2015 performance did not drop because the energy in the tank had been filled up to the brim by the collector before. So the collector is not a nice add-on but an essential part of the heat source. The tank is needed to buffer energy for colder periods; otherwise the system would operate like an air heat pump without any storage.

I am calling Data Kraken for help to give me more diagrams.

There are two kinds of energy balances:

1) From the volume of ice and tank temperature the energy still stored in the tank can be calculated. Our tank ‘contains’ about 2.300 kWh of energy when ‘full’. Stored energy changes …

• … because energy is extracted from the tank or released to it via the heat exchanger pipes traversing it.
• … and because heat is exchanged with the surrounding ground through the walls and the floor of the tank.

Thus the contribution of ground can be determined by:

Change of stored energy(Ice, Water) =
Energy over ribbed pipe heat exchanger + Energy exchanged with ground

2) On the other hand, three heat exchangers are serially connected in the brine circuit: The heat pump’s evaporator, the solar air collector, and the heat exchanger in the tank. .

Both of these energy balances are shown in this diagram (The direction of arrows indicates energy > 0):

The heat pump is using a combined heat source, made up of tank and collector, so …

Ambient Energy for Heat Pump = -(Collector Energy) + Tank Energy

The following diagrams show data for the season containing the Ice Storage Challenge:

From September to January more and more ambient energy is needed – but also the contribution of the collector increases! The longer the collector is on in parallel with the heat pump, the more energy can be harvested from air (as the temperature difference between air and brine is increased).

As long as there is no ice the temperature of the tank and the brine inlet temperature follow air temperature approximately. But if air temperature drops quickly (e.g. at the end of November 2014), the tank is still rather warm in relation to air and the collector cannot harvest much. Then the energy stored in the tank drops and energy starts to flow from ground to the tank.

On Jan 10 an anomalous peak in collector energy is visible: Warm winter storm Felix gave us a record harvest exceeding the energy needed by the heat pump! In addition to high ambient temperatures and convection (wind) the tank temperature remained low while energy was used for melting ice.

On February 1, we turned off the collector – and now the stored energy started to decline. Since the collector energy in February is zero, the energy transferred via the heat exchanger is equal to the ambient energy used by the heat pump. Ground provided for about 1/3 of the ambient energy. Near the end of the Ice Storage Challenge (mid of March) the contribution of ground was increasing while the contribution of latent energy became smaller and smaller: Ice hardly grew anymore, allegedly after the ice cube has ‘touched ground’.

Mid of March the collector was turned on again: Again (as during the Felix episode) harvest is high because the tank remains at 0°C. The energy stored in the tank is replenished quickly. Heat transfer with ground is rather small, and thus the heat exchanger energy is about equal to the change in energy stored.

At the beginning of May, we switched to summer mode: The collector is turned off (by the control system) to keep tank temperature at 8°C as long as possible. This temperature is a trade-off between optimizing heat pump performance and keeping some energy for passive cooling. The energy available for cooling is reduced by the slow flow of heat from ground to the tank.

# How Does It Work? (The Heat Pump System, That Is)

Over the holidays I stayed away from social media, read quantum physics textbooks instead, and The Chief Engineer and I mulled over the fundamental questions of life, the universe and everything. Such as: How to explain our heat pump system?

Many blog postings were actually answers to questions, and am consolidating all these answers to frequently asked questions again in a list of such answers. However, this list has grown quickly.

An astute reader suggested to create an ‘animation’ of the gradual evolution of the system’s state. As I learned from discussions, one major confusion was related to the role of the solar collector and the fact that you have to factor in the history of the heat source: This is true for every heat pump system that uses a heat source that can be ‘depleted’, in contrast to a flow of ground water at a constant temperature for example. With the latter, the ‘state’ of the system only depends on the current ambient temperature, and you can explain it in a way not too different from pontificating on a wood or gas boiler.

One thing you have to accept though is how a heat pump as such works: I have given up to go into thermodynamical details, and I also think that the refigerator analogy is not helpful. So for this pragmatic introduction a heat pump is just a device that generates heating energy as an output, the input energy being electrical energy and heat energy extracted from a rather cold heat source somewhere near the building. For 8kW heating power you need about 2kW electrical energy and 6kW ambient energy. The ratio of 8kW and 2kW is called the coefficient of performance.

What the typical intro to heat pumps in physics textbooks does not point out is that the ambient heat source actually has to be able to deliver that input energyduring a whole heating season. There is no such thing as the infinite reservoir of energy usually depicted as a large box. Actually, the worse the performance of a heat pump is – the ratio of output heat energy and input electrical energy, the smaller are the demands on the heat source. The Chief Engineer has coined the term The Heat Source Paradox for this!

The lower the temperature of the heat source, the smaller the coefficient of performance is: So if you run an air source heat pump in mid-winter (using a big ventilator) then less energy is extracted from that air source than a geothermal heat pump would extract from ground. But if you build a geothermal heat source that’s too small in relation to a building’s heating demands, you see the same effect: Ground freezes, source temperature decreases, performance decreases, and you need more electrical energy and less ambient energy.

I am harping on the role of the heat source as the whole point of our ‘innovation’ is our special heat source that has two components, both of them being essential: An unglazed solar / air collector and an underground water / ice tank plus the surrounding ground. The collector allows to replenish the energy stored in the tank quickly, even in winter: Air temperature just needs to be some degrees warmer than the cold brine. The tank is a buffer: When no energy is harvested by the collector at ambient temperatures below 0°C, water freezes and releases latent heat. So you can call that an air heat pump with a huge, silent and mainentance-free ‘absorber’ plus a buffer that provides energy for periods of frost and that allows for storing all the energy you don’t need immediately. Ground does provide some energy as well, and I am planning to post about my related simulations.(*) It can be visualized as an extension of the ice / water energy storage into the surroundings. But the active volume or area of ground is smaller than for geothermal systems as most of the ambient energy actually comes from the solar / air collector: The critical months in our climate are Dec-Jan-Feb: Before and after, the collector would be sufficient as the only heat source. In the three ‘ice months’ water is typically frozen in the tank, but even then the collector provides for 75-80% of the ambient energy needed to drive the heat pump.(*)

(*) Edit: This post written in 2017 show how much energy is stored / exchanged by each component. An overview of essential numbers is given here; emphasis on the volume of ice – which is compared to simulations here.

Components are off-the-shelf products, actually rather simple and cheap ones, such as the most stupid, non-smart brine-water heat pump. What is special is 1) the arrangement of the heat exchanger in the water tank and 2) the custom control logic, that is programming of the control unit.

So here is finally the series of images of the system’s state, shown in a gallery and with captions: You can scroll down to see the series embedded in the post, or click on the first image to see an enlarged view and then click through the slide-show.

More information on the system (technical data, sizing) and measurement data since 2012 can be found in this documentation – updated every few months.

Information for German readers: This post contains the German version of this slide-show.

# Ice Storage Challenge: High Score!

Released from ice are brook and river
By the quickening glance of the gracious Spring;
The colors of hope to the valley cling,
And weak old Winter himself must shiver,
Withdrawn to the mountains, a crownless king.

These are the first lines of the English version of a famous German poem on spring, from the drama Faust, by Johann Wolfgang von Goethe. Weird factoid about me: I was once inclined to study literature, rather than physics. But finally physics won, so this is a post about joyful toying with modeling heat transport in ice and water.

After 46 days we had a high score: The ice cube, generated by our heat pump, stopped growing at about 15m3. About 10mof water remained unfrozen. After the volume of ice had been in a steady state for a a while, we turned on the solar collector again to return to standard operations.

Where did the energy for the heat pump come from before?

The lid of the tank is insulated against ambient air, the solar collector was not operational, and no ice had been created: The remaining energy has to be provided by the 5th element that cannot be shut off: 1) water 2) ice, 3) ambient air, 4) solar radiation … 5) ground.

Normally ground supplies about 15 W per m2 surface area – deduced from monitoring the power transported with the brine flow and energy accounting for the tank. The active interface between tank and ground below frost depth is about 35 m2. This results in about 0,5 kW in total, thus just 12 kWh per day, much lower than the ~ 50 kWh ambient energy fed into the heat pump.

After much deliberation and playing with the heat transfer equation we came up with this description of the evolution of the ice cube:

Phase 1: Growth of ice into water.

• Ice starts to grow from the heat exchanger tubes into the remaining water. These tubes are installed in a meandering pattern, traversing the storage tank.
• At some point the thick layers of ice covering adjacent parts of the pipes touch each other. The surface of this solid ice cube is smaller than the interface between the meandering ice formations and water before. The power needed by the heat pump has to be pushed through a smaller surface – which is only possible if the temperature gradient within the ice gets larger. As the temperature at the ice-water interface has to be 0°C, the temperature at the heat exchanger has to decrease. This is exactly what we see from monitoring data – brine temperature drops well below 0°C.
• Side-effect: Due to the lower brine temperature the coefficient of performance  decreases slightly. So more of the total heating energy needs to be provided by the electrical input. We call this the heat source paradox: The worse performance is, the more you spare the energy stored in the heat source. Thanks to this self-protection mechanism, the energy in the tank will not suddenly drop to zero.

Evolution of the volume of ice, ambient temperature and brine temperature over time. The ice is now quickly melting again – in March the collector is already harvesting enough energy again for balancing heating demands!

Phase 2: Ice touching ground.

• As long as there is some water between ice and ground, the water temperature is 0°C. This is the temperature ground ‘sees’ and the temperature which is relevant for the low heat transport from ground to water.
• Ice touches some surfaces of the cuboid tank – the ones where the heat exchanger tubes are closest to the surface. Now ground is directly connected to ice with its temperatures < 0°C. The temperature gradient between ground and ice provides for a higher flow of energy. This is also indicated by the evolution of the temperature in the ground below the tank: While temperatures of undisturbed ground and the region below the tank had been aligned before, ground temperature beneath the tank still kept getting lower – although a few meters away from the tank ground is already warming up again.
• If enough heat is delivered by ground, no more heat is needed by freezing the remaining water in the tank. When ground temperature reaches zero, it can even freeze – which happens with geothermal systems, too. We might have extended the ice storage into ground.

The temperature sensor closest to the tank in 30 cm underneath. A few meters from the tank ground is already warming up again (following the standard ‘yearly temperature wave’), but below the tank the temperature is still getting lower – as highlighted by the blue rectangle.

Heat transport within ice is actually more efficient than transport in water: Ice has 4 times the heat conductivity of water, and 10 times the thermal diffusivity. The latter is a measure for the time a deposited ‘lump of heat’ will be spread in space:

So we have built a very efficient cold bridge between the heat exchanger and ground. Everything is consistent with the poetry of the differential equation of heat transfer.

I marvel at the intriguing and mathematically appealing physics in my backyard!

For the backyard (‘Office desk farming’). Happy Easter everybody!

# The Ice Storage Challenge

The more we enjoyed our spring-like winter, the more we were worried if we will ever see much ice in our underground water tank this heating season.

So we did what I had announced – we switched off the solar collector completely: Since February 1st our heat pump has been extracting heat energy from the tank only.

Ice started to grow from the heat exchanger pipes into the bulk of the tank, our re-purposed root cellar. These images have been taken at Day 9; showing the part of the pipes above water level:

The density of ice is about 10% lower than the density of water, and the supporting construction needs to make sure that this home-grown iceberg does not start to float. The visible tip of a floating iceberg would be equivalent to those 10% of the total volume. If the whole bulk of ice has to remain under water, the water level increases gradually. There is a tiny heater in the right place, to prevent the formation of a continuous slab of ice that would trap the not-yet-frozen water.

From the change in water level we determine the volume of ice:

The volume of the tank is about 27m3, so it can hold about 25m3 of ice. Now a bit more than half of it is frozen. About 0,3m3 of ice has been created per day, and the rate is slowing down as the weather gets milder and milder. We need less than 70kWh heating energy per day, for space heating and hot water.

The water level has now raised above the top part of the tubes shown in the previous image – and some interesting new structures have emerged which are not directly related to heat pump operations (Image taken at Day 25):

If it rains heavily, some water trickles down into the former cellar which still uses the original ‘ceiling’ built decades ago. The rain water hits the bulk of ice and is frozen as well.

Next part of this story: High Score!

_________________________________________

More data: Our documentation of measurement data now covers two full seasons (since autumn 2012) plus this season until end of January.