**Chapter 12** Digital Twin of Heat Exchanger

## *Miha Bobič*

## **Abstract**

How to create proper digital twin of plate heat exchange. The chapter will investigate alternatives for modeling and propose the best choices in modeling to trade-off between accuracy of the calculation vs. speed. Lump sum vs. distributed model, how to make plate model discrete, temperature distribution models, and details in heat exchange wall modeling will be considered. Last part would focus on experimental evaluation of the simulation using IR temperature vision measurements. Conclusions would take concrete examples of dynamic responses of heat exchange units, comparing numeric simulation with experimental results.

**Keywords:** dynamic behavior, heat exchanger, models with distributed parameters, digital twin

## **1. Introduction**

Heat exchangers are used for heat and cool transfer in different process industries, heating/cooling, and other applications. Plate heat exchangers are for majority of applications where pressure is not so high, the most cost-efficient solution. Brazed are used for smaller sizes, and more clean media, while gasketed are meant more for applications where occasional cleaning is needed.

The simplest control loop depicted in **Figure 1** is to control the outlet temperature on secondary side of the heat exchanger by controlling the flow on primary at given conditions of inlet temperatures on primary and secondary sides.

The system consists of a heat exchanger as the central component of the control system. The control can be performed using an electronic controller (usually due to the speed of response of the proportional-integral controller system) and the valve with a drive. In the return control loop, the measuring sensor for the temperature represents the guiding control variable.

The speed of process, seen from dynamics point of view, is a major challenge, as the time constants of the measuring system for the temperature of the heat carrier or hot water need to be small due to a relatively small control volume. Time constants may be less than few seconds.

## **2. Literature overview**

In the field of heat exchangers, we can highlight the contribution [1], which is an overview of the current literature on the state of heat exchangers. Aslam Bhutta et al. [2] and Zhang et al. [3] discuss the review of contributions dealing with the topic of

#### **Figure 1.** *Control loop schematics.*

numerical simulation of heat transmissions. In Ref. [4] the case of a numerical simulation of a plate heat exchanger is considered. Freund and Kabelac [5] also deals with modeling of dynamic responses of plate heat transmissions [6]. The topic of evaporation in heat exchangers is discussed. The paper [7] deals with the experimental determination of the operational factors of heat exchangers in district heating for heating and hot sanitary water. Sources [8, 9] present an experimental method of visualizing the flow in a flat heat exchanger with a microplate structure [10]. The authors are also involved in visualizing the flow and determining local heat transfer. The source [11] analyzes the local heat transfer conditions. While the source [12] addresses the impact of the final plates on the operation of the heat exchanger. The paper [13] focuses on a study of the sensitivity of the geometric properties of a heat exchanger to its efficiency. The paper [14] deals with the technical requirements and cost analysis of heat transmissions in low-temperature systems. A similar theme of low-temperature systems is also in the paper [14, 15], with the emphasis in publication [15] on the impact of low-temperature systems on heating and the preparation of hot water without accumulation and in the source [14], on the main changes in the structure of the exchangers and thus the dynamic properties of heat transfer for such systems. A comparison between the plate and the tube heat exchanger is made in Ref. [16]. The source [17] discusses experimental and simulation methods in the field of heat exchangers. In Ref. [5], a study is carried out on how to determine local heat transfer coefficients in heat exchangers [18]. It also addresses the sensitivity of the parameter analysis to the numerical model of the heat exchanger. The link between the flow rate and the heat transfer under the micro heat exchanger conditions is given in Ref. [19]. Laszczyk [20] provides a method of a simplified heat exchanger model. Improved heat transfer in heat exchangers deals with [21]. Heat exchangers shell and tube are treated in Refs. [22]. Sharifi et al. [23] present a dynamic simulation of the heat exchanger. General dimensioning and dynamic analysis can be found in Refs. [24, 25]. The more generic described heat transfer dynamics depend, of course, on the configuration and sample on the heat exchanger plate, the various examples described in Refs. [26–28]. Many studies use CFD to analyze current patterns, distribution of pressure field, and temperature field, such as Ref. [4]. Sharifi et al. [23] provide a general model of heat transfer. Semiempirical equations are important for determining the coefficient of heat transfer, such as in Refs. [29–31]. Carlomagno and Cardone [17] describe the method of noninvasive measurement of heat exchanger parameters. Aprea and Renno [32] describe the design of tests for the portable evaporator function, that is, similar functionality to a

