Simulative Investigation of Thermal Capacity Analysis Methods for Metallic Latent Thermal Energy Storage Systems

: Latent heat storage systems are a promising technology for storing and providing thermal energy with low volume, mass and cost requirements, especially when operated at high temperatures. Metallic phase change materials are particularly advantageous for high thermal input and output, which is especially important for mobile applications. When designing a storage system, it is essential to have precise knowledge about the potential storage capacity. However, the system’s storage capacity is typically calculated from material properties determined at lab scale, although systemic boundary conditions can have a considerable inﬂuence. Systemic inﬂuences can result from thermal and reactive interfaces or from the storage design. In order to consider these inﬂuences, we propose three calorimetric procedures to thermally analyse high-temperature metallic latent energy storage systems at an application scale. We examined the procedures in a transient simulation environment, monitoring the storage capacity of the system. The procedure, based on adiabatic conditions, shows the least deviation from the simulation input parameters, but is limited to the heating process of the storage. Discharging the storage can be represented by isoperibolic conditions with controlled heat exchange. The precision of the procedures depends on the evaluation routine, the calibration routine, the heat extraction rate and the thermal inertia of the test bench.


Introduction
As the electrification of vehicles increases, the challenge of ensuring thermal management without excess waste heat from combustion gains importance [1][2][3]. Utilising the traction battery to cover heating demand leads to reduced vehicle range in winter conditions. The driving range can be maintained under these conditions by adding thermal energy storage (TES) in addition to the electrochemical energy storage. Carefully designed TES can save mass, space, cost and rare raw materials compared to a larger battery [4,5]. Thermal energy can be stored latently in a phase transition of a phase change material (PCM) at constant temperature. The capacity of latent thermal energy storage (LTES) can be extended by using additional storage of thermal energy with sensible heating (including the heat stored, as materials change temperature). In LTES, thermal energy can be stored in a smaller temperature range with a higher energy density than when using sensible heat storage alone. Commonly used PCMs like salt or paraffin are of interest due to their high storage density, but they suffer from low thermal conductivity [6,7]. In vehicle applications, the storage system needs to be charged and discharged especially fast. Thus, PCMs with high thermal conductivity are preferred. Because of their high thermal conductivity compared to other kinds of PCMs, metallic PCMs (mPCMs) are particularly well suited to this application [8][9][10].
Storage with high thermal input and output can be termed thermal high-performance storage (THS). Kraft et al. first proposed a THS design for a battery electric subcompact car with a defined reference scenario for ambient temperatures up to −20 • C [5]. The design was based on the mPCM Al-12wt% Si, with an operating temperature up to 600 • C. Charging power of 11 kW and thermal output of 5.1 kW were suggested. The resulting system mass of the THS was 50 kg. The study showed that the designed THS could increase the vehicle's range in winter conditions by up to 26.7% for the chosen driving scenario. An even higher potential is expected for vehicles with higher heating demand, like buses [11].
For the adequate design of a THS system, thermophysical properties of the mPCM must be known precisely. Thermophysical properties provide information about the amount of energy the material can store and are typically determined via calorimetric methods. Experiments with small samples are often performed via differential scanning calorimetry (DSC). However, this method has some disadvantages, as the sample size is far from real application [12][13][14]. One essential disadvantage is the deviation of the ratio of surface area to volume of the storage material in a real application. Due to the high reactivity of some liquid mPCMs, the influence of the chosen crucible material for DSC measurements must not be underestimated. The severity of the reaction can depend on the surface-to-volume ratio of the mPCM sample [15,16]. A particular mechanism may become relevant when the reaction process involves diffusion-controlled processes. Those depend on the concentration gradient, which can vary for different sample dimensions. It is conceivable that a very small sample will reach a saturation limit of the solute, consuming less container material compared to a large sample.
Heat input and output in laboratory-scale experiments have different heat flux orientation and magnitude compared to the application scale. This may have an influence on the natural convection of liquid mPCM, which can be a relevant factor in damage to container material by erosion and abrasion. Another effect of a deviating heat flux situation is solidification of the mPCM, which leads to a characteristic microstructure with volume defects. It is not clear whether such phenomena could have a significant influence on damage processes or cycling stability, and consequently on the evolution of thermophysical properties.
One promising approach to assess thermophysical properties of larger PCM samples is the T-history method, first proposed by Yinping [12]. The basic principle is to log the temperature change of the PCM in a tube upon a sudden temperature change in the surroundings. The heat transfer between the PCM and the surroundings is determined from a twin setup sample with known properties. The setup was improved in order to reduce temperature gradients in the sample, which arise especially for low thermal conductivity PCMs such as salts, salt hydrates and paraffin. The T-history method has become widely accepted and used [14,17,18].
The T-history method was developed especially to characterise large PCM samples with volumes up to several millilitres with a rather simple setup. In this study, the desired test volume is in the range of several litres. A twin setup, as used in the T-history method, would result in a clumsy, space-consuming setup. It will always be necessary to build a second component completely identical to the test sample but with known properties and stability over the planned temperature program. In order to reduce the material and spatial effort, a calorimetric setup without a twin sample is favoured by the authors.
In the recent past, many studies related to lab-scale or application-scale thermal analysis of LTES systems were published. However, most of them focused on low-or mid-temperature non-metallic PCMs like paraffin or salt hydrates. Xu et al. [19] found deviations in the storage capacity determined by T-history and full-scale tests which were performed with sodium acetate trihydrate-based PCM in a temperature range up to 61 • C. Klimes et al. [20] reported on the influence of experimental techniques and conditions on the behaviour of PCMs. However, their review article only discusses non-metallic PCMs in the low or medium temperature range. Zauner et al. [21] characterised hybrid sensible-latent heat storage based on thermal oil and high-density polyethylene. Their large-scale test was modelled via Dymola and validated experimentally. Heat losses were not determined, probably due to their minor relevance in the examined temperature range up to 150 • C. Jansone et al. [22] found limits of the T-history method when determining the specific heat capacity of PCMs in the mid-temperature range. Waser et al. [23] developed a computational model for salt-hydrate based PCMs used in a low to medium temperature range. The model was validated experimentally and was found to be an important tool for properly designing subcomponents like heat exchangers.
Large-scale models and experiments are also necessary for high-temperature applications. However, in those temperature ranges, heat loss may play a larger role and should be considered carefully. Pirasaci et al. [24] developed a model to find critical design aspects for LTES using non-metallic PCMs with a melting temperature of around 500 • C. The model was validated experimentally and heat losses were considered through the energy balance of the heat transfer fluid. A full-scale thermal analysis of mPCM at high temperatures was performed by Rea et al. [25,26]. Heat flow rates for a storage system with 100 kg of Al-12wt%Si were calculated based on many assumptions. A final statement about the model's accuracy was not given. The dependency of heat flow rates and efficiency of the application scale (because of the surface-to-volume ratio) and the thermal cycling conditions was reported. They also mentioned some practical issues when designing a real scale container, such as corrosion issues. In their experimental setup, heat losses were not measured, which was explained by deviating losses of the test setup compared to a real system. An accurate assessment of losses of a real system was recommended.
To the authors' knowledge, none of the known thermal analysis techniques were adapted to application scale for mPCM at high temperatures for vehicle application with the adequate treatment of heat loss in the energy balance equation. In this study, we propose different system-calorimetric procedures for this purpose and analyse them via transient simulation for their accuracy to determine the storage capacity of metallic hightemperature THS. These novel procedures are the first to be applied to metallic LTES at application scale at high temperatures with proper heat loss consideration.

