Abstract

Plate heat exchangers (PHEs) play an important role in different technical fields, namely, in energetics, chemical industry, food industry, and others. To use PHE effectively, it is necessary to have correct data for pressure drop. Unfortunately, in open literature, a large difference among different authors occurs. In this work is shown that an essential portion of this difference lies in the choice of the typical length for the calculation of the friction coefficient. Care must be taken to consider the pressure drop of the distribution zone. A three-component model for hydraulic resistance of PHE in turbulent flow regime is proposed in this work. The proposed model shows good agreement with experimental data.

1. State of Research

Plate heat exchangers (PHE) are widely used in industry. A compact size due to high heat transfer coefficients is one major advantage over other heat exchanger types like shell and tube type. The special design of the flow channels results in early transition from laminar to turbulent flow and significant increase in heat transfer coefficients. Unfortunately, this comes along with an increase in pressure drop compared to flat channels, too. A reliable model for pressure drop calculation is required for a successful use of PHE.

The following model is traditionally applied to present experimental data of a PHE:where is the pressure drop in a channel of a PHE, is the average velocity in a channel, is the effective length of the PHE, is the equivalent diameter of a channel, and is the Fanning friction factor or Fanning factor in a channel.

In literature, the model according to equation (1) is widespread [119]. The Fanning friction factor is a function of Reynolds number and is presented in graphical or analytical form. If Fanning factors were found for a specific unit, this model can reproduce the specific unit’s pressure drop within 10% to 20% margin. Applied to a PHE with different geometries, the margin of deviation can be much higher. Fanning factors calculated with models of different authors can vary by up to 50–100%. Therefore, Fanning factors found for one plate type are mostly valid only for its specific design. Applications of equation (1) to various PHEs are stated in textbooks of Wang et al. [20] and Klemes et al. [21]. Some results cannot be explained well according to equation (1). The study [2] is a good example: it shows different correlations for refrigerants R 245fa and R 236a in turbulent flow regime for the same PHE. This means that equation (1) cannot always consider a strong influence of physical properties correctly.

Two citations of recent publications summarize the state of the art very well:(i)“Most of the previously reported empirical correlations are restricted to specific experimental geometry and flow conditions and cannot be generalized for all PHE applications” [9].(ii)“It can be seen from the experimental studies…, a general correlation applicable to all chevron angle even with the same chevron angle is not possible. A gasketed plate heat exchanger design continues to be proprietary in nature; the correlations obtained are applicable to this exchanger that is marketed” [12].

Several attempts of creating a generalized theory were made in the past. Especially, Dovic et al. [7], Martin [14], Muley and Manglik [16], Klemes et al. [21], and Arsenyeva et al. [22] deserve attention. The most physically reasonable approach is the three-component model proposed by Arsenyeva et al. [23].

This work is dedicated to study the influence of reference length L on the pressure drop of a PHE. At first, it is shown that usage of the model according to equation (1) results in Fanning factors that always depend on effective length L. Subsequently, the three-component model proposed by Arsenyeva et al. [23] is enhanced. It divides the channel in the two distribution zones—at the inlet and outlet of a channel—and the corrugation field in the middle of a channel. Each of the three zones of the channel has its own friction factor correlation. Furthermore, it is shown that the corrugation field correlation has universal characteristics and is independent of its length in turbulent regime.

2. Experimental Setup and Procedure

Experimental setup involves three PHEs (warm-side setup, test set, and cold-side set-up), arranged in two closed loops, as shown in Figure 1.

On the warm side, a pump circulates the warm fluid (water or ISO VG46 Oil). The warm fluid is heated by using a boiler. On the cold side, three air coolers maintain a constant temperature of the fluid (water or glycol). A fully automatic control unit adjusts the temperature and flow rates to comply with the set values. As test PHE were used four plate heat exchangers of FUNKE Wärmeaustauscher Apparatebau GmbH (Germany) [24]. One of these plates is shown in Figure 2.

The geometric parameters of a plate are illustrated in Figure 3, and the geometric values for plates P1–P4 are shown in Table 1.

