The influence of geometrical parameters on heat and mass transfer processes in rotary heat exchangers

This paper focuses on the influence of geometrical parameters of the rotary heat exchanger on its operation under high-speed rotor conditions. Proposed mathematical model of heat recovery wheel is based on a structure of the counter-flow heat exchanger model. After implementation of a numerical method on the modified α-model, the computer simulations were conducted. They allowed considering the distribution of the active heat and mass transfer zones (“dry,” “wet” and “frost”) on the matrix channel surface depending on the outdoor air conditions for different variants of rotor size. The obtained results indicate that the increase in the wheel’s depth leads to the increase in temperature effectiveness of heat recovery (from 0.667 to 0.814). Moreover, it was also established that the threshold temperature under which the frost is accumulated on the core surface also rises with the rotor depth (from − 10.7 to −9.6 °C). It was concluded that under certain outdoor air temperature conditions (t1i> –10.7 °C, for rotor depth equal to 0.40 m), the mass transfer rate of the condensed water vapor in the channels of the return air side is equal to the evaporation mass transfer rate in the channels of the supply air side, and hence, the temperature efficiency reaches the same value level as under “dry” operating conditions. It was established that only under frost accumulation conditions, the increased temperature effectiveness is observed.


Introduction
Nowadays, in the era of a big economic development, the significant increase in global energy consumption is observed. The expected development of the building sector, next to the industry and transport, has become responsible for this state. It is estimated that the level of buildings energy use constitutes about 40% of total primary energy consumption in developed countries [1]. Therefore, the legal regulations related to the reduction of energy demand have appeared. In modern buildings, energy is used mostly to ensure the required indoor air parameters [2]. Thermal insulation improvement allows reducing the static heat loss or heating gain in a building. However, due to the stringent indoor air quality requirements, more fresh air has to be supplied using mechanical ventilation systems. It causes an increase in energy demand for ventilation [3]. Energy saving can be achieved by using high-efficiency HVAC systems, including exhaust air heat/energy recovery exchangers. Different types of heat recovery devices such as: fixed-plate exchangers, rotary exchangers, coil energy recovery loops and heat-pipe heat exchangers were compared in ASHRAE Handbook [4]. However, only fixed-plate and rotary heat exchangers have become widespread and are being used by more and more customers. Thermal wheels constitute a group of the most often used heat recovery exchangers in the HVAC systems in Poland. One of their major advantages is high efficiency exceeding even 80% [5]. The reason for achieving such a high level is the potential to recover not only sensible, but also latent heat. In the case of water vapor condensation in the return airflow channels, this generated moisture can be assimilated by the supply air flowing through the same channel in the next part of the cycle. For similar reasons, many manufacturers produce rotors covered with hygroscopic material, which allows for a more intensive heat and mass exchange between two airflows. Such devices are called energy wheels.
In order to analyze the effectiveness of rotary heat exchanger, their construction and working principles should be presented. The rotary heat exchanger ( Fig. 1) consists of a rotor made of 0.06-0.10 mm thick aluminum foils [6], which constitute a storage mass. Alternately wrapped wavy and straight tapes (sheets) form a large number of small axially arranged channels through which both airstreams flow. The height of these channels varies from 1.4 to 2.7 mm [6,7]. As the wheel matrix rotates with specified speed between two air channels (supply and return), it transfers the heat from warmer to cooler airflow. The rotation is caused by the electric motor with belt pulley and v-belt connected to the rotor. The whole construction is fitted in a metal housing.
In winter, the warm and humid exhausted air passes through one half of the matrix, transmitting the majority of its heat to the storage mass, while the cold and dry outdoor air flows through the other half and receives heat accumulated in the previous part of the revolution. This cycle repeats as long as the wheel rotates.
Potential of air-to-air heat recovery exchangers has caused them to be the subject of many considerable and detailed studies. For many years, researchers have considered the correct operation of heat exchanger. There are a lot of research papers including analysis of heat transfer processes (under dry and phase transition conditions) [8][9][10][11][12]. In addition, many authors analyzed also the influence of different thermodynamic and geometrical parameters on heat recovery efficiency. In order to effectively predict the behavior of these heat exchangers under different climatic conditions, several numerical methods were developed. The researchers have also raised the issue of freezing problem.
Justo Alonso et al. [13] emphasized the need for using mechanical ventilation systems with heat recovery devices while designing Zero Energy Building in cold climate. They focused on comparison of the properties for different types of heat/energy recovery systems: recuperative and regenerative heat exchangers and also regenerative energy exchangers. Both heat and energy wheels are more susceptible to air leakage than recuperators since the latter ones have both air streams separated. At the same time, they concluded that heat wheels are characterized by lower risk of frosting than the flat-plate heat exchangers. The researchers also presented a variety of methods that allow determining the sensible and latent effectiveness of the rotary heat/energy exchanges. Fathieh et al. [14] analyzed the possibility to predict sensible effectiveness of air-to-air heat wheel on the basis of cyclic and single-step change transient experiments performed on the parallel-plate heat exchanger made of the same material. Their results indicated that this method may be successfully used to determine the sensible effectiveness of the parallel-flow rotors. In the case of counter-flow wheels, predicted values differed from the real efficiency, which made this method unreliable. In papers written by Abe et al. [15,16], the analytical model for estimating the effectiveness of energy wheel was also developed. The model uses the characteristics measured on the same non-rotating wheel operating under a step change in temperature and humidity. Then the predicted results of the sensible and latent effectiveness were compared with the experimental standard test data and the numerical simulation results achieving agreement within the uncertainty boundary. Another study [17,18] shows correlations which may be used to determine sensible, latent and total effectiveness of energy wheel on the basis of known operating conditions. The presented equations use the dimensionless group for heat and moisture transfer. Simulation data confirmed the reliability of these correlations. Seo et al. [19] described a simple model which allows receiving the effectiveness of heat wheels with a maximum error of 5% within a wide range of rotational speed. It was based on the governing equations of a periodic-flow heat exchanger. They indicated that this method is very useful in practical applications.
The researchers have also developed many theoretical models of heat exchangers to investigate heat and mass transfer processes in detail as well as the impact of different physical and thermophysical parameters on the heat recovery efficiency. Büyükalaca and Yilmaz discussed the influence of rotational speed on the effectiveness of rotary heat exchangers [20]. The authors focused especially on low-speed operating conditions. As a result of their numerical and experimental studies, an empirical equation to calculate the effectiveness at low values of rotational speed was obtained. The authors also proved that the empirical equation of Kays and London [21] is not valid in the case of low-speed region. Nizovtsev et al. [10] studied the processes of condensation and vaporization in the regenerative air-to-air heat exchanger with periodic change in the air flow direction. For this purpose, the physical and mathematical model was developed. The authors analyzed the influence of relative humidity of the rooms on heat and mass transfer processes. They paid particular attention to the three characteristic entire indoor RH ranges that differed in moisture transfer processes (RH < 30%, RH = 30 … 80%, RH > 80%) . In the first range of relative humidity, there is no mass transfer. The second range of indoor RH is characterized by balanced condensation and vaporization processes, while Heat recovery section with a rotary heat exchanger: 1metal housing, 2-rotor, 3-electric motor, 4-drive belt, 5return air side, 6-supply air side the moisture accumulation in the channels can only be observed in the third range of RH. The analysis of the influence of relative humidity is a very important issue of the operation of rotary heat exchanger. This makes it possible to determine the conditions responsible for freezing processes. Frost formation on the plate surface of the air-toair heat exchangers is a common problem leading to the decreasing performance of heat recovery systems operating in cold climates. In the past three decades, a number of researchers have sought to determine the operating conditions which frost can occur inside the heat/energy exchangers.
A lot of papers also contain the analysis of frost control strategies and defrosting methods. For example, in 1989 Holmberg [9] analyzed frosting limit in both nonhygroscopic and hygroscopic surfaces of energy wheels. In the case of non-hygroscopic rotary heat exchanger, the core was made of untreated aluminum. Hygroscopic exchanger was additionally covered by an oxide layer. The results showed that frosting limit of a hygroscopic wheel is even 10 °C lower than the limit of non-hygroscopic rotor. Laboratory tests also showed that it took twice more time to cause 50% increase in pressure drop while operating non-hygroscopic wheel. Ten years later, Bilodeau et al. [22] investigated freezing in the rotary heat and mass exchanger. Experimental part of their studies evidenced that glazed frost dominates over rough frost in these types of exchangers. They also developed the transient threedimensional model, which allowed indicating that the absolute humidity of the exhaust air is the most important parameter determining the frost formation. Additionally, they showed nonlinear variation of the AH (absolute humidity) and the temperature of the air streams flowing through the frost zone. The authors also made a conclusion that the traditional frost control strategies should be replaced by the new ones to improve the performance of the heat recovery systems. Nasr et al. [23] raised the issues of frost formation and freeze protection methods in air-to-air heat and energy exchangers. For this purpose, a detailed literature review was conducted. On its basis, it may be noticed that the problem of frost formation was studied mostly for flat-plate heat exchangers. Moreover, their paper also contains the comparison of the different frost control strategies and defrosting techniques including capital and operating costs. The authors added that the problem of frost formation had been a subject of studies for over 30 years; however, they stressed the need for more detailed researches, especially in the field of frost accumulation on surfaces of both rotary and flat-plate exchangers.
In recent years, several attempts have been made to determine the frosting limits. Anisimov et al. [11] performed the simulation based on the mathematical model to investigate both heat and mass transfer processes in a cross-flow exchanger operating under frost conditions. The influence of the relative humidity of return air on the temperature effectiveness has been proven-the increase in the effectiveness with increasing RH is a result of the effect of heat of phase transition. The authors also proposed safe operating conditions for different performances of the considered heat exchanger. Moreover, they noticed that the return air parameters with the value of dew point temperature equaled 0 °C are the variants with the most unfavorable threshold temperatures. Jedlikowski et al. [12] presented a modified α-model, which allows analyzing the heat and mass transfer processes inside a counter-flow heat exchanger. It was underlined that the formation of different heat and mass transfer regions depends mostly on the relative humidity of exhaust air and outdoor air temperature. The authors focused on the frost formation phenomenon and presented frost threshold temperatures for various operating conditions, on the basis of which the analysis of optimal bypass frost control strategies was conducted. Moreover, it was noticed that even the fully opened bypass damper is not a sufficient method under some sub-zero outside air temperature operating conditions. A continuation of the topic of frost protection methods for a counter-flow heat exchanger can be found in other paper [24]. Pacak et al. conducted the analysis of power demand for two commonly used techniques: bypassing and preheating the outdoor airflow. Their investigation was based on a theoretical model and examined for different return air parameters and values of heat exchanger effectiveness. The results revealed that preheating is characterized by higher heat recovery efficiency within the most typical return airflow conditions and lower total power demand. In the paper written by Liu et al. [25], the frost formation problem in both cross-flow heat and energy exchangers was presented. The aim of this research was therefore to determine the frosting limits and analyze the impact of operating parameters on these limits. Another important practical information was defining the sensible and latent effectiveness for the energy exchanger to avoid frosting. Liu et al. [26] constructed and tested a prototype of a novel air-to-air quasi-counter-flow membrane heat exchanger. The analysis concerned heat and moisture transfer for low operating temperatures. For this purpose, the analytical model was developed. The heat exchanger reached a high sensible and latent effectiveness. This study also showed that the effectiveness level does not depend on the outside air temperature as long as there is no condensation and frost. The third paper of the authors [27] comprised a model, which allowed defining the frost threshold temperature.
The literature review allows drawing the conclusion that recently the most widely studied devices are the fixed-plate heat exchangers (both counter-flow and crossflow) and the enthalpy rotary wheels. However, far too little attention has been paid to the non-hygroscopic wheels. In this paper, the influence of geometrical parameters on heat and mass transfer processes and effectiveness of rotary heat exchanger is presented.

Theoretical analysis of heat and mass transfer process
The possibilities of creating different areas of active heat and mass transfer inside the rotary heat exchanger require thorough analysis of its operation (Fig. 2). One of the most frequently considered cases of year-round operation of heat exchangers under high-speed rotor conditions is the variant of "dry" heat exchange. During such operation, only "dry" heat exchange takes place in the channels of heat exchanger as a result of which both sides of the airflow channels remain dry. This situation occurs mainly during the summer. The processes of heat exchange during the cold season (especially in winter) are completely different. Particularly noteworthy are the conditions of extremely low ambient air temperatures. The high temperature difference between the two airflows results in a significant reduction in the plates surface temperature of the heat exchanger. If the temperature of the plate surface is higher than the dew point, the "dry" heat exchange mode takes place X 2 = 0.00 … 0.34 . But, if the temperature of the plate surface is lower than the dew point, the condensation process begins, as a result of which an additional water film may be formed X 2 = 0.34 … 0.61 . However, if the plate surface temperature additionally drops below the freezing point, the condensed water will gradually freeze X 2 = 0.61 … 1.00 . It is worth noting that the water vapor condensation processes are accompanied by sensible and latent heat transfer. In this way, not only sensible but also latent heat can be recovered. An important problem connected with the operation of these heat exchangers is the danger of freezing the accumulated condensate, which can cause an increase in pressure drop and consequently block the air flow through the heat exchanger. Moreover, during the analysis of the active heat and mass transfer zones creation, it is also worth mentioning the possibility of formation of the so-called "transient" area between "wet" and "frost" areas X 2 = 0.61 … 0.73 . The presence of this area results from the fact that there are parts of surfaces covered with water film and frost layer. A detailed discussion of this phenomenon is presented in our earlier papers [11,12,24,28]. The possibility of creating different heat and mass transfer areas depends on the value of plate surface temperature in relation to the dew point and freezing point temperatures. The conditions for the formation of the three main zones are summarized in Table 1.
The convection heat transfer coefficients for "dry" heat exchange conditions and for some "wet" heat exchange cases can be calculated on the basis of the Nusselt number with uniform heat flux for both thermal entrance region and fully developed flow. Two air streams are considered to be laminar flow [28], due to the low airflow velocity (Re < 2000). In this case, the Nusselt number of air can be divided into two areas: the beginning undeveloped region (L ≤ l b ) and for the fully developed region (L > l b ). The length of the first undeveloped zone can be calculated from Eq. (1): For the first area, the Nusselt number can be determined by means of empirical Eqs. (2) or (3), depending on the boundary conditions [28]: For the fully developed air flow region, the Nusselt number is constant [29][30][31]. For the flat channel: In the case of sine passage [29][30][31][32][33][34][35][36][37]: In the case of triangular passage [29][30][31][32][33][34][35][36][37]: Calculations of local values of Nusselt numbers were carried out taking into account the variable height of the channel caused by the variable thickness of frost.

Description of the mathematical model
The mathematical model of rotary heat exchanger is based on the structure of the counter-flow heat exchanger model [12,24]. Depending on the value of rotor speed or solidto-airflow heat capacity rate ratio, the effectiveness of the heat recovery depends on the two characteristic operating states: high speed (non-oscillation mode) and low speed (oscillation mode). For the values of rotating speed lower than the critical speed n r < n critical r , the plate surface temperature will oscillate. A gradual decrease in rotor speed will therefore increase the amplitude of the plate surface temperature oscillation. Increasing the amplitude of the plate surface temperature oscillation reduces the driving force of the heat exchange processes, which in turn reduces the heat exchanger effectiveness [19,29]. Analyzing the influence of rotational speed on the heat transfer process, Kays  where o temperature effectiveness of the counter-flow heat exchanger, correction factor for rotational speed.
The correction factor for rotational speed can be determined by the following equation [21]: where W r solid-to-airflow heat capacity rate ratio.
The solid-to-airflow heat capacity rate ratio can be calculated by means of the formula: where W r solid heat capacity rate, W 1 supply airflow heat capacity rate.
A similar formula, in agreement with the data presented by Kays and London, has been proposed by Büyükalaca and Yilmaz [20,38]: It should be noted that in the case of periodic steady-state operation, the plate surface temperatures in the direction of air flow are similar to temperatures determined for the channels of the counter-flow heat exchanger.
However, an increase in rotational speed above the critical value n r > n critical r ≈ 5rev∕ min [29].
(W r > W r critical ≈ 5 [19,20,29]) has practically no effect on the temperature effectiveness of the heat exchanger. In this case, the local temperature of plate surface in each cross section of the rotor remains practically constant during rotor rotation t p1 ≈ t p2 ≈ t p . For this reason, the temperature profile of matrix in the direction of air flow would be similar to the temperature profile of the counter-flow heat exchanger with temperature effectiveness o .
Based on the main assumptions described above for the operation of the rotary heat exchanger under nonoscillating conditions, the basic equation for high rotor speeds n r > n critical r can be significantly simplified and presented as simple differential equations of the heat and mass balance developed for supply and return airflows with a first-order solution.
The proposed modified α-model of the rotary heat exchanger is considered according to the system of coordinates X 1 (outdoor airflow direction) and X 2 (return airflow direction) (Fig. 2). In order to simplify the structure of model, a number of basic assumptions were introduced, including those based on models available in the published papers [11,12,24,28]: • The heat exchanger operates in quasi-steady-state conditions. • The outdoor and return airflows in contact with the plate surface are treated as an ideal, homogeneous and incompressible gas. • The driving force of the mass exchange process is the moisture content gradient. • The heat and mass transfer inside the rotor matrix take place in the normal direction (α-model [29,30]). • There are no additional heat sources in the airflows. • Water vapor condensation may initially form on the part of plate surface located on the return air side (semi-circular geometry), and then as a result of rotor rotation, the condensate will be transferred to the outdoor air section. • The evaporated heat flux may in no case be greater than the heat of water vapor condensation. • The temperature of the airflows varies according to the directions of the coordinate system. • Heat losses to the ambient air are ignored.
On the basis of the above-mentioned assumptions, the equations of balances for the outdoor and return airflow and for the matrix of the rotary heat exchanger were developed.
Energy balance for the outdoor and return air streams under "dry" heat transfer conditions can be expressed as follows: In this case, only sensible heat transfer takes place inside the channels of rotary heat exchanger so that the moisture content of both air streams will be constant and unchanged x 1 = const and x 2 = const .
The set of equations for heat exchange processes under water vapor condensation conditions is much more complicated. In such case, the sensible and latent heat transfer may occur at the same time. For this reason, temperature changes will also be accompanied by changes in the moisture content of the air stream. Mass balance for the return and outdoor air streams under water vapor condensation conditions can be presented in the form: The component Δ̄e in Eq. (18) means the relative duration of the evaporation period of the condensed water film (or frost sublimation) in the outdoor air flow section. It is worth mentioning that this component refers to the condition from the model assumptions concerning the possibility of water evaporation (or frost sublimation) in an amount not exceeding the evaporating mass flow rate from the return air channel.
The relative duration of the evaporation period of the condensate can be calculated by means of a set of equations: Water vapor mass flux on the matrix surface for return m 2 and outdoor m 1 air channel necessary for determining the relative duration of the evaporation period of the condensate can be calculated as follows: It should be noted that if the water vapor mass flow rate in the return air channel exceeds the water vapor mass flow rate in the outdoor air channel, a gradual accumulation of water film or frost layer can be observed on the matrix surface. Such operating conditions of the heat exchanger can be very dangerous, especially when the outdoor air temperature is below freezing. In this case, there is a high risk of frost formation inside the heat exchanger passages, which reduces the effectiveness of the unit and increases the pressure drop in the airflow channels as a result of blockage by frost [11,12,28].
Energy balance for the outdoor and return air streams under water vapor condensation conditions can be presented as follows: • for the case of "wet" heat exchange • for the case of "frost" heat exchange The solution of the presented sets of simultaneous differential Eqs. (14)- (18) and (22)-(25) requires the introduction of inlet airflow conditions at the entrances to the outdoor and return airflow passages (Fig. 2) and boundary conditions for supply and exhaust air side matrix surface: The presented sets of equations describing the heat and mass exchange in the matrix of the rotary heat exchanger are nonlinear and cannot be solved using analytical methods. For this reason, it was decided to use numerical methods based on the modified Runge-Kutta method. This method has an adequate accuracy and stability, which has been confirmed by solving similar problems [11,12,24,28,[39][40][41]].

Mathematical model validation
In order to validate the mathematical model of the rotary heat exchanger, a test stand was built (Fig. 3) (Fig. 4). The results of the comparison are presented in the form of the changes in supply air temperature expressed as a function of the outdoor air temperature for the air flow rate of 500 m 3 /h. As one can see, deviations between the results of measurement and simulation are in the range − 10… + 14%. The visible slight differences are mainly due to non-uniform airflows distribution, system leakage and heat losses. The presence of frost was determined on the basis of measurements of pressure drops in the heat recovery system, causing an increase in resistance to flow through the exchanger. Good agreement of experimental and simulation results confirms the adequacy of the developed model, so it can be concluded that the mathematical model allows to effectively predict the efficiency of rotary heat exchanger operation.  are presented on the charts (Figs. 5,6,7,8,9,10,11,12,13,14) for each variant.

Results of simulations
On the basis of the chart (Fig. 5) it can be seen that within the considered range of outdoor air temperature different zones of heat and mass transfer may be formed on the rotor's surface: "dry," "wet" and "frost." As one can see, "dry" area dominates at higher values of t 1i ; however, as the outdoor air temperature drops, share of this area's surface decreases at the expense of the rest of the zones. An opposite trend is observed when "frost" area is analyzed-the colder air is supplied to the heat exchanger, the larger surface is covered by frost. Moreover, within a certain range of sub-zero outdoor air temperature t 1i = (−20 … − 10.7) • C , the increase in the wheel's heat recovery efficiency is noted, what can be explained as the effect of released heat of phase transition. It is also worth noting that "frost" area is additionally divided into two subzones, which distinguish the area with frost accumulation from the other one, where frost is molten as a result of the contact with the warmer air stream. It allows to indicate the important range of temperature t 1i = (−10.7 … − 5) • C , in which the whole formed frost melts down. For the outdoor  ε(t) t 1i , °C Accumulated "frost" area Molten "frost" area "Wet" area "Dry" area Heat exchanger effectiveness air temperature t 1i < −10.7 • C , the process of frost accumulation begins. At the same time, molten "frost" area gradually decreases. The single most striking issue to emerge from the analysis of Fig. 5 is the variability of the temperature effectiveness of heat recovery depending on the values of the outdoor air temperature. For the conditions of "dry" heat exchange t 1i = (3.3 … 5) • C , the lack of significant changes in the effectiveness of the exchanger seems obvious. On the other hand, with the decrease in outdoor air temperature t 1i < 3.3 • C , heat of phase transition is released, what can be observed in a form of the "wet" area. Despite the favorable phenomenon of condensation, the temperature effectiveness of the thermal wheel remains constant. The explanation of this result is the mass balance on the both sides of rotor M 1 = M 2 (Fig. 6). This condition stays unchanged up to the threshold temperature t 1i = −10.7 • C . On the left-hand side of this value, the process of frost accumulation begins due to the mass disbalance M 1 < M 2 . Therefore, the increase in the temperature effectiveness of the heat exchanger (Fig. 5) should be attributed to the excess mass rate of condensed water-vapor which occurs on the return air side M 2 (Fig. 6). This effect is extremely noticeable within the range of outdoor air temperature t 1i = (−20 … − 10.7) • C.
The exact variability in the effectiveness of rotary heat exchanger may be observed by analyzing both mass rates of condensed and evaporated water inside two airflow sections (Fig. 6). On the basis of this graph, Only within the third range of outdoor air temperature, the latent heat released during the process of condensation of water-vapor is partially taken during the evaporation of water in the adjacent channel and the rest of this heat is used for additional heating of the wheel matrix. As a result, the supply air stream temperature rises what leads directly to the increase in the temperature effectiveness of heat recovery.
On the basis of the presented charts (Figs. 5,6,7,8,9,10,11,12,13,14), one can notice that the function of the temperature effectiveness is shaped in a similar way for all presented variants of the rotor's depth, i.e., its value rises due to moisture accumulation (in the form of water film or frost layer). Another similarity is an increase in "dry" area with increasing outdoor air temperature. Moreover, the influence of the rotor's size on the heat and mass transfer process is also observed. The most noticeable relationship is an increase in the temperature effectiveness of heat recovery as the geometrical parameters of the wheel rise. The next conclusion is concerned with the formation of different heat and mass transfer zones. On the one hand, larger depth of the rotor causes an increase in "wet" area on the surface of the storage mass. On the other hand, for the lowest values of outdoor air temperature within the range of t 1i = (−20 … − 10) • C the total area covered with both "wet" and "frost" zones is greater for the smaller depth of the rotary heat exchanger. The influence of the wheel's depth is also visible analyzing the values of threshold outdoor air temperature, at which frost accumulation begins inside the channels. It turns out that there is a clear trend of increasing frosting limit temperature by raising the rotor's depth. In case of outdoor air temperature t 1i = −10 • C , while operating the wheel with a depth of 0.20 m, there is no accumulated "frost" area. For the same outdoor air temperature, the largest size of the rotor causes 10% of accumulated "frost" zone. Another important conclusion is that increasing depth of the wheel ensures the release of much larger latent heat.
The main limitation of operation of these types of regenerators is a possibility of frost building up on the inlet supply air side of the heat exchanger (Fig. 2). Unfortunately, such issues can only be analyzed by means of numerical simulations with the use of programs developed on the basis of mathematical models. Effective testing of rotary heat exchanger operation in freezing conditions is very dangerous because it can cause permanent damage or even destruction of unit. For this reason, future directions of research will be carried out in order to determine the inlet air parameters for safe operation of the heat exchanger. In the next stages, the possibilities of reducing the effectiveness of the device will be analyzed, which, as could be seen in the graphs (Figs. 5,6,7,8,9,10,11,12,13,14), allow for the elimination of frost accumulation zones, causing its melting.

Conclusions
In this paper, the influence of geometrical parameters of a rotary heat exchanger on heat and mass transfer processes as well as its temperature effectiveness is investigated.
For this purpose, the mathematical model of a heat recovery wheel based on a structure of the counter-flow heat exchanger model was developed. It can be used successfully when high-speed rotor conditions are considered. The next step was to implement a numerical method that could be used to solve developed heat and mass transfer balance equations, and afterward the computer simulations were carried out. As a result, the percentage distribution of different heat and mass transfer zones inside the wheel as a function of outdoor air temperature value for various rotor depths was presented. Furthermore, the analysis of the water-vapor mass transfer between supply and return air sides of the heat exchanger was conducted. It constitutes a fundamental Presented analysis allowed to receive the following conclusions: • It was indicated that moisture accumulation has a significant influence on the increase in temperature effectiveness of heat recovery. This property is a result of the released heat of phase transition, which in the cases under consideration means an increase in efficiency of 1.4%. • Rotary heat exchangers equipped with a larger rotor are characterized by a higher limit temperature, which determines the occurrence of condensation (in the case of rotor depth equal to 0.20 m → t 1i < 3.3 °C but for the rotor depth equal 0.40 m → t 1i < 4.5 °C). • The increase in rotor's depth leads to the increase in temperature effectiveness of the rotary heat exchanger (from 0.667 to 0.814) because of larger heat and mass transfer area. This effect is observed under both "dry" and "wet" heat transfer conditions. The geometrical parameters of the wheel have also the influence on the frost accumulation threshold temperature-the value of temperature increases with increasing of the rotor's depth (from -10.7 to -9.6 °C) (a deeper rotor causes a greater cooling of the return air stream, which in turn causes the surface temperature of the wheel matrix to fall below 0 °C and the dew point temperature to drop faster). • It was also shown that the total area of both "frost" and "wet" zones under conditions of the extremely low values of the outdoor air temperature t 1i = (−20 … − 10) • C is in the range (66…55)% and is smaller than for a smaller depth of rotor (71…56)%.
The received results will be the basis for further optimization studies, and they will allow evaluating the possibility of using a rotary heat exchanger depending on climatic conditions. Fig. 14 Two water-vapor mass transfer rates inside outdoor and return air heat exchanger sides expressed as a function of outdoor air temperature values (for return air parameters: t 2i = 20 • C , RH 2i = 40% and rotor depth 0.40 m)