Simulated Specimen
The specimen simulated in the characterisation procedure potentially includes subcomponents like mPCM, container and heat transfer components. These heat transfer components are located in a central cylindrical tube, as proposed by Kraft and Klein Altstedde [4]. This study was performed with a sample mass of m s = 5.688 kg. The chosen storage material was the eutectic alloy Al-12wt%Si, which was operated in a temperature range between 100 • C and 600 • C [4]. The storage material choice was made according to the storage concept proposed by Kraft and Klein Altstedde [4].

Methods
The aim of calorimetric characterisation of THS is to investigate the relevant thermophysical properties of the storage system. The properties are specific enthalpy of fusion h fus , specific heat capacity c p (ϑ) and melting temperature ϑ fus . The specific state variables are gained by normalising to the sample mass m s . The extrinsic values can be derived from the course of sample heat Q S over sample temperature ϑ S . Figure 1 shows a typical course of Q S for a sample with first-order phase transition and its context with the mentioned relevant thermophysical properties. Here, only systems with a sharp melting temperature-not a melting temperature range-are considered (e.g., pure metals or eutectics). In general, it is possible to define the sample, not only as the storage material, but as the full system, including further system components like the housing. Hence, the boundaries of the sample can be defined depending on the matter of interest.
In Figure 1, the ideal course of Q S is shown for a sample with a first-order phase transition. In Figure 2, it is qualitatively compared to a demonstrative "measured" Q Scurve. Deviations from the ideal course result from a temperature gradient in the sample and test setup, which in turn means that the entire sample volume is not passing through the phase transition at the same time. This effect becomes more dominant the bigger the sample is. However, this effect is not expected to be very severe in this case, as metallic samples with high thermal conductivity minimise the thermal gradient [27].  The heat capacity C p (ϑ1: ϑ2) in a temperature range from ϑ 1 to ϑ 2 can directly be deduced from the slope of Q s (ϑ): The specific heat capacity c p can be derived from the relation to the sample mass m s via c p = C p /m s . As specific heat capacity is a temperature-dependent quantity, it is recommended to compare c p values of different samples for similar temperature ranges. Enthalpy of fusion H fus and melting temperature ϑ fus are derived geometrically by the routine shown in Figure 3. This is a standardised routine described in DIN 51007. Q S is extrapolated linearly from the temperature ranges for C p determination closest to the melting temperature (range for lower temperature: a; range for higher temperature: b) to ϑ fus . This linear course of Q S before the phase transition step is called the initial base line. ϑ fus is defined as the intercept between the extrapolated initial base line and the inflection tangent of the step. The equation to determine H fus is: The specific enthalpy of fusion h fus is determined via h fus = H fus /m s . In order to identify the calorimetric principles which cover all operational boundary conditions, a literature review mainly based on calorimeter classification by Sarge et al. [28] was performed. The characterisation procedure should be able to test the charging and discharging conditions of the storage. Heat output and input should be controllable in order to investigate the influence of the cooling and heating rates. Principles that are not suitable are listed below:

