Numerical Analysis and Optimization of Solar-Assited Heat Pump Drying System with Waste Heat Recovery Based on TRNSYS

: In this paper, a new solar-assisted heat pump drying system with waste heat recovery and double water tanks (SCAHP) was established and the system optimized by TRNSYS and variable air volume experiment. The annual cumulative efﬁciency of the SCAHP ( COP ac ), buffer tank heating efﬁciency ( η BT ) and hot water storage tank heating efﬁciency ( η ST ) are all considered as optimization objectives, and this paper discusses the relationship between the three optimization objectives and studies the inﬂuence of hot air volume (q DR ), area of solar collector ( A SC ), inclination angle of solar collector ( I SC ) and volume of heat storage water tank ( V ST ) on system efﬁciency in the drying process. In order to explore the general rule of system optimization, numerical analysis and optimization were carried out in ﬁve typical cities in ﬁve climatic regions of China. The results show that the mixed variable air volume mode can increase the drying rate with less energy loss, and shorten the drying period by 28.6–33.3%. When the system surface body ratio (SBR) in Nanjing is between 3.1 and 4.1, the COP ac , η BT , and η ST can reach the maximum value simultaneously. It is estimated that the cost can be recovered in 5 years when the system conﬁguration parameters are optimized.


Introduction
Drying has a large number of applications in industry, agriculture and other fields, while drying is one of the most energy-consuming units [1]. In developed countries, about 10% of fuel is used for drying, and in China, drying energy consumption accounts for about 12% of the entire industrial energy consumption [2]. In some industries, the energy waste caused by drying is enormous, with Mujumdar [3] reporting that drying consumes up to 70% of the total energy in the manufacture of wood products, 50% of the total energy in the manufacture of finished textiles, and over 60% of the total energy in agricultural maize production. In order to develop sustainable energy sources, three important technological changes are required: energy conservation on the demand side, increased efficiency in energy production, and the replacement of fossil fuels with various renewable sources [4]. In this regard, the use of heat pump drying systems and solar heat collecting systems has improved energy efficiency and reduced consumption of fossil fuels.
The energy efficiency of traditional dryers is usually very low. Lawton J [5] points out that the efficiency of traditional dryers can only reach 35% at the highest. Queiroz R [6] found that heat pump drying saves 40% energy in comparison with using electrical heater dryers. The traditional dryer is not only inefficient, but also the quality of drying products is not up to the requirements. With the rapid development of industry and agriculture, vacuum drying, freeze drying, spray drying and other methods have been put forward one after another; however, these drying processes are highly energy-intensive and have the disadvantages of high economic costs or low drying rates, and the advantages of heat pump drying are further prominent. In some studies, heat pump drying is considered to provide better product quality under the premise of lower energy consumption [7][8][9][10][11]. The first batch of heat pump drying studies in the literature was completed by Hodgett [12] and Geeraert [13]. In recent decades, many scientists have studied how to improve the efficiency of heat pump systems from various aspects [14][15][16]. However, related studies show that heat pump technology still has problems such as using high-grade electric energy, the small drying scale of single heat pumps and high system maintenance costs [17,18].
The combination of multiple energy sources, especially the introduction of renewable energy, has become one of the current research hotspots. Solar energy has attracted much attention as the most readily available renewable energy source. Open solar drying is the earliest drying technology that has been applied [19], but its external drying conditions are difficult to be controlled, resulting in too-long drying times of materials or poor quality of dried products. The efficiency of solar direct drying of crops is between 20-40% [20].
To solve these problems, many scientists are trying to combine solar drying with heat pump technology.
Luo Huilong et al . [21] designed and constructed an air source heat pump-assisted solar energy drying system. Compared with conventional drying systems, the air source heat pump-assisted solar energy drying system has greater energy saving potential. Rahman et al . [22] economically optimized the evaporator and air collector of a solar-assisted heat pump drying system. Their economic analysis showed that using an air collector area of 1.25 m 2 and an evaporator collector area of 2 m 2 , the system had significant cost savings with a payback period of 4.4 years. In a recent work, Meltem [23] et al . designed and tested a new type of sustainable photovoltaic/heat-assisted heat pump drying. The results show that the average coefficient of performance of the heat pump can reach 4.18, according to solar radiation and ambient temperature, and an average electric efficiency of 12.27% and a maximum of 14.02% can be obtained.
However, in the drying system combined with solar energy and heat pump, there is often a mismatch between energy supply and demand and excessive waste heat waste. Energy storage technology can solve this problem well. The application of heat storage in the solar drying system is justified under the following three circumstances [24]: (1) the drying period can be prolonged by energy storage; (2) the excess energy at the peak of radiation can be stored to avoid excessive local drying; and (3) the temperature of dry air can be controlled to avoid damage to materials. Qiu et al . [25] used a solar-assisted heat pump drying system with heat recovery to dry 10 kg of carrots, and the results showed that the system saved about 40.5% of the energy. Wang et al . [26] studied the use of another form of solar-assisted heat pump drying system with heat recovery to dry 80 kg of mangoes, and the results showed that the efficiency of the system was improved by about 6%, and the waste heat recovery rate reached 41.7%.
On the basis of predecessors' research, this paper proposes a new solar-assisted heat pump drying system with waste heat recovery and double water tanks, where the waste heat recovery module can partially offset the variable air volume of energy loss. On the basis of experiment explores the influence of variable air volume on the drying process, and system efficiency, the whole-year numerical simulation is carried out based on TRNSYS, and the parameters and operation parameters of each component of the system are optimized. Figure 1 shows the schematic diagram of SCAHP. The system includes a solar collector (SC), air-water heat exchanger (AWEH), buffer tank (BT), hot water storage tank (ST), air-water heat pump (HP), water pump, fan and drying room (DR). The parameters of system components, operation parameters and corresponding variation ranges are shown in Table 1. Among them, the fan is a forced convection fan installed at the outlet of DR, so as to facilitate the variable air volume drying experiment and waste heat recovery.  All energy consumption of DR is provided by BT, the heat of BT comes from SC, HP, ST and waste heat recovery, the heat of ST comes from the SC and waste heat recovery, and the heat recovered from waste heat can be used in two ways: first, passing the waste hot air directly to the low temperature side of the HP through the fan; second, passing the waste hot air to the SC inlet side after being converted by AWEH.

