Next Article in Journal
A Comprehensive Review of Deep-Learning Applications to Power Quality Analysis
Previous Article in Journal
Investigation of Thermo-Hydraulic Performances of Artificial Ribs Mounted in a Rectangular Duct
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Scientific Approaches to Solving the Problem of Joint Processes of Bubble Boiling of Refrigerant and Its Movement in a Heat Pump Heat Exchanger

Department of Energy and Power Engineering, Institute of Engineering and Technology, South Ural State University, 76 Prospekt Lenina, 454080 Chelyabinsk, Russia
*
Author to whom correspondence should be addressed.
Energies 2023, 16(11), 4405; https://doi.org/10.3390/en16114405
Submission received: 22 April 2023 / Revised: 14 May 2023 / Accepted: 25 May 2023 / Published: 30 May 2023
(This article belongs to the Topic Heat Exchanger Design and Heat Pump Efficiency)

Abstract

:
The joint processes of bubble boiling and refrigerant movement in heat pump tubes are considered. A coil located on the back side of the photovoltaic panel (PVP) is used as an additional heating surface. A mathematical model of the theoretical determination of the temperature on the front surface of the PVP cooled by freon in the coil from the back side has been created. A feature of the numerical simulation of bubble boiling of refrigerant in a coil was the insufficiently detailed study of the process and the lack of references to the results obtained in earlier work. The authors have clarified the boundary conditions and assumptions for numerical simulation of the bubble boiling process of refrigerant R407C. The results were compared with the theoretically found values, and later used in the design of an energy complex consisting of a heat pump and a photovoltaic panel. The mathematical model of theoretical calculation of the PVP temperature and the methodology for constructing a numerical model of bubble boiling of refrigerant as part of the methodology for designing energy technology complexes based on the authors’ plan to present a unified methodology for energy technology complexes for desalination of seawater based on other types of renewable energy sources.

1. Introduction

To date, there are no universal methods for calculating the heat transfer coefficient or a universal description of the thermophysical picture of boiling in a liquid stream. However, many empirical methods have been created for various refrigerants under conditions of low-pressure values (up to 2–4 atm) and mass flow rates, and special formulas for calculating heat transfer to small channels (with a diameter of 1 mm and above) and to strictly specified pressure values have been derived [1,2,3,4,5].
Refrigerant R407C was chosen as the working medium for the study. Its thermodynamic properties are close to the common R22, while it does not contain chlorine, and is non-flammable and non-toxic.
Consider energy complexes that combine heat pumps and the use of solar energy. In the scientific paper [6], the authors developed a system based on solar energy for the production of fresh water and electricity. This system consists of a solar module, a solar-powered Rankine cycle, a heat storage subsystem, and a distillation device. Sea water used for desalination is heated by a saturated steam-water mixture from a steam turbine. Using the outlet liquid as a heat source for the distillation device eliminates the need for an external heat exchanger for condensation. The authors of [6] calculated the energy and exergy efficiency, power and fresh water production capacities for each element of the system. In addition, a parametric study was carried out to take into account environmental conditions.
In research paper [7] the authors analyzed the energy complex for four modules, including a solar heliostat and a receiver, a heat storage device, a Rankine cycle and a desalination system. The steam generator transfers heat to the Rankine cycle. The authors of the study [8] present the balance equations for each element separately and for the entire system as a whole. A mathematical model for the utilization of waste heat from the Rankine cycle has been developed. However, the authors make a number of assumptions: there are no chemical reactions, the potential and kinetic energy practically do not change, the energy costs for the pump are negligible, except for operation in the Rankine cycle, and heat losses occur only through the surface of the reservoir with molten salt and the solar collector.
In [9], the heat loss rate of a solar receiver was calculated. In the research paper [10], after calculating the temperature difference, the authors analyzed the heat losses through the tank with molten salt. The authors of [10] developed a mathematical model of heat accumulation in a tank. This is a rather complicated model, since the volume of the salt melt inside the tank is constantly changing. In [11], a mathematical model of the seawater desalination process was constructed. The study [12] considers the energy and exergy efficiency of the system as a whole. This can be achieved using the main inputs and outputs, while the useful outputs of the system are fresh water and electricity, and the only source of energy is solar energy.
The methodology considered by the authors is a scientific approach to solving the problems of calculating efficiency and optimization in the design of energy technology complexes that combine a heat pump (including a turbo expander for the implementation of the organic Rankine cycle), a photovoltaic panel, a desalination plant, and optionally a solar concentrator and a wind generator. In addition, the scientific approach does not deny the possibility of not only working on the Rankine cycle, since the methodological base does not change. A special feature of the study was the method of heating refrigerant—an additional heater in the form of a coil installed on the back of the photovoltaic panel, Figure 1a,b. Figure 1c also shows the basic scheme of operation of the energy technology complex implemented by the authors. The article considers the first part of the installation—a heat pump.
Let us explain the diagram in Figure 1c. The working substance is refrigerant R407c, which moves along the contour of the system. An evaporator works to convert a gaseous substance into a liquid at low temperature and pressure. The compressor then compresses the liquid refrigerant, increasing the pressure to 19 bar and the temperature to 85 °C. After passing through the compressor, the refrigerant enters the heat exchanger (condenser). Here it gives up its heat to the water, turning it back into a liquid form under pressure. Water enters the system through the filler neck and exits through the drain cock. After the condensed refrigerant passes through the expansion valve, the pressure drops significantly, causing some of the liquid to evaporate. This mixture of liquid and vapor goes to the evaporator. The liquid turns into vapor, evaporating under the action of the heat of the FEP or an additional heater. Due to this, boiling occurs in the evaporator, which consumes heat from the environment.