•
Isothermal calorimeters can only be operated at a defined temperature, thus it is not possible to test the whole operating temperature range of the specimen system. Depending on the heat path from the sample to the surroundings, the calorimeter signal can show a thermal lag [28,29]. Therefore, scanning calorimeters are not recommended for use with large sample sizes [30] and inhomogeneous samples, which is how the THS specimen system can be classified [31].
Suitable calorimetric principles are the isoperibol calorimeter with controlled heat exchange for both charging and discharging conditions and adiabatic calorimeters for charging conditions only. Three measurement procedures to realize these principles are proposed in this study. In order to keep the setup complexity as low as possible, procedures without a reference setup were used.
All proposed measuring procedures are based on solving a balancing equation of all heat amounts, resulting in the heat accumulation in sample Q S at a certain sample temperature ϑ S , referring to the chosen reference temperature ϑ Ref . ϑ S is logged over time t, making it possible to describe either Q S (ϑ S ) or Q S (t). Here, we chose to use the time-dependent values. The balancing equation, which can be used for every procedure, is: where Q H is the heat input of the electrical heater, Q C is the thermal energy in the components of the test bench, Q Loss is the heat loss to the surroundings and Q WF is the heat output by the working fluid. The difference between the procedures is in how the different contributions are determined. An overview is given in Table 1. m Isoperibolic with cooling correlation Zero Measure (ϑ C − ϑ C,0 ) Correlate Q Cool to ϑ ex Q C is always determined by measuring the difference of the components' instantaneous temperature ϑ C and initial temperature ϑ C,0 . With knowledge about the heat capacity of all test setup components C C , the heat put into the components is described by: The next sections describe the proposed measuring procedures in more detail.

Stepwise Adiabatic Procedure
In order to realise adiabatic conditions, the latent energy storage system is thermally isolated. Perfect adiabatic conditions, that is, where no heat is exchanged with the surroundings, cannot be achieved in reality. Thus, physical implementations are technically always isoperibolic, albeit with potentially very small heat loss [28,32].
Thermal energy is supplied to the thermally isolated sample via an electrical heater. If supply losses are neglected, heat input Q H can be determined by electrical heating power P H , which is supplied in the time interval between t i and t i+1 : Heat loss Q Loss cannot be directly measured but can be determined via calibration. To achieve this, the heater is power-controlled to hold the sample temperature constant for an equilibration time t eq until a steady state is reached. Under steady state conditions, the total heating rate equals the heat flow into the setup components and the heat loss rate to the surroundings . Q Loss . The heat loss rate is then correlated to the surface temperature via a second-degree polynomial. The heat loss to the surroundings Q Loss is finally received via integration over time. For this study, t eq was chosen to be 45 min for each temperature range. The temperature step between 500 • C and 600 • C, t eq , was chosen to be 60 min, as the phase change of the examined material Al-12wt%Si occurs in this temperature range.