Introduction of SCAHP
The SC and HP are connected in parallel, and the start and stop of the SC circuit is not directly related to the start and stop of the HP circuit. When the temperature of the BT is insufficient, the SC is first opened for heating, that is, pump 1 is turned on, and priority is given to the heat supply to the BT. The above can be expressed as that the hot water in the SC flows to the BT through DV 1, and then flows back to the SC through DV 3, CV 1, Pump 1, and CV 3 after coming out of the BT. When the temperature of the BT reaches the set value, but the temperature of the ST does not reach the set value, the hot water in the SC flows to the ST through DV 1, and passes through DV 2, CV 1, Pump 1 and CV 3 after coming out of the ST flow back to the SC. When the temperature of the BT and ST reach the set value, the SC circuit is closed. When the temperature in the BT is lower than a set value lower than the set temperature of the SC circuit, the HP and Pump 2 are turned on, and the hot water flows into the BT and returns to the HP until the BT reaches the set upper limit and closes the loop.
Waste heat recovery can be achieved by passing the waste heat to the inlet of the SC or the low temperature side of the HP through the fan. There are two points to note: 1, the carrier of waste heat recovery heat is air. When it leads to the low temperature side of the SC, it needs to be converted by a water-air heat exchanger. At the inlet of the SC, it is mixed with the return water of the original SC circuit, heated by the SC, and then led to the BT or ST. 2, When the heat recovered by the waste heat directly passes to the low temperature side of the HP, it may not be able to meet the flow demand of the low temperature side of the HP. At this time, the hot air from the waste heat recovery needs to be mixed with the air in the atmospheric environment, that is, pass through CV 4 to the low temperature side of the HP.
The heat supply from the ST to the BT is only related to the temperature difference between the two water tanks, and whether the BT reaches the set temperature is related. From the perspective of heat supply, the heat supply from the ST is in parallel with the heat supplied by the SC and HP. When the heating conditions are met, pump 6 is turned on, and the hot water flows from the ST to the BT and returns to the ST. It should be noted that the control logic of this system and the set temperature range of each loop will be described in detail in Section 2.2 Control Logic.  Figure 2a shows the control logic of SCAHP except for the waste heat recovery module. The figure also provides the temperature settings of the water tank in each heating mode mentioned in Section 2.1. According to the temperature value of the BT, the control mode of the system is divided into three modes.  Mode 1: If the average temperature of BT (T BT ) is lower than 50 • C, turn on the HP to heat the BT to 65 • C, and then turn off the HP. At the same time, it is judged whether the outlet temperature of SC (T SC-OUT ) or the outlet temperature of ST (T ST ) is greater than T BT . When the judgment value is 1, the heating circuit of the SC or ST is started to heat the BT to 55 • C. When T BT reaches 55 • C or T SC-OUT is less than T BT , judge whether T SC-OUT is greater than T ST . When the judgment value is 1, the SC transfers heat to the ST.