2. Materials and Methods

2.1. Accepted Assumptions

2.1.1. Vapor-Liquid Mixture

The condensate is throttled before being fed to the evaporator. At the same time, along with a decrease in pressure, the enthalpy constancy (1) can be set with sufficient accuracy.
i 1 = i 2
The throttling process in the P-i diagram for the selected refrigerant is shown in Figure 2 below in Process 3–4.
Given the necessary pressure drop and knowing the condensation and evaporation temperatures, it is possible to calculate the degree of dryness of steam at point 4, corresponding to the state of refrigerant at the inlet to the evaporator. The degree of dryness under the specified conditions is equal to x = 0.4. Thus, as a mandatory boundary condition, a fraction of 40% of the mass flow rate for steam is set.

2.1.2. Boiling

The calculation takes into account the flow rate in a pipe with a diameter of 5 mm. Developing the previous suggestion, a mix of bubble and foil boiling is assumed. When the liquid boils and evaporates, bubbles form on the walls, which quickly collapse and form a thin film of liquid when there is enough steam in the channel section.

2.1.3. Flow Mode

This flow mode corresponds to the form of a mixture in which steam forms the core of the flow, while the liquid phase flows in the shape of a foil on the surface of the wall. From the first assumption, it can be concluded that the vapor-liquid blend enters the evaporator in the emulsion mode, which is typical for the initial stage of vaporization. Due to the increase in steam consumption (up to a fraction of 0.4, as can be concluded from the diagram in Figure 2) bubbles tend to accumulate and stick together in the flow into larger ones.
It was assumed that the emulsion regime would prevail throughout the site until the formation of steam shells, that is, large steam plugs (Figure 3a). As a result of computer modeling, a similar picture was obtained in the initial formulation of the problem. Figure 3b shows the distribution of the proportion of liquid (red) and steam plugs (gradient to blue) in the tube section. Such a picture was achieved based on the considerations below.
Intensive evaporation of liquid along the boundaries of dry spots (centers of vaporization) is taken as the determining mechanism of bubble boiling. A calculation that takes this factor into account explains the extremely high rate of heat transfer during bubble boiling. When studying a physical model and generalizing experimental data, the analytical Formulas (2) and (3) are adopted for the calculation of heat transfer:
q b o i l . = 3.43 10 4 λ 2 Δ T s 3 ν σ T s 1 + h g h f Δ T s 2 2 R i T s 1 + 1 + 800 B + 400 B
B = h g h f ρ g ν 1.5 σ λ T s 0.5
The thermal transfer at boiling in the liquid flow is calculated as a combination of convective heat exchange and heat exchange at boiling with the predominance of the latter (4):
q = q c o n v . + q b o i l .
The final criterion similarity equation takes the form (5):
N u = R e P r ξ / 8 1 + 900 / R e + 12.7 ξ / 8 0.5 P r 0.67 1 .