Isoperibolic Procedure with Loss Correlation
In this procedure, the amount of thermal energy stored in sample Q S is determined by measuring the heat balance of the working fluid which is cooling the sample. This procedure belongs to the class of isoperibol calorimeters with controlled heat exchange. Therefore, the sample is again thermally isolated, but this time, the working fluid is conducted through the central tube of the sample system. At first, the sample is heated to a maximum temperature by an electrical heater. Then, the heater is shut off and the flow of the working fluid is started. The heat transferred from the sample to the working fluid Q WF is determined by measuring the temperature difference of the working fluid before passing sample ϑ in and after passing sample ϑ out . The mass flow rate . m of the working fluid is set constant and the specific heat capacity of the working fluid c p,WF is known. Q WF can be determined via the following equation: Q Loss is determined via calibration and correlation during the heating phase, as described for the stepwise adiabatic procedure. Thus, calibration and measurement can be performed concurrently. For this procedure, t eq = 30 min was chosen for calibration. The mass flow rate . m was set to 50 kg/h.

Isoperibolic Procedure with Cooling Correlation
This procedure also belongs to the class of isoperibol calorimeters with controlled heat exchange. Again, the sample is isolated thermally and a working fluid is conducted through the centre of the sample. This time, the heat transferred to the working fluid Q WF and the heat loss to the surroundings Q Loss are not determined separately, but are combined in the expression Q Cool = Q WF + Q Loss . Calibration is performed in order to correlate Q Cool to the temperature of the heat exchanger ϑ ex in the heat transfer zone of the sample to the working fluid. Therefore, the power of the electrical heater P H is logged for different steady state temperature steps during heating. According to the other calibration methods, the following balance is true: Q Cool . Cooling rate . Q Cool is then correlated to ϑ ex . It cannot be distinguished which amount of thermal energy is put out by the working fluid or by the loss to the surroundings. To keep the correlation as precise as possible, the working fluid's mass flow . m and the surrounding conditions need to be constant for both calibration and measurement. After the calibration and heating phase, the electrical heater is shut off, but the working fluid is kept flowing. Then, ϑ ex (t) is logged over the cooling phase. For this procedure, t eq = 30 min was chosen for calibration with active working fluid flow.

Experimental Setup
A test bench was designed to perform the proposed measurement procedures. In Figure 4, a schematic drawing of all parts relevant for the simulation is shown. The specimen for this study was a thermal energy storage consisting of an AISI 321 containment and the mPCM Al-12wt%Si. The container had the shape of a hollow cylinder with rounded edges and an empty mass of 3.3 kg, and was filled with m s = 5.688 kg Al-12wt%Si cast into it. The heating unit was a hollow cylindrical copper piece with bore holes, into which 8 electric cartridge heaters were inserted within the central tube of the containment. Each cartridge heater had a maximum electrical power of 500 W. Within the centre of the heating unit, air could flow through the specimen for discharging. In order to enable high heat transfer, a heat transfer structure with fins was present in the hollow shape of the copper cylinder. The components of the air duct through the centre of the storage were thermally decoupled from the rest of the setup by ceramic fittings. The thermal energy storage component was surrounded by two layers of microporous thermal insulation to reduce thermal losses. PID controllers controlled the air mass flow and the heating power.