Control Logic of SCAHP
Mode 2: When T BT is greater than 50 • C and less than 55 • C, it is judged whether T SC-OUT or T ST is greater than T BT . When the judgment value is 1, the SC or ST is activated to transfer heat to the BT until T BT gets to 55 • C, When T BT reaches 55 • C or T SC-OUT is less than T BT , judge whether T SC-OUT is greater than T ST . When the judgment value is 1, the SC is started to transfer heat to the ST, otherwise SCAHP is in shutdown state, and T BT gradually drops and enters mode 1.
Mode 3: When the temperature of the BT is greater than 55 • C, it is judged whether the outlet temperature of the SC is greater than the temperature of the ST. When the judgment value is 1, the SC will be started for heating, otherwise the system is in shutdown state, and the temperature of the BT will gradually drop into mode 2. Figure 2b shows the control logic of the system's waste heat recovery module. The first thing to note is that when the waste heat recovery fan is not turned on, the waste heat mode is turned off and the system does not enter this control module. When the fan is turned on, first determine whether the HP is turned on. When the HP is turned on, the hot air obtained in the waste heat recovery and ambient air will be mixed and sent to the low temperature side of the HP. If the HP is in a shutdown state, judge whether the SC circuit is open, and whether the temperature of the hot water exchanged with the hot air recovered by the waste heat in the heat exchanger is greater than the inlet temperature of the SC; when the two judgments both return the value 1, the heat recovered by the waste heat is led to the entrance of the SC, otherwise the waste heat is directly discharged into the surrounding environment.

TRNSYS Simulation Model of SCAHP
According to the SCAHP physical model, control logic and parameter settings described above, a system simulation model is established in TRNSYS 17. TRNSYS (Transient System Simulation Program) is an extremely flexible graphically based software environment used to simulate the behavior of transient systems, and most TRNSYS-based simulations focus on evaluating the performance of thermal and electrical energy systems. The model contains weather file, SC, HP, BT, ST, DR, water pump, fan, heat exchanger and control strategy. The main components of the TRNSYS simulation model and their descriptions are shown in Table 2. The on/off differential controller generates a control function which can have a value of 1 or 0.

Drying room Type 682
This model simply imposes a user-specified load (cooling = positive load, heating = negative load) on a flow stream and calculates the resultant outlet fluid conditions.
When performing model verification, the time step is set to 1 min, which is controlled at the same magnitude as the time step of the verification object. In the simulation analysis of the whole year, the step length is set to 1 h. The system control strategy is implemented by Type 2b and calculator. The shunt valve at the outlet of the BT and ST is replaced by a calculator. The TRNSYS simulation model diagram is shown in Figure 3.

Research Method
This research adopts a combination of experiment and simulation, taking 731 g of carrots as the experimental material, and exploring the influence of variable air volume on material drying, suitable variable air volume control modes, drying rate and energy consumption data of different variable air volume modes in the variable air volume experiment. Exploring the optimal configuration and optimal operating parameters of the system in the simulation, input of the real-time data of the DR energy consumption is obtained through the experiment through Type 682. Taking COP ac , η BT and η ST as optimization goals, orthogonal experiments are used to optimize A SC , I SC and V ST . Finally, the economic performance of the system is also verified.
This study takes five typical cities in China's five major climate regions as the research objects, namely Harbin, Tianjin, Nanjing, Guangzhou, and Kunming, to explore the general laws of the system proposed in this paper. In the simulation calculation, the meteorological parameters are based on the typical yearly parameters of five cities. The typical meteorological year (TMY) is based on the monthly average of the past 10 years. The parameters of the five cities are shown in Table 3. Taking Nanjing City as the main analysis object, the analysis and optimization results of other cities will be summarized in the table and given uniformly. It should be pointed out that in our model only the difference in ambient temperature and solar radiation between different cities is considered.

