energies
Article Improvement of a Nusselt-Based Simulation Model for Heat Transfer in Rotary Heat Exchangers
Eloy Melian 1,* , Harald Klein 2 and Nikolaus Thißen 1
1 Institute of Industrial Ecology, Pforzheim University of Applied Sciences, 75175 Pforzheim, Germany; [email protected] 2 Department of Mechanical Engineering, Technical University of Munich, 85748 München, Germany; [email protected] * Correspondence: [email protected]; Tel.: +49-72-3128-6407
Abstract: In the last 50 years, the technology of rotary heat exchangers has not changed considerably. A reliable simulation can help improve the design of this technology. In this work, a simulation for rotary heat exchangers was developed and validated with multiple experimental data. This simulation takes an innovative approach based on locally calculated heat transfer coefficients and considers the entry region effect. This approach proved to be accurate since the average difference between the experimental results and the proposed model with a constant heat boundary condition is 0.1% and the maximum absolute deviation 1%. Experimental, as well as simulation results, indicate that lower empty tube gas velocity (1 m/s) and higher rotational speed (12 rpm) improve thermal efficiency compared to commonly used operating conditions. Additionally, a new model for predicting the local internal Nusselt number for sine ducts in the rotor channels is proposed, which considers the entry region effect.
Keywords: rotary heat exchanger; heat recovery wheel; thermal wheel; local Nusselt number; thermal efficiency; temperature effectiveness; heat recovery; heat recovery system; regenerator; heating, ventilation and air conditioning (HVAC)
Citation: Melian, E.; Klein, H.; Thißen, N. Improvement of a Nusselt-Based Simulation Model for Heat Transfer 1. Introduction in Rotary Heat Exchangers. Energies 2021, 14, 10. https://dx.doi.org/10.3390/ Rotary heat exchangers are already among the most efficient [1,2] and compact [3] en14010010 technologies for air-to-air heat recovery. Over the last decade, heat transfer efficiency requirements in the European Union (EU) for air-conditioning systems have been steadily Received: 8 November 2020 increasing. In the EU after 2020 it will be legally required in non-residential buildings that Accepted: 17 December 2020 the thermal efficiency must be greater than 85% [4]. The thermal efficiency η refers to the . Published: 22 December 2020 relationship between the actual heat flow transferred by the heat exchanger Qtrans f erred and the maximal heat flow that is theoretically possible in an adiabatic and infinite long heat Publisher’s Note: MDPI stays neu- . exchanger Q as defined by different authors [5,6] tral with regard to jurisdictional claims max in published maps and institutional . affiliations. Qtrans f erred η . (1) ≡ Qmax This definition of efficiency does not include electricity and mechanical power losses Copyright: © 2020 by the authors. Li- due to ventilators in the system. These loses are out of the scope of this work and a possible censee MDPI, Basel, Switzerland. This point for research in further work. article is an open access article dis- The relevance of these legislative efficiency measures relies on the amount of installed tributed under the terms and condi- devices. According to a study by Kaup and Kampeis [7], in Germany alone more than tions of the Creative Commons Attri- 60,000 ventilation systems were delivered in 2009 and approximately 60% of them had a bution (CC BY) license heat recovery unit. By 2012, deliveries had doubled and more than 80% of the systems (https://creativecommons.org/licenses/ by/4.0/). were equipped with a heat recovery unit. This development is due to higher energy
Energies 2021, 14, 10. https://dx.doi.org/10.3390/en14010010 https://www.mdpi.com/journal/energies Energies 2021, 14, 10 2 of 26
costs, since ventilation represents between 35% and 38% of energy loss in buildings [7]. These numbers do not consider the demand for industrial purposes and yet show how significant a thermal efficiency improvement on heat recovery systems can influence energy demand. Among the different types of heat recovery system, rotary heat exchangers can provide thermal efficiency greater than 80%, making them the most efficient technology known to date [8,9]. This technology was originally developed and patented in 1922 by the Swedish engineer Frederik Ljungström [10]. In the following decades, thermodynamic, aerodynamic and structural engineering principles were examined to further develop rotary heat exchangers [11]. Rotary heat exchanger innovations stagnated during the 1960s. Since then, for example the rotor length in air ventilation systems has been a fixed 20 cm as a standard design parameter. If a change in these design or operational parameters could increase the thermal efficiency by 5–10%, this would already meet the European target [5], making it possible for rotary heat exchangers to comply with the regulations in the future. This would open the possibility of saving more energy, and decrease costs and CO2 emissions with regard to the current rotary heat exchangers. The purpose of this work is to develop and validate a simulation of rotary heat exchangers that can predict their thermal efficiency by a given set of operational and design parameters. With this simulation, a better design of rotary heat exchangers can be achieved in the future. Previous work has tried to fulfil this task; the accuracy of these state-of-the-art simulations was in the best cases limited (maximum absolute deviation of 3.5% points for Leong [12] or 3.71% in average for Özdemir [13]). To achieve the goal of this work, a high-precision simulation of the heat transfer for rotary heat exchangers, with the following innovative strategies are implemented in this work: (1) To improve the accuracy of the simulation results of rotary heat exchangers by using a locally calculated Nusselt number model in order to determine the heat transfer coefficients between the gas streams and the solid material. To date, many studies have used constant heat transfer coefficients along the tube, where the heat transfer takes place [13–15]. (2) Since there is no literature on local Nusselt numbers for the developing region in sine ducts, this effect is usually neglected in other works. In this work, this knowledge gap is addressed. (3) For validating the developed simulation model, multiple experimental data will be obtained from a pilot plant and from a publicly available online database from Eurovent [16]. This will ensure that the developed simulation model is validated with data from different rotary heat exchangers. Previous simulations were only validated with one rotary heat exchanger [2,12,15,17–22] or a single rotor from a literature source [13,23–25]. (4) Previous authors have used constant thermophysical properties [12,13,15,19,20,23,26]. Here, for improving simulation accuracy, thermophysical properties are calculated locally and are temperature dependent. (5) The state-of-the-art simulations for a rotary heat exchanger use either the ε-NTU method [2,12,14,15,17,20–24,27–32] or computer fluid dynamics [13,18,25,26]. Finite difference methods, with locally calculated Nusselt number models and entry-region con- siderations, as done in this work, have until now been non-existent in literature. In the current context, the development of a more accurate simulation model for rotary heat exchangers is highly relevant because on the one hand, it predicts the extent of the impact of each parameter on the operation parameters, in particular the heat recovery efficiency of the system. On the other hand, such a model can serve as a basis for optimization. Without this accuracy, a further optimization would be futile. This is of highest significance because it can satisfy both sustainable manufacturing of this sys- tem and compliance with regulations regarding energy recovery. These requirements have become more and more strict over the last decades and probably will do so in future too. Energies 2021, 14, 10 3 of 26
Energies 2020, 13, x FOR PEER REVIEW 3 of 26
2.2. Materials Materials and and Methods Methods ForFor practical practical reasons, reasons, rotary rotary heat heat exchangers, exchangers, which which are are also also known known in the in literaturethe literature as “thermalas “thermal wheels”, wheels”, “heat “heat recovery recovery wheels”, wheels”, “Kyoto “Kyoto wheels” wheels” or simply or simply “rotors” “rotors” will in will most in casesmost hereafter cases hereafter be referred be referred to as “rotor” to as “rotor” or “thermal or “thermal wheel”. wheel”. ThermalThermal wheels wheels are are placed placed inside inside housing housing that th isat dividedis divided into into two two sectors, sectors, as shown as shown in Figurein Figure1. The 1. thermalThe thermal wheel wheel rotates rotates continuously continuo fromusly one from sector one tosector the other. to the In other. one sector, In one hotsector, process hot airprocess (1.1) flowsair (1.1) into flows thewheel, into the the wheel, wheel the metal wheel material metal is material in contact is in with contact the hotwith process the hot air, process which transfersair, which heat transfers to the wheel,heat to thethewheel wheel, heats the wheel up and heats stores up the and energy stores providedthe energy by theprovided process by air. the The process process air. air The flows process continuously air flows through continuously the wheel through leaving the thewheel wheel leaving at a lower the wheel temperature at a lower (exhaust temperatur air 1.2).e Since(exhaust the wheelair 1.2). is Since continuously the wheel rotating, is con- thetinuously heated partrotating, of the the wheel heated rotates part intoof the the wheel other rotates sector whereinto the it other meets sector the colder where ambient it meets airthe (2.1). colder The ambient heated solidair (2.1). material The heated transfers solid the material stored heat transfers to the ambientthe stored air. heat The to ambient the am- airbient leaves air. the The wheel ambient at a air higher leaves temperature the wheel (supplyat a higher air 2.2).temperature The now (supply cold wheel air 2.2). section The continues rotating to the first sector where it again heats up. Since the thermal wheel is now cold wheel section continues rotating to the first sector where it again heats up. Since rotating continuously from one sector into the other, one half of the wheel is continuously the thermal wheel is rotating continuously from one sector into the other, one half of the being heated and storing energy from the process air, while the other half is continuously wheel is continuously being heated and storing energy from the process air, while the cooling and returning the stored energy to the ambient air. other half is continuously cooling and returning the stored energy to the ambient air.
FigureFigure 1. 1.Rotor Rotor working working principle principle and and stream stream names. names. . Q trans f erred is the is the energy energy that that the the rotor rotor provides provides to to the the supply supply air, air, which which is is based based on on thethe enthalpy enthalpy change change from from the the ambient ambient air air (2.1) (2.1) to to the the supply supply air air flow flow (2.2). (2.2). Its Its value value is is identicalidentical to to the the energy energy decrease decrease of of the the process process air air (1.1) (1.1) to to the the exhaust exhaust air air (1.2). (1.2). If If the the heat heat exchangerexchanger is is considered considered without without heat heat losses, losses, the the energy energy balance balance for for the the fluid fluid streams streams can can bebe expressed expressed as: as: . . . Qtrans f erred = m1.1 (h1.1 h1.2) = m2.1 (h2.2 h2.1) (2) = . · ∙ (ℎ− . −ℎ . ) = · . ∙ (−ℎ . −ℎ . ) (2) . Q is defined as the energy that would be transferred if the heat exchanger were max is defined as the energy. that would be transferred if the heat exchanger were adiabaticadiabatic and and infinitely infinitely long. long.Q max is limitedis limited by by the the stream stream with with the the lower lower mass mass flow: flow: . . ( . ∙ (ℎ . −ℎ . ) =min m2.1 (h1.1 h2.1) (3) Q = min . · ∙ (ℎ− −ℎ ) (3) max m . (h . h ) . 1.1· 1.1 − 2.1 Since both mass flows in this paper are assumed to be equal on each side and in nar- rowSince temperature both mass ranges, flows inthe this specific paper heat are assumedcapacity tocan be be equal considered on each sideconstant. and in The narrow equa- temperaturetion is commonly ranges, presented the specific as: heat capacity can be considered constant. The equation is commonly presented as: h ℎ . h−ℎ . ϑ . ϑ− . η = =2.2 − 2.1 = =2.2 − 2.1 (4)(4) h ℎ . h−ℎ . ϑ . ϑ− . 1.1 − 2.1 1.1 − 2.1 TheThe right right side side of Equationof Equation (4) is(4) what is what is often is often called called “temperature “temperature effectiveness” effectiveness” by Shah by andShah Sekuli´c[ and Sekuli6] orć “temperature [6] or “temperature efficiency” efficiency” by Eurovent by Eurovent [16] and for[16] comparison and for comparison reasons thisreasons will bethis the will most be importantthe most important parameter para to bemeter further to be used further in this used work. in this Hereafter work. Here- it is onlyafter referred it is only to referred as “temperature to as “temperature effectiveness” effectiveness”η. .
Energies 2021, 14, 10 4 of 26
Energies 2020, 13, x FOR PEER REVIEW 4 of 26
2.1.2.1. Research Research Tools Tools TheThe methodological methodological approach approach in this in paperthis pa usesper simulations uses simulations based onbased thermodynamic on thermody- equationsnamic equations as stated as in stated Section in 2.1.1 Section and 2.1.1 a pilot and plant a pilot set-up plant as set-up depicted as depicted in Section in 2.1.2Section. These2.1.2. are These be the are main be the tools main used tools for used obtaining for obtaining and analyzing and analyzing data afterwards. data afterwards.
2.1.1.2.1.1. Simulation Simulation InIn this this section, section, the the simulation simulation tool tool and and methods methods are are explained: explained: this this comprises comprises the the energyenergy balance balance (Section (Section 2.1.1.1 2.1.1.1),), assumptions assumptions and and boundary boundary conditions conditions (Section (Section 2.1.1.2 2.1.1.2),), thethe simulation simulation concept concept and and numeric, numeric, as as well we asll theas the thermodynamic thermodynamic and and fluid fluid dynamic dynamic considerations.considerations.
2.1.1.1.2.1.1.1. Energy Energy Balance Balance TheThe base base model model for for the the heat heat transfer transfer of rotorsof rotors has has already already been been explained explained in detailin detail in in existingexisting literature literature [30 [30,33].,33]. It isIt basedis based on on the the energy energy balances balances between between a metal a metal channel channel and and thethe gas gas that that flows flows through through it. Thisit. This model model assumes: assumes: • One-dimensionalOne-dimensional heat heat conductivity conductivity in in axialaxial directiondirection z for for thethe gas and and solid solid phase phase (see • (seeFigure Figure 22 for for the the zdirection).direction). • TheThe solid solid temperature temperature is homogeneous is homogeneous in the in the radial radial direction directionr as depicted as depictedon Figure on Figure2 • (infinite2 (infinite radial radial thermal thermal conductivity). conductivity). Please Please note note that that the the radial radial direction direction is perpen-is perpen- diculardicular to theto thez coordinate. coordinate. • AtAt each each given given cross-section cross-sectio insiden inside the the duct duct along alongz the gas the isgas perfectly is perfectly mixed mixed and hasanda has • homogeneousa homogeneous temperature temperature distribution. distribution. • TheThe solid solid phase phase is homogeneousis homogeneous in in thermophysical thermophysical properties properties at at the the start. start. • • No back mixing in the gas flow takes place. • No back mixing in the gas flow takes place. • The system is considered to have no heat losses. • The system is considered to have no heat losses. • Potential and kinetic energy effects are negligible. • Potential and kinetic energy effects are negligible. • Only one-dimensional gas flow in the direction z is considered. • Only one-dimensional gas flow in the direction is considered. • The gas is incompressible and at constant pressure. • The gas is incompressible and at constant pressure. • Enthalpy is only a function of temperature. • Enthalpy is only a function of temperature. . . • No mass transfer takes place between gas and solid, (Please note that if m1.1 = m2.1, • No mass transfer takes place between gas and solid, (Please note that if . = . , thenthen thermal thermal effi efficiencyciency and and temperature temperature effectiveness effectiveness are are equal). equal). • No adsorption, absorption or condensation is considered. • No adsorption, absorption or condensation is considered.
: Rotor length
Internal flow channels
coordinate
Rotor foil layers coordinate
Figure 2. z and r coordinates, empty tube gas velocity va and internal gas velocity vin and rotor length l (transverse view). Figure 2. and coordinates, empty tube gas velocity and internal gas velocity and rotor length (transverse view). The energy balances result in the following partial parabolic differential equation system which can be found easily in literature [34]: The energy balances result in the following partial parabolic2 differential equation dϑG dϑG d ϑG Gas phase : ε ρsystemG cp, G which= canα beav found(ϑS easilyϑG) inε literatureρG vin c p[34]:, G + ε λG (5) · · · dt · · − − · · · dz · · dz2 Gas phase: ⋅ ⋅ ⋅ = ⋅ ⋅ ( − ) − ⋅ ⋅ ⋅ + ⋅ ⋅ (5) , ,
Energies 2020, 13, x FOR PEER REVIEW 5 of 26 Energies 2021, 14, 10 5 of 26