Solar Thermal Technologies for Low-Carbon Industrial Processes: Dynamic Modelling, Techno-Economic Analysis and Environmental Impact

: Solar thermal technologies are already available on the market, and they are robust and relatively cheap. Unfortunately, solar heat is seldom used in the industrial processes, and the main obstacle of solar heat diffusion is often the lack of adequate predictive modelling of solar plant integration that identifies its energy potential, economic feasibility, and environmental benefits. In this paper, we aim to investigate and evaluate the possibility of supplying solar heat to the pasta-drying process located in the northeast of the Italian Alps (“Felicetti”). The methodology proposed is structured with the combination of several software, namely, PVGIS ® , Matlab ® , Dymola ® . The methodology developed is tested, considering solar thermal energy as the primary source, in different geographical contexts.


Introduction
Thirty billion tons of carbon dioxide are emitted yearly in the atmosphere of the Earth, and this is causing global warming and drastic climate changes. Average global temperatures have risen 1 °C above pre-industrial levels and, at current rates of warming, are projected to reach 1.5 °C within two decades [1]. Therefore, it is urgent to plan small-and large-scale decarbonization interventions. The use of alternative eco-sustainable and carbon-free technologies for the generation of electrical and thermal energy is now a compulsory and conscious choice [2]. Solar thermal technologies are consolidated at economically attractive prices. The total carbon-free solar thermal energy should be a cornerstone of the current energy transition. The industrial sector uses more energy than any other end-use sector [3], and in Europe, it covers the 22.8% of total energy demand [3], with 71% [4,5] accounting for heat demand. Solar heat has enormous potential for several industrial processes but is still mostly unused [6,7]. This potential depends firstly on the amount of solar radiation available. This is an unpredictable variable depending on several factors, such as the location of the place, cloudiness, air mass etc. Moreover, a careful analysis and industrial process characterization in terms of thermal energy demand and process temperature profile is crucial to also achieve adequate integration of solar process heat. Worldwide, 152 operating plants provide solar process heat for industrial purposes, with more than 140,600 m 2 of solar collectors and power of about 100 MWth [8]. Most of the industrial heat is required at low to medium temperatures, with temperature below 250 °C [9] in processes such as sterilization, drying, hydrolysis, distillation, washing, cleaning, evaporation and polymerization [10]. Several applications for the integration of solar heat in industrial processes have been developed: applications of a solar drying system for agro-food industries [11,12] such as pasta [13,14] and tomatoes [15]; for the manufacturing of clay bricks and tiles [16]; and for the meat [17] and dairy industries [18]. In the IEA SHC Task 49 [8], two types of integration are identified: supply level and process level [8]. The supply level integrations differ according to the heat transfer fluid (steam or liquid) used and according to the application, whereas the process level integrations differ according to the way of exchanging heat.
Unfortunately, solar energy is a discontinuous source of energy, and for this reason, it is not easy to match with industrial processes that require energy on demand. However, as reported by many authors [13,16,18], this problem can be overcome by using dynamic predictive models to study and analyze the potential of solar heat integration. Many research works are based on modelling of the solar plant using Transys ® as software, and solar plants located in highly irradiated places such as Sicily, Indonesia, and Morocco. In this work, the authors propose a methodology approach based on the synchronized use of different simulation tools, modelling and comparing the potential of the integration of solar heat in the new pasta factory headquarters in Molina (Trento Area), and a hypothetical one in Palermo. A 100% renewable solution was found considering the land use and performance of the three investigated solar collector technologies. Finally, the use of Dymola ® has made it possible to verify the flexibility and dynamic capacity of the solar field adaptation to intermittent variations of solar radiation. An adequate control system was, therefore, necessary to govern the proper operation of the system. The considered industrial process is a real case study of the pasta-drying process. In the first part of the manuscript, the utilized methodology scheme that shows the different steps and information used to build the model is reported. In the second part, there is a detailed description of each step.