Total Heat Required for Drying Process
The heat consumed by carrots during the drying process can be divided into three parts: The heat required for preheating carrots, the heat required for evaporation of water from carrots, and the heat loss of the entire system during the drying process. According to the first law of thermodynamics: The heat required for the carrot preheating stage is: where c r is the specific heat of the carrot under constant pressure (kJ/kg·k), m r0 is the initial mass of the carrot (kg), and T r0 is the carrot temperature before preheating ( • C), equivalent to ambient temperature, and T r1 is the carrot temperature after preheating ( • C), equivalent to the temperature after the constant temperature of the DR stabilizes. The heat required by carrots to evaporate water is: where γ is the latent heat of vaporization of the moisture in the carrot (kJ/kg), and m r2 is the mass of the carrot after the drying (kg). The heat loss during the drying process of carrot is:

Solar Collector
The total amount of radiation during the working hours of SC is: where A c is the area of the solar collector m 2 , I t is the instantaneous solar radiation value incident on the SC, W/m 2 , and t 1 is the total open time of the SC (h). The energy used by the SC to provide the DR is: where G is the total energy received by the SC (kJ), η sc is the conversion efficiency of SC, η BT is the heating efficiency of the BT, η ST is the heating efficiency of the ST, η AW HE is the heat exchange efficiency of the air-water heat exchanger, the air-water heat exchange efficiency of the system is set to a fixed value, 0.65, and Q SC−ST is the total energy obtained by the ST from the SC (kJ).
The total heating share (SF) of the system occupied by the heat supply of solar collectors can be expressed as:

Tank
The heat supplied by the SC to the BT is: The heat supplied by the SC to the ST is: where m sc is the mass flow of water passing through the solar collector (kg/s), c sc is the specific heat of water under constant pressure kJ/kg·k, and T scout is the outlet temperature of the SC ( • C), T scin is the inlet temperature of the SC ( • C), t 2 is the time it takes for SC to supply heat to the BT (h), t 3 is the time it takes for the SC to supply heat to the ST (h). The heat supplied by the HP to the BT is: where m HP is the mass flow of water passing through HP (kg/s), c HP is the specific heat of water under constant pressure (kJ/kg·k), and T HPout is the outlet temperature on the condenser side of the HP ( • C), T HPin is the inlet temperature on the condenser side of the HP ( • C), t 4 is the total open time of the HP (h). The heat supplied by the ST to the BT is: m HB c HB (T HBout − T HBin )dt (11) where m HB is the mass flow rate of water flowing from the ST to the BT (kg/s), c HB is the specific heat of water under constant pressure (kJ/kg·k), and T HBin is the outlet temperature of the ST flowing to the BT ( • C), T HBout is the outlet temperature of the BT flowing to the ST ( • C), t 5 is the time it takes for the ST to supply heat to the BT (h). The heat obtained by the BT is: The heat supplied by BT to DR is: where m BD is the mass flow of water from the BT to the air-water heat exchanger connected to the DR (kg/s), and c BD is the specific heat of water under constant pressure (kJ/kg·k), T BDin is the outlet water temperature of the BT flowing to the air-water heat exchanger connected to the DR ( • C), T BDout is the temperature entering the BT from the air-water heat exchanger connected to the DR ( • C), and t 6 is the time it takes for the BT to supply heat to the DR (h).
The heating efficiency of the SC is: The heating efficiency of the BT is: The heating efficiency of the ST is:

Heat Pump
The simulation system uses TRNSYS's own Type941 air-water heat pump. For detailed performance parameters of this heat pump, please refer to its own mathematical model.
The instantaneous COP and annual cumulative COP of HP are: where P HP is the instantaneous power of the HP (W), and t 7 is the annual cumulative number of hours, taking 8760 h.

System
The total power consumed by SCAHP is: where P pump is the sum of instantaneous powers of all pumps in the system (W), and P f an is the sum of instantaneous power of all fans in the system (W). The benefits of SCAHP are: SCAHP instantaneous COP and annual cumulative COP are:

Material Drying Characteristics Analysis
The specific moisture extraction rate can be expressed as (kg/kW h): The dry basis moisture content of the material can be expressed as: where m ri is the mass of carrot at time i (kg). The drying rate can be expressed as: where M t+∆t is the moisture content at time t + ∆t (g moisture/g dry matter), and M t is the moisture content at time t (g moisture/g dry matter).

System Economic Analysis
This paper uses the "life cycle" method and the "payback period" method to conduct an economic analysis of SCAHP.

"Life Cycle" Method
Life cycle cost refers to the sum of all expenditures in the entire life cycle of a product, including the acquisition of raw materials, product use costs, etc. It refers to the sum of the production cost of the enterprise and the user's use cost.
The cumulative net value of the dryer life, the cost savings during the life cycle: where V npt is the net present value of the dryer in a period of time calculated by year (t), and the expression is: where i is the inflation rate, taking 0.08, n is the long-term investment interest rate, taking 0.1, D is the number of days the dryer is used per year, and S d is the money saved per day in a period of time per year (RMB/d). It can be expressed as: where C s is the selling price of dried carrots (RMB/d), and C ds is the manufacturing cost per kilogram of dried carrots (RMB/kg). It can be expressed as: where C f d is the price of fresh carrots required to produce one kilogram of dried carrots (RMB/kg), and C d is the cost of using a dryer to produce one kilogram of dried carrots (RMB/kg). It can be expressed as: where C f is the price per kilogram of fresh carrots (RMB/kg).
where C a is the cost of a dryer (RMB). It can be expressed as: where C ac is the annual investment cost of the dryer (RMB), C m equipment annual maintenance cost (RMB), and C r is the equipment annual operating cost (RMB). It can be expressed as: where C cc is the investment cost of the dryer (RMB), and j is the life of the dryer (year). The annual operating cost of the system is: where C e is the electricity price per kWh (RMB/kW h).

Payback Cycle
The payback period (N) can be expressed as: where V 1pt is the net present value of the dryer in the first year (RMB).

Simulation Model Verification
This study uses the experimental results of Qiu [25] to verify the simulation model. The system parameters are consistent with Qiu. The meteorological parameters adopt the typical meteorological year parameters of Nanjing City. The DR inlet temperature and the DR outlet temperature are the verification objects. As shown in Figure 4, the simulated DR inlet temperature (T in−sim ) has an error of 0.07-6.25% compared with the experimental DR inlet temperature (T in−exp ), and the average error is 2.36%. Compared with the experimental DR outlet temperature (T out−exp ), the error of simulated DR outlet temperature (T out−sim ) is 0.04-10.57%, and the average error is 3.29%. Compared with the model proposed by Qiu [25], the model proposed in this paper has an additional air-water heat exchanger that connects the BT and DR. The heat exchanger can be regarded as an efficiency device in essence. As shown in Figure 5, when the heat exchanger efficiency is set to 1, the inlet and outlet temperature of the DR is exactly the same as when the heat exchanger is not installed. Set the heat exchanger efficiency to 0.8, 0.6, 0.4, and get the temperature lines of T in−0.8 , T out−0.8 parallel to T in−1 , T out−1 . Based on the above comparison and verification results, the simulation system is correct.

