Issue 
Renew. Energy Environ. Sustain.
Volume 6, 2021
Achieving Zero Carbon Emission by 2030



Article Number  27  
Number of page(s)  7  
DOI  https://doi.org/10.1051/rees/2021026  
Published online  16 August 2021 
Research Article
Datadriven Modelling of Thermal Solid Sorption Storage Systems
University of Applied Sciences Upper Austria, Energy Research Group ASIC, Wels, Austria
^{*} email: harald.kirchsteiger@fhwels.at
Received:
28
June
2021
Received in final form:
19
July
2021
Accepted:
19
July
2021
An approach fofsr datadriven modelling of open sorption storage systems using zeolite as storage material is presented. The overall dynamic simulation model has the inflow air stream mass flow and absolute humidity as inputs and computes the outflow air temperature. The model is subdivided into several components, where dynamic state space and process model identification techniques are applied. A comparison of the proposed modelling technique with simulated data from a validated model based on first principles shows that a reasonable accuracy − for a model application in temperature control systems design − can be obtained. It was found that using the proposed strategy, only a limited number of experiments are required, thus saving experimental time. Moreover, the computational requirements for a simulation using the proposed model are greatly reduced compared to a simulation model where differential equations discretised in time and space must be solved.
© H. Kirchsteiger and N. DaborerPrado, Published by EDP Sciences, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
As per December 2020, EU decided to cut down greenhousegas emissions by 55% (with respect to 1990) by the year 2030 [1]. Among the sectors responsible for a large share of the emissions are industry, electrical power generation, transportation, agriculture, and heating and cooling in commercial and residential buildings. In order to reach the ambitious goals, a significant emission reduction through technological innovation and political incentives in all the individual sectors is essential.
In the heating and cooling sector for commercial and residential buildings, technologies for renewable generation are available and well established. However, one major drawback is the intermittency and seasonal gap between high production in summer and high demand in wintertime of solardriven energy production. An overproduction during periods of increased irradiation may not be required and is therefore wasted, while the undersupply must be covered by potentially nonrenewable energy sources. Therefore, several techniques for longterm solar energy storage were analysed in the literature [2]. Among those techniques proposed, thermochemical storages are a very promising alternative [3] since they can significantly reduce the required space for the storage because of increased energy density of the material itself.
Among thermochemical storages, sorption storage systems were gaining significant attention lately [4]. Many different combinations of two materials can form a sorption pair (for example silica gel and water vapor) and many different process concepts were proposed in the literature, see for example [5]. Arguably the most straightforward design of a sorption storage system is an open storage (in contrast to a closed storage under vacuum) as it was demonstrated for example in [6,7]. The basic and simplified working principle of sorption storages using the material combination water and zeolite is illustrated in Figure 1. The storage is charged by supplying heat to the sorption material which is getting dried in this process. This operating mode is called desorption. In a dry state, the sorption material can keep the stored energy without losses unless it is getting in contact with watervapor. To discharge the storage, a controlled amount of water vapor is supplied to the material which in turn releases high amounts of heat. This operating mode is called adsorption. The whole process is repeatable all over again, however, thermomechanical stability of the zeolite might become an issue to consider [8]. Since there are no heat losses when kept in a dry state, seasonal storage becomes feasible [9].
Besides the properties of the material itself, the precise quantities of heat released per second depend mainly on the following three properties of the air stream which is guided through a storage compartment filled with the material: mass flow, inflow temperature and inflow humidity. The relation among those quantities is not a static one because it is well known that sorption materials have a limited storage capacity for water and moreover the uptake rate is nonlinear [10]. In order to use sorption storages in the best possible way in future renewable energy systems, mathematical models which numerically describe these relations are of great relevance, especially for modern modelbased control algorithms and energy management systems.
While the literature on models for the adsorption process of single grains of material is extensive, the scientific literature on dynamical models of whole sorption storage systems is rather scarce. Models based on physical first principles were developed in [11–13]. A neural network model was identified in [14]. As with all physically motivated models, parameters have to be adapted to the particular system under observation through measurements.
The approach for modelling presented in this paper is starting straight away from measured data to obtain an input/output relation of temperature and humidity of the air flow through the storage. Instead of relying on material properties, all parameters are estimated from data and are not interpretable as physical quantities.
Fig. 1 Basic working principle of solid sorption storage. 
2 Methodology
Ultimately, the goal of this study is to derive a mathematical relation between the temperature T_{IN} and humidity h_{IN} of an inflowing airstream into a fixed bed sorption storage schematically shown in Figure 2, and the temperature T_{OUT} and humidity h_{OUT} of the outflowing airstream. Another variable having a major influence on the output quantities is the mass flow of the air. We are considering zeolite 13X as storage material in this study.
It is worthwhile to mention that from a systems theory point of view this system is a switching system (the dynamics are different in adsorption and desorption mode) and it is nonlinear (one nonlinearity being a saturation of the water uptake of storage material, which depends in level on the input variables).
The general strategy of the datadriven model development is subdivided into 4 stages:

Obtain informative measurement data on a single adsorption experiment.

Combine several experiments and develop a submodel for the zeolite water saturation.

Develop a submodel for the output temperature.

Combine the submodels to an overall simulation model.
In the study presented here, measurement data was obtained from a physically motivated simulation model [12] which was already validated against a realworld storage setup. In this way, more information on the required amount of data and on experimental conditions can be learned without the need of expensive and timeconsuming lab measurements since the very beginning. An important sideresult of the investigations presented here are guidelines for the exact procedure and setup of lab experiments which will make it possible to save considerable amounts of experimental time.
Fig. 2 Fixed bed reactor shown without additional appliances (right) and schematic view (left). 
2.1 Measurement data
The three main influence factors on the storage behaviour which can be easily adjusted and measured are the temperature, absolute humidity, and mass flow of the incoming air stream. It is assumed the temperature T_{IN} is always 20 °C, which is in agreement with the intended use of the storage as a heat generator for domestic buildings, driven by inhouse ambient air. Additionally, a modest shift of this temperature would mainly result in a corresponding shift of the outflow temperature if the storage were properly insulated.
A single experiment conducted started with a welldefined initial water saturation of the storage material x_{S} (0) given in kg water per kg zeolite dry mass. Two randomly chosen constant values inside of the ranges shown in Table 1 were attributed to the quantities and h_{IN}. The experiment was stopped when the gradient of the water saturation was smaller than a preset threshold value d_{xs}. A total of N_{exp} experiments were conducted.
Symbols and values used in the simulation.
2.2 Water saturation submodel
The water saturation submodel maps the input quantities and h_{IN} to the intermediate quantity x_{S}. The general model structure assumption is a first order differential equation, which is given in state space form as(1) where the coefficients a, b_{1}, b_{2}, c, d_{1}, d_{2} are to be estimated from data and the state x_{k} is of no particular interest. The observed nonlinearity between the inflow humidity and the maximum storage material saturation was compensated with a piecewise affine input transformation of the form(2)if h_{IN} is within a region defined by h_{IN,i−1} < h_{IN} ≤ h_{IN,i+1}. The coefficients α_{i} and β_{i} are again estimated from data. For the estimation process, an increasing number of experiments was processed, with the aim to find the minimum number of required experiments for a reasonable model quality. The identification of the state space model was done with the prediction error method implemented in the System Identification Toolbox [15].
2.3 Output temperature submodel
The model structure assumption for the temperature submodel is a parametervarying first order process model with two inputs of the form(3)where K_{1}, τ_{1}, K_{2}, τ_{2} are the unknown parameters to be estimated from data and s is the Laplace variable. Note that except for K_{1}, all parameters depend on the mass flow and are thus not constant parameters, making this model parameter varying. The estimation strategy for those parameters was a twostep procedure. First, individual process models were estimated on every single set of experimental data giving N_{exp} different estimates for the parameters K_{1}, τ_{1}, K_{2}, τ_{2}. While K_{1} was found to be constant, the other parameters showed a strong dependency on the mass flow. Thus, second, polynomial dependencies of those parameters on the mass flow were estimated: (4)
The input quantity u_{1} is derived from a linear regression of the inflow humidity with regards to the observed maximum outflow temperature T_{OUT,max} (5)where γ_{1}, δ_{1} are again parameters to be estimated, and the input quantity u_{2} is the time derivative of the previously modeled storage material water saturation(6)
2.4 Overall simulation model
The overall model is simply formed by combining all the submodels and polynomial expressions defined above. Graphically, this is represented in Figure 3. For comparison purposes and evaluation, this model was compared to the simulation model presented in [12].
A flowchart showing the process of model estimation is given in Figure 4. The linear regressions were solved using standard LeastSquares approaches while the state space model and process model were estimated using sophisticated tools available from the System Identification Toolbox [15]. For the branching criterion, the average FIT value (see Eq. (7)) was computed on all available experiments, not only on the ones which were used for model parameter estimation.
Fig. 3 Overall model structure. 
Fig. 4 Flowchart of the model estimation process. 
3 Data
Measurement data was obtained from the simulation model [12] using the following set of parameters defined in Table 1. A total of N_{exp} = 30 individual experiments were recorded, where each single one has a random value of and h_{IN} within the range as defined in the table. The recorded data also include information on the storage material saturation x_{S}, which can be obtained by gravimetric analysis in a lab.
4 Results
In this section graphical results, as well as numerical results using the FITvalue (7) as qualitative criterion how well an estimate coincides with the signal y to be modeled with mean value are presented. A FIT value of 100% indicates perfect agreement between model and measurement while a value of 0% means the model is only able to correctly reproduce the mean value of the measurement. Also, negative values are possible.(7)
4.1 Results on the water saturation submodel
The first N = 10 experiments were used to determine the piecewise affine relation in equation (2). It was found that two linear pieces are sufficient to describe the relation, as it is show in Figure 5.
For the estimation of the state space model, sequentially an increasing number N of experiments was used. Clearly, using only one dataset N = 1 results in a good model fit for this particular case, but on the remaining 29 datasets the model performance will be not acceptable. In order to find a reasonable number of required datasets for estimation, N was sequentially increased, until the average fit value (7) on all the experiments used for identification did not change significantly anymore.
It was found that this is already the case for five experiments, where an average fit value of 88.99% is obtained for the datasets used for estimation and an average fit of 89.49% for the 25 independent datasets. Note that, in general, fit values on an independent validation dataset not used for model estimation are less than fit values on identification datasets. Using just one dataset for estimation, the fit is 93.86% compared to an average of 85.86 on the validation datasets.
4.2 Results of the output temperature submodel
The linear regression model equation (5) was estimated to γ_{1} = 3296, δ_{1} = 20.35. While the polynomial dependencies of the estimated model parameters of the linear parameter varying model structure equation (3) cannot be seen in Figure 6, they can be appreciated in Figure 7. It was found that a polynomial of third degree fits the data well. Note that the parameter K_{1} is identified to be 1, regardless of the specific dataset used. Also note that once all the polynomials are fixed, model structure (3) is uniquely defined by the input air mass flow.
Fig. 7 Polynomial dependencies of the process model parameters in equation (3). Stars indicate measured datapoints, dashed line is the estimated polynomial of 3rd order. 
4.3 Results of the overall simulation model
All the estimated models were combined together (see Fig. 3) and simulated in parallel to the physically motivated simulation model from [12]. Once again, the input values for and h_{IN} were chosen randomly within the intervals previously defined in Table 1. The final result of the output temperature can be seen in Figure 8 and shows a good agreement between the proposed databased model and the reference signal.
Fig. 8 Comparison between databased model and simulation data as reference: zeolite saturation x_{s} (left) and outflow air temperature T_{OUT} (right). 
4.4 Remarks
It has to be mentioned that extrapolation is critical for databased models. The presented model is therefore only valid within the limits of the mass flow and humidity as defined in Table 1. It has to be mentioned also that the presented model is only considering the adsorption process. In the next step, the desorption process has to be analysed in detail. Also, the methodology is to be applied on the realworld experimental setup.
5 Conclusion
An entirely databased modelling approach for solid sorption based thermal storage systems using zeolite as storage material was presented. The approach was tested against data from a previously validated simulation model. The endresults demonstrate a good agreement between the model output quantity which is the produced outflow temperature. Another finding is that only 5 experimental datasets are required to have sufficiently informative data that covers the entire region of interest. Due to its low computational requirements, the data driven model is well suited for implementation on a programmable logic controller used to control the operation of sorption storages. Moreover, the proposed model enables modelbased control design approaches, for example model predictive control, which can significantly enhance the overall system performance. Compared to standard PID (proportional integral derivative) type approaches, they offer the possibility to inherently consider saturations and constraints and optimize the cyclic operation of the storage.
Acknowledgments
This project is financed by research subsidies granted by the government of Upper Austria.
References
 European Commission, Stepping up Europe's 2030 Climate Ambition − Investigating in a ClimateNeutral Future for the Benefit of our People, Brussels, 17.09.2020 [Google Scholar]
 P. Pinel, C.A. Cruickshank, I. BeausoleilMorrison, A. Wills, A review of available methods for seasonal storage of solar thermal energy in residential applications, Renew. Sustain. Energy Rev. 15, 3341–3359 (2011) [CrossRef] [Google Scholar]
 G. Krese, R. Kozelj, V. Butala, U Stritih, Thermochemical seasonal solar energy storage for heating and cooling of buildings, Energy Build. 164, 239–253 (2018) [CrossRef] [Google Scholar]
 D. Aydin, S.P. Casey, S. Riffat, The latest advancements on thermochemical heat storage systems, Renew. Sustain. Energy Rev. 41, 356–367 (2015) [CrossRef] [Google Scholar]
 Y. Zhang, R. Wang, Sorption thermal energy storage: concept, process, applications and perspectives, Energy Storage Mater. 27, 352–369 (2020) [CrossRef] [Google Scholar]
 B. Zettl, G. Englmair, G. Steinmaurer, Development of a revolving drum reactor for opensorption heat storage processes, Appl. Therm. Eng. 70, 42–49 (2014) [CrossRef] [Google Scholar]
 B. Zettl, H. Kirchsteiger, An open sorption heat storage application, in Proc. Int. Sustainable Energy Conference, Graz, Austria, 2018, pp. 605–611 [Google Scholar]
 F. Fischer, W. Lutz, J.C. Buhl, E. Laevemann, Insights into the hydrothermal stability of zeolite 13X, Micropor. Mesopor. Mater. 262, 258–268 (2018) [CrossRef] [Google Scholar]
 B. Fumey, R. Weber, L. Baldini, Sorption based longterm energy storage − process classification and analysis of performance limitations: a review, Renew. Sustain. Energy Rev. 111, 57–74 (2019) [CrossRef] [Google Scholar]
 D.D. Do, Adsorption Analysis: Equilibria and Kinetics (Imperial College Press, London, 1998) [CrossRef] [Google Scholar]
 G. Engel, S. Asenbeck, R. Köll, H. Kerskes, W. Wagner, W. van Helden, Simulation of a seasonal, solardriven sorption storage heating system, J. Energy Storage 13, 40–47 (2017) [CrossRef] [Google Scholar]
 N. DaborerPrado, Modeling and Simulation of an Innovative Domestic Sorption Storage System, M.S. Thesis, University of Applied Sciences Upper Austria, Wels, Austria, 2019 [Google Scholar]
 B. Mette, H. Kerskes, H. Drück, Experimental and numerical investigations of different reactor concepts for thermochemical energy storage, Energy Proc. 57, 2380–2389 (2014) [CrossRef] [Google Scholar]
 L. Scapino, H.A. Zondag, J. Diriken, C.C.M. Rindt, J. van Bael, A. Sciacovelli, Modeling the performance of a sorption thermal energy storage reactor using artificial neural networks, Appl. Energy 253, 113525 (2019) [CrossRef] [Google Scholar]
 The MathWorks, R2019a, System Identification Toolbox, Natick, Massachusetts, US [Google Scholar]
Cite this article as: Harald Kirchsteiger, Nayrana DaborerPrado, Datadriven Modelling of Thermal Solid Sorption Storage Systems, Renew. Energy Environ. Sustain. 6, 27 (2021)
All Tables
All Figures
Fig. 1 Basic working principle of solid sorption storage. 

In the text 
Fig. 2 Fixed bed reactor shown without additional appliances (right) and schematic view (left). 

In the text 
Fig. 3 Overall model structure. 

In the text 
Fig. 4 Flowchart of the model estimation process. 

In the text 
Fig. 5 Piecewise affine input transformation equation (2). Stars indicate measured datapoints. 

In the text 
Fig. 6 Individually (using one dataset at a time) estimated model parameters of equation (3). 

In the text 
Fig. 7 Polynomial dependencies of the process model parameters in equation (3). Stars indicate measured datapoints, dashed line is the estimated polynomial of 3rd order. 

In the text 
Fig. 8 Comparison between databased model and simulation data as reference: zeolite saturation x_{s} (left) and outflow air temperature T_{OUT} (right). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.