Methodology
The industrial process considered is a real case study of a pasta drying process at the Felicetti Pasta Factory. The production site is historically located in Predazzo, while a new plant will be built in 2021 in the location of Molina di Fiemme. These two locations are situated in the northeast Italian Alps, and they are characterized by a complex topography of ground, which affects the amount of solar radiation available at ground level mainly in the winter months when the solar elevation is lower. In this study, we want to propose a new methodology whose scheme is shown in Figure 1, which compares the potential of solar integration and the relative land use in a low-irradiated and complex topography location such as Molina and in a high-irradiated and simple topography location such as Palermo in the south of Italy.
The two locations were studied separately, considering the different geographical coordinates. The time series of a typical meteorological year (TMY) from 2007 to 2016 by PVGIS [19] are used for solar radiation and air temperature in Molina and in Palermo. The data collections used are the hourly series of DNI (direct normal irradiance), GHI (global horizontal irradiance), and environmental air temperature. Subsequently, three solar collector technologies were considered using data related to the efficiency, size and incidence angle datasheets of three commercial companies, such as Soltigua [20] for parabolic through collectors (PTC), Industrial Solar [21] for linear Fresnel reflectors (LFR) and Kloben [22] for compound parabolic collectors (CPC). In the coefficients used for calculating efficiency are shown.  In Table 2, the three solar thermal collectors are presented in terms of sun-tracking technology and land use, which is defined as the ratio of total net area to the gross area used. For the PTC and LFR technology to work correctly, usually, 1-D sun tracking with north-south orientation is used, whereas the CPC is stationary technology that can be tilted in order to optimize the thermal energy production. Table 2. Technology and land use of the three solar thermal collectors used in examination.