Variable Air Volume Experiment
In this paper, an electric heating thermostat was modified. On the basis of retaining the original drying room air outlet, a vent was opened, as shown in Figure 6, and a fan added to this vent to perform forced convection for variable air volume experiments and heat recovery. At this time, the original air outlet has become an air inlet, as shown in Figure 7.  In this experiment, 731 g of fresh carrots were dried in two drying modes: the forced convection fan was turned off mode (TOFF) and the forced convection fan was turned on mode (TONF). The flow rate of the fan is 39 kg/h. The carrots are equally divided into two plates, and the internal dimensions of the incubator are 0.45 m × 0.45 m × 0.4 m. The initial moisture content of drying is 9.15 g moisture/g dry matter, and the final moisture content is 0.25 g moisture/g dry matter. Therefore, when the weight of dry carrot reaches 0.09 kg, the drying process ends. The equipment used in the experiment is shown in Table 4. In the experiment, the weight of carrot was measured every half an hour, and the moisture content and drying rate of carrot in the two modes were shown in Figures 8 and 9.  It can be seen from Figure 8 that the drying rate of carrot in the two modes decreases with the increase of time. The drying rate of the TONF mode is 3.5 h faster than that of the TOFF mode, saving 33.3% of the drying time. This is due to the forced convection in the DR, and the humidity is maintained at a low level during the drying process. Natural convection prevents the drying moisture during the drying process from being discharged in time and stays in the DR, making the air humidity close to saturation and inhibiting the drying rate.
It can be seen from Figure 9 that starting from the 4th hour, the drying rate of TONF mode (the forced convection fan was turned on mode) is lower than that of TOFF mode (the forced convection fan was turned off mode). This is because the moisture content of carrots in TONF mode has been reduced to 1.03 g moisture/g dry matter at this time, and the moisture content of the carrots in TOFF mode corresponds to 4.29 g moisture/g dry matter, the latter is much higher than the former, and the lower moisture content of carrots reduces the drying rate.
It can be seen from Figure 10 that the drying energy consumption in TONF mode (the forced convection fan was turned on mode) is 28.26% higher than that in TOFF mode (the forced convection fan was turned off mode). When TONF mode is switched to TOFF mode in the 4th hour, it can save about 20% of energy compared to TONF mode, and only about 8% higher than TOFF mode. In combination with Figures 8-10, we recommend using the mixed mode from TONF to TOFF in the 4th hour (TOCF mode). The drying rate is faster and the energy consumption is low. The following discussion of the article will be based on the TOCF model.