2.2. Theoretical Calculation

2.2.1. The Calculation with Assumptions Is Divided into Two Stages

1. Determining the value of heat entering the system, taking into account losses to the environment;
2. Determination of the refrigerant temperature at the outlet of the evaporator. The first part is presented in the format of a flowchart in Figure 4 according to the Gauss methodology (with input of initial data, forward and reverse, and output of results). As a result, the temperature on the surface of the PVP is calculated with an error of less than 0.1%. The results with an ambient temperature of 20 °C and a wind of 1 m/s are summarized in Table 1.
A 3D diagram is constructed based on the calculation of the temperature of the PVP wall depending on two input arguments: wind speed and ambient temperature. The result is shown in Figure 5:
The result obtained is used to calculate the temperature on the surface of the evaporator tube (6):
t ev . wall = t PEC   wall   q i = 1 n R i
The temperature on the evaporator wall is a boundary condition for simulating the refrigerant boiling process. For further theoretical calculation, it is a test value when calculating the thermal transmission ratio.
Theoretically, the approximate value of the thermal transmission ratio is calculated according to the Labuntsov formula:
α = 0.075 1 + 10 ρ ρ ρ 2 / 3 λ 2 ν σ T s 1 / 3 q 2 / 3
According to the reference books, the saturation temperature is determined to be TH under pressure p at point 4 in the diagram according to Figure 2. The physical properties of saturated liquid refrigerant are searched for at the obtained temperature according to reference data [14]. The Reynolds number is calculated at a given flow velocity ω (8):
R e = ω d ρ μ
The Nusselt number is calculated according to the Petukhov formula (it implies a developed turbulent regime that will arise due to the presence of turns of the evaporator and high flow velocity of the refrigerant) (9):
N u = 0.021 R e 0.8 P r 0.43 ε
Next, the convective thermal transmission and boiling heat transfer are compared by comparing the corresponding coefficients. The calculation of the thermal transmission ratio for the forced flow process in the tube is determined as follows (10):
α c o n v . = N u λ d
The thermal transmission ratio for refrigerant boiling is computed according to the equation corresponding to the geometric and thermophysical similarity. Its approximate calculation is performed by Formula (7) under the condition of a boiling pressure above 1 atm.
IF αboil/αconv. > 3, then the thermal transmission ratio during boiling in the pipe is assumed to be equal to αboil, otherwise—αconv.
As a verification condition for the correctness of the selection of the criterion equation, the surface temperature of the evaporator tube is calculated to be higher. The calculation is correct if the temperature value in this part coincides with the previous one (11):
t e v . w a l l = T s 273.15 + q α e v . + q 2 λ w a l l d ln d d + 2 δ
To determine the refrigerant temperature at the evaporator outlet, the amount of heat to bring the steam-water mixture to the state of saturated steam is determined, as well as for its overheating (12):
Q b o i l . = ω ρ π d 2 4 r 1 x
Evaporation section length (13):
l e v . = Q b o i l . q π d
The remaining section is a refrigerant superheater (14):
l s . h = l t l e v .
The refrigerant outlet temperature is (15):
t f . o u t = 4 q l s . h ω ρ d c p + T s 273.15

2.2.2. Comparison of Data from Theoretical Calculation and Analytical Calculation for Temperature

Figure 4 and Formula (15) compared the obtained data. The discrepancy was about 0.1–0.2 °C. Comparison of Data from Theoretical Calculation and Analytical Calculation are presented in Table 2.

2.3. Differential Equations When Solving by Numerical Method

2.3.1. Energy Equations

To apply solutions by the numerical method, the formulas are used:
ρ u i t + x j ρ u i ρ u j + P x i = x τ i j + τ i j R + S i
ρ E t + V ρ E + ρ = k e f Δ T j h j J j + τ e f V + S h

2.3.2. The Transfer Equation