Type of Collectors
Technology Land Use (%) PTC It should be noted that the parabolic tracking of PTC collectors requires more space between rows than what is necessary for LFR collectors. In the specific case of CPC, on the other hand, the proper distance was calculated in order to avoid self-shading phenomena considering an optimal tilt angle of 37 degrees. Moreover, the three technologies present different values of longitudinal (Kl) and transversal (Kt) incident angle modifiers. These are calculated as polynomial functions of longitudinal ( ) and transversal ( ) components of angle of incidence, interpolating the values present in the solar collector datasheet. The components of the angle of incidence are defined as follows [13]: θ = arctan sin γ tan θ ; θ = arctan cos γ tan θ cos θ In which is the azimuth angle and θ is the zenith angle, which is the complementary angle of the sun elevation angle (α ) and depends on a number of hourly parameters, such as the hourly angle h, the latitude and the declination angle , as shown in the following equations: α = arcsin sin ϕ δ + ϕ δ ℎ ; γ = arcsin cos δ cos ℎ cos α The industrial process characterization is based on the drying process of pasta that requires a large amount of thermal energy. Super-heated water at 135 °C is used to guarantee adequate temperature and humidity conditions in four different drying tunnels. Each tunnel is used for a specific shape of pasta. The thermal energy is necessary to progressively reduce the fraction of water contained in the pasta until it reaches a rate below 12.50% [23], which is imposed by law (DPR 9. 2.2001 No. 187). The production process is continuous, 24 h a day. The average net production of pasta is 2730 kg/h, leading to a total annual production of 19,500 tons. The total thermal demand of the company is 6.69 GWh/year, while the total electrical demand is 3.81 GWh/year. The thermal demand is satisfied by two gas boilers (each of 2000 kW, efficiency: 90%) and a gas cogenerator (450 kWth and 460 kWel). With 4.73 GWh/year, the two boilers cover about 70% of the company's heat requirements, as shown in Figure 2. In contrast, the remaining part is covered by a combined heat and power system (CHP), which operates continuously, generating simultaneously 1.96 GWhth/year of heat and 3. 12 GWhel/year of electricity, thus covering about 82 % of the electricity demand. With an hourly analysis performed in MATLAB, the thermal power was calculated as follows: in which is the total net surface of the different solar field; is the product of and ; is the environmental temperature of the respective location obtained from TMY via PVGIS; is the hourly solar radiation consisting of DNI for tracking 1D technologies and GHI for CPC; and , and are the coefficients used for the calculation of efficiency. In particular, the PTC technology is independent from , which assumes a unit value. Therefore, it performs better in terms of overall efficiency.
Finally, MATLAB ® was also used for the adequate design of the heat exchanger (HEX) and the solar field hydro-dynamic parameters that are the necessary inputs for the subsequent dynamic modeling in Dymola. The idea was to integrate solar heat into the pressurized water heating loop of the conventional thermal energy system. The industrial water loop used in this assessment consists of an operating temperature of 135 °C and a return temperature of 125 °C. The objective of the simulation was to see how this integration works, with the interface of the solar loop consisting of the use of a commercial Therminol SP [24] diathermic oil, and with the water loop positioned through a suitable tubular heat exchanger. The design of the HEX was based on the calculation of a series of geometrical and fluid thermo-dynamic parameters, as shown in Table 3, considering the turbulence conditions concerning the number of Re and Pr exposed by the Dittus-Boelter correlation [25] and considering a Therminol SP loop based on a temperature of 185-155 °C. The thermodynamic properties of the heat thermal fluid (HTF), such as specific heat or viscosity, were considered as polynomials in the function of the fluid temperature [25,26]. In addition, a ΔT of 30 °C (185-155 °C) allows the design flow rate to be calculated in a completely turbulent regime. Therefore, a design flow rate of approximately 5.16 kg/s was calculated. Subsequently, considering a thickness 71% 29% Gas boilers CHP appropriate to the commercial pipes under consideration, the number of tubes (nt) inside the tubular heat exchanger was calculated, as demonstrated in Figure 3. The internal convection coefficient ℎ (Therminol SP side) is calculated according to the internal diameter with equation ( = ), while the external convection coefficient, ℎ (water side), is calculated in the function of the hydraulic diameter. Therefore, considering a thermal conductivity of λ steel of 15 W/mK, the overall transmission coefficients, referred to as the internal (Ui) and external(Ue) surface as a function of the respective convection coefficients, are calculated. Finally, the geometrical parameters, such as the internal surface (Ai), external surface (Ae) and length (L) of the HEX, were found according to a design of thermal power, considering a 90% efficiency of heat exchange. The volumes of the exchange mediums, Therminol SP (Vsf) and water (Vwf), and the volume of the metal were also found to be dependent on the geometric parameters. The total metal mass of the HEX was calculated considering a constant density of steel equal to 7980 kg/m3 (ANSI B36). The heat exchanger geometry and dimensions are highly influenced by the design of the exchange temperature. It can be noted that a tubular heat exchanger of 8.56 m has a weight of 175 kg. In order to minimize the space occupied, these types of heat exchangers have the possibility to possess coil-type configurations. These parameters will be used as input in the dynamic simulation in Dymola ® , whose simulation scheme is presented in Figure 4.
In the case of the Felicetti pasta factory, solar heat is integrated in the water loop that supplies the gas boiler, which is characterized by a moss flow rate of 15 kg/s and an input and output temperature of 135 and 125 °C, respectively. With the counter flow tubular heat exchanger (HEX), it is possible to distinguish the solar loop in which the Therminol SP circulates (on the left) and the water loop (on the right) that feeds the industrial process. Localized pressure drops are considered in each component of the system as constant, while the value of distributed pressure drops inside the pipes is proportional to the mass flow rate. Considering the limited available space, a small volume of the tank proportional to the flow rate of about 8 m 3 was considered. A series of PIDs that regulate the pump and the closing-opening of the three valves located at the solar field is set to a temperature of 185 °C at the tank outlet.

