Next Article in Journal
Fracture Modeling of Shale Oil and Gas Reservoirs in Texas
Next Article in Special Issue
Heat and Flow Characteristics of Aerofoil-Shaped Fins on a Curved Target Surface in a Confined Channel for an Impinging Jet Array
Previous Article in Journal
Measures and Prescriptions to Reduce Stray Current in the Design of New Track Corridors
Previous Article in Special Issue
A Direct Numerical Simulation Assessment of Turbulent Burning Velocity Parametrizations for Non-Unity Lewis Numbers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determination of Heat Transfer Correlations for Fluids Flowing through Plate Heat Exchangers Needed for Online Monitoring of District Heat Exchanger Fouling

1
Department of Energy, Cracow University of Technology, 31-155 Cracow, Poland
2
Department of Thermal Processes, Air Protection, and Waste Utilisation, Cracow University of Technology, 31-155 Cracow, Poland
*
Author to whom correspondence should be addressed.
Energies 2023, 16(17), 6264; https://doi.org/10.3390/en16176264
Submission received: 17 July 2023 / Revised: 17 August 2023 / Accepted: 22 August 2023 / Published: 28 August 2023

Abstract

:
This article deals with the problem of estimating the degree of fouling of plate heat exchangers (PHEs) used in district heating substations (where the working medium is water). A method for calculating the thermal resistance of fouling is proposed based on a comparison of the thermal resistance of a fouled and clean heat exchanger. The thermal resistance of the heat exchanger for both fouled and clean apparatuses is determined as the inverse of their overall heat transfer coefficient. In the method, the heat transfer coefficients necessary to determine the overall heat transfer coefficient of the clean exchanger are calculated using a modified Wilson method. Moreover, the heat transfer coefficients on the clean heat exchanger plates’ cold water side are determined based on experimental tests. The computational algorithm presented in this paper will make it possible to develop software to monitor and thus optimise the operation of district heating substations.

1. Introduction

PHEs are used in many areas of industry. They are commonly used in heating; refrigeration; district heating; and various industries, such as the food processing, brewing, pharmaceuticals, oil, gas, and chemical industries. These heat exchangers are composed of a number of plates with a thickness of approximately 0.3 (Brazed PHE) or 0.4 mm (Plate-and-Frame PHE) [1]. Depending on the area of use, heat exchangers differ in the way the plates are connected (brazed, welded, gasketed, and semi-welded) and the material of the plates (stainless steel or more sophisticated alloys, e.g., alloy 254 SMO, alloy C-276, titanium, palladium-stabilised titanium, Incoloy, and Hastelloy) [1]. Their popularity is due to their small size, but, at the same time, high thermal efficiency. For example, compared with shell-and-tube heat exchangers, the total heat transfer coefficient of the PHE is about twice as high [2]. The large value of the overall heat transfer coefficient is caused by the embossing on the plates, due to which the fluid flowing between the plates is disturbed, increasing the intensity of heat transfer. The embossing on the plates also increases their stiffness, which allows PHEs to operate at high pressures. The high heat output of PHEs is due not only to the higher values of overall heat transfer coefficients, but also to the fact that the embossing on the plates increases the surface area of heat transfer of these apparatuses.
A significant problem in the operation of PHEs is the efficient assessment of the fouling condition of the apparatus. They should be subjected to regular cleaning; otherwise, fouling will cause a decrease in the heat output of the heat exchanger, as well as an increase in the pressure drop of fluid. Contaminants are deposited in heat exchangers regardless of the fluid physical state [3,4,5]. The literature indicates that in some industries (e.g., the dairy industry [5,6,7,8]), the occurrence of contaminants in this type of apparatus is a serious problem.
Although fouling in heat exchangers is a problem, as mentioned earlier, in the case of PHEs, due to disturbed fluid flow, there is partial self-cleaning of the apparatuses. This means that PHEs do not foul as quickly as, for example, shell-and-tube heat exchangers. In addition, the deposition of fouling can be reduced by using the appropriate shape of plate embossments or appropriate coatings [3,4,5,9,10] or by choosing the optimal fluid velocity [11,12].
The optimisation of the cleaning process for PHEs at district heating substations could be aided by either the capability to determine the rate at which the layer of deposits builds up on the plates or by monitoring the thermal resistance of the heat exchangers. One of the methods proposed in the literature is to predict the formation of a deposit layer using numerical models of PHEs [7,8,11]. All numerical models to predict the formation of deposits in PHEs, described in [7,8,11], also considered their partial washout by the flowing fluid. In [11], the authors have numerically analysed the formation of a layer of fouling in a plate heat exchanger using 3D Realisable κ-ε with non-equilibrium wall functions modelling and the von Karman analogy for three different apparatus (with different geometries) and taking into consideration the influence of geometry on fluid dynamics. The analyses were supported by experimental studies. The accuracy of the proposed solution is satisfactory, but its limitation is that the research can only be performed for turbulent flow. In addition, more complex models based on Computational Fluid Mechanics (CFD) tend to have long calculation times. Most articles investigating fouling in PHEs are related to the dairy industry, and articles [7,8] are examples of this. Both articles take a different approach to modelling the formation of the contamination layer from the previous one, based on the principle of conservation of mass and energy. In paper [7], a 2D dynamic model of a PHE is proposed, allowing simulation with arbitrary cold and hot fluid flow configurations and with different fouling and cleaning kinetics. On the other hand, the paper [8] discusses the problem of milk pasteurisation in a High-Temperature-Pasteurisation-Process plate pasteuriser. A model of the entire plant, including PHEs, was made. The study adopted a simple 1D model of the rate of fouling of the PHEs with protein with minimum calculation time and investigated the effect of the thickness of the resulting fouling layer on the milk pasteurisation temperature. Both approaches presented are not applicable to PHEs in district heating substations because of the different mechanisms of fouling formation due to the different types of working fluid.
Regardless of the increasing development of CFD models, the best and simplest method of evaluating contamination in PHEs is to determine their thermal resistance via experiments. An example of just such an approach is presented in [13]. The authors of this paper aimed to determine the thermal resistance of a plate heat exchanger included in a district heating substation. Like the method developed and presented in this paper, the technique described by S.B. Genica et al. in [13] compares the thermal resistance of a fouled and clean heat exchanger. However, in their study, based on the four plate heat exchanger models analysed, the thermal resistance of the fouling was determined as a function of shear stress. In addition, the overall heat transfer coefficient for a clean heat exchanger was determined based on its manufacturer’s design software. The need to use this software is a limitation of this method. A literature review indicates that there is still a need to develop a quick and simple method to optimise the PHEs cleaning process.
This paper proposes a method for determining the thermal resistance of fouled PHEs using experimental studies. It consists of determining the thermal resistance of a clean and a fouled heat exchanger. If the difference between these resistances exceeds a preset value, the heat exchanger will require cleaning. The thermal resistance of a fouled heat exchanger is determined by measurements from the transformed Peclet equation. At the same time, the thermal resistance of a clean heat exchanger is determined as the sum of the thermal resistances of heat transfer and heat conduction, where the thermal transfer resistances are counted from correlations on the Nusselt number of the cold and hot fluid. These, in turn, are calculated from the proposed modified iterative Wilson method. This paper presents the results of calculations using a new modification of the Wilson method on a constructed experimental test stand.