heat exchanger. The use of MATLAB Simulink is shown in Ref. [33]. Integrated thermofluid dynamics and modeling reliability are shown in Ref. [34]. The literature also lists several approaches to modeling heat exchangers for different purposes, for example in Ref. [35], a model mode that is object-oriented using the Modelica programming language and provides a model based on the control volume method. Using Laplace's transformation, dynamic analysis of plate heat transfer is given in Ref. [36]. While in the Ref. [37] heat exchanger is described with the hyperbolic first order differential equations. Furthermore, in Ref. [38], a simplified mathematical model of a heat exchanger based on the final volume method is also used in the interpolation scheme. Similar is used in Ref. [39]. Further more modelling of the heat exchangers in the system and elements are the system are described in Refs. [40–44].

## **3. Basics of modeling**

Seen from dynamic behavior point of view heat exchanger is a system with distributed parameters, meaning that the temperature continually changes throughout its length, thus traditional approach of direct Laplace transformation becomes too complex to solve, as we end up in double Laplace transformation, causing numeric problems with solution.

If we want to reduce the model to lump sum, then discretization is needed. Of course, the question is how many sections we should use to get good approximation with minimum numeric calculation load. The answer is in what kind of constitutional equation we use for temperature distribution in every cell. One has three options, one is if cells represent perfect mixing cell, and the assumption is that temperature immediately drops to exit temperature, and the second option is the opposite where the assumption is that mixing is imperfect, so the temperature is kept on inlet value until it leaves the cell. One can intuitively feel that those options are not realistic. The third option is to use distribution, it can be linear or to get closer to reality logarithmic. Experiments show that linear distribution is the best choice between number of cells, calculation effort, and minimized calculation error.

The second important question is how to model the thin metal wall between two media streams. One can model it as a conductive element of just ignoring the conduction. Given the ratio between heat convection resistance factor and heat conduction in a thin metal wall, it can be very well assumed that this is low resistance to heat compared to heat transfer for media to wall.

It is also important that proper material and turbulence models are used, thus in next chapter, we will investigate properties of water.

#### **3.1 Basic properties**

The material properties depend, of course, on different parameters. In the example shown, the impact of pressure and temperature is greatest, but we are particularly interested in variables:


The density gives us the dependence of the mass change per volume unit. In general, we can write down the dependence:

$$d\rho = \left(\frac{\partial \rho}{\partial T}\right)\_p dT + \left(\frac{\partial \rho}{\partial p}\right)\_T dT \tag{1}$$

When the coefficient of temperature extension is applied, the following shall be obtained:

$$\frac{d\rho}{\rho} = -\beta dT + \frac{1}{\varepsilon} dp \tag{2}$$

The eq. (2) specifies that the density depends on the temperature and the suitability of the substance. In this case, the heat carrier will be used as a heat carrier for district heating, and for this purpose, the constitution equation is by source [45]:

$$\rho = \left(1000 + 0,067\left(1 - e^{-0,016T}\right) - 0,001T\right) \text{za } 0^\circ \text{C} < T < 50 \,^\circ \text{C} \tag{3}$$

$$\rho = \left( \mathbf{1.067} - e^{-0.023(T - 50)} - \mathbf{0.001T} \right) \text{za } \mathbf{50}^{\circ}\text{C} < T < 200^{\circ}\text{C} \tag{4}$$

Viscosity can be obtained from the viscous tensor as:

$$
\pi = \eta \frac{dv\_{\infty}}{dy} \tag{5}
$$

The viscosity is therefore greater, the smaller the shear deformity. For water, the kinematic viscosity equation can be summarized as [45]:

$$\nu = \mathbf{1}, \mathbf{7}5e^{-0.029T}\mathbf{10}^{-6} \qquad \text{for } \mathbf{0}^{\circ}\text{C} < \mathbf{T} < \mathbf{20}^{\circ}\text{C} \tag{6}$$

$$\nu = 0.98e^{-0.02(T-20)} \text{10}^{-6} \text{ for } 20^{\circ}\text{C} < \text{T} < 50^{\circ}\text{C} \tag{7}$$

$$
\nu = 0.55 \left( \frac{50}{T} \right)^{0.9} 10^{-6} \qquad \text{for } 50 \, ^\circ \text{C} < T < 200 \, ^\circ \text{C} \tag{8}
$$

Specific heat is a greatness that determines how much energy we need to produce to heat the substance by 1 K.

$$c\_p(T) = \frac{1}{m} \frac{dh}{dT} \tag{9}$$

or it can be written as an approximation equation as:

$$c\_p = -7 \cdot 10^{-7} T^3 + 8 \cdot 10^{-5} T^2 - 0,0031T + 4,2161 \text{ for } 0^\circ \text{C} < \text{T} < 50^\circ \text{C} \tag{10}$$

$$c\_p = -6 \cdot 10^{-11} T^5 + 2 \cdot 10^{-8} T^4 - 3 \cdot 10^{-6} T^3 + 0,0002 T^2 - 0,\tag{11}$$

$$10069T + 4,269 \text{ For } 51^{\circ}\text{C} < T < 100^{\circ}\text{C}$$

## **4. Heat exchanger**

Heat exchanger is the heart of the thermal substation, as we control heat transfer with changing thermal and hydraulic parameters. Despite the different types of heat transferors, in the practice of district heating, plate heat exchangers have established themselves, namely, for smaller dimensions and relatively good media quality (clean

**Figure 2.** *Schematics of heat exchanger.*

and noncorrosive), dedicated laptops, and all other applications are covered by gasketed heat exchangers.

In general, all plate laptops have a sample on plates shaped like a fish bone, the only other ones on the market are Danfoss's microplate heat exchangers, which have a pointshaped sample. The sample on the heat transfer plates is intended to increase the heat transfer area and increase turbulence in order to increase the heat transfer coefficient.

Modern heat exchangers also have an input distribution section and an output collection section to increase the efficiency of the panels, and the task of these two sections is to redistribute the medium to flow evenly throughout the plate and thus increase the efficiency of the surface.

There are quite a few dynamic models of heat transfers in the literature [1]. The simplest assumes a first-order system with a possible delay, but such a heat exchanger model is inappropriate because the validity of the mathematical model is too limited. More interesting are the models where the dual Laplace transformation is used, but they need to be linearized, which in turn is in too limited force range for our case. Therefore, the heat exchanger will be considered as a system with distributed parameters which means that the heat exchanger will be considered as several consecutive systems with no distributed parameters. The directions of the currents in the heat exchanger are schematically shown in **Figure 2**.

### **4.1 General assumptions and limitations of the mathematical model**

When the heating medium enters the heat exchanger, the more distant plates (for example, 100 to 200), represent a greater hydraulic resistance against the water flow. Therefore, for the first time, water flows over nearby plates faster and the second time there is a certain delay due to the flow of water. As this will be negligible for our case (about 20 panels), we will not take this effect into account.

In general, the heating medium has a higher temperature than the surrounding area, so the heat transfer takes place not only from the warmer medium to the cooler medium but also in the surrounding area. Due to the relatively good thermal insulation, we will neglect these losers. The translation of heat over steel plates exists, but due to the thickness of these (about 0,2 mm) and the relatively high coefficient of heat translation, this will be ignored. We will also pre-bet that the heat exchanger is clean without the suspended limestone on the walls. The most unrealistic assumption will be the uniformity of the distribution of the heating medium by plate. Assuming that the mixing section is relatively short toward the entire surface, we can also take this assumption into account.

In the heat exchanger model, we considered that the heat exchanger is a nonlinear system with distributed parameters. Therefore, in the model, the heat exchanger will be divided into several smaller subunits, in which linearity will be assumed. The heat exchanger has different dynamic responses according to the input parameters, whether it is the temperature or the mass flow rate of the heat carrier. It is important to determine the dynamic properties of the heat exchanger, however, to respond to the temperature change.

#### **4.2 Static characteristics**

The static characteristic of the heat exchanger is the ratio of mass flow rate to the heat ousted by the energy equation:

$$P = \dot{m}c\_p \Delta T \tag{12}$$

If only this is considered and when the efficiency of the heat exchanger is introduced, which determines the ratio between actual thermal power and maximum thermal power:

$$
\varepsilon = \frac{P\_{dcf}}{P\_{maks}} \tag{13}
$$

and for primary side

$$
\varepsilon\_h = \frac{T\_{11} - T\_{12}}{T\_{11} - T\_{21}} \tag{14}
$$

#### **4.3 Dynamic characteristics**

The mathematical model of the heat exchanger is basically the problem of switching heat from a hot medium over the wall to a cold medium. This can be described using the energy equation for a cold and hot medium:

$$\rho\_i \mathfrak{e}\_{pi} \frac{\partial^2 T\_i}{\partial \mathbf{x} \partial t} = \nabla \cdot (\lambda\_i \nabla T\_i) - \underline{w\_i} \rho\_i \mathfrak{e}\_{pi} \nabla T\_i + r\_i(-\Delta I\_i) + \nu\_{\text{Ti}} \left(\rho\_i T\_i \frac{\partial p\_i}{\partial t} + \underline{w\_i} T\_i \nabla p\_i\right) \tag{15}$$

Just as in the case of a temperature sensor, there are no chemical reactions and generated heat due to the mechanical part, which is why the last two articles of the eq. (15) are off. However, since the temperature gradient and the translation of heat over the wall are ignored, the equation can be simplified by consideration in:

$$m\_i c\_{pi} \frac{dT\_i}{dt} = \dot{m}\_i c\_{pi} \frac{dT\_i}{d\mathbf{x}} \pm a\_i A\_i (T\_{hi} - T\_{ci}) \tag{16}$$

Therefore, assuming one dimension model with complete mixing in the cell we obtain a model based on the energy balance of the heat exchanger:

• heat exchanger hot side:

$$m\_h \mathfrak{c}\_{ph} \frac{dT\_h}{dt} = \dot{m}\_h \mathfrak{c}\_{ph} (T\_{hi} - T\_{ho}) - a\_h A\_h (T\_h - T\_w) \tag{17}$$

• cool side of the heat exchanger:

$$m\_c c\_{pc} \frac{dT\_c}{dt} = \dot{m}\_c c\_{pc} (T\_{ci} - T\_{co}) + a\_c A\_c (T\_w - T\_c) \tag{18}$$

• and the wall of heat exchanger:

$$m\_w c\_{pw} \frac{dT\_w}{dt} = a\_h A\_h (T\_h - T\_w) - a\_c A\_c (T\_w - T\_c) \tag{19}$$

However, since the heat exchanger system is a system with distributed parameters, the analytical solutions of the above model are impossible. In the trace of this, the heat exchanger will be divided into several units (**Figure 3**). For each unit, the balance sheet equations can be recorded:

**Figure 3.** *(a) Heat exchanger section, (b) one cell, (c) resistance scheme.*

• heat exchanger hot side:

$$m\_h c\_{ph} \frac{dT\_{hi+1}}{dt} = \dot{m}\_h c\_{ph} (T\_{hi} - T\_{hi+1}) - a\_h A\_h \left(\frac{T\_{hi} + T\_{hi+1}}{2} - T\_{ui}\right) \tag{20}$$

• cool side of the heat exchanger:

$$m\_{\varepsilon}c\_{p\varepsilon}\frac{dT\_{ci-1}}{dt} = \dot{m}\_{\varepsilon}c\_{p\varepsilon}(T\_{ci-1} - T\_{ci}) + a\_{\varepsilon}A\_{\varepsilon}\left(T\_{wi} - \frac{T\_{ci-1} + T\_{ci}}{2}\right) \tag{21}$$

• and the wall of heat exchanger:

$$m\_w c\_{pw} \frac{dT\_{wi}}{dt} = a\_h A\_h \left(\frac{T\_{hi} + T\_{hi+1}}{2} - T\_{wi}\right) - a\_c A\_c \left(T\_w - \frac{T\_{ci-1} + T\_{ci}}{2}\right) \tag{22}$$

An important assumption is what the temperature distribution in each cell will be. In general, it can be assumed to jump (which is impossible in the heat exchanger, but if the sections are infinitesimally narrow, the assumption is good enough for the apportionment, while greatly simplifying the calculation), which means that the medium has an exit temperature from the cell immediately upon entering the cell. Such a model is very easy to use but the assumption is valid only if mixing is practically ideal. A more realistic model is where the temperature changes linearly along the cell. Thus, the model is numerically still sustainable, but to us, results that are close to real are why this model is used in the above equations. The third model considers the real logarithmic distribution of the temperature in the cell (**Figure 4**).

By rearranging eqs. (20)–(22) and inserting variables:

$$
\tau\_{h/c\_i} = \frac{m\_{h/c\_i}}{\dot{m}\_{h/c\_i}} \tag{23}
$$

$$\delta\_{w\_i} = \frac{m\_{w\_i} c\_{p\_{h/\epsilon\_i}}}{A\_i} \tag{24}$$

$$\gamma\_{h/c\_i} = \frac{\mathcal{A}\_i \alpha\_{h/c\_i}}{m\_{h/c\_i} \mathcal{C}\_{p\_{h/c\_i}}} \tag{25}$$

One can get a matrix for the i-th cell:

**Figure 4.** *Temperature distribution in the cell.*

… *<sup>T</sup>*\_ *hi*�<sup>1</sup> *T*\_ *hi* … *T*\_ *wi* … *T*\_ *ci <sup>T</sup>*\_ *ci*�<sup>1</sup> … ………… … …… …… …………… ………… … …… …… …………… … <sup>0</sup> <sup>1</sup> *τhi* � <sup>1</sup> <sup>þ</sup> *<sup>γ</sup>hi τhi* … <sup>0</sup> *<sup>γ</sup>hi τhi* ………… 0 ………… … …… …… …………… … <sup>0</sup> *<sup>γ</sup>hi δwi* … <sup>0</sup> � *<sup>γ</sup>hi* <sup>þ</sup> *<sup>γ</sup>ci δwi* … <sup>0</sup> *<sup>γ</sup>ci δwi* … 0 ………… … …… …… …………… … … <sup>0</sup> � <sup>1</sup> *τci* þ *γci τci* … <sup>0</sup> � *<sup>γ</sup>ci τci* ………… ………… … …… …… …………… ………… … …… …… …………… � … *Thi*�<sup>1</sup> *Thi* … *Twi* … *Tci Tci*�<sup>1</sup> … þ *τh*1 0 0 0 0 0 0 0 0 <sup>1</sup> *τc*1 0 0 0 0 0 0 � *Thin Tcin* ½ � (26)

or:

$$\dot{\underline{T}}\_{H/w/c} = \underline{\mathbf{A}} \left( \tau\_{H/c}, \gamma\_{H/c}, \delta\_w \right) \cdot \underline{T}\_{H/w/c} + \underline{\mathbf{B}} \left( \tau\_{H/c\_1} \right) \cdot \underline{T}\_{H/c\_w} \tag{27}$$

#### **4.4 Laplace transformation**

The above-mentioned equation system can be mapped to:

$$T\_{hi+1} = \frac{\left(\frac{1}{\tau\_h} - \frac{1}{\delta\_h}\right)}{s + \left(\frac{1}{\tau\_h} + \frac{1}{\delta\_h}\right)} T\_{hi} + \frac{\frac{2}{\delta\_h}}{s + \left(\frac{1}{\tau\_h} + \frac{1}{\delta\_h}\right)} T\_{wi} = G\_{h1} T\_{hi} + G\_{h2} T\_{wi} \tag{28}$$

$$T\_{c\bar{t}-1} = \frac{\left(\frac{1}{\tau\_{\epsilon}} - \frac{1}{\delta\_{\epsilon}}\right)}{s + \left(\frac{1}{\tau\_{\epsilon}} + \frac{1}{\delta\_{\epsilon}}\right)} T\_{c\bar{t}} + \frac{\frac{2}{\delta\_{\epsilon}}}{s + \left(\frac{1}{\tau\_{\epsilon}} + \frac{1}{\delta\_{\epsilon}}\right)} T\_{wi} = \mathbf{G}\_{\epsilon 1} T\_{c\bar{t}} + \mathbf{G}\_{\epsilon 2} T\_{wi} \tag{29}$$

$$T\_{wi} = \sum\_{j=1}^{n} \frac{\chi\_h}{s + \Im(\chi\_h + \chi\_c)} T\_{hj} + \sum\_{j=1}^{n} \frac{\chi\_c}{s + \Im(\chi\_h + \chi\_c)} T\_{cj} \tag{30}$$

to obtain a simplified record:

$$T\_{cn} = \prod\_{j=1}^{n} G\_h T\_{hj} + \prod\_{j=1}^{n} G\_c T\_{cj} \tag{31}$$

From what we see, we get an nth-order system that depends on changes in input temperatures as well as material, geometric, and flow parameters.

#### **4.5 Constitution equations**

The material properties are generally described in the literature [45] and will not be further stated, but the determination of the heat transfer coefficient is more important. When using no dimensional numbers and links between them:

$$Nu = CPr^m \operatorname{Re}^n \tag{32}$$

And for our case, Reynolds' number of

$$Re = \frac{\dot{m}}{\sigma \nu} \tag{33}$$

Nusselt number:

$$Nu = \frac{aA}{a\lambda} \tag{34}$$

The coefficients depend on the shape of the channel and the Reynolds number. For our example, the coefficients are:

$$\mathbf{C} = \mathbf{0}, \mathbf{2}; \mathbf{m} = \mathbf{0}, \mathbf{67}; \mathbf{n} = \mathbf{0}, \mathbf{4} \tag{35}$$

Valid for Reynolds numbers:

$$Re > 100\tag{36}$$

## **5. Experimental model**

Although most of the models of heat exchangers are on the market with a fish bone sample, it was still selected to test the heat exchanger with a microplate sample. The reason is that heat exchangers with a microplate are more efficient because the flow between the plates is more subdued and the distribution by plate is generally better, which gives us about a 20–30% lower pressure drop or a much better heat transfer. The comparison between the plates is shown in **Figure 5**.

The test was performed on a heat exchanger with three plates in different configurations. Four cases were tasted with different flow configurations and speeds as per **Table 1**. Although the test was carried out on only three channels, the multi plates' behavior is similar. The geometry of the heat exchanger is shown in **Figure 6**.

For this purpose, a device was designed to allow for a rapid switch between hot and cold water. The heat exchanger was flushed with cold water and, at the moment of disturbance, a stream of hot water was flushed in. The switching was made with magnetic valves. Despite the fast-switching magnetic valves, the disturbance was not stepped but was approximately similar to the second-order response. To record this phenomenon, the entry signal was recorded, simulated and compared to the real state. Comparisons are shown in **Figures 7**–**10**.

**Figure 5.** *Comparison between (b) microplate and (a) fishbone heat exchanger.*


**Table 1.** *Test design.*

In **Figure 7**�**10**, input disturbances are simulated, while in **Figure 11**, deviations are depicted.

The discrepancies between the approximated and the measured disturbance are shown in **Figure 11**.

A FLIR A600 IR thermal camera was used to measure the dynamic temperature change. Flow measurements were carried out through weighing. Since the thermal camera can take a large number of images per second of 4000 images per second, it was assumed that the images are static at 0.1 s, which virtually eliminates the dynamic error of measurement. The temperature measurement error is +/�0,4°C.

In addition to the temperatures on the outside of the plate, the thermal camera also visualized the temperature field. The change in the temperature field is shown in **Figures 12**–**14** as a comparison between type A and other types.

The results show us that in all cases the distribution of temperature is uneven, which is probably due to the flow conditions on the panel itself. Despite the use of microchannel technology, the apparent distribution of water across the heat exchanger plate is still uneven, or not ideal for heat transfer.

**Figure 6.** *Heat exchanger geometry.*

**Figure 7.** *Disturbance A.*

**Figure 8.** *Disturbance B.*

**Figure 9.** *Disturbance C.*

**Figure 10.** *Disturbance D.*

**Figure 11.** *Disturbance error between approximation and real one.*

**Figure 12.** *Case A and B.*

The main observations are:


#### **Figure 14.** *Case A and D.*

The results of the dynamic time response tests of the heat exchanger are shown in **Figure 15**.

One can conclude from **Figures 12** and **13**, that case B is not relevant for temperature front propagation. Comparing then cases A, C, and D one can see that speed of water is main contributor to temperature front propagation. As in A-D cases hot is on outside border, then there is very little difference between A and C, while D with half speed, is noticeable slower than case A.

## **6. Numeric experiment**

The solution to the eq. (27) considering the geometry of the heat exchanger in **Figure 16** (Heat exchanger was divided into 10 segments — **Table 2** indicates a shift from the starting position) and the specifics of the experiment (input disturbance does not step)) is shown in **Figure 17**�**20**. The equation system was solved with the MATLAB Simulink software package using the RK4 numerical method.


**Table 2.**

*Segments position.*

**Figure 15.** *Experimental results.*

**Figure 16.** *Heat exchanger sections.*

**Figure 17**�**20** show simulation results split into sections. One can see the development of temperature as a function of time. The last 10th curve represents result of heat exchanger outlet.

**Figure 17.** *Simulation A.*

**Figure 18.** *Simulation B.*

**Figure 19.** *Simulation C.*

**Figure 20.** *Simulation D.*

## **7. Comparison between numerical calculation and tests results obtained and discussion**

A comparison between the measured and simulated results is shown in **Figure 21**�**24**.

The deviation between the measured and simulated results is relatively small and is less than 10%. The only difference is case D in the initial state because the heat exchanger was flushed with cold water and was not at ambient air temperature, as assumed in the simulation. Therefore, the initial situation is different. The deviations are shown in **Figure 25**.

## **8. Conclusion**

Obtaining digital twin from plate heat exchanger is relatively easy when behavior of the heat transfer coefficients is known. It is necessary to sectionize heat exchanger to get lump sum model, which is easy to resolve. A number of sections determine the accuracy, but also the distribution within the cell. The model can be further simplified and improved.

**Figure 21.** *Comparison A.*

**Figure 22.** *Comparison B.*

**Figure 23.** *Comparison C.*

**Figure 24.** *Comparison D.*

The mathematical model can be simplified by ignoring the intermediate wall, complete cell mixing, and independence of material properties from temperature.

However, it can be improved by increasing the number of cells to improve the approximation of the real state or by considering the logarithmic distribution of the temperature in the cell.

## **Conflict of interest**

"The authors declare no conflict of interest."

## **Nomenclature**



## **Author details**

Miha Bobič Danfoss Trata d.o.o, Ljubljana, Slovenia

\*Address all correspondence to: miha.bobic@danfoss.com

© 2022 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## **References**

[1] Abu-Khader MM. Plate heat exchangers: Recent advances. Renewable an Sustainable Energy Reviews. 2012;**16**: 1883-1891

[2] Aslam Bhutta MM, Hayat N, Bashir MH, Khan AR, Ahmad KN, Khan S. CFD applications in various heat exchangers design: A review. Applied Thermal Engineering. 2012;**32**:1-12

[3] Zhang J, Zhu X, Mondejar ME, Haglind F. A review of heat transfer enhancement techniques in plate exchangers. Renewable and Sustainable Energy Reviews. 2019;**101**:305-328

[4] Gullapalli VS, Sundén B. CFD simulation of heat transfer and pressure drop in compact brazed plate heat exchangers. Heat Transfer Engineering. 2014;**35**(4):358-366

[5] Freund S, Kabelac S. Investigation of local heat transfer coefficients in plate heat exchangers with temperature oscillation IR thermography and CFD. International Journal of Heat and Mass Transfer. 2010;**53**(19–20): 3764-3781

[6] Ayub ZH, Khan TS, Salam S, Nawaz K, Ayub AH, Khan MS. Literature survey and universal evaporation correlation for plate type heat exchangers. International Journal of Refrigeration. 2019;**99**: 408-418

[7] Guo Y, Wang F, Jia M, Zhang S. Modeling of plate heat exchanger based on sensitivity analysis and model updating. Chemical Engineering Research and Design. 2018;**138**:418-432

[8] Bobič M, Ogorevc T, Mazej M, Petkovšek J, Pribošek J, Sevčnikar J, et al. Vizualizacija delovanja ploščnega prenosnika toplote, Akademija

strojništva 2014 in 4. In: Mednarodna konferenca strojnih inženirjev. Ljubljana, Slovenija: ZSIS; 2014

[9] Pribošek J, Bobič M, Golobič I, Diaci J. Correcting the periodic optical distortion for particle-tracking velocimetry in corrugated-plate heat exchangers. Strojniški Vestnik. 2016;**62**(1):3-10

[10] Jin S, Hrnjak P. A new method to simultaneously measure local heat transfer and visualize flow boiling in plate heat exchanger. International Journal of Heat and Mass Transfer. 2017; **113**:635-646

[11] Kabelac S. Analysis of local heat transfer in plate heat exchangers for flow pattern optimization. In: 14th International Heat Transfer Conference. Vol. 4. ASME; 2010. pp. 375-384

[12] Jin S, Hrnjak P. Effect of end plates on heat transfer of plate heat exchanger. International Journal of Heat and Mass Transfer. 2017;**108**:740-748

[13] Clarke DD, Vasquez VR, Whiting WB, Greiner M. Sensitivity and uncertainty analysis of heat exchanger designs to physical properties estimation. Applied Thermal Engineering. 2001;**21**:993-1017

[14] Thorsen JE, Iversen J. Impact of lowering dt for heat exchangers used in district heating systems. In: DHC 13 the 13th International Symposium on District Heating and Cooling. Copenhagen, Denmark; 3rd-4th 2012

[15] Thorsen JE, Gudmundsson O, Brand M. Performance Specifications for Heat Exchanger for District Heating Substations of the Future, the 14th International Symposium on District Heating and Cooling, 7th–9th Sep.

Stockholm, Sweden: Danish Technical University; 2014

[16] Ansari MS, Mishra P, Gaikwad K, Awachat PN. Heat transfer analysis of corrugated plate heat exchanger over coil type heat exchanger: A review. International Journal of Engineering Research and Technology. 2014;**3**(3): 446-449

[17] Carlomagno GM, Cardone G. Infrared thermography for convective heat transfer measurements. Experiments in Fluids. 2010;**49**:1187-1218

[18] Gustafsson J, Delsing J, Van Deventer J. Improved district heating substation efficiency with a new control strategy. Applied Energy. 2010;**87**:1996- 2004

[19] Huang C, Cai W, Wang Y, Liu Y, Li Q, Li B. Review on the characteristics of flow and heat transfer in printed circuit heat exchangers. Applied Thermal Engineering. 2019;**153**:190-205

[20] Laszczyk P. Simplified modeling of liquid-liquid heat exchangers for use in control systems. Applied Thermal Engineering. 2017;**119**:140-155

[21] Piper M, Zibart A, Djakow E, Springer R, Homberg W, Kenig EY. Heat transfer enhancement in pillow-plate heat exchangers with dimpled surfaces: A numerical study. Applied Thermal Engineering. 2019;**153**:142-146

[22] Qian Z, Wang Q, Cheng J. Analysis of heat and resistance performance of fin-and-tube heat exchanger with rectangle-winglet vortex generator. International Journal of Heat and Mass Transfer. 2018;**124**:1198-1211

[23] Sharifi F, Glokar Narandji MR, Mehravaran K. Dynamic simulation of plate heat exchangers. International

Communications in Heat and Mass Transfer. 1995;**22**:213-225