Results
Matlab ® was used to process the hourly data of solar radiation and climatic temperature (estimating the heat production and the components of IAM for the three different solar collectors) to size the heat exchangers (geometry, heat transfer coefficients and metal mass) and to evaluate the mass flow rate working conditions. In the second step, Dymola ® was used to monitor temperature transients and simulate the activity of the solar field in a predictive mode. The available space of the pasta factory is of 2000 m 2 . The yearly thermal energy produced by the solar field with the three technology in Molina and in Predazzo and the space required to totally cover the heat demand of the pasta factory using only solar energy (100% RES) are shown in Table 4. To satisfy 100% of the thermal demand by solar heat, it is necessary to produce at least 6.69 GWh/year to be combined with an adequate thermal storage to satisfy the daily and seasonal variability. The study of the thermal storage solution is not an aim of this work.  Table 4 shows that in order to cover 100% of the thermal energy demand of industrial processes, 27,710 m 2 of available space is necessary in Molina and 13,310 m 2 in Palermo. It is useful to notice that the LFR collectors are the technology that produces more thermal energy in one year in both locations. In the available space of the Felicetti pasta factory (2000 m 2 ), the LFR collectors are able to provide 0.341 GWh/year of solar heat in Molina and 0.73 GWh/year in Palermo.
In general, there is an increase in production in Palermo when compared to Molina of about 53% for PTC, 52% for LFR and 41% for CPC technology. Moreover, considering the same location, the LFR can produce 22% more than PTC and 53% more than CPC.
In Figure 5, it is possible to observe the yearly thermal energy production trend for the three solar collector technologies in the function of the gross area used: 2000 m 2 , 13,000 m 2 , and 27,000 m 2 . Almost 52% more available gross area is required in Molina, and with the available space of 27,710 in Palermo, it is possible to produce thermal energy of about 14 GWh/year. Finally, a typical day of 21st June was simulated in both locations, using Dymola. Then, a timestep simulation of one minute was carried out in order to assess the transient behavior of the solar field. The solar field is composed of 18 PTC collectors. The number was calculated considering the available surface of the pasta factory of 2000 m 2 . In Figure 6, the inlet and outlet temperature of the solar field and outlet temperature of tank considering the HTF loop (Thermino SP) are shown. It is possible to see the effect of a cloud passing between 9 a.m. and 10 a.m. and at around 1 p.m., during which time the decrease in DNI causes a change in the pump flow rate and recirculation conditions. As a consequence, the inlet and outlet solar field temperatures of HTF drop at a 136 °C, while the HTF temperature at the outlet of the tank decreases only to 163 °C. It is possible to notice that the maximum temperature reached by the HTF Therminol SP is under 250 °C and is in accordance with the operating conditions to ensure a certain durability, as expressed in the related datasheet of the commercial product. The black line is the trend of the HTF at the tank outlet. Differently from the temperature trend of the solar field input (blue) and output (orange), the variations due to flow rate tuning and PID operation are reduced due to the thermal inertia induced by a tank that acts as a buffer to variations and to the DNI intermittences. Figure 7 shows the trend of thermal power of the solar field made up of the same numbers of PTC collectors. In orange is the thermal power related to the location of Palermo, while in blue is that related to the location of Molina. Furthermore, the difference of DNI between Palermo and Molina is also shown. Thus, it is possible to see that when this difference is equal to zero, as observed at around 7 a.m. (424 min), the thermal powers have the same value (about 440 kW).

Discussion
The main objective of this paper was to evaluate the integration of solar thermal power into a generic industrial process. As an example, the pasta-drying process of the Felicetti pasta factory was chosen as a case study. A flexible methodology was adopted, in which to make the analysis, it is possible to choose any geographical location and any industrial process. The methodology was based on the joint use of different programs and datasheets and data that provide information on geographical and process coordinates, such as operating temperature, working cycle and flow rate. It has been shown that solar thermal integration works differently with the different solar technologies considered, even in geographically complex contexts, and can play a key role in the massive decarbonization of industrial processes. In fact, the carbon-free solution, found with the integration of 100% of solar thermal power, makes it possible to avoid an emission of about 1500 tons of CO2/year. In order to decarbonize the industrial drying process of pasta completely with solar heat, the needed available space in Molina is 27,710 m 2 , while in Palermo, it is 13,300 m 2 . Moreover, in the pasta factory's available space of 2000 m 2 , in Molina, 99 tons of CO2/year are avoided with LFR technology, while in Palermo, 208 tons of CO2/year are avoided. Finally, through the dynamic simulation performed in Dymola, it was possible to verify the behavior of the solar field with sudden cloud shading. The clouds determine the sudden decrease in the DNI, which makes it necessary to use a control system on the pump flow rate and recirculation conditions. It was possible to verify that despite the small size of the HTF tank, it is possible to reach and to maintain the temperature required by the industrial process by appropriately tuning the flow rate and recirculation. Furthermore, it can be noted that the geographical location and topography have a significant influence on DNI and thermal production. In particular in Figure 6, it is possible to observe that with the decrease in the DNI to 140 w/m2, it is possible to bring the amount of thermal power produced by the solar field to zero. In regard to the one-minute simulation on 21st June in Palermo, there is a peak thermal power of 536 KW with a solar fraction of 83%, while in Molina, there is a peak of thermal power of 439 Kw with a solar fraction of 67%.