2. Method for Evaluating Fouling of the Heat Exchanger

In the method, the degree of fouling of the PHE will be estimated by comparing the thermal resistance of the fouled R f o u l e d and clean R c l e a n heat exchanger:
Δ R = R f o u l e d R c l e a n
The thermal resistance of a contaminated heat exchanger will be calculated as the reciprocal of the overall heat transfer coefficient k f o u l e d determined using the measured data:
R f o u l e d = 1 k f o u l e d
In the case of a clean heat exchanger, the thermal resistance will be determined using the following equation:
R c l e a n = 1 k c l e a n = 1 α c + δ p λ p + 1 α h
where k c l e a n is the overall heat transfer coefficient of a clean heat exchanger in W/(m2K); α h and α c are the heat transfer coefficients of the heat exchanger’s hot and cold sides in W/(m2K), respectively; δ p is the thickness of the plate in m; and λ p is the thermal conductivity coefficient of the heat exchanger plate material in W/(mK). The heat transfer coefficients α h and α c will be calculated from the Nusselt numbers for hot fluid N u h and cold fluid N u c , respectively. The new method for determining constants in the correlations per Nusselt number is described in Section 4.
Cleaning of the heat exchanger will be necessary when the following inequality is met:
Δ R R l i m
where R l i m is the limit value of thermal resistance in (m2K)/W.

3. Determination of the Overall Heat Transfer Coefficient in the Heat Exchanger of a District Heating Substation

To determine the overall heat transfer coefficient k of the heat exchanger, the formula obtained from the Peclet equation can be used [14]:
Q ˙ = F k A Δ T m
where Q ˙ is the heat exchanger heat power in W, F is the correction factor that depends on the heat exchanger type ( F = 1 for parallel-flow and counter-flow heat exchangers [15]), k is overall heat transfer coefficient in W/(m2K), A is the total heat transfer surface area in m2, and Δ T m is the logarithmic mean temperature difference in K.
The heat output of the heat exchanger can be calculated from the energy balance equation for the hot fluid flowing through the apparatus, as Q ˙ h , based on the known mass flow rate m ˙ h and its average specific heat c ¯ w h , as well as the fluid temperatures at inlet T h 1 and outlet T h 2 .
Similarly, the heat flow rate in the heat exchanger can be calculated from the energy balance equation for the cold fluid, as Q ˙ c , based on the known mass flow rate m ˙ c , average specific heat c ¯ w c , and temperatures of the cold fluid at the inlet T c 1 and outlet T c 2 of the apparatus.
If there is no heat loss from the heat exchanger to the surroundings, then Q ˙ h = Q ˙ c . In practice, even though thermal insulation is used on the heat exchangers, heat losses occur and the heat flow rate value Q ˙ h is higher than the value Q ˙ c . In this case, the average value of the heat flow rate Q ˙ m that is transferred in the apparatus from hot to cold fluid can be calculated as [16]
Q ˙ m = Q ˙ h + Q ˙ c 2
Based on the average heat flow rate Q ˙ m , the heat transfer surface area A and the mean logarithmic temperature difference Δ T m , the overall heat transfer coefficient k for the fouled heat exchangers can be obtained using the transformed Equation (5):
k = Q ˙ m F A Δ T m