Generally used in models with one or two differential equations, the transfer equation can be presented in this form:
ρ Φ t + ρ u j ¯ Φ x j = P D + x j ( μ + Γ Φ ) Φ x j + A
The interpretation of the parameters for each type of equations is presented in Table 3.
The turbulence model is used in the calculations of transition k-kl-ω.

2.4. Modeling in a Software Environment ANSYS Fluent

Standard k-ε (+2E) is a two-parameter RANS model, which is very economical in terms of computational resources and is used for a wide class of flows. Wall functions are used and the model is suitable only for steady turbulent flows. The disadvantage is that in many problems there is a strong eddy diffusion. It should be noted that this model is able to take into account the effect of compressibility and natural convection, which is useful in the calculation of industrial cooling systems. However, this model is not suitable for solving the problem of refrigerant boiling, since it gives a large error and low convergence of results during boiling.
The authors came to the conclusion that transitional flow models should be used (the models only model the transition from laminar to turbulent flow; the wall function is used):
Transition k-kl-ω (+3E)
The k-kl-ω (+3E) transition turbulence model is a widely used turbulence model that combines the k-ε and k-ω models together with the Baldwin–Lomax model. This combination of models provides a reliable method for predicting flow behavior in transients (which is useful in calculating boiling due to its erratic nature and difficult prediction).
The k-kl-ω (+3E) transition model has a number of advantages over other turbulence models. It is able to accurately predict the onset of the transition from laminar to turbulent flow, which is a challenge for many turbulence models.
Transition k-kl-ω (+3E) is a three-parameter model, suitable for solving heat transfer problems in elements that create turbulence (for example, on a turbine blade).
However, the k-kl-ω (+3E) transition model is not without its drawbacks. The model is relatively complex and requires significant computational resources to solve. In addition, it is more prone to numerical instability than other simpler turbulence models, potentially making it difficult to use in some problems. Finally, the model relies on the use of a transition factor to another model, which can lead to inaccuracies and errors if the factor is set incorrectly.
Border conditions
  • 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.
Initial data
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.
Computer modeling of the flow process of refrigerant R407C in the channel of a brass evaporator with an external diameter of 10 mm and a wall width of 1 mm with a dispersed-ring flow is carried out. Figure 6a shows a geometric model.
The search for unknowns in the equations is carried out using structured and body-adapted prismatic grids (Figure 6b), which allows solving the problem when describing complicated geometry (sweep method for cylindrical bodies; at least 2 million prismatic and tetrahedral elements). A partition in the wall region (a layer of at least three cells) was also performed, which allows using the boundary layer theory for correct calculation of the flow at the walls of the tube, as well as correct calculation of the redistribution of the vapor and liquid phase.
The simulation includes the condition of fluid flow in the gravity field with the formation of vapor bubbles. The group of problems of the effect of the heating wall on the evaporated liquid is solved with satisfactory accuracy by a non-stationary formulation using the Transition k-kl-ω turbulence model. This is a three-parameter model of the RANS group, suitable for a narrow range of tasks, which allows modeling the transformation from laminar to turbulent flow, in which wall functions are used. The main solvable equations embedded in the mathematical model are the energy equation and the turbulence and turbulence-dissipation model equations.
The lifting force (causing the mixing of the flow) occurs due to convection, described by the Boussinesq approximation. This approximation consists in solving the system of Navier–Stokes equations and the energy equation using a special dependence of density on temperature only when calculating the volumetric force.
u ˜ t + u ˜ u ˜ = 1 ρ 0 p + ν u ˜ β θ g
θ t + u ˜ θ = χ Δ θ
div u ˜ = 0
ρ ( T ) = ρ 0 ( 1 β θ )
The approximation allows us to take into account the gravity field, that is, the aspiration of steam bubbles in the direction opposite to gravity, respectively, to take into account the orientation of the tube in space.

3. Results