Every device could have one of the three chevron patterns: HH (with chevron angle φ = 63°), HL (with chevron angle φ = 63°/27°), and LL (with chevron angle φ = 27°).

All four plates are similar, and only the length of the corrugated field varies. In Figure 4, the test set with test PHE is shown.

Platinum resistance thermometers Pt1000 (T1–T4) are used to measure the fluid temperatures at the inlet/outlet of the test PHE with the measurement error ±(0.15 K + 0.002 t). To improve the precision of the sensors, we calibrate all temperature sensors simultaneously at the area of working temperature. The deviation of measured temperature of all sensors is about 0.05–0.1 K. The flow rates in the test section are measured by volumetric inductive flow meters of type KROHNE OPTIFLUX 4300 (MF1 and MF2) with a measurement error of ±1% for measuring range 3–500 m³/h. Static mixers are used to reduce local temperature variations. Static mixers (2 and 4) at outlets are basically used, mixers (1 and 3) at inlets if necessary. The pressure losses of the fluid flowing through the test section is measured with a differential pressure transmitter model DPT10 (PD1 and PD2) with the measurement error of ±0.2%. The two devices are used for measuring ranges 0–500 mbar and 0–3 bar. Pressure at the inlet of PHE is measured with the pressure transmitter (P1 and P2) with the measurement error of ±0.2%.

The values are measured as close as possible to the test PHE, taking care to use calibrated measurement pipes (length/diameter ≥ 10). The test PHE was configurated as a single-pass U-type counterflow device. The corrugation profile of the plate is given in Figure 5. It is close to sinusoidal.

All data were obtained under steady-state conditions with close energy balance (less than ±5% difference between the measured performance of warm and cold side). A constant inlet temperature (≈60°C) for the hot fluid was maintained throughout each set of experiments. The inlet temperature of the cold fluid was constant (≈40°C), too. All sensors and measuring equipment are calibrated accordingly to standard ISO 9001 : 2008.

3. Data Reduction

The primary measurements of each fluid consist of the flow rate, inlet and outlet temperatures, pressure drop between the inlet and outlet of the PHE, and inlet pressure for every side. The channel flow Reynolds number is based on the equivalent diameter De:

All fluid properties were calculated at the bulk-mean temperature given by

From the measured overall pressure drop across the PHE , the pressure drop in the channel can be obtained from the following equation:

For the port, pressure drop according to [17] is used:

The pressure loss in inlet and outlet pipe was estimated on the basis of the smooth tube friction factor, and local pressure losses on the sudden area changes. The pressure loss considers the distance of pressure sensors to PHE.

The Fanning friction factor was calculated as

The choice of effective length of the PHE L is discussed later.

The experimental uncertainties were determined by the procedure outlined by [25]. The maximum errors in the primary measurements of flow rate , pressure drop Δp, and temperature T were ±1%, ±0.2%, and ±(0.15 K + 0.002 T), respectively. Based on these measurement uncertainties and the propagation of error analysis, the uncertainties in Reynolds number Re and Fanning friction factor f were 5.7% and 9.2%, respectively.

4. Experimental Results

For the valid presentation of the data for the pressure losses in the channel of PHEs, the choice of the characteristic length L of the PHEs is important. In common use is the vertical distance between holes in the plate Lport [14, 16]:

Some authors use the relation between the area of the plate and its width [3, 5]:

The authors [10] consider the choice:

For comparison, all our data and data of other authors have been processed with a characteristic length according to equation (7). Results for plate P2 are presented in Figure 6. For each chevron pattern, three zones can be identified graphically [16]:(1)Laminar zone without vortices, Re < 100(2)Vortices zone, 100 < Re < 500(3)Turbulent zone, Re > 200–500

By increase in chevron angle φ, the hydraulic resistance in the channel increases and transition to turbulence is displaced to smaller Reynolds numbers.

Figure 7 illustrates experimental data in the turbulent zone for all four PHE P1–P4 with chevron pattern HH and data of other authors.