4. Determination of the Nusselt Number Correlation for Cold and Hot Fluid in a Clean Heat Exchanger

The method proposed in the article is experimental and modifies the popular Wilson graphical method. Thus, the laboratory studies to identify the cold fluid Nusselt number correlation are based on measuring cold and hot fluid parameters at a constant hot flow rate and for various cold fluid flow rate values. In turn, to determine the hot fluid’s Nusselt number correlation, laboratory testing in the heat exchanger must be carried out at a constant hot fluid flow rate and a variable flow rate of the hot fluid.
The developed method is based on iterative calculations. The single-iteration calculation algorithm presented in this chapter is for a single measurement point. Performing a series of iterations for several measurement points based on the described algorithm allows us to determine the Nusselt number correlation of the cold fluid N u c .
At the beginning of the iteration, the Nusselt number N u c is calculated for the selected N measurement points using a modified Colburn correlation [17]:
N u c = C c 1 R e c 0.8 P r c 0.33 + C c 2
in which the constant C c 1 replaces the value of 0.023 and an additional constant C c 2 is inserted. In Equation (8), R e c and P r c denote the Reynolds and the Prandtl numbers of the cold fluid, respectively.
Based on the calculated Nusselt number N u c , the heat transfer coefficient at the surface washed by the cold fluid α c is calculated using the following equation:
α c = N u c λ c d h y d r , c
where d h y d r , c is the hydraulic diameter of the channel the cold fluid flows through in m and λ c is the cold fluid thermal conductivity coefficient in W/(mK).
Next, using the heat transfer coefficient on the cold side of the heat exchanger α c , the heat transfer coefficient on the surface washed by the hot fluid α h is determined using the transformed Equation (3) for the overall thermal resistance of the PHE:
1 α h = 1 k 1 α c δ p λ p
Thereafter, based on the computed hot-side heat transfer coefficient α h , we can compute the hot fluid Nusselt number N u h using the following equation:
N u h = α h   d h y d r , h λ h
where d h y d r , h is the hydraulic diameter of the channel the hot fluid flows through in m and λ h is the hot fluid thermal conductivity coefficient in W/(mK).
The calculated values of the Nusselt number N u h for N measurement points allow for the determination of the constants C h 1 and C h 2 via approximation with a linear function of the following form:
N u h = C h 1 R e h 0.8 P r h 0.33 + C h 2
using the linear regression method.
In Equation (12), R e h and P r h are the Reynolds and the Prandtl numbers of hot fluid, respectively.
With the determined coefficients C h 1 and C h 2 (slope and vertical intercept of the function, respectively), the Nusselt number N u h is recalculated using Equation (12).
In the next step of the algorithm, the hot-side heat transfer coefficient α h is computed:
α h = N u h   λ h d h y d r , h
Based on the calculated heat transfer coefficient on the hot side of the heat exchanger α h , the cold-side heat transfer coefficient α c can be calculated as follows:
1 α c = 1 k 1 α h δ p λ p
The Nusselt number of the cold fluid N u c is then determined as follows:
N u c = α c   d h y d r , c λ c
From the calculated Nusselt number N u c and the Prandtl P r c and Reynolds R e c numbers of the cold fluid at the end of the iterative loop, the constants C c 1 and C c 2 are determined for N measurement points by approximating the data with a linear function (8) using the linear regression method.
The iterations are continued as long as the following conditions are met:
C c 1 j C c 1 j 1 C c 1 j ε c 1
and
C c 2 j C c 2 j 1 C c 2 j ε c 2
where j is the number of the iteration and ε c 1 and ε c 2 have small values (less than 0.001 for ε c 1 and 0.01 for ε c 2 ) .
The calculation algorithm for determining the constants C h 1 and C h 2 in the Nusselt number correlation for hot fluid can be written analogously to the algorithm for determining the constants C c 1 and C c 2 in the Nusselt number correlation for cold fluid presented above.
In the single-iteration loop, an apostrophe “ ” has been used in some symbols, i.e., in the heat transfer coefficient α and the Nusselt number N u for cold and hot fluid. This was performed as such because, in the single-iteration loop, these values appear and are calculated twice. The rule of thumb was adopted, with the symbol “ ” denoting the heat transfer coefficients calculated from the formula for the overall heat transfer resistance between the cold and hot fluid in the heat exchanger (Equations (10) and (14)) and the Nusselt numbers calculated from them (Equations (11) and (15)).
Plate heat exchangers in district heating substations may differ in design, i.e., in the number of plates or the type of corrugations. The proposed computational algorithm will allow the determination of constants in the correlations on the Nusselt numbers for cold and hot fluid flowing through a specific clean apparatus model. This, in turn, will make it possible to determine its total thermal resistance for any operation conditions.