Figure 7 displays the main results of numerical modeling in post-processing. Figure 7a shows that due to the instability of the process, convergence proceeds slowly. However, Figure 7b–d shows the totals for temperatures and fractions of the liquid and vapor phases of refrigerant. The results are indeed close to theoretical and experimental data in accuracy.
It should be noted that there are few research papers devoted to the joint processes of heat transfer, motion and boiling, both in the Russian Federation and in the world today. Basically, scientific laboratories use experimental data obtained earlier on stands and on-site objects. The complexity of research in the field of bubble boiling consists in the transience of processes (high intensity), and a large volume of newly obtained mixtures (for example, the same R407c is a non-azeotropic mixture of three refrigerants). In other words, there are a lot of new artificially produced mixtures, and full-scale work on numerical modeling and verification of at least heat transfer coefficients according to previously obtained formulas and ratios has not been carried out recently. The authors in their works try to partially compensate for the lack of “fresh” data that could be used by both design engineers and research engineers.

Comparison of the Received Data

In the course of the research, the authors were able to compare the temperatures on the back surface of the PVP: obtained experimentally, theoretically, and by numerical modeling. The discrepancy of the data turned out to be acceptable for engineering techniques. The closest results were obtained with theoretical calculation and numerical modeling. The difference was 0.1 °C. At the same time, when compared with experimental data on thermal sensors, the discrepancy is already more significant, approximately 0.3 °C, if we take into account the error of measuring instruments—0.5 °C. However, even with such values of temperature deviations, it is necessary to understand the order of the measured values—from −10 to + 90 °C. That is, with a measurement range of 100 ° C, the error turns out to be insignificant.
Thus, it can be said with confidence that the method of determining the temperature of the PVP and the data introduced into the boundary conditions during numerical modeling open up new opportunities for designers of energy technology complexes based on renewable energy sources, as new scientific approaches to temperature prediction appear, and consequently, more accurate determination of the design features of equipment, in particular, heat exchangers.
Numerical modeling and Formula (15) compared the obtained data. The discrepancy was about 0.1–0.5 °C.
If we compare the results in Table 2 and Table 4, we can conclude that numerical simulation gives results that are closer to the theoretical calculation data, although, for the analytical calculation, the convergence turned out to be quite good.

4. Discussion

4.1. Limitations of the Study

Regarding the limitations of the flow regime under study, it should be noted that after optimizing the mathematical model of heat transfer to the evaporator tube and a number of numerical calculations, it was found that the flow regime in the evaporator is dispersed-annular. Steam in the emulsion mixture is redistributed so that the main phase (liquid with 60% of the flow rate) moves along the pipe surface and in the form of small drops in steam, that is, in an annular mode.
The authors cannot guarantee the applicability of the obtained results for other ratios of the vapor and liquid phases.

4.2. Prospects for the Use of a Heat Pump in an Energy Complex for Seawater Desalination

The goal of creating plants with a thermal process for obtaining fresh water from sea water has led to the development of innovative concepts and design schemes with high efficiency. Particular attention is paid to the issues of optimizing operating modes, replacing expensive structural materials, pipe surfaces with other heater configurations, and reducing scale formation by an increase in the initial temperature of heating demineralized water in order to reduce the consumption of thermal and electrical energy. This goal can be achieved through the utilization of low-grade and waste heat. The cost of water produced in the desalination process is largely determined by the energy source used in the construction of large distillation plants [6,7,8,9,10,11,12].
In the future, the authors will develop a methodology for calculating the exergy fluxes of a desalination plant. An organic Rankine cycle with a turbo expander is used instead of a heat pump for the advanced scheme of the energy complex as a whole, and additional components for the use of renewable energy sources, such as a wind turbine and a solar concentrator, can be implemented. This approach allows for an integrated approach to energy conversion, using multiple sources for maximum efficiency.
First of all, new achievements in the field of low-temperature technology will be considered; for example, the authors of articles [15,16] explore the properties of refrigerant R407c. At the same time, such installations, which are energy technological complexes, must be environmentally friendly. The authors currently have scientific developments in this area. For example, in the article [17], the basics of constructing mathematical models for trigeneration and multigeneration complexes are considered.

5. Conclusions