Simulation
For a transient simulative investigation of the measurement procedures, the Dymola simulation environment and the Modelica language were used to implement appropriate models. Sources for these models were the Modelica Standard Library v. 3.2.2 and Modelon Base Library v. 2.6 from Modelon AB, as well as in-house models, e.g., for the mPCM. As a solver, Esdirk45a (explicit singly diagonally implicit Runge-Kutta) [33] was used in this work.
For investigation of temperature gradients within components, several cells were assembled, where each cell consisted of a capacity and 4 heat conductors. This enabled heat conduction in 2 spatial dimensions. The discretisation level of each component was chosen in relation to the expected temperature gradient. The discretised simulation model is schematically displayed in Figure 5. Due to the rotational symmetry of most components, this approach results in a rotationally symmetric 2-dimensional model. For some test bench components, like the outer insulation layer, rotational symmetry is not true, which leads to a simplification of these parts. Due to their large lateral dimension compared to their small thickness, one-dimensional heat conduction through the outer insulation layer to the surroundings was also assumed. The PID controllers for air mass flow control and heating power control were also included in the simulation model. In order to detect heat flows of interest authentically, virtual heat flow sensors were established in the simulation model. The copper cylinder with the dual function of heater and heat transfer body was geometrically simplified to a hollow cylinder. The geometry of the heat storage component itself was simplified to a cylindrical shape with a rectangular cross-section.
For the investigation of measurement procedures, different measured variables had to be logged, as outlined in Section 2.2. In the simulation model, only sensor positions that are accessible in the real setup were used. All sensors were assumed to operate ideally error-free. An overview of the measured variables and their corresponding sensor positions within the setup is shown in Table 2. The indicated positions can also be found in Figure 5. Table 2. Temperature sensors used to measure relevant temperatures for procedure evaluation (for cylindrical components, logarithmic average was used). In italics are positions of sensors as marked in Figure 5.

Stepwise Adiabatic Procedure
The course of sample temperature ϑ S over time t for the stepwise adiabatic procedure is shown in Figure 6a. The temperature plateaus show the six equilibration steps where the heat loss to the surroundings was determined by measuring the heating power. The melting process of the storage material can be observed at an additional temperature plateau at 576.5 • C. The course of Q S (ϑ) resulting from Equations (4) and (5) with temperature and power outputs of the simulation is shown in Figure 6b as well as the course of Q S (ϑ) calculated from the simulation input values for c p , ϑ fus and h fus . The slope of Q S (ϑ), which corresponds to c p , is observed to be lower in the simulation output than the input property. Thus, the overall amount of energy stored in the studied temperature range is underestimated.

Isoperibolic Procedure with Loss Correlation
The course of sample temperature ϑ S over time t for the isoperibolic procedure with loss correlation is shown in Figure 7a. The heating process for calibration and the cooling process for measurement can be seen. The course of Q s resulting from simulation and the input course of Q s are shown in Figure 7b. The simulated course of ϑ s again shows the equilibration steps for heat loss determination and reveals the solidification process with an additional plateau at 576.5 • C. Upon further cooling, the temperature decreased more slowly and approximated room temperature asymptotically. The simulated course of Q s revealed that c p (ϑ) was not constant, although this was assumed for the input value. This interesting observation is elaborated further in the discussion. Again, the overall energy amount stored in the studied temperature range is underestimated. In contrast to the stepwise adiabatic procedure, this time Q s resulted from a cooling half-cycle. This can be seen at the sharp edge at higher energy amounts, where solidification starts, and the diffuse transition to sensible cooling, as a solidification of the whole sample, takes some time.

Isoperibolic Procedure with Cooling Correlation
The course of sample temperature ϑ S over time t for the isoperibolic procedure with cooling correlation is shown in Figure 8a. The heating process for calibration and the cooling process for measurement can be seen. The course of Q s resulting from simulation and the input course of Q s are shown in Figure 8b. The simulation of ϑ S shows a similar course as in the procedure with loss correlation. The simulated course of Q S was nearly linear in the sensible range. However, the deviation of c p in the solid range was visibly larger compared to the other procedures, and the energy stored in the studied temperature range was more greatly underestimated.

Deviation of Simulation Input Values
In order to evaluate the three proposed measuring procedures, simulated measured values were compared to the simulation input values for c p (ϑ), ϑ fus , h fus and the overall stored energy in sample Q S between 100 • C and 600 • C. The input values were determined via differential scanning calorimetry. An overview of the thermophysical properties generated by the different procedures in comparison to the input values is shown in Figure 9, with indication of the absolute deviation (primary ordinate) and relative deviation E (secondary ordinate). The isoperibolic procedure with loss correlation shows the lowest error regarding the determination of c p , but the largest deviation for h fus . However, the isoperibolic procedure with cooling correlation shows contrasting behaviour. The assignment of stored energy to the sensible or the latent regime depends on the choice of extrapolation limits and the width of the now diffuse phase transition (compared with Figure 3). The stored heat in the whole temperature range Q s 100:600 is a better indication for the accuracy of the procedure, as it is not as affected by these evaluation artefacts. Regarding this quantity, the stepwise adiabatic procedure has the lowest deviation from the simulation input value. However, this procedure is only suitable for analysing heating processes, not half-cycles with active cooling.