5. Description of Laboratory Stands for the Verification of the Modified Wilson Method

The research discussed in this article will be conducted at two laboratory stands. One is a district heating substation belonging to the Municipal Heat Engineering Company (MPEC) in Krakow, which has existed for many years in one of the Cracow University of Technology buildings (Figure 1). The other is a new laboratory stand, a heat substation located in the laboratory of the Department of Energy (DE) of the Cracow University of Technology (Figure 2 and Figure 3a).
The heating substation in the DE laboratory is a parallel substation supplying two heating systems: a domestic hot water system (with PHE 1 in Figure 2) and a central heating system (with PHE 2 in Figure 2). The substation is supplied with hot water from an individual heat source, an electric boiler of the company ELTERM of 15 kW. The hot water reservoir used during the tests consists of two storage tanks with a capacity of 800 litres. A Danfoss heat exchanger model XB12M-1-16 (composed of 16 plates) is installed in the central heating system. In comparison, a Hexonic heat exchanger model LC110-30-2 (consisting of 30 plates) is installed in the domestic hot water system (Figure 3a). A drawing of a single heat exchanger plate with dimensions is shown in Figure 3b, while Figure 3c shows the type of embossing on the plates (chevron). Both heat exchangers are used as counter-flows.
The flow rate of cold and hot fluid in the system is measured using axial turbine flow meters ranging from 4 to 160 L/min (measuring accuracy: ±3% of measured value). Water temperatures in the installation are measured using type K (NiCr-Ni) sheathed thermocouple sensors with an outer diameter of 1.5 mm class 1. Figure 2 shows the location of sensors for measuring the temperature and flow rate. The location of sensors in the hot water system with PHE 1 are as follows: Th1,1 and Th2,1—hot water temperature measurement at the inlet and outlet, Tc1,1 and Tc2,1—cold water temperature measurement at the inlet and outlet, and mc1 and mh1—cold and hot water flow rate measurement. The location of sensors in the central heating system with PHE 2 are as follows: Th1,2 and Th2,2—measure the hot water temperature at the inlet and outlet, Tc1,2 and Tc2,2—measure the cold water temperature at the inlet and outlet, and mc2 and mh2—measure the cold and hot water flow rate. Compact pressure sensors with piezo-resistive measuring cells with temperature compensation and a pressure membrane made of special steel, whose basic error is ±0.1% of the final value, are used to measure the pressure in the installation. In contrast, the pressure drop in the heat exchangers is measured using differential pressure transmitters whose error for the basic measuring range is less than 0.05%. All sensors used for the measurements are from Ahlborn. All sensor measurement data are recorded by Ahlborn’s ALMEMO logger. The ALMEMO system enables real-time visualisation of the data collected from the sensors, data logging and the generation of reports in tables and graphs. Flow rate control in both systems is achieved via control valves. The research described in detail in this article will be carried out at the second stand.

6. Computation of the Cold Fluid Heat Transfer Coefficient in the Clean Heat Exchanger Using Measurement Data