When optimizing the mathematical model of heat transfer to the evaporator tube and a number of numerical calculations, it was found that the flow regime in the evaporator is dispersed-annular. The steam in the emulsion mixture is redistributed so that the main phase (liquid with a fraction of 60% of the flow rate) moves along the surface of the tube in the form of small droplets in steam—in an annular mode (Figure 2).
As the two-phase flow moves, there is a change to the opposite extreme case, when steam flows along the wall, and the liquid moves in the form of small droplets—in the dispersed mode (Figure 2). The wet steam reaches the saturation state, and then the dry steam overheats (process 4–1 in Figure 1).
It is also revealed that in order to obtain the correct refrigerant heat transfer coefficient for a coil type evaporator, it is necessary to derive a criterion similarity equation based on experimental data.
The results of the calculation of the evaporator with the turbulence model Transition k-kl-ω showed a sufficiently acceptable convergence for engineering computations with the calculation results and empirical data. The discrepancy between the values of the obtained temperatures ranged from 0.1 to 0.5 °C in the measured temperature range 100 °C.

Author Contributions

Conceptualization, K.O., S.A., A.K., Y.B., S.K. and A.O.; Methodology, K.O., S.A., A.K., Y.B., S.K. and A.O.; Validation, K.O.; Formal analysis, K.O.; Investigation, K.O.; Data curation, K.O.; Writing—original draft, K.O.; Writing—review & editing, K.O.; Visualization, K.O.; Supervision, K.O.; Project administration, K.O. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by the Russian Science Foundation grant No. 22-19-20011, https://rscf.ru/en/project/22-19-20011/ (accessed on 18 April 2023).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

i1enthalpy of condensate, kJ/kg;
i2enthalpy of vaporization, kJ/kg;
qboil.boiling heat flux density, kW/m2;
λthermal conductivity of liquid Freon, kW/(m·K);
ΔTstemperature difference between wall and liquid, K;
hgenthalpy of gas, kJ/kg;
hlenthalpy of liquid, kJ/kg;
Bcorrection factor of free convection heat transfer coefficient, non-dimensional;
νkinematic viscosity, m2/s;
σsurface tension coefficient, N/m;
Rigas constant for freon vapor, kJ/(kg·K);
ρggas density, kg/m3;
qconv.heat flux density at convective heat transfer kW/m2;
NuNusselt number, non-dimensional;
ReReynolds number, non-dimensional;
PrPrandtl number, non-dimensional;
ξcorrection factor that takes into account the features of the flow in a given system, non-dimensional;
tev. wallevaporator wall temperature, °C;
tPEC wallPVP surface temperature, °C;
qheat flow through PVP, kW/m2;
Rthermal 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;
Tssaturation temperature, K;
ppressure, Pa;
μdynamic viscosity, Pa·s;
ωflow velocity, m/s;
dtube 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);
λwallthermal conductivity of the wall, kW/(m·K);
δwall thickness, m;
Qboil.heat of saturated steam, kW;
xdegree of dryness, non-dimensional;
rheat of vaporization, kJ/kg;
lev.evaporation section length, m;
ls.hrefrigerant superheater section, m;
lttotal length of the PVP section, m;
tf.outfreon outlet temperature, °C;
cpfreon steam heat capacity, kJ/(kg·K);
( ρ E ) t partial differential of unstable heat transfer, kJ;
V ( ρ E + ρ ) a term that takes into account thermal conductivity, kJ;
k e f Δ T a term that takes into account the diffusion of mixtures, kJ;
j h j J j sum taking into account viscous dissipation, kJ;
S h a term that takes into account the enthalpy of inflows/drains, kJ;
u ˜ current velocity, m/s;
ρ0liquid density at some equilibrium temperature T0, kg/m3;
Θdeviation of temperature from equilibrium state (Θ = TT0), K;
χcoefficient of thermal diffusivity, m2/s;
g acceleration of gravity, m/s2;
βvolume expansion coefficient, K−1.

References

  1. 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]
  2. 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]
  3. 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]
  4. 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]
  5. 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]
  6. 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]
  7. 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]
  8. 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]
  9. 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]
  10. 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]
  11. 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]
  12. 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]
  13. Khrustalev, B.M. Heat and Mass Transfer; BNTU: Minsk, Belarus, 2007; 606p. [Google Scholar]
  14. 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]
  15. 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]
  16. 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]
  17. 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]