System Optimization
A solar collector (SC) can effectively improve the overall performance of SCAHP. However, due to high initial investment and different climate conditions in various regions, there are certain difficulties in system design. At the same time, as a key component of the system coupling, a hot water storage tank (ST) also has some difficulties in design. There are few researches on the optimization of heat pump drying system configuration. This research tries to give some suggestions.
As shown in Figure 11, the annual cumulative efficiency of the system (COP ac ) shows an increasing trend with the increase in the area of solar collector (A SC ). This is because as the A SC increases, the proportion of solar heating in the overall heating of the system increases, and the proportion of air-water heat pump (HP) heating is reduced, the total energy consumption of drying room (DR) remains unchanged, the power consumption decreases, and the COP ac increases accordingly. In this configuration, COP ac is 3.571~4.805. As the volume of heat storage water tank (V ST ) increases, COP ac gradually slows down as the A SC increases. When the V ST is 0.02 ∼ 0.1 m 3 , as the A SC increases, the COP ac with a larger V ST system will always surpass the system with a smaller V ST at a certain point. When the V ST is 0.12 ∼ 0.4 m 3 , as the A SC increases, the COP ac shows an overall downward trend. This is due to the large V ST , the temperature rise of ST is slow, and the heat in the ST cannot be effectively and timely passed to the BT. As a result, on the one hand, the heat collection pump is idling, which increases energy consumption. On the other hand, it cannot effectively reduce the heat supplied by the heat pump and the solar collector cannot be used to the maximum extent, which leads to a decrease in the COP ac . Figure 11 also reveals that when the A SC is 0.1 ∼ 0.25 m 2 , the COP ac shows a decreasing trend as the V ST increases, and when the A SC is 0.3 ∼ 0.8 m 2 , the COP ac increases first and then decreases as the V ST increases; this is because the smaller V ST cannot fully release the heat supply potential of the larger A SC , and as the A SC increases, the V ST corresponding to the highest value of the system efficiency also increases. When A SC is 0.4, 0.5, and 0.6 m 2 correspond to ST with volumes of 0.04, 0.08, and 0.1 m 3 , respectively, so that the COP ac reaches the corresponding maximum value. In addition, it can be seen from the figure that when the A SC is 0.1 ∼ 0.4 m 2 , the distance between the lines is relatively sparse. Within this range, as the A SC increases, the COP ac increases rapidly. When the A SC is 0.45 ∼ 0.8 m 2 , the distance between the lines is relatively dense. In this range, as the A SC increases, the system efficiency increases slowly, which shows that when the A SC is increased to 0.45 m 2 , the system has been mainly heated by the SC, and the HP has played a smaller role. It is recommended to use solar collectors with an area less than 0.45 m 2 in this system. In summary, the coupling of the SC and ST is very important. Different A SC require designing a ST of appropriate volume. In order to find this appropriate value, we introduce the surface-to-body ratio (SBR): SBR = A op / V op + V BT , Among them, A op and V op respectively represent the A SC and the V ST corresponding to the optimal COP ac . V BT is the volume of BT, taking 0.05 m 3 . It can be calculated that the SBR value is concentrated between 3.182~4.091 under the premise of the recommended value.
As shown in Figure 12, when the V ST is 0.02 ∼ 0.34 m 3 , the η ST increases first and then decreases as the A SC increases. When the V ST is 0.36 ∼ 0.4 m 3 , the η ST shows an increasing trend as the A SC increases. As the V ST increases, the A SC corresponding to the maximum η ST becomes larger. This is because the larger the V ST , the slower the temperature rise, and a larger A SC is required to provide sufficient heat.  Figure 12 also reveals that when the A SC is 0.1 ∼ 0.55 m 2 , the η ST shows a decreasing trend as the V ST increases, when the A SC is 0.6 ∼ 0.8 m 2 , the η ST will first increase and then decrease as the V ST increases. This is because the V ST is too small compared to the A SC , and it cannot store enough solar collector heat. As a result, the heat generated by the SC cannot pass into the BT in time, and the HP needs to be turned on to provide heat to the BT, and the heat generation potential of the SC is not fully utilized. When the A SC is 0.45 ∼ 0.8 m 2 , the distance between the lines is relatively dense. In this range, as the V ST increases, the η ST decreases slowly. This shows that when the A SC increases to 0.45 m 2 , it is enough to meet the energy demand of the system. The decrease in the η ST is mainly caused by the increase in V ST . This shows that more heat from the SC is wasted, which confirms the analysis result in Figure 11.
To sum up, take the η ST as the goal, and based on the effective value of this calculation, the best SBR can be calculated from 1.923 to 4.444, and the A SC corresponding to the maximum η ST is greater than 0.3 m 2 . It can be seen that the η ST does not match the overall efficiency of the system. It is recommended to choose a small V ST after selecting the A SC , and increase the η ST to avoid too much heat from the solar collector being wasted. Combining increasing the A SC and V ST will also increase the initial investment, and it is recommended not to take the η ST as the first choice for system design.
As shown in Figure 13, the η BT shows a trend of increasing first and then decreasing as the A SC increases; this is because when the A SC is relatively small, the heat of the BT comes from the HP.The HP heats the BT to 65 • C, and during the process of lowering the temperature of the BT, the heat loss is large, which causes the η BT to be low. As the A SC increases, the heat received by the BT from the SC increases, and the BT will no longer rise to 65 • C frequently, thereby reducing heat loss and improving η BT . With the further increase in the A SC , enough heat can be stored in the ST. When the system shuts down at night, the heat from the ST will be transferred to the BT, and then it will interact with the surrounding atmosphere in the BT. The heat exchange causes waste, and the η BT is gradually reduced. This problem will be solved in a non-stop system.  Figure 13 also reveals that when the A SC is 0.1 ∼ 0.45 m 2 , the η BT shows a basically unchanged trend with the increase of V ST , during the whole process, the η BT changes at most by 1.49% and at least 0.44%, which can be ignored. This is because the SC cannot provide sufficient heat for the BT to lose heat at night. When the A SC is 0.5 ∼ 0.8 m 2 , the η BT decreases first and then increases as the V ST increases, this is because SC can provide sufficient heat at this time, the ST has a small volume and a high temperature. When it is shut down at night, the ST can provide more heat to the BT, resulting in greater heat loss, as the V ST increases, the temperature of the ST decreases the amount of heat transfer is reduced, the heat loss is also reduced, and the heating efficiency of the BT is gradually improved. In summary, the system configuration required for the maximum efficiency of the BT is different from the system configuration required for the maximum efficiency of the ST.
As shown in Figure 14, the heat provided by the SC shows an increasing trend as the A SC increases. When the A SC is 0.1 ∼ 0.65 m 2 , the heat supply of the SC shows an increasing trend as the V ST increases, when the A SC is 0.7 ∼ 0.8 m 2 , the heat supply of SC will decrease first and then increase as the V ST increases. This is because when the A SC is small, as the V ST increases, it can absorb more heat from the SC, and when the A SC is larger, the smaller V ST quickly heats up, which can make better use of the heat recovered from the waste heat, causing the SC to decrease the heat first. The heat provided by the HP shows a decreasing trend as the A SC increases, when the A SC is 0.1 ∼ 0.25 m 2 , the heat supply of the HP shows an increasing trend as the V ST increases, when the A SC is 0.3 ∼ 0.8 m 2 , the heat supply of the HP decreases first and then increases as the V ST increases; this is because when the A SC is 0.1 ∼ 0.25 m 2 , the system mainly relies on the HP heating. As the A SC increases, the system mainly uses the SC for heating, the heat supplied by the HP decreases in the initial stage. Later, as the V ST increases, the η ST decreases, and the BT needs the HP to provide more heat, which is consistent with the analysis in Figure 12. When the A SC is 0.15 ∼ 0.35 m 2 , the heat supply ratio of the SC is 34.14~68.03%, which is reasonable. Based on the above analysis, water tank efficiency and heat supply analysis, we suggest that the A SC of this system is 0.15 ∼ 0.35 m 2 , the SBR is between 3.2~4.1, and the V ST is 0.02 ∼ 0.06 m 3 . The recommended values for other climate zones are shown in Table 5. The following analysis in this paper will be based on Nanjing, with a A SC of 0.3 m 2 and a V ST of 0.04 m 3 . As shown in Figure 15a, the difference in heat supplied by the SC with different I SC in 12 months ranges from 4.83 to 58.27%. October is the month with the smallest difference and June is the month with the largest difference, and the annual average difference is 23.2%. It can be seen that the design of the I SC plays a key role. The annual heat provided by the SC increases first and then decreases with the increase of the I SC . When the I SC is 30 • , the total heat provided by the SC throughout the year is the largest. If the SC is designed to have no tilt angle change throughout the year, it is recommended to select a tilt angle of 30 • . As shown in Figure 15b, when Nanjing is in January to February and November to December, and the I SC is 50 • , the most heat is collected. From March to October, the solar collector can collect the most heat when the I SC is 15 • . Therefore, we recommend changing the I SC twice a year. In March, change the I SC to 50 • , and in November Change the I SC to 15 • until March of the next year. The configuration parameters of I SC in other cities are shown in Table 5.