The proposed algorithm, which is a modification of the Wilson method, was verified by carrying out test computations. They used measured data from the heat exchanger in the domestic hot water system of the laboratory stand described in the previous section.
Measurements carried out during the laboratory tests included measurements of water temperature at the inlet and outlet of the heat exchanger, as well as flow rates of cold and hot water. The data sampling time during the measurements was 1 s. Cold and hot water physical properties were determined for the average temperature values at the inlet and outlet of the heat exchanger, separately for each fluid.
The domestic hot water system under study uses a Hexonic LC110-30-2 PHE, and its heat transfer area A = 3.3   m 2 . The correction factor F = 1 , because the analysed heat exchanger is a counter-flow one.
The measurements were carried out so that the flow rate of cold water was constant (about 40 L/min.) while the flow rate of hot water was changed (in the range from about 5 to 40 L/min). From the archived database, N = 8 measurement points were selected. At the beginning of the test calculations, the mean heat flow rates Q ˙ m and overall heat transfer coefficients k of the PHE under test were determined for a selected number of measurement points ( i = 1 , ,   N ). The results of these calculations are shown in Table 1. Table 2 lists the average values of the temperatures T m , h and T m , c and the water properties, i.e., densities ρ h and ρ c , kinematic viscosities ν h and ν c , Prandtl numbers P r h and P r c and heat transfer coefficients λ h and λ c , for hot and cold water, respectively.
In the next step, the cold water Nusselt numbers N u c , i were calculated according to the modified Colburn correlation expressed by Equation (8). The initial values of the constants C c 1 0 and C c 2 0 in the first iterative loop for calculating the Nusselt number values N u c , i were taken as C c 1 0 = 0.05 and C c 2 0 = 0 . The water’s physical properties were computed for the average temperature using [18]. The calculated values of the Nusselt number for the selected measurement points N u c , i are shown in Table 3. In Table 3 et seq., i denotes the number of measurement points in the sequence.
Then, the cold-side heat transfer coefficients α c , i of the heat exchanger were calculated using Equation (9), and later, the hot-side heat transfer coefficients α h , i were calculated using Equation (10). In order to calculate α h , i , it is necessary to know the thickness of the plates in the heat exchanger, which is δ p = 0.001   m , and the thermal conductivity of the material from which they are made (stainless steel), which is λ p = 16   W / mK . In the next step, based on Equation (11), the values of the Nusselt number for hot water N u h , i were determined, which were then approximated by the linear function (12), which allowed us to determine the constants C h 1 1 and C h 2 1 , the values of which are −0.23977 and 40.00581, respectively.
The results of the subsequent calculations listed above needed to calculate the constants C h 1 1 and C h 2 1 are shown in Table 4. Table 4 also contains the values of the Reynolds R e h , i and Prandtl P r h , i numbers of hot water necessary to perform the approximation with a linear function.
In Table 4, for the first measurement point, the heat transfer coefficient α h , i and the Nusselt number N u h , i have negative values. This only occurs in the first iteration and is due to the fact that when determining α h , i from Equation (10), the sum of the heat transfer resistance on the cold side and heat conduction was larger than the total heat transfer resistance of the heat exchanger. Therefore, α h , i in this case has a negative value. The Nusselt number N u h , i was determined from the α h , i of Equation (11), so its value is also negative. The problem was eliminated in subsequent iterations, due to the increasingly better fit of the constants in the Nusselt number correlations.
After determining the constants C h 1 1 and C h 2 1 , the Nusselt numbers N u h , i from Equation (12) were again determined, and from them, the heat transfer coefficients α h , i from Equation (13) were calculated, followed by the heat transfer coefficients α c , i from Equation (14). The Nusselt number values for cold water N u c , i were then computed from Equation (15). The results of the calculations are presented in Table 5.
The first iteration ends with the approximation of the calculated values of the Nusselt numbers N u c , i by Equation (8) using the linear regression method, resulting in the constants C c 1 1 and C c 2 1 , which are 0.01988 and 2.55663, respectively. The coefficient of determination of the approximation using the linear function (8) performed is r 2 = 0.9314202 . The approximation using the linear function (8) is also shown graphically in Figure 4. The computed constants C c 1 1 and C c 2 1 are the initial values for the next iteration loop.
The computational algorithm for determining the Nusselt number correlation for a cold fluid is shown as a block diagram in Figure 5.
Table 6 shows the calculated Nusselt numbers N u c , i and heat transfer coefficients α c , i for cold water in five successive iteration steps.
The computational results of the successive iterations presented in Table 6 show that the proposed iterative method converges. The values of the constants C c 1 j and C c 2 j determined in subsequent iterations and the coefficients of determination r 2 are summarised in Table 7 and shown graphically in Figure 6a,b.
The modified Wilson method described in this article enables the determination of constants in the correlation of the Nusselt number for cold fluid flowing through a heat exchanger. Both conditions (16) and (17) were satisfied after 17 iterations. The final form of the correlation for the Nusselt number for cold water, when the hot water flow rate is about 40 L/min, is as follows:
N u c = 0.02493 R e c 0.8 P r c 0.33 + 2.01   for   30   R e c 190
However, it can also be used to determine the constants in the correlation of the Nusselt number of a hot fluid, but this will require subsequent measurements in which the hot fluid flow is constant and the cold fluid flow is variable.

7. Conclusions

This paper presents a method for determining the thermal resistance of fouling of PHEs in district heating substations and a method for determining the correlation of the Nusselt numbers of cold and hot fluid flowing through a clean heat exchanger (and, indirectly, the thermal resistance of the apparatus). The computational test of the second method was carried out for the measurement data obtained on a laboratory stand specially built for this purpose, which is a district heating substation of the domestic hot water system. Several iterations were carried out, and the method gave good results because they converged. As the final result of the presented research, a correlation with the Nusselt number for cold water flowing through the heat exchanger was obtained. It will allow us to determine the heat transfer coefficient on the cold side of the heat exchanger for each of its operating conditions.
The developed method will allow us to assess the degree of contamination of the PHEs. The method is fast and simple to use in practice. It will enable us to monitor the operation of PHEs in district heating substations and optimise heat exchanger cleaning intervals.

Author Contributions