In the turbulent zone, it is obvious that data depend on plate length Lport: the higher the length is, the higher the Fanning friction factor is. The shortest plate P1 has the smallest friction factor. The dependence of friction factor fch(Lport) on the relative length of a plate Lport/Bch is confirmed by comparison of our experimental data with literature, too. Our data for the long plates P3 (Lport/Bch = 3.69) and P4 (Lport/Bch = 4.73) agree well with the data of Bond [4], Heavner [12], and Thonon [19], who investigated industrial PHEs. Data for our shorter plates P1 (Lport/Bch = 2.06) and P2 (Lport/Bch = 2.86) are considerably lower.

The same characteristic behavior can be observed for data of other authors. Fanning factor f is the highest in the publication of Focke et al. [9] who used a plate with Lport/Bch = 4.0. Data for a plate with Lport/Bch = 3.57 from Shaji and Das [18] are lower and agree well with our data for plate P3 (Lport/Bch = 3.69). Fanning factors of Lee et al. [13] for a plate with Lport/Bch = 2.94 are approximately the same as our data of plate P2 (Lport/B = 2.86). Muley and Maglik [16] investigated a plate with Lport/Bch = 2.4 and received considerably lower Fanning factors compared to all other mentioned authors.

In Figure 8, experimental data for plates P1–P4 with chevron pattern LL are compared to data of other authors. For corrugation angle φ = 27° (LL), the dependence on the length Lport changes in comparison with HH. For LL chevron pattern, Fanning factors decrease by increase in length Lport. Therefore, in Figure 8, the values of the shortest plate P1 are the biggest ones and values for longest plate P4 are the smallest ones.

As before, our experimental data are consistent with the data of Bond [4], Heavner [12], and Thonon [19]. Generally, variation between the correlations has increased. A relatively high influence of distribution zone on channel pressure drop for LL pattern seems to be a probable cause of this behavior.

Data of Focke et al. [9] (L/B = 4.0) are now in the range of the lowest data. Fanning factors of Shaji and Das [18] for a plate with Lport/Bch = 2.4 are the second highest now. The proportionality of the Fanning factor and plate length is not as expected for data of Shaji and Das [18]. This may be caused by different geometries of the distribution zones.

5. Three-Component Model

The typical approach (equation (6)) to make a connection between hydraulic resistance and characteristic length of the plate leads to length-dependent Fanning friction factors. The resulting correlations are of opposite behavior for low and high (LL and HH) chevron angles. It is caused by a different character of the flow in all three characteristic zones—inlet and outlet distribution zone and corrugated field. Due to this, the idea of Arsenyeva et al. [23] to divide a plate in three areas and describe the friction factor independently for each area seems reasonable.

A main challenge of a three-component model is the complexity of measurements for all components. A calculation-based approach to determine the different parts of the pressure drop is an alternative to measurements. This is possible if experimental data for one series of plates (equal corrugation type, distribution zone, and width) with different lengths of the corrugated field are available.

From experimental data, the pressure drop of a channel ∆pch is received. At first, data for a low corrugation angle (LL) are used. The turbulent flow region is initially considered because of very limited amount of data in laminar region. A dimensionless form of the channel pressure drop is introduced:

The channel pressure drop ∆pch consist of three components in this model:with pressure drop of the corrugated field ∆pch and pressure drop of the inlet ∆pinlet and outlet ∆poutlet distribution zone. Each component is described as follows:

The characteristic velocity of the channel is used for all three zones. This makes the data handling easier, and comparability of results between distribution zones and corrugated field is improved. A main disadvantage is that velocity is not the physical velocity in distribution zones. The cross section in distribution zones differs from the one of the corrugated field and therefore velocity, too. As a result, friction factor finlet has only computational meaning.

From equations (10)–(13), we get

For PHEs P1–P4, four functions are received from experimental data:

Characteristic values for Reynolds number in the range Re = 1000,…, 10000 are chosen and called . Values are according to equations (14) and (15), linear depending from . Using known mathematical procedure, we receive linear approximation for four points for all characteristic values of . The resulting y-intercept of this linear function is while is the slope. The obtained values are presented in Table 2.

By approximating the data from Table 2, we obtain

The deviation of equations (16) and (17) is about 5% from data of Table 2. With equation (14), dimensionless pressure drop Fch becomes the form

The dimensionless pressure drop can be rearranged to Fanning friction factor of the channel (equations (6) and (14)). We get

According to equation (19), Fanning factor of the channel fch is related to characteristic length Lcorr. Earlier in Figures 68, the Fanning factor was based on characteristic length Lport. By using equations (16), (17), and (19) for plates P1–P4 with LL corrugated field, we receive

Calculated values from equation (20) are drawn as lines in Figure 9, and they are labeled with “Theory.” As can be seen, calculated data agree well with experimental data. The deviation is within ±5%. It must be emphasized that both equations (16) and (17) are only valid for specific geometry of distribution zone and corrugated field of plates P1–P4.

Distribution zones are equal for all patterns of corrugated fields. Due to this fact, equation (16) is valid for HH and HL channels, too. Based on experimental data for pressure drop ∆pch of HH and HL channels, the Fanning factor can be calculated from equations (11)–(13) for corresponding channel type.

To receive fcorr (HH), at first, dimensionless pressure drop is calculated for the HH corrugated field:

This allows us to calculate the Fanning friction factor of the channel for characteristic length Lcorr:

With known from equation (16), we get

Values calculated from equation (23) result in a function depending on the Reynolds number. Following approximation is obtained from experimental data of the HH channel type:

The Fanning factor of the channel from equations (16) and (24) is

Calculated values of from equation (25)—labeled as “Theory”—and experimental data according to equation (22) are compared in Figure 10.

For the HH channel type, data show only a slight dependency from the length of the corrugated field. This is due to a relatively low hydraulic resistance of distribution zones compared to HH corrugated field.

All calculations according to equations (21)–(25) are performed for HL channel type, too. We get

Calculated values of from equation (27)—labeled as “Theory”—and experimental data are compared in Figure 11.

Experimental data agree with theoretical values within ±10%, and these are compatible with the uncertainty of the measurements.

The three-component model gives the opportunity to use data sets for twelve configurations (four plate lengths with three corrugation angles each) and describes with only four correlations for friction factors: , , , and . It is important to note that these correlations are length independent.

In Figure 12, all experimental data are compared with theoretical values for of the three-component model. For this purpose, experimental data for all channels were recalculated to experimental data for the corrugated field according to the following equation:

It can be seen that data from experiments of all four plates correspond to the theoretical calculations. Experimentally is confirmed that the factor does not depend on the length of the corrugated field .

5.1. Dependence of Fanning Friction Factor from Plate Length

Based on the proposed three-component model, a simple physical explanation of observed opposite behavior of Fanning friction factor for HH and LL channels can be obtained. It is assumed that friction factor is linear dependent from distribution zone length:where is a unique relative pressure drop coefficient for each series of plate with equal geometry of distribution zone. It describes the specific hydraulic resistance of this zone. We derive from equation (19)

We introduce the relative length of the corrugated fieldand get

By increase in plate length Apl, both Lport and Lcorr increase; however, length of distribution zone remains unchanged:

With equation (31), we get,

As a result, if Lport increases, X approaches one (X ⟶ 1) and from equation (32),

For plates with very short corrugated field, X approaches zero (X ⟶ 0) and

The characteristic dependence of from relative length of corrugated field X is shown in Figure 13. As a general rule, the Fanning friction factor of HH corrugated field is bigger than relative pressure drop of the inlet , while the friction factor of LL corrugated fields is smaller:

For example, for PHE P1–P4 at Re = 1000, the friction factor is calculated as follows:

According to equation (35), an increase in plate length for the HH channel type should result in an increase in friction factor that approaches . LL corrugated field has a relatively low pressure drop according to equation (37); therefore, with increase in plate length, should decrease and approach . Our experimental data confirm this behavior.

5.2. Pressure Drop Fraction in Inlet and Outlet Distribution Zones

It is important to know what part of overall channel pressure drop corresponds to distribution zone. We introduce the inlet pressure drop ratio :

With equations (10) and (12), can be written as

Equation (41) shows that always . The higher is, the higher the portion of pressure drop in distribution zones is. For plates P1–P4, Figures 1416 show calculated values for for all three channel types.

The inlet pressure drop ratio is highest in shortest plate P1 with LL corrugation angle. It is between 66% and 77%. Plates P2–P4 show values from 39% up to 61% as it can be seen in Figure 14. Plates with HH corrugated field shows values from 12% to 17% for shortest plate P1 (Figure 16). Longer plates have an inlet pressure drop ratio from 4% up to 11%. The results of intermediate corrugation angle HL lay in-between HH and LL channels. Values for from 28%–40% for plate P1 and 11%–29% for plates P2–P4 are calculated.

The calculated values of inlet pressure drop ratio show that the distribution zones play an important role in channel pressure drop for chevron angles of 30°–45° and lower.

5.3. Comparison of the Three-Component Model to Existing Data

Our approximation for the pressure drop factor of distribution zone (equation (16)) is received by investigation of plates with particular distribution zone. Arsenyeva et al. [23] show the following equation:with friction coefficient fDZ calculated by

This equation is also obtained for a particular plate design. Therefore, a full quantitative match with our data can not be expected. Nevertheless, a comparison is of interest.

By comparing equations (12) and (42), it has to be noted that

For Re = 2700, we get from equation (43), fDZ = 38 and according to equation (44), finlet = 38/4 = 9.5 and from equation (16), finlet = 5.525 + 7500/2700 = 8.3.

The difference between equations (43) and (16) is about 13%; this is within the range of experimental uncertainties. Figure 17 shows calculated values of equations (43), (44), and (16). In the range Re = 1000, …, 40000, values from both equations agree well.

6. Conclusions

(i)Pilot study of pressure losses of plate heat exchangers is conducted. The calculated Fanning friction factor is consistent with data from open literature.(ii)Experimentally, the dependence of Fanning friction factor on distance between ports Lport is determined. In case of corrugation pattern HH (chevron angle φ > 60°), Fanning friction factor increases with increase in length Lport; in case of corrugation pattern LL (φ < 30°) in contrast, friction factor decreases with increase in length Lport. This dependence is confirmed by the data published in open literature.(iii)For pressure losses in PHE, the enhanced three-component model is offered. This model accounts for the hydraulic resistance of distribution zones and the correlations for corrugated field. It corresponds to experimental data with an error up to 10%.(iv)On the basis of three-component model is explained why at chevron angles φ > 60° Fanning friction factor increases with increase in length Lport, and in case of φ < 30° in contrast, friction factor decreases with increase in length Lport.(v)The portion of hydraulic resistance of distribution zone in general resistance of channel is calculated. This portion can be considerably (up to 70%) in case of corrugation angles φ < 30°. For corrugation angles of φ > 60°, it can reach 10–15%.(vi)Comparison of the received dependence of hydraulic resistance of a distribution zone on the Reynolds number with the known data from literature is carried out. Good correspondence is obtained.

Nomenclature

A:Area of PHE (m2)
B:Width (m)
D:Port diameter (m)
De:Equivalent diameter = 2t (m)
Dh:Hydraulic diameter = 2t/FK (m)
F:Fanning friction factor, dimensionless
FK:Surface enlargement factor, dimensionless
L:Length (m)
Re:Reynolds number, , dimensionless
t:Wave amplitude (m)
T:Temperature, °C
IR:Inlet pressure drop ratio, dimensionless
w:Velocity (m/s)
X:Relative length of corrugated field, dimensionless
Δp:Pressure losses, Pa
Δwall:Wall thickness (m)
Λ:Wavelength of corrugation (m)
λ:Heat conduction coefficient (W/mK)
µ:Dynamic viscosity (kg/(m·s))
ρ:Density (kg/m³)
φ:Chevron corrugation angle relative to flow direction, degree
corr:Corrugation
ch:Channel
inlet:Inlet
mean:Mean
outlet:Outlet
pipe:Pipeline
pl:Plate
port:Port distance
proj:Projection
wall:At wall temperature.

Data Availability

The measurement data used to support the findings of this study have not been made available because these data are business secret of FUNKE Wärmeaustauscher Apparatebau GmbH.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors would like to thank our employer FUNKE Wärmeaustauscher Apparatebau GmbH for providing the opportunity to publish this article.