1. Introduction
Approximately 8% of the worldwide electric energy consumption is caused by the food industry for refrigeration and maintaining the cold chain [
1]. Much of this is attributable to refrigerated transport by road, whose cooling units are generally exposed to harsher conditions and have much poorer energy efficiency than stationary refrigeration units [
2]. Therefore, for both economic and environmental reasons, substantial efforts are devoted to reducing the energy consumption of refrigerated vehicles [
2,
3,
4]. However, the energy efficiency of refrigerated vehicles is still compromised by improper temperature control concepts, such as commonly used heuristic controllers. This is particularly evident for last-mile transport, where frequent door openings and the resulting heat ingress from the outside stress temperature control due to the high cooling demand required to cool down the cooling chamber after the door is closed [
5].
Already in the design stage, the correct dimensioning of the individual components of the refrigerated vehicle is of great importance. For example, Maiorino et al. [
6] studied the optimal design of the battery and photovoltaic panels of a hybrid refrigerated van. However, particular attention should be paid to the sizing of the cooling system, as this is one of the most critical components in terms of energy consumption and meeting the vehicle’s cooling needs. Generally, refrigeration units are oversized for steady-state operation due to sizing for peak cooling loads [
7,
8], leading to increased energy consumption, wear of the compressor, and acquisition costs. In the literature, the sizing of the cooling equipment in stationary conditions is broadly discussed, e.g., for HVAC systems in buildings [
9,
10] and a refrigeration system in an industrial plant [
11]. Although already some simulation-based design methods for optimizing refrigeration systems exist [
12,
13,
14], little has been done regarding the sizing of the cooling equipment of small-scale refrigerated vehicles with the consideration of door openings.
Therefore, this work analyzes the sizing of a secondary loop refrigeration unit for a small-scale refrigerated vehicle together with a mixed-integer model predictive controller (MPC). For this purpose, closed-loop simulations of a refrigerated vehicle with differently sized thermal energy storage (TES) capacities of the secondary loop are conducted. In the simulations, the system is disturbed by multiple door openings. Furthermore, the same simulations are performed using a rule-based proportional–integral (PI) controller, commonly used in industry for temperature control of refrigerated vehicles [
15]. The two controllers and the differently sized secondary loops are compared based on their energy consumption and time required to bring the cooling chamber temperature back to the reference after door openings.
The use of TES has already found application in many fields [
16], and it has been shown that they can significantly increase the efficiency of energy systems, e.g., photovoltaic panels [
17]. TES systems are also very promising in the cold chain industry [
18]. Cold TES can either use latent heat mediums like phase change materials (PCM) or sensible heat mediums. Storage systems with latent heat mediums generally have a much higher energy density than systems with sensible heat mediums, while TES systems using sensible heat mediums are more straightforward to realize, require less maintenance, and are less expensive [
16], making them more suitable for small-scale refrigerated vehicles. In [
19], the parameters of a PCM-based TES in a refrigerator cycle are analyzed. The results show that a suitably sized TES can significantly improve the cooling system’s efficiency compared with a system without a TES. However, the analysis only included the cooling unit with the TES, without any cooling chamber. Jeong et al. [
20] proposed a cooling system that stores low-temperature liquid refrigerant by actively controlling the refrigerant mass flow rate for the evaporator. Experimental results show that this allows the cooling capacity to be temporarily increased for peak cooling demand, thereby allowing a smaller cooling unit design.
Refrigeration systems with TES are especially beneficial for small-scale refrigerated vehicles, which often face door openings, resulting in high peak cooling demand [
5]. Systems for road refrigeration can solely be based on TES, which are charged at the base, or be combined with a conventional cooling unit [
3]. Mousazade et al. [
21] studied the thermal performance of a cooling system with only a TES for a 6-ton refrigerated truck. In their experimental evaluation of three different PCM materials, the cooling chamber of a stationary truck could be kept solely by the TES at a constant temperature for a maximum of 5.1 h. However, cooling systems solely relying on a TES have a relatively high mass and practically no setpoint change capability. The storage capacity of PCM-based TES operated in parallel with a conventional cooling unit is analyzed in [
14]. The influence of the sizing is given by means of closed-loop simulations with a heuristic control concept, although door openings were not taken into account. Shafiei et al. [
22] used a similar cooling system architecture for their control concept for a refrigerated truck, which can achieve considerable energy savings by future load predictions and utilizing the TES. This paper also shows that only a suitable control structure can unleash the full potential of such combined systems, but their control concept lacks the inclusion of door openings.
Due to the high global warming potential of HFC refrigerants, increasingly more environmentally friendly liquids are being used as refrigerants, such as propane, ammonia, or hydrocarbon, which have the major drawback of being flammable. As a result, alternative refrigeration technologies, such as secondary loop refrigeration systems [
23], are becoming more popular for physically separating hazardous refrigerants from the cooling chamber by adding an additional loop to the primary cooling cycle. Both single-phase (e.g., glycol mixtures and hybrid nanofluids [
24]) and two-phase fluids (e.g., CO
2 [
25], paraffin emulsion [
26], and salt hydrates [
27]) are used as secondary refrigerants. Advantages of this architecture are reduced primary refrigerant charge and leakage and the ability to use the secondary loop as a TES [
23]. Storage management of the TES allows a more flexible operation of the refrigeration unit without limiting the cooling of the cooling chamber, which means that the refrigeration unit can be operated generally more energy efficiently.
Due to the aforementioned reasons, TES systems with sensitive heat mediums are better suited for small-scale refrigerated vehicles. Therefore, this work considers a secondary loop cooling system with a glycol mixture as the secondary refrigerant, which is always kept above its freezing point. The modeling of a refrigerated vehicle with such a cooling system is described by Fallmann et al. [
28]. In their work, a dynamic low-order model was estimated by a gray-box modeling approach using measurement data from a small-scale refrigerated vehicle with door openings. Based on that work, a control concept that takes advantage of the storage capability of the secondary loop was proposed and evaluated on a test bed in [
29]. The experimental results show significant performance benefits regarding both energy consumption and peak cooling capacity for cooling after door openings. However, in the literature, the sizing of the secondary loop storage capacity was not studied until now.
The remainder of the paper is structured as follows: First,
Section 2 describes the model of the small-scale refrigerated vehicle with a secondary loop cooling unit. Next,
Section 3 explains the variation of the secondary loop storage capacity in the model. The two controllers to regulate the temperature inside the cooling chamber are introduced in
Section 4.
Section 5 summarizes the parameters of the simulation study, followed by the results in
Section 6. The paper concludes with a discussion and conclusion in
Section 7 and
Section 8, respectively.
2. Model Description
For modeling the cooling chamber and the cooling unit of the small-scale refrigerated vehicle, the model proposed by Fallmann et al. [
28] was adopted since it shows good agreement with measurement data from a real-world refrigerated vehicle. The only modification to the original model was the removal of the heater acting as a disturbance heat flow in the cooling chamber, as this disturbance is irrelevant to this work. In the following, a brief overview of the model is given, and the interested reader is referred to [
28], where the modeling of the system and its validation is explained elaborately. Note that the thermodynamic quantities in this work correspond to ISO 80000-5:2019 [
30].
Figure 1 shows the system scheme with the cooling unit and cooling chamber, as well as important modeling variables. The cooling unit is a commercially available secondary loop refrigeration unit [
31], which can be divided into a cooling loop and a storage loop. The cooling unit can be switched on and off by
, activating/deactivating the compressor of the cooling loop, the condenser fan, and the glycol pump of the storage loop. The cooling loop is a standard vapor compression refrigeration cycle with propane as refrigerant, powered by a compressor rotating with the speed
. The storage loop, containing glycol, can store thermal energy supplied by the cooling loop and release thermal energy to the air in the cooling chamber through an air chiller. The glycol temperatures at the inflowing and outflowing position of the air chiller are given by
and
, respectively. A fan mounted to the air chiller evokes a heat flow,
, either due to natural or forced convection between the secondary loop and the air inside the cooling chambers, depending on the status of the fan
. The interior of the cooling chamber, with the lumped air temperature
, is separated from the environment by insulated walls and a door. A second-order system describes the heat transfer characteristics of the insulated walls with the two wall temperatures
and
. Door openings are indicated by
, which entails a heat flow between the air inside the cooling chamber and the ambient air (ambient air temperature
). Furthermore, the model describes the power consumption of the compressor
, the condenser fan
, the glycol pump
, the air chiller fan
, and the total consumption of all components
.
The mathematical description of the dynamic model relies on first principles, which is explained comprehensively in [
28]. Given that the model comprises both continuous and binary variables, it is a hybrid model [
32]. Hence, a model formulation with a mode selector and a switched affine system is chosen. This approach allows to separate binary and continuous variables. The mode selector maps the three binary variables of the model on one of the eight modes
, see
Table 1.
Depending on the mode, an associated affine model is selected in the switched affine system. This model comprises solely continuous variables and can be written in state-space formulation as
where
denotes the continuous time. This state-space model describes the system dynamics by transforming the state vector
according to
with the continuous input
and disturbance
given by
and
, respectively, into the model outputs
according to
The transformation is defined by the system matrix
, input vector
, system disturbance vector
, affine system vector
, output matrix
, feedthrough vector
, output disturbance vector
, and affine output vector
. The model matrices of the state space systems and their parameters were adopted from [
28], where they were experimentally identified and validated. Their values are given in
Appendix A.
5. Simulation Setup
A simulation-based approach was used to study the refrigeration system’s performance with differently sized secondary loop storage capacities, a methodology frequently used in the literature to design other cooling systems [
12,
13,
14]. Matlab/Simulink [
41] was used as the simulation environment, where the continuous-time model described in
Section 2 and the discrete-time controllers are simulated in a closed loop. The parameters of the simulation study and the individual simulations are given in
Table 4.
A total of 19 differently sized storage loops were examined by varying the scaling factors between and 4. Each storage loop size is simulated seven times with different disturbances due to door openings, for both the MPC and the rule-based PI controller. This results in 266 individual simulations.
Each simulation has a length of
. Within this period, the system is disturbed by four door openings, lasting between 1 and 4 min. Because the timings of these switching disturbances strongly influence the results, the starting time of each of the door openings is randomly shifted by
, normally distributed random numbers in the interval (−125 s, 125 s). These random door opening shifts are generated seven times for each of the simulations with the same secondary loop sizing and controller (values listed in
Appendix C). In this manner, the controllers and the differently sized storage loops are evaluated with the same seven different door opening timings.
Logistics information systems are widely used in the industry for route planning for refrigerated vehicles. These systems also provide information about the time of door openings, when cargo is loaded or unloaded. In this work, it is assumed that these timings are perfectly accurate. Hence, the door openings in the simulation are precisely embedded in the MPC predictions of disturbances. Furthermore, the profile of the temperature window is derived from those predicted door openings. The temperature window is inactive during the duration of the door openings and the period T2TWmax after the door is closed. Further, it is assumed that the full state measurement or state reconstruction is available without uncertainty or measurement noise for this simulation study.
Each simulation starts at a steady state, where the cooling chamber temperature equals the reference temperature
, and the cooling unit and fan are active with a closed door. The ambient temperature is 22 °C for all simulations, which corresponds to the year-round temperatures in Central Europe for which the vehicle under consideration was configured. Furthermore, the adopted model [
28] was identified at similar ambient temperatures and thus is most accurate in this outside temperature range.
6. Results
In this section, the results of the simulation study are presented. First, the MPC and the rule-based PI controller are compared using two simulations with equally sized secondary loops. Then, the results of all simulations are statistically evaluated, and the optimal storage loop sizing is determined.
6.1. Controller Evaluation
Figure 4 shows the results of two exemplary simulations with the same secondary loop sizing.
The two uppermost graphs show the simulated temperatures of the cooling chamber and the thermal energy storage. Below, the power consumption of the compressor, the condenser fan and glycol pump, and the air chiller fan are displayed, allowing to conclude on the respective control variables selected by the controllers. Both controllers hold the cooling chamber temperature during steady-state operation without door openings inside the temperature window, marked in blue. However, after the door openings, indicated by the orange background shading, the MPC can cool down the cooling chamber significantly faster than the rule-based PI controller. Before the door openings, the MPC increases the compressor speed but does not activate the air chiller fan. As a result, thermal energy is stored in the secondary loop, lowering the glycol temperature. When the door is closed again, the MPC activates the air chiller fan and cools the cooling chamber quickly due to the large temperature difference between the glycol and the air inside the cooling chamber. The rule-based PI controller lacks information about future door openings and thus does not react in advance.
A difference between the two controllers is also evident in the energy consumption,
, which can be derived from the power consumption according to
In the simulation shown in
Figure 4, the MPC consumes 693 Wh, and the rule-based PI controller 795 Wh. The power consumption graphs also indicate that the most energy consumption is attributable to door openings. Stationary operation with a closed door requires significantly less power due to the good insulation of the walls. Furthermore, the computational time of the closed-loop simulation with the MPC is significantly longer than with the rule-based PI controller due to the complexity of the optimization task. In fact, the single simulation with the MPC takes 62.7 min, while with the rule-based PI controller, it only takes 3.24 s.
6.2. Simulation Study
For the statistical evaluation of the simulations, two performance parameters are introduced to improve the comparability of the results. The energy consumption of each simulation is related by
to the reference energy consumption according to
where
is the mean energy consumption of the seven simulations with the MPC and
.
The relative time to reach the target temperature after door openings
is defined according to
where
is the time to reach the upper bound of the temperature window
after each door opening.
Additionally, the scaling factor of the storage loop was generalized to make the results of this work applicable to cooling chambers with different dimensions by defining the dimensionless relative thermal storage capacity of the secondary loop,
, given by
where
is the scaled storage capacity of the glycol loop, and
the thermal capacity of the air inside the cooling chamber with the model parameter
(value, see
Appendix A).
Figure 5 shows the two performance parameters for each simulation depending on the relative capacity of the storage loop. A line and background shading highlight the mean and standard deviation of the simulations for each relative thermal storage capacity.
For both controllers, the energy consumption decreases as the storage capacity of the secondary loop increases. However, the energy consumption of the rule-based PI controller is generally around 20% higher compared with the MPC. Further, the relative time for reaching the temperature window increases with increasing the storage capacity of the glycol loop. As shown in the individual simulations above, the rule-based PI controller takes about twice as long to cool down the cooling chamber after a door opening compared with the MPC. Furthermore, at very small thermal storage capacities of the secondary loop, the power consumption and the time to reach the temperature window increase significantly with the MPC. In fact, the energy consumption is around 20% higher compared with the slightly larger storage capacities with the same controller.
The amplitude statistic of the glycol temperature for different storage capacities is shown in
Figure 6.
With the MPC, the glycol temperature is much lower at small thermal storage capacities and spans a broader temperature range. Moreover, the system is more affected by the minimum glycol temperature at low thermal storage capacities. However, with the rule-based PI controller, the glycol temperature remains relatively constant over the entire range and has a small variance. This also shows that this conventional temperature controller cannot utilize the thermal energy storage due to its simple heuristic control law in which the cooling unit and fan are operated in parallel.
To facilitate the sizing of the secondary loop storage capacity, the parameter
is introduced, which gives the controller performance according to
where
and
are the means of
and
for the different secondary loop storage capacities of the individual simulations, and
is a weighting factor between 0 and 1. The weighting factor can either prioritize low energy consumption or fast cooling after door openings in the secondary loop storage capacity design.
Figure 7 shows
for three different weighting factors and indicates the optimal sizing of the secondary loop with the MPC.
With , the design exclusively focuses on minimal energy consumption, while for , only the time to reach the temperature window after door openings is considered. When both targets should be equally satisfied, with , a clear optimum between and is observable for the MPC. The rule-based PI controller performs much worse than the MPC, and no obvious optimum exists. Therefore, the design can either focus on fast cooling after door openings with high energy consumption or vice versa.
7. Discussion
The thermal storage capacity of the storage loop must be chosen appropriately, as both an over- and undersized secondary loop can hinder the overall performance of the cooling system. Employing an MPC, an undersized storage loop increases energy consumption and reduces the temperature control capability due to the limitations of the minimum glycol temperature, which only allows an insufficient amount of thermal energy to be stored in the secondary loop for the given cooling chamber. With an oversized secondary loop, the pass-through of the control input to the cooling chamber temperature gets too slow. Due to the finite horizon of the MPC, caused by limited computational resources, the time before door openings is too short to cool down the storage loop sufficiently to achieve fast cooling after door openings. On the other hand, the rule-based PI controller is unsuitable for utilizing the thermal storage of the secondary loop and thus has significantly higher energy consumption and takes longer to cool down the cooling chamber after door openings, which is in line with [
22,
29]. From the results, it is also evident that a smaller storage capacity leads to greater agility of the system since the pass-through of the controlled variables on the temperature inside the cooling chamber is faster, as also stated in [
23]. This allows the controller to respond quickly to incorrectly predicted door openings, model errors, or other disturbances.
Overall, this work shows the importance of considering the control strategy besides the vehicle operating conditions, such as ambient temperature, door openings, or other exogenous factors, when sizing the secondary loop storage capacity for the implementation in real-world refrigerated vehicles. Furthermore, this work recommends using predictive controllers for temperature control of secondary loop cooling systems, as only those are suitable to utilize the thermal energy storage of the secondary loop by predicting the future required cooling capacity.
However, the proposed sizing concept for the secondary loop storage capacity comes at the expense of extensive closed-loop simulations. Especially with the MPC, the computational time is substantial due to the complexity of the optimization task, even though, for example, limited horizon length and move blocking were applied. Although it is not a problem for the offline simulations performed, both controllers are real-time capable and thus suitable for application in actual refrigerated vehicles, which was also experimentally shown for similar control schemes in [
29,
44].
Since the vehicle under consideration is intended for use in Central Europe, simulations with an ambient temperature of 22 °C, characteristic for this region of the world, are considered in this work. Nevertheless, the proposed methodology offers the possibility to consider different ambient temperatures to adapt the sizing of the secondary loop thermal capacity to the specific operating conditions under which the vehicle shall be operated. However, for ambient temperatures below the freezing point, it is recommended to extend the model of the refrigerated vehicle to account for ice formation and mass infiltration during door openings as discussed in [
28].
Furthermore, when designing a refrigeration system, one should consider economic factors, as was done in [
11,
12,
45,
46]. The costs of the refrigeration system over its entire lifetime are mainly driven by its energy consumption when considering a fixed system architecture (secondary loop cooling system with a sensible heat medium). A change in the thermal storage capacity of the secondary loop and the associated change in the amount of glycol only marginally affects the costs. Hence, energy consumption was the only economic factor considered when optimizing the secondary loop storage capacity.
By expanding the model with a submodel for the cargo as described by [
47,
48], the thermal capacity of the cargo could be considered in the optimization of the secondary loop sizing. However, as is common for small-scale refrigerated vehicles, the frequent change of different types of cargo, whose parameters are often unknown, makes this approach impractical. Therefore, an empty cooling chamber was assumed in this work, which is the most conservative approach, as there is no additional thermal storage within the cooling chamber.
The proposed sizing concept can be utilized to design secondary loop refrigeration systems, increasing their efficiency and cooling performance. In addition, optimizing the sizing of the secondary loop can allow for a smaller cooling unit design, which further entails increased efficiency and lifetime, and reduces the acquisition costs of the cooling system. In general, combining a secondary loop refrigeration system with a suitable control scheme can be a simple alternative to more complex and costly cooling systems with PCM-based TES [
14,
22].
Future work includes the consideration of the vehicle’s powertrain when optimizing the size of the refrigeration system’s thermal energy storage. Since in most refrigerated vehicles, the energy for the cooling system is provided directly by the vehicle’s powertrain, optimizing the thermal energy storage of the refrigeration system can enable the use of load-shifting strategies for the powertrain to increase the vehicle’s energy efficiency further. Another task of interest is optimizing the cooling power of the primary loop in parallel to the size of the thermal energy storage.