Conceptualisation, J.T.; methodology, J.T.; software, T.R. and M.J.; validation, T.R. and M.J.; formal analysis, T.R.; investigation, T.R., M.J. and T.S.; data curation, T.R. and M.J.; writing—original draft preparation, T.R. and M.J.; writing—review and editing, T.R. and M.J.; visualisation, T.R.; supervision, M.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All data are available in the main text. Additional datasets related to this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arsenyeva, O.; Tovazhnyanskyy, L.; Kapustenko, P.; Klemeš, J.J.; Varbanov, P.S. Review of Developments in Plate Heat Exchanger Heat Transfer Enhancement for Single-Phase Applications in Process Industries. Energies 2023, 16, 4976. [Google Scholar] [CrossRef]
  2. Joniec, W. Wymienniki płytowe. Rynek Instal. 2010, 3, 74–77. [Google Scholar]
  3. Arsenyeva, O.P.; Crittenden, B.; Yang, M.; Kapustenko, P.O. Accounting for the thermal resistance of cooling water fouling in plate heat exchangers. Appl. Therm. Eng. 2013, 61, 53–59. [Google Scholar] [CrossRef]
  4. Bell, I.H.; Groll, E.A. Air-side particulate fouling of microchannel heat exchangers: Experimental comparison of air-side pressure drop and heat transfer with plate-fin heat exchanger. Appl. Therm. Eng. 2011, 31, 742–749. [Google Scholar] [CrossRef]
  5. Ahn, H.S.; Kim, K.M.; Lim, S.T.; Lee, C.H.; Han, S.W.; Choi, H.; Koo, S.; Kim, N.; Jerng, D.-W.; Wongwises, S. Anti-fouling performance of chevron plate heat exchanger by the Surface Modification. Int. J. Heat Mass Tran. 2019, 144, 118634. [Google Scholar] [CrossRef]
  6. Kananeh, A.B.; Peschel, J. Fouling in Plate Heat Exchangers: Some Practical Experience. In Heat Exchangers—Basics Design Applications; Mitrovic, J., Ed.; IntechOpen: Rijeka, Croatia, 2019; pp. 533–550. [Google Scholar]
  7. Thonon, B.; Grandgeorge, S.; Jallut, C. Effect of Geometry and Flow Conditions on Particulate Fouling in Plate Heat Exchangers. Heat Transf. Eng. 1999, 20, 12–24. [Google Scholar]
  8. Li, W.; Li, H.; Li, G.; Yao, S. Numerical and experimental analysis of composite fouling in corrugated plate heat exchangers. Int. J. Heat Mass Tran. 2013, 63, 351–360. [Google Scholar] [CrossRef]
  9. Lee, E.; Jeon, J.; Kang, H.; Kim, Y. Thermal resistance in corrugated plate heat exchangers under crystallisation fouling of calcium sulfate (CaSO4). Int. J. Heat Mass Tran. 2014, 78, 908–916. [Google Scholar] [CrossRef]
  10. Blanpain-Avet, P.; André, C.; Azevedo-Scudeller, L.; Croguennec, T.; Jimenez, M.; Bellayer, S.; Six, T.; Martins, G.A.S.; Delaplace, G. Effect of the phosphate/calcium molar ratio on fouling deposits generated by the processing of a whey protein isolate in a plate heat exchanger. Food Bioprod. Process. 2020, 121, 154–165. [Google Scholar] [CrossRef]
  11. Sharma, A.; Macchietto, S. Fouling and cleaning of plate heat exchangers: Dairy application. Food Bioprod. Process. 2021, 126, 32–41. [Google Scholar] [CrossRef]
  12. Indumathy, M.; Sobana, S.; Panda, R.C. Modelling of fouling in a plate heat exchanger with high temperature pasteurisation process. Appl. Therm. Eng. 2021, 189, 116674. [Google Scholar] [CrossRef]
  13. Genic, S.B.; Jacimovic, B.M.; Mandic, D.; Petrovic, D. Experimental determination of fouling factor on plate heat exchangers in district heating system. Energy Build. 2012, 50, 204–211. [Google Scholar] [CrossRef]
  14. Taler, D. Numerical Modelling and Experimental Testing of Heat Exchangers; Springer Nature: Cham, Switzerland, 2019. [Google Scholar]
  15. Srivastva, U.; Malhotra, R.K.; Kumar, K.R.; Kaushik, S.C. Comparative assessment of classical heat resistance and Wilson plot test techniques used for determining convective heat transfer coefficient of thermo fluids. Therm. Sci. Eng. Prog. 2019, 11, 111–124. [Google Scholar] [CrossRef]
  16. Hewitt, G.F. (Ed.) Handbook of Heat Exchanger Design; Begell House: New York, NY, USA, 1992. [Google Scholar]
  17. Seiwert, J. Independent measurement of condensation and vaporisation heat transfer coefficients: An alternative to the Wilson plot method suitable for multiphase exchangers. Int. J. Heat Mass Transf. 2019, 136, 1171–1185. [Google Scholar] [CrossRef]
  18. Verein Deutscher Ingenieure. VDI Heat Atlas; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
