Scientific Approaches to Solving the Problem of Joint Processes of Bubble Boiling of Refrigerant and Its Movement in a Heat Pump Heat Exchanger
Abstract
:1. Introduction
2. Materials and Methods
2.1. Accepted Assumptions
2.1.1. Vapor-Liquid Mixture
2.1.2. Boiling
2.1.3. Flow Mode
2.2. Theoretical Calculation
2.2.1. The Calculation with Assumptions Is Divided into Two Stages
2.2.2. Comparison of Data from Theoretical Calculation and Analytical Calculation for Temperature
2.3. Differential Equations When Solving by Numerical Method
2.3.1. Energy Equations
2.3.2. The Transfer Equation
2.4. Modeling in a Software Environment ANSYS Fluent
- From the experience of modeling the emulsion flow, the condition of a two-phase flow at the inlet is set (previously, only one liquid phase was set at the inlet, and not a mixture);
- The extreme pressures in the cycle were determined (3.5 and 19 bar, respectively);
- The saturation temperature is set from the theoretical calculation;
- Based on points 2 and 3, a certain degree of dryness of the vapor-liquid mixture at the inlet to the evaporator is set.
- −
- Element size a: 133.16 × 133.16 × 3 mm, diagonal: 205 mm;
- −
- Tube inner diameter: 3 mm, tube outer diameter: 5 mm;
- −
- Temperature on the front surface of the element: 60 °C;
- −
- Ambient air temperature: 30 °C.
3. Results
Comparison of the Received Data
4. Discussion
4.1. Limitations of the Study
4.2. Prospects for the Use of a Heat Pump in an Energy Complex for Seawater Desalination
5. Conclusions
Author Contributions
Funding
Data Availability Statement
Conflicts of Interest
Nomenclature
i1 | enthalpy of condensate, kJ/kg; |
i2 | enthalpy of vaporization, kJ/kg; |
qboil. | boiling heat flux density, kW/m2; |
λ | thermal conductivity of liquid Freon, kW/(m·K); |
ΔTs | temperature difference between wall and liquid, K; |
hg | enthalpy of gas, kJ/kg; |
hl | enthalpy of liquid, kJ/kg; |
B | correction factor of free convection heat transfer coefficient, non-dimensional; |
ν | kinematic viscosity, m2/s; |
σ | surface tension coefficient, N/m; |
Ri | gas constant for freon vapor, kJ/(kg·K); |
ρg | gas density, kg/m3; |
qconv. | heat flux density at convective heat transfer kW/m2; |
Nu | Nusselt number, non-dimensional; |
Re | Reynolds number, non-dimensional; |
Pr | Prandtl number, non-dimensional; |
ξ | correction factor that takes into account the features of the flow in a given system, non-dimensional; |
tev. wall | evaporator wall temperature, °C; |
tPEC wall | PVP surface temperature, °C; |
q | heat flow through PVP, kW/m2; |
R | thermal resistance of layers (protective glass, silicon wafers, etc.), (m2·K)/kW; |
α | heat transfer coefficient, kW/(m2·K); |
Ρ″ | saturated steam density, kg/m3; |
ρ | density of a saturated liquid, kg/m3; |
Ts | saturation temperature, K; |
p | pressure, Pa; |
μ | dynamic viscosity, Pa·s; |
ω | flow velocity, m/s; |
d | tube inner diameter, m; |
ε | correction factor taking into account the heating of the liquid from the wall, non-dimensional; |
αconv. | convective heat transfer coefficient, kW/(m2·K); |
αboil. | heat transfer coefficient for freon boiling, kW/(m2·K); |
λwall | thermal conductivity of the wall, kW/(m·K); |
δ | wall thickness, m; |
Qboil. | heat of saturated steam, kW; |
x | degree of dryness, non-dimensional; |
r | heat of vaporization, kJ/kg; |
lev. | evaporation section length, m; |
ls.h | refrigerant superheater section, m; |
lt | total length of the PVP section, m; |
tf.out | freon outlet temperature, °C; |
cp″ | freon steam heat capacity, kJ/(kg·K); |
partial differential of unstable heat transfer, kJ; | |
a term that takes into account thermal conductivity, kJ; | |
a term that takes into account the diffusion of mixtures, kJ; | |
sum taking into account viscous dissipation, kJ; | |
a term that takes into account the enthalpy of inflows/drains, kJ; | |
current velocity, m/s; | |
ρ0 | liquid density at some equilibrium temperature T0, kg/m3; |
Θ | deviation of temperature from equilibrium state (Θ = T − T0), K; |
χ | coefficient of thermal diffusivity, m2/s; |
acceleration of gravity, m/s2; | |
β | volume expansion coefficient, K−1. |
References
- Romsy, T.; Zacha, P. CFD simulation of upward subcooled boiling flow of freon R12. Acta Polytech. CTU Proc. 2016, 4, 73–79. [Google Scholar] [CrossRef]
- Milman, O.O.; Ananyev, P.A.; Korlyakova, M.O.; Miloserdov, V.O. Experimental studies of non-stationary thermo-hydraulic processes at freon R113 boiling. J. Phys. Conf. Ser. 2019, 1382, 012114. [Google Scholar] [CrossRef]
- Kuzmin, A.Y.; Bukin, A.V. Experimental study of heat transfer during boiling on a smooth tube under conditions of free convection of alternative refrigerants R407c and R410a. South Russ. Ecol. Dev. 2010, 5, 121–124. [Google Scholar] [CrossRef]
- Aleksandrov, A.A.; Orlov, K.A.; Ochkov, V.F. Thermophysical Properties of Working Substances of Heat Power Industry: Internet Reference Book; MPEI Publishing House: Moscow, Russia, 2009. [Google Scholar]
- Osintsev, K.V.; Alyukov, S.V. Experimental Investigation into the Exergy Loss of a Ground Heat Pump and its Optimization Based on Approximation of Piecewise Linear Functions. J. Eng. Phys. Thermophys. 2022, 95, 9–19. [Google Scholar] [CrossRef]
- Demir, M.E.; Dincer, I. Development and analysis of a new integrated solar energy system with thermal storage for fresh water and power production. Int. J. Energy Res. 2017, 42, 2864–2874. [Google Scholar] [CrossRef]
- Demir, M.E.; Dincer, I. Development of an integrated hybrid solar thermal power system with thermoelectric generator for desalination and power production. Desalination 2017, 404, 59–71. [Google Scholar] [CrossRef]
- Saidur, R.; Elcevvadi, E.; Mekhilef, S.; Safari, A.; Mohammed, H. An overview of different distillation methods for small scale applications. Renew. Sustain. Energy Rev. 2011, 15, 4756–4764. [Google Scholar] [CrossRef]
- Nafey, A.S.; Mohamad, M.; El-Helaby, S.; Sharaf, M. Theoretical and experimental study of a small unit for solar desalination using flashing process. Energy Convers. Manag. 2007, 48, 528–538. [Google Scholar] [CrossRef]
- Dincer, I.; Rosen, M.A. Chapter 2—Exergy and Energy Analyses, 2nd, ed.; Dincer, I., Rosen, M.A., Eds.; Elsevier: Amsterdam, The Netherlands, 2013; pp. 21–30. [Google Scholar] [CrossRef]
- Nasrabadi, A.M.; Korpeh, M. Techno-economic analysis and optimization of a proposed solar-wind-driven multigeneration system; case study of Iran. Int. J. Hydrogen Energy 2023, 48, 13343–13361. [Google Scholar] [CrossRef]
- Yapicioglu, A.; Dincer, I. A newly developed renewable energy driven multigeneration system with hot silica sand storage for power, hydrogen, freshwater and cooling production. Sustain. Energy Technol. Assess. 2023, 55, 102938. [Google Scholar] [CrossRef]
- Khrustalev, B.M. Heat and Mass Transfer; BNTU: Minsk, Belarus, 2007; 606p. [Google Scholar]
- Zalepugin, D.Y.; Tilkunova, N.A.; Chernyshova, I.V.; Vlasov, M.I. Application of Sub- and Supercritical Freons in Xenogenic Bone Matrix Processing. Russ. J. Phys. Chem. B 2017, 11, 1051–1055. [Google Scholar] [CrossRef]
- Fang, Y.; Wu, M.; Guo, Z.; Ni, H.; Han, X.; Chen, G. Evaluation on Cycle Performance of R161 as a Drop-in Replacement for R407C in Small-Scale Air Conditioning Systems. J. Therm. Sci. 2022, 31, 2068–2076. [Google Scholar] [CrossRef]
- Deb, S.; Mahesh, K.P.; Das, M.; Das, D.C.; Pal, S.; Das, R.; Das, A.K. Flow boiling heat transfer characteristics over horizontal smooth and microfin tubes: An empirical investigation utilizing R407c. Int. J. Therm. Sci. 2023, 188, 108239. [Google Scholar] [CrossRef]
- Alabugin, A.; Osintsev, K.; Aliukov, S.; Almetova, Z.; Bolkov, Y. Mathematical Foundations for Modeling a Zero-Carbon Electric Power System in Terms of Sustainability. Mathematics 2023, 11, 2180. [Google Scholar] [CrossRef]
Title | Designation | Calculation Formula or Method of Determination | Units of Measurement | Result |
---|---|---|---|---|
Heat flow | Qsol | Given | W/m2 | 700 |
Thermal conductivity of silicon | λsil | Given | W/(m2·K) | 148 |
Absorption capacity of silicon | Asil | Given | - | 0.2 |
Degree of silicon blackness | εsil | Given | - | 0.71 |
Silicon thickness | δsil | Given | m | 0.001 |
Stefan-Boltzmann coefficient | σ | Given | W·m−2·K−4 | 5.67 × 10−8 |
Thermal conductivity of glass | λgl | Given | W/(m2·K) | 0.96 |
Absorption capacity of glass | Agl | Given | - | 0.1 |
The degree of blackness of the glass | εgl | Given | - | 0.93 |
Glass thickness | δgl | Given | m | 0.002 |
Length and width of the PVP | l | Given | m | 0.15675 |
PVP Square | F | m2 | 0.025 | |
Reflectivity of silicon | Rsil | Given | - | 0.33 |
Heat absorbed by the PVP | Q | W | 11.524 | |
Heat absorbed by silicon | Qabs.sil | W | 3.44 | |
Heat passed through silicon | Qpast sil. | W | 8.084 | |
Heat absorbed by glass | Qabs.gl. | W | 1.72 | |
Heat passed through the glass | Qpast gl. | W | 6.364 | |
Air velocity | ω | Given | m2/s | 1 |
Thermal conductivity of air | λ | Given | W/(m2·K) | 2.4 × 10−2 |
Air viscosity | νair | Given | m2/s | 14.61 × 10−6 |
Prandtl number for air | Pr | Given | - | 0.67 |
Ambient air temperature | tair | Selected | °C | 20 |
Ambient air temperature | Tair | K | 293 | |
Outside temperature of the PVP | tsur1 | Selected | °C | 29.12 |
Outside temperature of the PVP | Tsur1 | K | 302.12 | |
Reynolds number | Re1, Re2 | - | 1.073 × 104 | |
Nusselt number | Nu1, Nu2 | - | 60.727 | |
Heat transfer coefficient | α1, α2 | W/(m2·K) | 9.298 | |
Temperature of the inner surface of the PVP | tsur2 | °C | 29.1 | |
Losses from above | qabove | W | 2.083 | |
Silicon radiation losses | qsil.rad | W | 0.951 | |
Losses from below | qund. | W | 2.083 | |
Glass radiation losses | qgl.rad | W | 6.399 | |
Heat losses in the system | Qlosses | W | 11.516 | |
Error rate | ΔQ | % | 0.067 |
Wind Speed, m/s | 1 | 5 | 10 |
---|---|---|---|
The temperature of the outer wall according to the theoretical calculation, °C | 40.74 | 37.20 | 34.30 |
External wall temperature according to analytical calculation, °C | 40.57 | 37.06 | 34.21 |
Ф | P | D | Γϕ | |
---|---|---|---|---|
Kinetic energy | k | or | ||
Kinetic energy dissipation rate | ε | |||
Specific rate of dissipation | ω |
Wind Speed, m/s | 1 | 5 | 10 |
---|---|---|---|
Numerical modeling calculation, °C | 40.23 | 36.94 | 34.10 |
External wall temperature according to analytical calculation, °C | 40.57 | 37.06 | 34.21 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Osintsev, K.; Aliukov, S.; Kovalev, A.; Bolkov, Y.; Kuskarbekova, S.; Olinichenko, A. Scientific Approaches to Solving the Problem of Joint Processes of Bubble Boiling of Refrigerant and Its Movement in a Heat Pump Heat Exchanger. Energies 2023, 16, 4405. https://doi.org/10.3390/en16114405
Osintsev K, Aliukov S, Kovalev A, Bolkov Y, Kuskarbekova S, Olinichenko A. Scientific Approaches to Solving the Problem of Joint Processes of Bubble Boiling of Refrigerant and Its Movement in a Heat Pump Heat Exchanger. Energies. 2023; 16(11):4405. https://doi.org/10.3390/en16114405
Chicago/Turabian StyleOsintsev, Konstantin, Sergei Aliukov, Anton Kovalev, Yaroslav Bolkov, Sulpan Kuskarbekova, and Alyona Olinichenko. 2023. "Scientific Approaches to Solving the Problem of Joint Processes of Bubble Boiling of Refrigerant and Its Movement in a Heat Pump Heat Exchanger" Energies 16, no. 11: 4405. https://doi.org/10.3390/en16114405
APA StyleOsintsev, K., Aliukov, S., Kovalev, A., Bolkov, Y., Kuskarbekova, S., & Olinichenko, A. (2023). Scientific Approaches to Solving the Problem of Joint Processes of Bubble Boiling of Refrigerant and Its Movement in a Heat Pump Heat Exchanger. Energies, 16(11), 4405. https://doi.org/10.3390/en16114405