1. Introduction
The need to reduce the environmental impact of industries and prepare industries for the intermittent behavior of the increasingly renewable future-energy-system necessitates disruptive approaches for transitioning the industry to become increasingly flexible [
1]. The typical energy-flexibility approaches considered in industries include the introduction of carbon pricing and time-of-use-based electricity pricing [
2]. The industrial facilities can reduce the overall operation costs while aiding in electricity-system balancing through the flexible operation of the industrial units in combination with time-of-use-based electricity pricing. The facility flexibility can be derived from the production process units within the facility, providing an instantaneous response. The flexibility can also be derived through flexible control of the production environment, which can be curtailed based on the thermal inertia of the system.
Due to the electricity consumption and thermal inertia, refrigeration and cooling processes and ventilation units have been identified as end-use applications with significant flexibility potential [
3]. One of the industrial sectors with a substantial use of refrigeration and cooling processes is the meat industry [
4]. An example of a meat production facility is the production of canned meat that can be used for long-term storage. However, there is still hesitance among industrial consumers to provide energy flexibility, due to uncertainties surrounding energy potential, production implications, and economic influence. For instance, in Denmark, the meat industry constitutes 2% of the annual electricity consumption, and can hence provide substantial flexibility. Within the Danish meat industry, specific steps of the meat-processing stage are subject to legislation that hinders flexible consumption by ensuring food safety and spoilage. Therefore, the primary refrigeration and cooling for flexible operation are found in the meat-production facilities that produce consumer products from the meat cuttings.
However, process cooling for food production is an energy-intensive industry with complex interactions and restrictions that complicates the ability to utilize energy flexibility, due to unforeseen consequences to production. Meanwhile, introducing energy flexibility to industrial facilities has been coupled with several barriers that hinder participation in energy-flexibility programs, e.g., a lack of understanding of the consequences flexible operation may induce on production [
5,
6]. These barriers within the industrial-process operations are complicated, due to complex system-interactions [
7]. Furthermore, research has mainly focused on single-process evaluation of flexibility potentials that do not adequately consider up- and downstream propagation of production decisions [
8]. Moreover, little literature has focused on the energy-flexibility potential in the meat-canning process. Methods for assessing the potential flexibility in individual facilities which can enable the active participation of process-cooling facilities in the electricity system are essential, but not yet well discussed in the literature.
Therefore, this paper introduces a simulation-based assessment method for investigating energy flexibility in process cooling with a case study of a Danish process-cooling facility for canned-meat food production. Multi-method simulation is used in this paper: multi-agent-based simulation to investigate individual entities within the process-cooling system and the system’s behavior; discrete-event simulation to explore the entire process-cooling flow; and system dynamics to capture thermophysical properties of the refrigeration unit and states of the refrigerated environment. Furthermore, a multi-method simulation library is developed in this paper with the OptQuest optimization engine for accessing energy-flexibility potential in the case study.
The rest of this paper is structured as follows: firstly, a literature review is presented, examining the current literature on energy flexibility in process-cooling facilities. Secondly, the applied methods, the agent-based modeling and multi-objective optimization approach, are introduced. The simulation-development section outlines the simulation architecture and agent communication. Subsequently, the individua- agent logics are outlined, covering the production environment, meat product, process-flow agents, and the refrigeration unit. Together with the refrigeration unit, the schedule-optimization explanation is provided. The scenario-design section introduces the scenarios conducted to verify and validate the designed agents. Based on the developed simulation-solution, three scenarios are presented; the scenarios are presented to verify and test the developed solution. The case-study section outlines the case study used for testing the scenarios with the specific agent-population sizes. Subsequently, in the results section, the resulting operation based on the case study is presented, emphasizing the refrigeration-cycle performance for three scenarios. Lastly, the results and limitations are discussed, before the conclusion section.
2. Literature Review
Industries have high demand-response potential due to high baseloads, less coordination effort, and energy-management systems already in place [
9]. For instance, in Denmark industries consume 26.3% of the Danish end-user electricity consumption [
10], and the amount of flexible electricity consumption within Danish manufacturing companies is equivalent to 380 MW [
11].
Implicit demand-response entails the industrial company reacting to market-price signals, e.g., from Nord Pool Spot Market, and thereby adapting the electricity-consumption behavior to lower operating-costs. In practical terms, this incites the industrial companies to allocate their bulk production of goods to hours where the day-ahead electricity price is low, to benefit from arbitrage. In addition to implicit demand-response, it is also possible for the industrial company to participate in explicit demand-response. Explicit demand-response is dispatchable and tradable flexibility [
12]. A third-party, e.g., an aggregator, often facilitates explicit demand-response. Explicit demand-response is often divided into programs based on the capacity and response times of the actors.
Refrigeration and cooling represent most of the electricity consumption in the food-processing industry where, on average, approximately 50% of the electricity consumption is used for cooling and refrigeration [
4]. The food processing can be subdivided into processing and production. The processing stage uses the most electricity. Some food, e.g., meat, requires fast freezing-speed to avoid spoilage [
3]. Therefore, the electricity used for freezing tunnels is considered inflexible. The electricity usage related to the storage cooling that keeps the food under acceptable conditions between process steps can be considered for energy-load deviations and flexibility.
There is large energy-flexibility potential for ventilation, refrigeration, and cooling. For instance, according to [
3], cooling processes could be used for load shifting for multiple hours and for around 50%, equivalent to 1513 TJ, of electricity consumption. However, process-production systems often involve complex systems and processes that are challenging to model and describe, and often operate based on extensive tacit knowledge [
7]. For instance, ref. [
13] examines the meat-production process, and presents the initial processing of different meat types (cattle, pork, poultry) and the energy usage associated with production across European countries. Similarly, in [
14], the saving potential in the food industry for Latvia and Kazakhstan are compared, comparing the potentials for shifting for cooling and ventilation.
Process cooling has previously been examined for its ability to provide flexibility in the electricity market. Ref. [
15] examines the potential for demand response through fan-speed variations in the drying chambers and cooling interruption in the drying chamber, and shows a 6.65% reduction in total cost of electricity. Ref. [
16] evaluates the effects of demand response in a refrigerated room, and shows that utilizing demand response once per day yielded a core-temperature increase of the stored product of 1.1 °C after three days. Ref. [
17] reviews the economic and environmental impact of the meat industry’s participation in energy flexibility, and shows that the meat industry can save between 3 and 5% on annual CO
2 emissions by participating in demand response and 5–6% on the overall electricity cost. Ref. [
18] examines a frozen-food manufacturing plant for load balancing in the processes of cooking, deep freezing, and cold warehouse, and finds that the facility can provide balancing for the local power system. Ref. [
19] examines the potential for HVAC systems to provide demand response in industrial facilities with high thermal-inertia, and finds a potential saving of 15.23% to 17.33%.
Furthermore, ref. [
20] proposes a multi-criteria evaluation framework for mapping various industrial processes to the feasibility of electricity-market participation. The study considers a cooling process for food storage, which has fewer barriers in capacity-based reserve markets. However, the study also remarks on the need to consider the facility’s flexibility and the interactions of the industrial processes with other processes along the production chain. Ref. [
21] finds that an industrial-scale bakery with a glycol coolant-system (that cools the bread-dough mixers and the fermentation room) can respond to electricity system changes.
Process cooling has also been considered extensively in combination with various types of thermal-energy-storage technologies. In [
22], slurry-ice thermal storage is proposed for the process cooling of cheese production, to aid in reducing the cost of energy costs. In [
23] a site-specific feasibility study is performed to investigate the potential for aquifer thermal energy storage. In [
24], an HVAC chiller is examined in conjunction with an ice-storage system in a supermarket, and is investigated for the potential to provide ancillary services.
However, industries’ active participation in implicit demand-response schemes is associated with various barriers. In the process-cooling facilities, the primary barriers and concerns regarding implicit demand-response participation are linked to the quality constraints of the product. As described in [
8], there is a need for to quantify the implications for a production process in terms of up- and downstream propagation that may be incurred due to a flexible operation. In [
5], barriers to demand-response programs are examined in which the most frequent barriers are related to technical uncertainty in terms of damage to the product; however, economic considerations such as idle labor and unpredictability are mentioned. Similar findings are presented in [
6], with the addition of research that can identify portfolios of demand-response measures that can be executed with minimal impact on the production.
Furthermore, the demand-response potential for a process-cooling facility involving food will often be subject to food-temperature restrictions. Ref. [
25] introduces a supermarket-food-temperature estimation approach for a portion of ground beef. The estimation approach is a time-constant approach for determining the temperature increase over time, and the result shows that the food can reach its maximum allowable temperature within 150 min. The shift in refrigeration load under consideration of the food temperature has also been discussed in [
26,
27].
Moreover, there are often limitations from a legal standpoint regarding how much temperatures can vary in the cooling process. The cooling processes are often subject to some leeway regarding temperature setpoints. Ref. [
28] finds that varying temperatures do not affect the food products significantly, enabling the potential for flexible operation.
Several studies have used optimization to effectively utilize energy flexibility. For instance, ref. [
29] applies a linear program for optimizing a food-dehydration process using a drum dryer. Ref. [
30] applies a mixed-integer linear program across seven production lines to optimize the food-production scheduling with the consideration of multiple factors for each production line, e.g., energy consumption, productivity, and labor costs. Ref. [
31] uses a mixed-integer nonlinear program for optimal chiller-loading to improve the operation of a large-scale process-cooling facility, and an 8.57% energy saving is achieved. The study also considers variable electricity price and shows a potential economic saving of approximately 42.2%. Using an optimization-based framework examines the potential for chiller plants to provide demand response through load curtailment [
32]. Ref. [
33] investigates the potential for a large constant-speed centrifugal chiller to provide grid-frequency regulation, and finds 5–7.5% regulation-capacity potential by using a model-predictive-control strategy.
However, as the optimization methods do not consider the internal-production-process flow and optimization, start conditions may change before an optimal solution has been found in a dynamic production system. Simulation methods can solve the above issues and can capture the inherent uncertainty in production processes [
8,
34]. Some previous research has applied agent-based simulation to the food industry. For instance, [
35] uses agent-based simulation to examine the beef-cattle supply chain. Ref. [
36] introduces agent-based modeling in the meat industry as a control layer in production. Ref. [
37] utilizes agent-based modeling for brewery-fermentation optimization. Furthermore, to efficiently examine the potential for implicit demand-response in a canned-meat production facility, the applied methods should be able to address all aspects of the facility sufficiently.
Based on the reviewed literature, it was found that limited literature has discussed the energy-flexibility potential in the meat-canning process. Furthermore the literature has emphasized single-process evaluation of flexibility potentials that do not adequately consider up- and downstream propagation of production decisions. Therefore, the novelty value of this paper is to introduce a simulation-based assessment method for investigating energy flexibility in process cooling, using a case study of a Danish process-cooling facility for canned-meat food production. The importance of the findings proposed in this paper is their contribution to the enablement of increasingly flexible behavior in industries using process cooling.
4. Simulation Development
The developed simulation represents the chosen canned-meat production, which can be used for examining the potential for energy flexibility. Canned-meat production consists of several facility-specific process steps. While in production, the meat is contained in various cold-storage areas collectively controlled by one or multiple refrigeration units. Facility resources can be required for operating the processes, e.g., personnel for operating machinery or forklifts for transporting goods.
An overview of the simulation architecture is shown in
Figure 1. The top-level process-cooling agent represents the given facility. Each facility may have multiple rooms or storage areas with various temperature-requirements and profiles. It is important to note that the primary production-flow is modeled and matched with the given cooling-environment.
The process cooling is a placeholder for all the environments in the different rooms across the entirety of the facility. Several processes and conveyor blocks are associated with each environment that receives its temperature data. The user controls the associated processes and conveyor blocks. The environment is not required to contain any blocks. All the environments are connected to the centrally controlled refrigeration unit. The indoor temperature is determined based on the refrigeration unit’s electricity consumption and current outdoor temperature. As seen from
Figure 1, there is a continuous flow of products between processes and conveyors, and the processes and conveyors are subject to the environment in which they reside. Thereby, the environment will continuously notify its observers of any state changes. Using an observer pattern enables a dynamic number of observers within the environment subjects. Similarly, the environments observe the refrigeration unit and the internal-state changes in response to the refrigeration-unit input are determined by the individual environments. The communication of the developed simulation follows the sequence diagram, which is shown in
Figure 2, and the refrigeration unit determines the state of the refrigerated rooms. The state of the refrigerated room subsequently propagates into the agents contained in the environment.
In
Figure 2, the adjustment of the refrigeration unit is started from the top-level industrial facility. The refrigeration adjustment can happen at specified intervals. Once the signal for refrigeration-unit control is sent, the refrigeration unit will collect the temperatures of the rooms and the current outdoor temperature, which will determine the refrigeration power-consumption for the forthcoming hour. The refrigeration power is published to the cooling environments, which adjust their temperature based on the supplied cooling power and the current outdoor temperature. Once the temperature adjustment has been completed, the alarm setpoint is checked, and if the temperature has been above the alarm setpoint for more than the specified alarm-delay-time, the specific cooling-environment will trigger an alarm at the top-level facility agent.
Subsequently, the state of the cooling environment is published to all the environment listeners. The environment listeners consist mainly of process-type agents that will pass the state information to agents currently contained in the process, i.e., the meat-product containers and meat products. The meat agent will adjust its temperature based on the environment state received, which in turn updates the thermophysical properties of the meat for future state-changes.
As shown in
Figure 2, the agent-based system follows an observer-based design pattern. The use of an observer design-pattern enables system flexibility and reusability. An event-driven system furthermore complements the domain-specific use of observer responses to a subject, e.g., cooling-environment state changes are reflected in all agents currently residing in the specific environment that implements the required listener-interface.
An overview of the individual agents developed for the process-cooling facility and their functionality can be seen in
Table 1:
The primary parameters that are modified during simulation are those of the operating schedule of the refrigeration unit, which propagate to the cooling environments. The state of the cooling environments will determine the temperature of the meat. The temperature estimation of the meat is, in this case, dependent on the composition. There are four primary agents in the simulation: production environment, meat product, process agent, and refrigeration unit.
4.1. Production Environment
The production-environment agent was created as a placeholder for the facility’s environment in a specific room. A central parameter to be estimated for the operation of the process cooling is the room temperature. The temperature is considered a main decision parameter, as the correct temperature ensures the quality of the products. The literature’s general approach to cooling-load estimation is centered around Newton’s cooling law [
51]. Cengel states that the temperatures evolve according to an exponential curve [
52]
where
T(
t) is the temperature at the next time step,
Ti is the internal temperature at the current time step, and
T∞ is the surrounding temperature. The temperature is assumed to be uniform, and the reciprocal time constant is defined as:
where
is the overall heat-transfer-coefficient,
is the surface area,
is the density,
the volume, and
is the specific heat under constant pressure [
52]. The reciprocal-time-constant parameter,
, can be estimated based on the temperature segments found in collected data, similar to the approach presented by Herten et al. [
51]. The reciprocal time constant is used to ensure that the temperature development does not exceed the observed maximum hourly change.
The production-environment agent adheres to the company’s temperature setpoints and restrictions. The state-chart representation of alarm behavior is shown in
Figure 3.
The production-environment agent will initially start with the alarm off. It will start the alarm timer if the temperature is above the setpoint. If the temperature is not reduced before the alarm timer runs out, the alarm is triggered, and a warning message will be issued. Alternatively, the temperature can become too low, damaging the products and directly triggering an alarm. In terms of optimization, an alarm trigger corresponds to an infeasible run. The temperature within the compartment is updated continuously. The compartment-environment temperature is determined based on Equations (1) and (2).
4.2. Meat Product
Limited information was obtained about the meat within the production flow. However, based on an example provided by [
3] and the ASHRAE chapter on the food-temperature thermal properties, modeling of the individual meat products could be performed [
53,
54]. The temperature modeling of the meat was carried out using Newton’s law of cooling, previously presented in (1). As surface temperature is used for the rules on the hygiene of food of animal origin, a uniform-temperature estimation is assumed to be sufficient [
55]. The thermophysical properties of meat are highly dependent on the product’s composition. It is assumed that the pork meat adheres to the composition found in the shoulder of the pig. The composition is outlined in
Table 2.
The meat’s thermal properties were regularly evaluated based on the above composition and the thermal-property model for food components presented in [
5]. The density of the meat is calculated according to (3):
The porosity, ε, of the meat is zero, as no granularity is present, and
and
are the mass fractions and density of the meat constituent, respectively. In a similar way to the density, the specific heat of the meat is modeled according to its constituents, as seen in (4).
The heat transfer coefficient of the meat was modeled based on the experimental findings presented by [
56], which were extended to be described using the regression model seen in (5):
where
describes the air velocity. It is assumed that the meat was always kept above its initial freezing-point, for modeling the thermal properties. Falling below the initial freezing-point will lead to significant changes in the thermal properties and potential damage to the product. The initial freezing-point for a pork shoulder is −2.2 °C, significantly below the compartments’ storage temperature. Within the production-flow simulation, the meat-product agent adheres to the production environment associated with the process in which it currently resides.
4.3. Process-Flow Agents
The process agent refers to the fundamental process logic inherited by any specific process agent used for the simulation. The process agents represent the subprocesses shown in
Figure 4, which the product agent will pass through during its lifetime in production. The flowcharts used for specific agents within the simulation can be seen in
Figure 4.
The meat-product agent will be inserted into the process block through the enter block. After entering the block, the product will go through the process steps in the specific agent before exiting the next step in production or finishing the production. All the process blocks were modeled based on the above concept. The cold storage had to be slightly altered to include a rack system for storage. For raw-material storage, the products are transferred to a storage rack, where they will stay for a specified time before being picked automatically and continuing production. The rack-store delay may be case-specific, as the meat is retrieved based on its usage in a required recipe. The current implementation delays all the products for the same duration of time. The meat grinder will wait for a specified amount of meat agent to enter the process; once the amount is satisfied, the products are batched into one new meat agent which is delayed, based on the grinder setting. The filling machine will take two inputs of a container, i.e., a can and a containable entity, i.e., meat, and will add the meat to the can, which can be delayed for a set period. The hydrostatic cooker will take the filled-can agents and delay them for a period corresponding to the cooking time.
4.4. Refrigeration Unit
The refrigeration-unit agent represents the facility’s control and operation of the refrigeration unit. The state-chart logic used for the agent can be seen in
Figure 5.
The refrigeration unit responds to an inputted operation schedule. The initial state of the refrigeration unit is chosen based on the initial schedule, through the branch shown in
Figure 4. The ON or OFF states do not represent the refrigeration unit being entirely on or off. The states correspond to an increase or decrease in power consumption, as it is currently enforced that the refrigeration unit should run all the time, and the capacity is adjusted. The refrigeration unit’s primary power-consumption is attributed to the compressors that adjust the inputted work, based on the refrigeration load. As modeling the correct behavior of the refrigeration unit poses a significant challenge, with underlying assumptions that can be challenging to address, a data-driven modeling approach is applied.
Symbolic regression has been identified as a promising solution for obtaining explanatory models without imposing any a priori assumptions [
57]. Symbolic regression has seen development using both genetic programming and simulated annealing; however, as shown in [
58], symbolic regression using simulated annealing generally outperforms symbolic regression based on genetic programming. The performance was further underlined in another study comparing state-of-the-art symbolic-regression tools; the study showed that the software TuringBot, based on simulated annealing, outperformed its counterparts, Eureqa and AI Feynman, discovering target equations faster and with fewer data points [
59]. Symbolic regression was enabled through TuringBot to establish the refrigeration-unit electricity consumption under various conditions. Therefore, TuringBot was used to develop a refrigeration-unit power function, relating the individual room-cooling-load to the refrigeration-unit electricity consumption.
Multi-Objective Optimization of Refrigeration-Operation Schedule
Enabled by the OptQuest optimization engine, multi-objective optimization was used to select the refrigeration-unit-operation schedule. The objectives used for the multi-objective optimization are stated below (6):
The minimization is subject to operational constraints of the refrigeration unit. Furthermore, the problem is constrained by the number of cooling-system alarms being equal to zero. The optimization was performed in incremental steps for the simulation. Each simulation step corresponded to the optimization horizon. A sequence diagram providing an overview of the multi-objective optimization utilization can be seen in
Figure 6.
As seen from the sequence diagram presented in
Figure 6, the optimization is started by generating all possible operation schedules. Using the developed simulation, the operational schedule performance is tested under the outlined constraints. Once an optimized schedule has been identified, it is used to for simulating the forthcoming 24 h in the simulation before the optimization procedure is repeated. As the simulation-optimization engine does not support schedule optimization, an alternative approach is used for selecting the appropriate schedule. A discrete schedule contains a finite number of possible combinations and a finite number of possible states. Therefore, the possible schedules can be constructed as arrays of states from the current timestamp,
n, until the desired planning horizon,
m, i.e., the number of future timesteps considered.
The possible schedules can be collected in an array with a length corresponding to the number of possible permutations. The possible permutations are derived from the number of possible states, s.
Each possible schedule is linked to the resulting operation’s associated electricity price, carbon dioxide emission, and facility operation. The simulation-optimization engine can subsequently use an index list referring to all possible permutations, where each index holds a unique combination of states. Based on the index list, the simulation environment can be used as a testbed for examining the impact of a specific operation schedule.
Assuming that the scheduling is carried out on an hourly basis and that the state of the refrigeration cycle is operated in an ON/OFF manner, inferring two states, the total number of operational possibilities can be expressed as (8) [
60]
where
m represents the scheduling horizon and (8) represents the total number of permutations with repetition allowed. It is essential to include repetitions, as the system’s current state will influence its future states, due to changing temperatures. Each operational possibility can be expressed as a binary number with a number of bits corresponding to the scheduling horizon. In this case, zeros correspond to OFF, and a one in the system is ON. Generating all possible schedules hence presents a matrix of size
, in which each row represents a possible operational schedule. Based on the possible operational schedules, the objective functions are optimized under the previously outlined constraints. The scheduling is carried out on an hourly basis, in order to coincide with the electricity spot-prices which are supplied on an hourly basis. As the scheduling is carried out hourly, it is hence assumed that the refrigeration operation is varied on an hourly basis. The room-temperature development could be simulated in response to each binary string’s composition. As the optimization is carried out room-wise, each room’s preferred schedule may vary, based on the parameters and constraints; i.e., one room may increase in temperature more rapidly than another. Therefore, the optimization was carried out using the temperature development for all rooms, and only solutions that adhere to all the individual room constraints were considered for use. The approach was deemed necessary because the refrigeration unit is centrally controlled, and any operational change will affect all connected rooms.
The schedule was optimized using the CO
2 emissions and electricity spot-price. The Danish transmission system operator, Energinet, supplies a forecast of the electricity spot-price for the next 24 h and for the forthcoming 9 h for CO
2 emissions [
61]. Note that for the forecasted CO
2 emission, Energinet ensures a minimum of 9 h of prognosis, but can include more hours. The energy-system forecasts were combined with a weather forecast to correspond to the scheduling horizon. The electricity price and CO
2 emissions associated with the operation were calculated as shown in (9) and (10).
where
, is the electricity consumption in the hour
i, and
, is the cost of electricity in the hour
i. Similarly, the CO
2 emissions are calculated as:
where
, refers to the amount of CO
2 emission per consumed electricity-unit.
6. Case Study
The case study process-cooling is located in Denmark, primarily produces canned meat, and has a small production line for sausages. In this case study, the primary production line of canned meat is examined, and some aspects of the temperature development for other facility parts are also included. The production facility is divided into several rooms that represent a specified space. The specific rooms may be physically interconnected or separate. An overview of the facility layout is shown in
Figure 7.
As shown in
Figure 7, meat cuttings are first received in the raw-materials section, room 1071. Subsequently, the meat is transported to the curing station; the meat is placed in containers, each of approximately 30 kg, with a distribution of cuttings matching the specific recipe. After being placed into containers, the containers are transported by conveyor to the storage room, and stay there for a prolonged time, in room 1060. Once ready, the containers adhering to the same recipe are transported by conveyor to a selected mincing-and-mixing machine. Each mincing-and-mixing machine marks the beginning of a factory line. The facility has a total of 7 lines for producing canned meat. After mincing and mixing, the minced meat is filled into cans and sealed. Once the cans are sealed, they are transported to the cooking station, where the meat is cooked within the can. The cooking station is the final step before the cans are ready to be packed and shipped to customers. The primary canned-meat production process can be represented by the flowchart seen in
Figure 8.
Consumer requirements, food-safety law, and regulations partially specify the room-temperature requirements. The temperature of each room within the production flow is monitored. The restrictions and setpoints of the different rooms are shown in
Table 3. The temperature setpoint is the temperature at which the room should be maintained. The alarm setpoint is the maximum temperature allowed to trigger the alarm delay. If the temperature stays above the alarm setpoint for the alarm-delay duration, the alarm will trigger, leading to food spoilage and economic loss for the company.
A central refrigeration-unit refrigerates the facility to maintain the required temperature within the individual rooms. The refrigeration unit comprises several components, and it is operated using five compressors and three condensers with individual fans and pumps. Each condenser operates with two fans and a pump that sprays water for increasing heat-transfer. The refrigeration cycle uses ammonia and has three tanks to store cooled ammonia temporarily. The cooled ammonia is utilized through heat exchangers; the heat is exchanged into glycol, used as the medium for refrigerating the facility. The refrigeration unit’s power consumption can be divided based on the components shown in
Table 4.
The refrigeration unit is continuously operating, and the electricity consumption is regulated based on the cooling demand. The compressors are turned on hierarchically once the current compressor(s) is insufficient to meet the cooling demand.
The examined facility provided several datasets for production and refrigeration; an overview of the available data can be seen in
Table 5. Besides the facility-specific data obtained from the company, several external-data sources were integrated into the simulation. The data were retrieved using the associated APIs. All the data were extracted based on the Unix timestamp of the simulation or the correlated UTC timestamp, to ensure the correct model time and data match. Currently, the Danish Meteorological Institute API is only used to retrieve the temperature; however, it is also possible to extract other parameters, e.g., wind speed, wind direction, and humidity, which can also be influential. All the Danish transmission system APIs were sorted based on data from DK1. If the simulation model is used for future projects in other physical locations than DK1, this should be adjusted. Note that the CO
2 emissions in DK1 and DK2 are not differentiated.
8. Discussion
Three scenarios cover the baseline operation and flexible operation based on electricity price and carbon dioxide emissions, respectively. The baseline scenario verifies the system’s ability to reflect the case study. Simulating during 5–12 October 2020 shows consistency between the historical and simulated electricity-consumption, with an absolute error of 0.04 kWh. The flexible-operation scenarios show that the refrigeration load can be shifted without violating room-temperature constraints, and could save the company approximately 32% on electricity costs and 30% on associated carbon dioxide emissions. Similarly, previous research examining energy-cost saving has reported potential from 5% to upwards of 70% [
15,
37,
62,
63].
As in results shown in
Table 8,
Table 9 and
Table 10, the facility can adjust its load to take advantage of the electricity-market conditions and forecasts. It is possible to reduce the facility’s electricity costs and CO
2 emissions. The current price optimization is the implicit demand-response that reacts to the electricity-market-price signals. However, electricity-market participation through explicit demand-response is expected to provide a more significant monetary benefit. Currently, the energy flexibility is fully based on the refrigeration unit within the facility. However, some flexibility may also be available in the operations through adequate planning. A comparison of the three scenarios and the actual consumption is shown in
Table 11. The simulated and actual consumption are comparable, indicating that the developed simulation-model can capture the underlying system. The savings can be realized effectively, and there are differences depending on the emphasis, i.e., wanting to save more on the electricity price or CO
2, even though both operations induce savings.
It should be noted that the Danish government has agreed to introduce carbon taxation from 2030 on 750 DKK/tonne of emitted CO
2 [
64]. As shown in
Table 11, accounting for the taxation in the overall evaluation of the results would effectively make the operating cost of the CO
2-based operation marginally cheaper, indicating that industrial facilities should consider this when considering a flexible operation.
The relationship between electricity price and carbon intensity has been shown to follow a positive relationship [
65,
66]. Therefore, optimizing toward either electricity cost or carbon intensity is likely to entail cost savings, which can also be observed from the results in
Table 11.
The observed savings resulting from the multi-objective optimization should be considered in response to the facility not utilizing their underlying leeway allowing temperature intervals. The difference in energy consumption results from operating the refrigeration unit closest to the temperature limits. The difference in CO
2 emissions is a result of the variability of the emissions associated with electricity consumption at various hours. Therefore, savings can be observed as a combination of the refrigeration unit operating less and becoming flexible. The saving should also be considered in response to uncertain room-temperature behavior. Generally, including sensor data in the production-flow model could increase the resolution of the simulation. Integrating IoT devices within the production-flow simulation allows new opportunities to be tested in a risk-free environment. Furthermore, this would allow AI-driven solutions to further the operation and achieve optimal refrigeration-load scheduling. Currently, the temperature estimation of the compartments is carried out based on the data analysis. As [
51] noted, adjacent cold stores influence each other in terms of refrigeration load, and this could hence be included in further investigation. In the process-cooling field, the application of reciprocal time constants, e.g., [
25,
51], are used to investigate the cooling-room-response over time. However, approaches using deep-learning artificial-neural-network models for predicting the room temperature may provide improved approximations over time [
67]. Other rooms and food-temperature-estimation approaches include model-based estimations using Kalman Filtering, as presented in [
68]. In this paper, reciprocal time constants were used for estimating the temperature response of the individual rooms, based on data provided for the case study.
Previous research has shown the ability to use various machine-learning approaches to examine the curtailment period for energy flexibility. Hoang et al. [
67] presented LSTM RNN models for predicting power and temperature variations. Similarly, Herten et al. [
51] introduced least-squares support-vector regression and Gaussian processes for accurately predicting the time to reach a specified boundary-condition. This paper introduces the use of a symbolic-regression method, which enables the possibility of visual inspection of the provided function and ease of integrating the function into the simulation environment. In addition, this paper considers the underlying production-process flow in combination with load curtailment for implicit demand-response.
This paper utilizes agent-based modeling to promote the reusability of the agents across the domain, similar to [
42]. As shown in the sequence diagram presented in
Figure 2, the agent communication follows an event-driven observer pattern. Through the generic-agent interfaces, the agent behavior in response to, e.g., the refrigeration cooling-load, can be altered to utilize other methods or customized to the specific case study.
Most previous research applies optimization approaches to address the production process and energy cost [
29,
30,
69]. However, using optimization alone is challenging for capturing underlying uncertainties in production processes, verification issues, and long runtimes [
8,
34]. The multi-method-simulation and optimization approach in this paper can capture underlying uncertainties in production processes and improve the performance of the energy-flexibility-potential assessment.
9. Conclusions
Process cooling for food production is an energy-intensive industry with complex interactions and restrictions that complicates the ability to utilize energy flexibility, due to unforeseen consequences in production. Therefore, methods for assessing the potential flexibility in individual facilities to enable the active participation of process-cooling facilities in the electricity system are essential but not yet well discussed in the literature. Therefore, this paper applies multi-method simulation and multi-objective optimization for investigating energy flexibility in process cooling with a case study of a Danish process-cooling facility for canned-meat food production.
The multi-method simulation with multi-agent based discrete-event and system dynamic-simulations have been developed for the canned meat production-flow. The developed simulation-model is centered around a data-driven approach with the integration of external-data sources. Furthermore, a simulation architecture is proposed, and the architecture provides a foundation for simulating production-flows in varying process-cooling facilities. Several generic agents are designed to interact in the simulation environment. The core agents are a production-environment agent representing the physical environment within the different rooms; a meat-product agent that was processed throughout the facility; a process agent which is used to represent the different processes that the product goes through; and a refrigeration-unit agent to control the cooling load of the production environments, centrally. Using the developed simulation, multi-objective optimization was used to select a refrigeration-unit load schedule to minimize overall electricity cost and CO2 emissions.
Three scenarios are designed and tested in the paper. The baseline scenario shows that the current operation could be captured using the developed agents. Building on the verified simulation-model, the decoupled operation could be examined for flexibility potentials using multi-objective optimization. The results show potential for the examined process-cooling facility to adjust its load based on electricity-market signals, to reduce the overall electricity cost and CO2 emissions. For one week (in October 2020), 32% of electricity costs and 822 kg of CO2 can be saved.
Furthermore, this paper proposes a domain-specific multi-method simulation and optimization approach for accessing energy-flexibility potential in industrial processes. The application to the process-cooling industry shows that this approach can not only decode the uncertainties in production processes, but also provide a holistic perspective on energy-flexibility potential. A simulation library is developed, and is able to represent a generic production-flow of the canned-food process cooling. This library allows building case-based simulations with pre-built facility components. Furthermore, various external factors, e.g., changed weather conditions, external electricity-market signals, etc., are also included in the library.
The food-process-cooling industry can use the proposed method to examine the energy-flexibility potential in the production-flows, ensuring the energy flexibility does not compromise the quality or food-safety requirements. Furthermore, the multi-method simulation and optimization can be modified to represent multiple facilities for examining the flexibility potential across an entire electricity end-user segment. Moreover, the method can be used by the distribution- and transmission-system operators to identify flexible and available loads within their electricity system.
This paper only considers implicit demand-response as a means of energy flexibility, due to the increased resistance towards explicit demand-response in the industry, due to the potential loss of control. The economic potential in explicit demand-response may be higher and recommended for examination in future research. Furthermore, the seasonal variance might influence the maximum curtailment in the facility, but only one week’s data is obtained in this paper. Therefore, a longer period should be examined to obtain a holistic understanding of the energy-flexibility potential. Furthermore, it would be beneficial to implement a sensitivity analysis to examine the impact of different cooling systems. Moreover, the general agent-based production-process flow could be extended to other food industries, e.g., breweries, to examine if any flexibility potential could be available during the food-processing stages. Furthermore, machine-learning methods might be useful for improving agent behaviors; the temperature development especially within the rooms could be improved. The optimization could be improved by expanding the optimization criteria to consider health and quality indicators, as well as additional key performance-indicators in the food industry. Lastly, the analysis should be extended to consider the entire supply chain, considering all stages for the manufacturers and consumers.