Figure 1. General view of the tested MPEC district heating substation.
Figure 1. General view of the tested MPEC district heating substation.
Energies 16 06264 g001
Figure 2. Scheme of the laboratory stand: red line—hot water, blue line—cold water.
Figure 2. Scheme of the laboratory stand: red line—hot water, blue line—cold water.
Energies 16 06264 g002
Figure 3. The heat substation located in the laboratory of the DE: (a) general view of the laboratory stand, (b) view of the single plate of Hexonic heat exchanger model LC110-30-2 with dimensions, and (c) view of a single panel with embossing included.
Figure 3. The heat substation located in the laboratory of the DE: (a) general view of the laboratory stand, (b) view of the single plate of Hexonic heat exchanger model LC110-30-2 with dimensions, and (c) view of a single panel with embossing included.
Energies 16 06264 g003
Figure 4. Graphical illustration of the linear regression method (red line) using calculated values of the Nusselt number N u c , i (symbols).
Figure 4. Graphical illustration of the linear regression method (red line) using calculated values of the Nusselt number N u c , i (symbols).
Energies 16 06264 g004
Figure 5. Block diagram of the computational algorithm for determining the Nusselt number correlation for a cold fluid.
Figure 5. Block diagram of the computational algorithm for determining the Nusselt number correlation for a cold fluid.
Energies 16 06264 g005
Figure 6. Results of subsequent iterations: (a) values of the constants C c 1 j ; (b) values of the constants C c 2 j .
Figure 6. Results of subsequent iterations: (a) values of the constants C c 1 j ; (b) values of the constants C c 2 j .
Energies 16 06264 g006
Table 1. Selected measurement points from tests carried out on 12 May 2023 and computed values for the mean heat flow rate Q ˙ m and the overall heat transfer coefficient k .
Table 1. Selected measurement points from tests carried out on 12 May 2023 and computed values for the mean heat flow rate Q ˙ m and the overall heat transfer coefficient k .
Time
T h 1
T h 2
T c 1
T c 2
V ˙ h
V ˙ c
c ¯ w h
c ¯ w c
Q ˙ h
Q ˙ c
Q ˙ m
Δ T m
k
°C°C°C°CL/minL/minJ/(kgK)J/(kgK) WWWKW/(m2K)
11:46:5063.757.412.663.640.14.94184.364178.6017,319.217,279.217,299.27.32715.94
11:49:4962.950.212.261.240.110.24182.604178.6134,969.734,575.634,772.711.68901.89
11:54:3460.643.312.056.339.715.04180.944178.7647,245.946,010.846,628.313.601038.80
11:57:4459.039.011.851.739.920.04180.094179.0754,956.655,301.155,128.815.131104.21
12:00:1056.934.911.646.539.824.94179.404179.6160,371.060,278.860,324.915.991143.08
12:01:2456.332.511.643.139.830.04179.144180.0765,346.765,588.765,467.716.761183.96
12:03:0954.829.611.440.039.734.84178.844180.6169,072.569,118.169,095.316.441273.49
12:03:5554.028.011.737.439.739.74178.734181.0571,296.370,883.371,089.816.451309.60
Table 2. Average temperature values and physical properties of hot and cold water for selected measuring points from tests carried out on 12 May 2023.
Table 2. Average temperature values and physical properties of hot and cold water for selected measuring points from tests carried out on 12 May 2023.
Time T m , h T m , c ρ h ρ c ν h ν c P r h P r c λ h λ c
°C°Ckg/m3kg/m3×10−6 m2/s×10−6 m2/s--W/(mK)W/(mK)
11:46:5060.638.1983.0992.80.47010.67752.964.480.6510.626
11:49:4956.636.7985.0993.30.49800.69593.164.610.6470.624
11:54:3452.034.2987.2994.20.53400.73173.424.880.6430.621
11:57:4449.031.8988.5995.00.55970.76863.605.160.6390.618
12:00:1045.929.1989.8995.80.58920.81403.825.510.6360.614
12:01:2444.427.4990.4996.20.60450.84503.935.750.6340.611
12:03:0942.225.7991.3996.70.62830.87704.116.000.6310.609
12:03:5541.024.6991.8997.00.64200.90064.216.180.6300.607
Table 3. List of selected measuring points from 12 May 2023 and calculated overall heat transfer coefficients k i and Nusselt numbers N u c , i of cold water ( i = 1 , ,   8 ).
Table 3. List of selected measuring points from 12 May 2023 and calculated overall heat transfer coefficients k i and Nusselt numbers N u c , i of cold water ( i = 1 , ,   8 ).
iTime T h 1 , i T h 2 , i T c 2 , i T c 1 , i m ˙ h , i m ˙ c , i Δ T m , i k i N u c , i
°C°C°C°Ckg/skg/s°CW/(m2K)-
111:46:5063.757.412.663.60.6570.0817.32715.941.284
211:49:4962.950.212.261.20.6580.16911.68901.892.282
311:54:3460.643.312.056.30.6530.24913.601038.803.041
411:57:4459.039.011.851.70.6570.33215.131104.213.748
512:00:1056.934.911.646.50.6570.41315.991143.084.359
612:01:2456.332.511.643.10.6570.49816.761183.964.980
712:03:0954.829.611.440.00.6560.57816.441273.495.520
812:03:5554.028.011.737.40.6560.66016.451309.606.064
Table 4. Results of first iteration calculations for i = 1 , ,   8 measurement points to determine the constants C h 1 1 and C h 2 1 .
Table 4. Results of first iteration calculations for i = 1 , ,   8 measurement points to determine the constants C h 1 1 and C h 2 1 .
i k i N u c , i α c , i α h , i N u h , i R e h , i P r h , i
W/(m2K)-W/(m2K)W/(m2K)---
17161.284651−4002−7.58344.422.96
29022.2821154763014.55325.113.16
310393.041153050629.72300.153.42
411043.748187638307.39287.823.60
511434.359216833116.43272.723.82
611844.980246630525.94265.813.93
712745.520272232636.38255.104.11
813106.064298231576.18249.654.21
Table 5. Results of first iteration calculations for i = 1 , , 8 measurement points to determine the constants C c 1 1 and C c 2 1 .
Table 5. Results of first iteration calculations for i = 1 , , 8 measurement points to determine the constants C c 1 1 and C c 2 1 .
i k i N u h , i α h , i α c , i N u c , i R e c , i P r c , i
W/(m2K)-W/(m2K)W/(m2K)---
17163.273172714162.7931.154.48
29024.171218818513.6663.134.61
310395.504286619913.9688.284.88
411046.054313621044.20112.075.16
511436.857353320824.19131.745.51
611847.217370821574.36152.905.75
712747.818400023604.78170.896.00
813108.112414124354.95189.856.18
Table 6. Calculated values of Nusselt numbers N u c , i and heat transfer coefficients α c , i in five subsequent iterations.
Table 6. Calculated values of Nusselt numbers N u c , i and heat transfer coefficients α c , i in five subsequent iterations.
Iteration Number (j)
(1)(2)(3)(4)(5)
i N u c , i α c , i N u c , i α c , i N u c , i α c , i N u c , i α c , i N u c , i α c , i
-W/(m2K)-W/(m2K)-W/(m2K)-W/(m2K)-W/(m2K)
11.284651.393.0671556.283.0121528.462.9621503.042.9171480.23
22.2821154.343.4641752.463.4261733.043.3901715.183.3581699.00
33.0411529.813.7661894.623.7401881.613.7161869.503.6941858.38
43.7481875.734.0472025.234.0332018.134.0192011.334.0072004.87
54.3592167.464.2902133.154.2862131.094.2812128.844.2762126.37
64.9802465.994.5372246.624.5432249.644.5482251.954.5512253.49
75.5202722.054.7512343.254.7662350.644.7792356.894.7892361.89
86.0642981.924.9682442.904.9922454.695.0132464.885.0302473.35
Table 7. Values of C c 1 j and C c 2 j and the coefficients of determination r 2 for successive iterations.
Table 7. Values of C c 1 j and C c 2 j and the coefficients of determination r 2 for successive iterations.
j C c 1 j C c 2 j r 2
10.019882.560.93142
20.020712.480.93501
30.021452.410.91899
40.022102.350.94114
50.022662.300.94367
60.023132.250.94587
70.023532.210.94776
80.023862.170.94938
90.024122.140.95076
100.024342.120.95194
110.024512.090.95294
120.024642.070.95381
130.024742.060.95455
140.024822.040.95519
150.024882.030.95574
160.024912.020.95623
170.024932.010.95666
180.024942.000.95704
190.024941.990.95739
200.024931.980.95770
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