[24] Wang L, Sundén B, Manglik RM. Plate Heat Exchangers: Design, Applications and Performance. Southampton: WIT Press; 2007

[25] Whalley R, Ebrahimi KM. Heat exchanger dynamic analysis. Applied Mathematica Modeling. 2018;**62**:38-50

[26] Arsenyeva OP, Kapustenko PO, Tovazhnyanskyy LL, Khavin G. The influence of plate corrugations geometry on plate heat exchanger performance in specified process conditions. Energy. 2016;**57**:201-207

[27] Dović D, Palm B, Svaić S. Generalised correlations for predicting heat transfer and pressure drop in plate heat exchanger channels of arbitrary geometry. International Journal of Heat and Mass Transfer. 2009;**52**(19–20):4553-4563

[28] Muley A, Manglik RM. Experimental study of turbulent flow heat transfer and pressure drop in a plate heat exchanger with Chevron plates. Journal of Heat Transfer. 1999;**121**(1):110-117

[29] Arsenyeva OP, Tovazhnyanskyy LL, Kaputenko PO, Demirskiy OV. Generalized semi empirical correlation for heat transfer in channels of plate heat exchanger. Applied Thermal Engineering. 2014;**70**(2):1208-1215

[30] Khan TS, Khan MS, Chyu MC, Ayub ZH. Experimental investigation of single-phase convective heat transfer coefficient in a corrugated plate heat exchanger for multiple plate configurations. Applied Thermal Engineering. 2010;**30**(8–9):1058-1065

[31] Sidebotham G. Nusselt number correlations. In: Heat Transfer Modeling. Springer-Verlag; 2015

[32] Aprea C, Renno C. Experimental analysis of a transfer function for an aircooled evaporator. Applied Thermal Engineering. 2001;**21**:481-493

[33] Maddali RK. Modeling ordinary differential equations in Mathlab Simulink. Indian Journal of Computer Science and Engineering. 2012;**3**(3): 406-410

[34] Badami M, Fonti A, Carpignano A, Grosso D. Design of district heating networks through an integrated thermofluid dynamics and reliability modelling approach. Energy. 2018;**144**(1):826-838

[35] Brand M, Thorsen JE, Svendsen S. Numerical modelling and experimental measurements for a low- temperature district heating substation for instantaneous preparation of DHW with respect to service pipes. Energy. 2012;**41** (1):392-400

[36] Das SK, Roetzel W. Dynamic analysis of plate heat exchangers with dispersion in both fluids. International Journal of Heat and Mass Transfer. 1995; **38**(6):1127-1140

[37] Hamze S, Witrant E, Bresch-Pietri D, Fauvel C. Estimating heat-transport and time-delays in a heat exchanger. In: 2018 IEEE Conference on Control Technology and Applications (CCTA). Copenhagen: Danish Technical University; 2018. pp. 1514-1519

[38] Michel A, Kugi A. Accurate loworder dynamic model of a compact plate heat exchanger. International Journal of Heat and Mass Transfer. 2013;**61**(33): 323-331

[39] Michel A, Kugi A. Model based control of compact heat exchangers independent of the heat transfer behavior II. Journal of Process Control. 2014;**50**(11–12):1859-1868

[40] Bobič M. Development of Measuring System of Coupled Pressure-Temperature Control System[Thesis]. Ljubljana, Slovenia: University of Ljubljana, Faculty for Mechanical Engineering; 2021

[41] Bobič M, Gjerek B, Golobič I, Bajsić I. Dynamic behaviour of a plate heat exchanger: Influence of temperature disturbances and flow configurations. International Journal of Heat and Mass Transfer. 2020;**163**:120439

[42] Bobič M. Dinamika preizkuševališča za razvoj in simuliranje delovanja toplotne podpostaje [magistrsko delo], mentor: Bajsić I., Ljubljana. Ljubljana, Magistrska dela: Fakulteta za strojništvo; 2000. p. 1079

[43] Bobič M., Muhič D., Zupnčič M., Sedmak I., Golobič I.: Dinamika ploščnega prenosnika toplote ogrevalnega sistema, Akademija strojništva 2017: inženirstvo - za kakovostnejše življenje, 6 mednarodna konferenca strojnih inženirjev 2017, Ljubljana, oktober 2017. Zveza strojnih inženirjev Slovenije, 2017. Vol. 6, No. 3/ 4, 78, Svet strojništva.

[44] Miao Q, You S, Zheng W, Zheng X, Zhang H, Wang Y. A Grey-box dynamic model of plate heat exchangers used in an urban heating system. Energies. 1398; **10**(9):2019

[45] Chapman AJ. Fundamentals of heat transfer. MacMilian Publishing Company; 1987

Section 4