Economic Analysis
The price of material is an important factor affecting the economic analysis of the drying system. According to the heat conservation relationship and calculation, using the life cycle method, the total amount of dried carrot in the system throughout the year is 45.99 kg. Table 6 shows the annual net present value of Nanjing and the cumulative net present value of the system life. The life cycle method is used to determine the dry payback period of carrot as 5 years. The cumulative net present value is 13,654 yuan respectively. Table 6. Annual revenue diagram of Nanjing drying system.

Year
Annual Net Present Worth Cumulative Net Present Worth

Conclusions
The energy analysis of SCAHP is carried out. The energy coupling relationship of the system in different climate regions is established, and the optimization model of the system is discussed. It can be concluded: (1) It is recommended to use the mixed variable air volume mode (TOCF mode) for drying, which can save about 20% of the drying time; (2) The system configuration to achieve the maximum system efficiency is not consistent with the system configuration to achieve the maximum water tank efficiency. It is suggested to choose a smaller V ST after selecting A SC to increase the heating efficiency of the ST and improve the system energy utilization rate. (3) It is suggested that the SBR of SCAHP in Nanjing should be between 3.182 and 4.091, and I SC should be changed twice a year. In March, I SC should be changed to 50 • , and in November, I SC should be changed to 15 • until the following March. The configuration parameter recommendations for the other cities are shown in Table 5. (4) Taking Nanjing as an example to analyze the economy of SCAHP, it can be concluded that the recovery cycle is 5 years, and the system is economically feasible.
Author Contributions: Conceptualization, Z.X. and Y.G.; methodology, Z.X.; software, Z.X. and C.Y.; validation, Z.X. and C.Y.; formal analysis, Y.Y. and Y.L.; investigation, Y.Y. and Y.L.; writing-original draft preparation, Z.X.; writing-review and editing, Z.X. and Y.G. All authors have read and agreed to the published version of the manuscript.