Romanowicz, T.; Taler, J.; Jaremkiewicz, M.; Sobota, T. Determination of Heat Transfer Correlations for Fluids Flowing through Plate Heat Exchangers Needed for Online Monitoring of District Heat Exchanger Fouling. Energies 2023, 16, 6264. https://doi.org/10.3390/en16176264

AMA Style

Romanowicz T, Taler J, Jaremkiewicz M, Sobota T. Determination of Heat Transfer Correlations for Fluids Flowing through Plate Heat Exchangers Needed for Online Monitoring of District Heat Exchanger Fouling. Energies. 2023; 16(17):6264. https://doi.org/10.3390/en16176264

Chicago/Turabian Style

Romanowicz, Tomasz, Jan Taler, Magdalena Jaremkiewicz, and Tomasz Sobota. 2023. "Determination of Heat Transfer Correlations for Fluids Flowing through Plate Heat Exchangers Needed for Online Monitoring of District Heat Exchanger Fouling" Energies 16, no. 17: 6264. https://doi.org/10.3390/en16176264

APA Style

Romanowicz, T., Taler, J., Jaremkiewicz, M., & Sobota, T. (2023). Determination of Heat Transfer Correlations for Fluids Flowing through Plate Heat Exchangers Needed for Online Monitoring of District Heat Exchanger Fouling. Energies, 16(17), 6264. https://doi.org/10.3390/en16176264

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