Discussion
In the simulation, all sensors were assumed to operate ideally error-free. The results of this study therefore only reflect the deviations resulting from the courses and setups of the different procedures themselves and numerical inaccuracy. In a real test setup, additional uncertainties of measurement from the sensors have to be considered.
Phase change characteristics: The melting step resulting from simulation is not as sharp as in the input. The reason could be a temperature gradient in the storage material upon melting. For the procedures where the cooling half-cycle is evaluated, the liquidus temperature can be determined more precisely and the solidus temperature is very diffuse. For the stepwise adiabatic procedure, only the heating half-cycle is evaluated, therefore the solidus temperature can be observed well but the liquidus temperature is diffuse.
Correlation temperature: The correlation temperature used for calibration has a significant systematic influence on the output. Thermal inertia in the system can lead to an incorrect estimation of heat loss if an unsuitable correlation temperature is chosen. Depending on the position of the correlation temperature, large thermal inertia can lead to the situation that no thermal equilibrium between correlation and sample temperature is reached after equilibration time during heating. Furthermore, a completely different thermal gradient can arise during cooling, so that the correlation function is no longer transferable from the heating to the cooling half-cycle.
Equilibration time: A further significant influence on the quality of calibration is the time chosen to reach steady state conditions during calibration. All procedures underestimated c p , which could result from a temperature gradient in the sample setup, meaning that steady state was not reached in the chosen time t eq . Thus, Q Loss was overestimated and c p underestimated. The deviation of c p for t eq = 45 min was smaller than for a tested alternative with t eq = 30 min, as the state between the temperature steps was closer to steady state conditions. This suggests that increasing t eq leads to increased accuracy in the determination of c p and especially heat losses. Increasing the time step allows large thermal inertias to equilibrate, but this increases experiment time. An optimum between correlation temperature and equilibration time should be determined when commissioning physical experiments. In the simulation, thermal inertia only results from bulk thermal resistance, whereas in the experiment, thermal contact resistance also needs to be considered. Thermal contacts are strongly geometry-dependent and are affected by thermal expansion and the properties of layers or coatings at interfaces. For storage systems with poor thermal contact, the accuracy of the procedures can be reduced if the equilibration time is not increased.
Heat loss: Another important aspect that influences the accuracy of the result is the proportion of heat lost in the total energy balance Equation (3). For the stepwise adiabatic procedure, deviations between simulation input and output values increased with increasing temperature as the ratio between lost thermal energy Q Loss and stored thermal energy Q S increased. For the isoperibolic procedure with loss correlation, the biggest deviation between simulated and input c p occurred at low temperatures, where the cooling rate is slowest. The slower the cooling, the faster Q WF decreases compared to Q Loss and the larger the proportion of Q Loss in the total balance. The accuracy of this procedure might be increased by increasing the portion of Q WF in the total balance with an increased mass flow rate of coolant. Similar behaviour was observed for the isoperibolic procedure with cooling correlation. However, in this procedure, the working fluid mass flow cannot be increased without limit, as the thermal output is not allowed to be larger than the thermal input realisable by the installed electrical heaters.
Numerical inaccuracy: Potential numerical inaccuracy can depend on the chosen time increment used for simulation and the residual error of the solver. These did not have considerable influence.

Conclusions
Three measuring procedures for the thermophysical characterisation of THS utilising mPCM at the system scale were derived and described. A test setup was designed and the procedures were examined through simulation. All sensors were assumed to operate ideally error-free so that deviations resulted only from the course of each procedure and potentially numerical inaccuracy. The most suitable procedure to determine the overall stored heat in the specimen system was the stepwise adiabatic procedure. However, this procedure is limited to heating processes and cannot be used to test half-cycles with active cooling. The accuracy of isoperibolic procedures involving either loss or cooling correlation depends on several parameters, including the extrapolation limits, the correlation temperature, the equilibration time, the mass flow of the working fluid and the thermal contact between thermal subsystems. The procedures could be optimised by adjusting these parameters through a sensitivity analysis. For a final evaluation of the proposed procedures, it is necessary to additionally assess errors that arise from the sensors and devices. The procedures outlined in this work provide a first step for commissioning an experimental apparatus. This device may be used to validate the absolute outputs of the simulation and check the reproducibility of the proposed procedures.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.