Figure 1. Designed experimental stand: (a) side view (overall dimensions of the stand: height—1706 mm, width—610 mm, depth—642 mm, angle of inclination of the photovoltaic panel—55 degrees, thickness of the photovoltaic panel redesigned for a thickness of 5 mm); (b) isometry; (c) basic diagram of the energy technological complex: Designations are as follows: 1—low-temperature heat source, 2—evaporator, 3—photovoltaic panel, 4—additional (backup) heater, 5—compressor, 6—condenser, 7—expander (capillary tube), 8—heated water, 9—valve, 10—main desalination tank, 11—steam compressor, 12—distillate heat exchanger-condenser, 13—desalinated water outlet, 14—salt water outlet.
Figure 1. Designed experimental stand: (a) side view (overall dimensions of the stand: height—1706 mm, width—610 mm, depth—642 mm, angle of inclination of the photovoltaic panel—55 degrees, thickness of the photovoltaic panel redesigned for a thickness of 5 mm); (b) isometry; (c) basic diagram of the energy technological complex: Designations are as follows: 1—low-temperature heat source, 2—evaporator, 3—photovoltaic panel, 4—additional (backup) heater, 5—compressor, 6—condenser, 7—expander (capillary tube), 8—heated water, 9—valve, 10—main desalination tank, 11—steam compressor, 12—distillate heat exchanger-condenser, 13—desalinated water outlet, 14—salt water outlet.
Energies 16 04405 g001
Figure 2. Heat pump cycle in the P-i diagram. Symbols and units of measurement: red lines—temperature, °C, blue lines—entropy, kJ/(kg·K), green lines—specific volume, m3/kg, horizontally—enthalpy in kJ/kg, vertically—pressure in bar.
Figure 2. Heat pump cycle in the P-i diagram. Symbols and units of measurement: red lines—temperature, °C, blue lines—entropy, kJ/(kg·K), green lines—specific volume, m3/kg, horizontally—enthalpy in kJ/kg, vertically—pressure in bar.
Energies 16 04405 g002
Figure 3. Phase transitions: (a) Emulsion flow mode of the vapor-liquid mixture; (b) Proportion of the liquid phase in the tube cross section, non-dimensional; (c) Ring flow mode (steam inside the liquid); (d) Dispersed flow mode (drops inside the steam). (a,b,d) are borrowed from the literature [13], (c) is the author’s numerical simulation (not described in the literature).
Figure 3. Phase transitions: (a) Emulsion flow mode of the vapor-liquid mixture; (b) Proportion of the liquid phase in the tube cross section, non-dimensional; (c) Ring flow mode (steam inside the liquid); (d) Dispersed flow mode (drops inside the steam). (a,b,d) are borrowed from the literature [13], (c) is the author’s numerical simulation (not described in the literature).
Energies 16 04405 g003
Figure 4. Block diagram of the calculation by the Gauss method. Designations of physical quantities are given in the nomenclature.
Figure 4. Block diagram of the calculation by the Gauss method. Designations of physical quantities are given in the nomenclature.
Energies 16 04405 g004
Figure 5. Dependence of solar panel wall temperature on ambient temperature and wind speed: the figure was developed by the authors based on the results of theoretical calculation in a specialized dynamic simulation program Math Lab Simulink.
Figure 5. Dependence of solar panel wall temperature on ambient temperature and wind speed: the figure was developed by the authors based on the results of theoretical calculation in a specialized dynamic simulation program Math Lab Simulink.
Energies 16 04405 g005
Figure 6. Model building steps: (a) Geometry of the evaporator model; (b) Mesh structure of the evaporator tube.
Figure 6. Model building steps: (a) Geometry of the evaporator model; (b) Mesh structure of the evaporator tube.
Energies 16 04405 g006
Figure 7. Post-processing: (a) Fluctuations of the discrepancy calculations of the dispersion-ring mode, non-dimensional; (b) Convergence of the degree of dryness at the evaporator output port, non-dimensional; (c) Mixing of the vapor-liquid mixture (Transition k-kl-ω), non-dimensional; (d) Temperature field (Transition k-kl-ω), °C; (e) Share of vapor in the volume (Transition k-kl-ω), non-dimensional.
Figure 7. Post-processing: (a) Fluctuations of the discrepancy calculations of the dispersion-ring mode, non-dimensional; (b) Convergence of the degree of dryness at the evaporator output port, non-dimensional; (c) Mixing of the vapor-liquid mixture (Transition k-kl-ω), non-dimensional; (d) Temperature field (Transition k-kl-ω), °C; (e) Share of vapor in the volume (Transition k-kl-ω), non-dimensional.
Energies 16 04405 g007
Table 1. Calculation data.
Table 1. Calculation data.
TitleDesignationCalculation Formula or Method of DeterminationUnits of MeasurementResult
Heat flowQsolGivenW/m2700
Thermal conductivity of siliconλsilGivenW/(m2·K)148
Absorption capacity of siliconAsilGiven-0.2
Degree of silicon blacknessεsilGiven-0.71
Silicon thicknessδsilGivenm0.001
Stefan-Boltzmann coefficientσGivenW·m−2·K−45.67 × 10−8
Thermal conductivity of glassλglGivenW/(m2·K)0.96
Absorption capacity of glassAglGiven-0.1
The degree of blackness of the glassεglGiven-0.93
Glass thicknessδglGivenm0.002
Length and width of the PVPlGivenm0.15675
PVP SquareF l 2 m20.025
Reflectivity of siliconRsilGiven-0.33
Heat absorbed by the PVPQ Q s o l F 1 R s i l W11.524
Heat absorbed by siliconQabs.sil Q s o l F A s i l W3.44
Heat passed through siliconQpast sil. Q s o l F 1 R s i l A s i l W8.084
Heat absorbed by glassQabs.gl. Q s o l F A g l W1.72
Heat passed through the glassQpast gl. Q s o l F 1 R s i l A s i l A g l W6.364
Air velocityωGivenm2/s1
Thermal conductivity of airλGivenW/(m2·K)2.4 × 10−2
Air viscosityνairGivenm2/s14.61 × 10−6
Prandtl number for airPrGiven-0.67
Ambient air temperaturetairSelected°C20
Ambient air temperatureTair t a i r + 273 K293
Outside temperature of the PVPtsur1Selected°C29.12
Outside temperature of the PVPTsur1 t sur 1 + 273 K302.12
Reynolds numberRe1, Re2 ω 1 ν a i r -1.073 × 104
Nusselt numberNu1, Nu2 0.67 R e 0.5 P r 1 / 3 -60.727
Heat transfer coefficientα1, α2 N u λ 1 W/(m2·K)9.298
Temperature of the inner surface of the PVPtsur2 t sur 1 Q s o l F 1 R s i l A g l A g l × δ s i l λ s i l + δ g l λ g l °C29.1
Losses from aboveqabove α 1 F T sur 1 T a i r W2.083
Silicon radiation lossesqsil.rad ε s i l σ F T sur 1 4 T a i r 4 W0.951
Losses from belowqund. α 2 F T sur 2 T a i r W2.083
Glass radiation lossesqgl.rad ε g l σ F T sur 2 + 273 4 T a i r 4 + Q s o l F A s i l + Q c o л F A g l W6.399
Heat losses in the systemQlosses q a b o v e + q s i l . r a d + q u n d . + q g l . r a d W11.516
Error rateΔQ Q Q l o s s e s Q 100 %0.067
Table 2. Data Comparison.
Table 2. Data Comparison.
Wind Speed, m/s1510
The temperature of the outer wall according to the theoretical calculation, °C40.7437.2034.30
External wall temperature according to analytical calculation, °C40.5737.0634.21
Table 3. Explanation of parameters.
Table 3. Explanation of parameters.
ФPDΓϕ
Kinetic energyk τ y u ¯ i x j p ε or β p k ω μ t σ k
Kinetic energy dissipation rateε c ε 1 ε k τ i j u ¯ i x j c ε 2 p ε 2 k u t σ e
Specific rate of dissipationω α ω k τ i j u ¯ i x j β p ω 2 σ ω u t
Table 4. Data Comparison (compare the results in Table 2 and Table 4).
Table 4. Data Comparison (compare the results in Table 2 and Table 4).
Wind Speed, m/s1510
Numerical modeling calculation, °C40.2336.9434.10
External wall temperature according to analytical calculation, °C40.5737.0634.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.

Share and Cite

MDPI and ACS Style

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

AMA Style

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 Style

Osintsev, 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 Style

Osintsev, 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

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop