Boundary condition equations.

## 1. Introduction

Heat exchangers are devices used to transfer energy between two fluids at different temperatures. They improve energy efficiency, because the energy already within the system can be transferred to another part of the process, instead of just being pumped out and wasted. In the new era of sustainability, the growing urgency to save energy and reduce overall environmental impacts has placed greater emphasis on the use of heat exchangers with better thermal efficiency. In this new scenario, the plate heat exchanger can play an important role.

A plate heat exchanger is a compact type of heat exchanger that uses a series of thin plates to transfer heat between two fluids. There are four main types of PHE: gasketed, brazed, welded, and semi-welded. The plate-and-frame or gasketed plate heat exchanger essentially consists of a pack of thin rectangular plates sealed around the edges by gaskets and held together in a frame (Figure 1). Plate heat exchangers were first introduced in 1923 for milk pasteurization applications, but are now used in many applications in the chemical, petroleum, HVAC, refrigeration, dairy, pharmaceutical, beverage, liquid food and health care sectors. This is due to the unique advantages of PHEs, such as flexible thermal design (plates can be simply added or removed to meet different heat duty or processing requirements), ease of cleaning to maintain strict hygiene conditions, good temperature control (necessary in cryogenic applications), and better heat transfer performance.

## 2. Mechanical characteristics

A PHE consists of a pack of thin rectangular plates with portholes, through which two fluid streams flow, where heat transfer takes place. Other components are a frame plate (fixed plate), a pressure plate (movable plate), upper and lower bars and screws for compressing the pack of plates (Figure 2). An individual plate heat exchanger can hold up to 700 plates. When the package of plates is compressed, the holes in the corners of the plates form continuous tunnels or manifolds through which fluids pass, traversing the plate pack and exiting the equipment. The spaces between the thin heat exchanger plates form narrow channels that are alternately traversed by hot and cold fluids, and provide little resistance to heat transfer.

### 2.1. Thermal plates and gaskets

The most important and most expensive part of a PHE is its thermal plates, which are made of metal, metal alloy, or even special graphite materials, depending on the application. Stainless steel, titanium, nickel, aluminum, incoloy, hastelloy, monel, and tantalum are some examples commonly found in industrial applications. The plates may be flat, but in most applications have corrugations that exert a strong influence on the thermal-hydraulic performance of the device. Some of the main types of plates are shown in Figure 3, although the majority of modern PHEs employ chevron plate types. The channels formed between adjacent plates impose a swirling motion to the fluids, as can be seen in Figure 4. The chevron angle is reversed in adjacent sheets, so that when the plates are tightened, the corrugations provide numerous points of contact that support the equipment. The sealing of the plates is achieved by gaskets fitted at their ends. The gaskets are typically molded elastomers, selected based on their fluid compatibility and conditions of temperature and pressure. Multi-pass arrangements can be implemented, depending on the arrangement of the gaskets between the plates. Butyl or nitrile rubbers are the materials generally used in the manufacture of the gaskets.

### 2.2. Design characteristics

This section presents some of the main advantages and disadvantages of a PHE, compared to shell-and-tube heat exchangers.

**Advantages**

*Flexibility:* Simple disassembly enables the adaptation of PHEs to new process requirements by simply adding or removing plates, or rearranging the number of passes. Moreover, the variety of patterns of plate corrugations available, together with the possibility of using combinations of them in the same PHE, means that various conformations of the unit can be tested during optimization procedures.

*Good temperature control:* Due to the narrow channels formed between adjacent plates, only a small volume of fluid is contained in a PHE. The device therefore responds rapidly to changes in process conditions, with short lag times, so that the temperatures are readily controllable. This is important when high temperatures must be avoided. Furthermore, the shape of the channels reduces the possibility of stagnant zones (dead space) and areas of overheating.

*Low manufacturing cost:* As the plates are only pressed (or glued) together, rather than welded, PHE production can be relatively inexpensive. Special materials may be used to manufacture the plates in order to make them more resistant to corrosion and/or chemical reactions.

*Efficient heat transfer:* The corrugations of the plates and the small hydraulic diameter enhance the formation of turbulent flow, so that high rates of heat transfer can be obtained for the fluids. Consequently, up to 90% of the heat can be recovered, compared to only 50% in the case of shell-and-tube heat exchangers.

*Compactness:* The high thermal effectiveness of PHEs means that they have a very small footprint. For the same area of heat transfer, PHEs can often occupy 80% less floor space (sometimes 10 times less), compared to shell-and-tube heat exchangers (Figure 5).

*Reduced fouling*: Reduced fouling results from the combination of high turbulence and a short fluid residence time. The scale factors for PHEs can be up to ten times lower than for shell-and-tube heat exchangers.

*Ease of inspection and cleaning*: Since the PHE components can be separated, it is possible to clean and inspect all the parts that are exposed to fluids. This feature is essential in the food processing and pharmaceutical industries.

*Easy leak detection*: The gaskets have vents (Figure 6) that prevent fluids from mixing in the case of a failure, which also facilitate locating leaks.

**Drawbacks**

*Temperature and pressure limitations*: An important limitation of PHEs is related to the plate gaskets. Pressures and temperatures exceeding 25 atm and 160 °C, respectively, are not tolerated because they can cause the standard gaskets to leak. However, gaskets made of special materials can withstand temperatures up to 400 °C, and it is possible to weld or braze the plates to each other in order to operate under more severe conditions. This would have the additional advantages of increasing the operational limits, as well as the possibility of working with corrosive fluids, because it would eliminate the need for gaskets. However, the PHE would lose its major advantages of flexibility and ease of cleaning, and the equipment would become more expensive.

*High pressure drop*: Because of the corrugated plates and the small flow space between them, the pressure drop due to friction is high, which increases pumping costs. The pressure drop can be reduced by increasing the number of passages per pass and splitting the flow into a greater number of channels. This diminishes the flow velocity within the channel, hence reducing the friction factor. However, the convective heat transfer coefficient is also reduced, decreasing the effectiveness of the heat exchanger.

*Phase change*: In special cases, PHEs can be used in condensation or evaporation operations, but are not recommended for gases and vapors due to the limited space within the channels and pressure limitations.

*Types of fluids*: The processing of fluids that are highly viscous or contain fibrous material is not recommended because of the high associated pressure drop and flow distribution problems within the PHE. Compatibility between the fluid and the gasket material should also be considered. Highly flammable or toxic fluids must be avoided due to the possibility of leakage.

*Leakage:* Friction between the metal plates can cause wear and the formation of small holes that are difficult to locate. As a precaution, it is advisable to pressurize the process fluid so that there is less risk of contamination in the event of leakage from a plate.

### 2.3. Arrangement of a plate heat exchanger

The simplest types of arrangements of plate heat exchangers are those in which both fluids make just one pass, so there is no change in direction of the streams. These are known as 1-1 single-pass arrangements, and there are two types: countercurrent and concurrent. A great advantage of the single-pass arrangement is that the fluid inlets and outlets can be installed in the fixed plate, making it easy to open the equipment for maintenance and cleaning, without disturbing the pipework. This is the most widely used single-pass design, known as the U-arrangement. There is also a single-pass Z-arrangement, where there is input and output of fluids through both end plates (Figure 7).

Countercurrent flow, where the streams flow in opposite directions, is usually preferred due the achievement of higher thermal efficiency, compared to concurrent flow, where the streams flow in the same direction. Multi-pass arrangements can also be employed to enhance the heat transfer or flow velocity of the streams, and are usually required when there is a substantial difference between the flow rates of the streams (Figure 8).

There are five parameters that can be used to characterize the PHE configuration [6]:

Number of channels (

Number of passes (

Hot fluid location (

Feed connection (

The plates of a PHE can provide vertical or diagonal flow, depending on the arrangement of the gaskets. For vertical flow, the inlet and outlet of a given stream are located on the same side of the heat exchanger, whereas for diagonal flow they are on opposite sides. Assembly of the plate pack involves alternating between the “A” and “B” plates for the respective flows. Mounting of the plate pack in vertical flow mode only requires an appropriate gasket configuration, because the A and B arrangements are equivalent (they are rotated by 180°, as shown in Figure 10a). This is not possible in the case of diagonal flow, which requires both types of mounting plate (Figure 10b). To identify each type of flow, Gut (2003) considered the binary parameter

## 3. Mathematical modeling

Due to the large number of plate types and pass arrangements, there are many possible configurations of a particular PHE design. As a result, a number of mathematical modeling approaches have been proposed for the calculation of performance. Two different modeling approaches are described below.

### 3.1. Model 1

A mathematical model was developed to simulate the general configuration of a PHE operating under steady state conditions, characterized using six different parameters [6]. In this model, the parameters considered are the number of channels, the number of passes for each side, the fluid locations, the feed connection locations, and the type of channel flow. The following assumptions are made:

The PHE operates at steady state;

The main flow is divided equally among the channels that make up each pass;

The velocity profile in the channels is flat (plug flow);

Perfect mixture in the end of each pass;

There are no heat losses to the environment;

There are no phase change;

There is no heat transfer in the direction of flow, either in the fluids or in the plates, so heat transfer only occurs in the direction perpendicular to the flow;

The physical properties of the fluids remain constant throughout the process.

The last assumption listed above implies an overall heat transfer coefficient

Applying the energy conservation law to a given volume of control of a generic channel

where

and

This system of linear differential equations can be written in the matrix form as:

where

The boundary conditions, which are dependent on the PHE configuration, can be divided into three different categories:

*Fluid inlet temperature*: In the channels of the first pass, the fluid inlet temperature is the fluid feed temperature.

*Change of pass temperature*: The temperature at the beginning of the channels of a particular pass is equal to the arithmetic average of the temperatures in the channels of the previous pass.

*Fluid outlet temperature*: The outlet temperature of the fluid is the arithmetic average of the outlet temperatures of the channels of the last pass.

The analytical solution is given by Eq. (10), where

Application of Eq. (10) in the boundary condition equations for the fluid inlet and change of pass enables the creation of a linear system of

**Example:** Creation of the linear system of

In order to illustrate the generation of the linear system, a PHE containing 7 thermal plates (or 8 channels), with the cold fluid making two passes and the hot fluid making one pass, is shown in Figure 12.

Applying Eq. (10), the following analytical solution can be achieved:

Using the boundary condition equations (7) and (8) for all the channels of the PHE under investigation, the equations presented in Table 1 are generated.

Cold fluid | Hot fluid |

These equations can be written in the following way:

To achieve the matrix form is Eq. (10) is applied to the linear system:

where

where

### 3.2. Model 2

The assumption is made that any multi-pass PHE with a sufficiently large number of plates (so that end effects and inter-pass plates can be neglected) can be reduced to an arrangement consisting of assemblies of single-pass PHEs [11]. This enables the development of closed-form equations for effectiveness, as a function of the ratio between the heat capacities of the fluids and the number of transfer units, for the arrangements 1-1, 2-1, 2-2, 3-1, 3-2, 3-3, 4-1, 4-2, 4-3, and 4-4 (Table 2). In other words, most multi-pass plate heat exchangers can be represented by simple combinations of pure countercurrent and concurrent exchangers, so that a multi-pass PHE is therefore equivalent to combinations of smaller single-pass exchangers (Figure 13).

The assumptions considered are the same as in the first mathematical model. The derived formulas are only valid for PHEs with numbers of thermal plates sufficiently large that the end effects can be neglected. This condition can be satisfied, depending on the required degree of accuracy. For example, a minimum of 19 plates is recommended for an inaccuracy of up to 2.5% [12]. Elsewhere, a minimum of 40 thermal plates was used [11, 13]. In the formulas,

Formulas | Arrangements | ||

Arrangement 111where: | |||

Arrangement 211where: | |||

Arrangement 221where: | |||

Arrangement 311where: | |||

Arrangement 321where: | |||

Arrangement 331where: | |||

Arrangement 411where: | |||

Arrangement 421where: | |||

Arrangement 431 where: | |||

Arrangement 441where: |

## 4. Design of a plate heat exchanger

### 4.1. Basic equations for the design of a plate heat exchanger

The methodology employed for the design of a PHE is the same as for the design of a tubular heat exchanger. The equations given in the present chapter are appropriate for the chevron type plates that are used in most industrial applications.

#### 4.1.1. Parameters of a chevron plate

The main dimensions of a *chevron* plate are shown in Figure 14. The corrugation angle, *β*, usually varies between extremes of 25° and 65° and is largely responsible for the pressure drop and heat transfer in the channels.

The corrugations must be taken into account in calculating the total heat transfer area of a plate (effective heat transfer area):

where

The enlargement factor of the plate is the ratio between the plate effective heat transfer area,

For the effective heat transfer area, the hydraulic diameter of the channel is given by the equivalent diameter,

where

#### 4.1.2. Heat transfer in the plates

The heat transfer area is expressed as the global design equation:

where

where

where

The convective heat transfer coefficient,

#### 4.1.3. Design methods

There are two main approaches used in the design of PHEs, namely the log-mean temperature difference and the thermal effectiveness methods. For the first method, the rate of heat transfer is given by:

where

Where

The correction factor is a function of the heat exchanger configuration and the dimensionless parameters

when the PHE has only one thermal plate, so that only two channels are formed by the end plates and the thermal plate, with each stream flowing through one channel;

when the number of thermal plates is sufficiently large that the edge effect can be neglected.

The adimensional parameters

The second method provides a definition of heat exchanger effectiveness in terms of the ratio between the actual heat transfer and the maximum possible heat transfer, as shown in Eq. (25):

The actual heat transfer can be achieved by an energy balance:

Thermodynamically,

Using Eqs. (26), (27) and (28), the PHE effectiveness can be calculated as the ratio of temperatures:

The effectiveness depends on the PHE configuration, the heat capacity rate ratio (

#### 4.1.4. Pressure drop in a plate heat exchanger

The pressure drop is an important parameter that needs to be considered in the design and optimization of a plate heat exchanger. In any process, it should be kept as close as possible to the design value, with a tolerance range established according to the available pumping power. In a PHE, the pressure drop is the sum of three contributions:

Pressure drop across the channels of the corrugated plates.

Pressure drop due to the elevation change (due to gravity).

Pressure drop associated with the distribution ducts.

The pressure drop in the manifolds and ports should be kept as low as possible, because it is a waste of energy, has no influence on the heat transfer process, and can decrease the uniformity of the flow distribution in the channels. It is recommended to keep this loss lower than 10% of the available pressure drop, although in some cases it can exceed 30% [3].

where

The values for

#### 4.1.5. Experimental heat transfer and friction correlations for the chevron plate PHE

Due to the wide range of plate designs, there are various parameters and correlations available for calculations of heat transfer and pressure drop. Despite extensive research, there is still no generalized model. There are only certain specific correlations for features such as flow patterns, parameters of the plates, and fluid viscosity, with each correlation being limited to its application range. In this chapter, the correlation described in [14] was used.

where

In Reynolds number equation,

The constants

### 4.2. Optimization

Any industrial process, whether at the project level or at the operational level, has aspects that can be enhanced. In general, the optimization of an industrial process aims to increase profits and/or minimize costs. Heat exchangers are designed for different applications, so there can be multiple optimization criteria, such as minimum initial and operational costs, minimum volume or area of heat transfer, and minimum weight (important for space applications).

Heat transfer | Pressure drop | |||||

0.718 | 0.349 | 50.000 | 1.000 | |||

> 10 | 0.348 | 0.663 | 10 - 100 | 19.400 | 0.589 | |

> 100 | 2.990 | 0.183 | ||||

45 ° | < 10 | 0.718 | 0.349 | < 15 | 47.000 | 1.000 |

10 - 100 | 0.400 | 0.598 | 15 - 300 | 18.290 | 0.652 | |

> 100 | 0.300 | 0.663 | > 300 | 1.441 | 0.206 | |

50 ° | < 20 | 0.630 | 0.333 | < 20 | 34.000 | 1.000 |

20 - 300 | 0.291 | 0.591 | 20 - 300 | 11.250 | 0.631 | |

> 300 | 0.130 | 0.732 | > 300 | 0.772 | 0.161 | |

60 ° | < 20 | 0.562 | 0.326 | < 40 | 24.000 | 1.000 |

20 - 200 | 0.306 | 0.529 | 40 - 00 | 3.240 | 0.457 | |

> 400 | 0.108 | 0.703 | > 400 | 0.760 | 0.215 | |

< 20 | 0.562 | 0.326 | < 50 | 24.000 | 1.000 | |

20 - 500 | 0.331 | 0.503 | 50 - 500 | 2.800 | 0.451 | |

> 500 | 0.087 | 0.718 | > 500 | 0.639 | 0.213 |

The optimization problem is formulated in such a way that the best combination of the parameters of a given PHE minimizes the number of plates. The optimization method used is based on screening [15], where for a given type of plate, the number of thermal plates is the objective function that has to be minimized. In order to avoid impossible or non-optimal solutions, certain inequality constraints are employed. An algorithm has been proposed in a screening method that uses MATLAB for optimization of a PHE, considering the plate type as the optimization variable [16]. For each type of plate, local optimal configurations are found (if they exist) that employ the fewest plates. Comparison of all the local optima then gives a global optimum.

## 5. Formulation of the optimization problem

Minimize:

Subject to:

If closed-form model is considered, as the closed-form equations are limited for some number of passes, there are two more constraints:

Depending on the equipment model, the number of plates can vary between 3 and 700. The first constraint (38) is imposed according to the PHE capacity. Constraints (39) and (40) can also be imposed, depending on the available pumping power. The velocity constraints are usually imposed in order to avoid dead spaces or air bubbles inside the set of plates. In practice, velocities less than 0.1 m/s are not used [5].

The optimization problem is solved by successively evaluating the constraints, reducing the number of configurations until the optimal set (OS) is found (if it exists). The screening process begins with the identification of an initial set (IS) of possible configurations, considering the channel limits. A reduced set (RS) is generated by considering the velocity and pressure drop constraints. The constraint of thermal effectiveness is then applied to the RS, in increasing order of the number of channels. Configurations with the smallest number of channels form the local optima set. The global optimum can therefore be found by comparing all the local optima. It is important to point out that the global optimum configuration may have a larger total heat transfer area. However, it is usually more economical to use a smaller number of large plates than a greater number of small plates [17]. The optimization algorithm is described in Table 4.

In Step 5, both methods can be used. The model using algebraic equations has the limitation of only being applicable to PHEs that are sufficiently large not to be affected by end channels and channels between adjacent passes. Industrial PHEs generally possess more than 40 thermal plates, although the limitation of the number of passes can still be a drawback. The major advantage of the model using differential equations is its general applicability to any configuration, without having to derive a specific closed-form equation for each configuration.

Steps | Mathematical relations | Comments |

Step 1: Input data for each type of plate. | PHE dimensions, fluids physical proprieties, mass flow rate and inlet temperature of both streams, constraints. | |

Step 2: Verification of the number of plates constraint. The initial set of configurations (IS) is determined. | The vector | |

Step 2.1: For each element of the vector | They are integer divisors of the number of channels of the corresponding side. If one is using closed-form equations, the number of passes constraint must be considered. | |

Step 3: Verification of the hydraulic constraints. The reduced set of configurations (RS) is determined. | Fluid velocity and pressure drop constrains are verified. | |

Step 3.1 Take | Cold fluid flows in the side I and the hot fluid in the side II. | |

Step 3.1.1 Verification of the velocity constraint for the fluid in side I. | Cold fluid velocity | |

Step 3.1.2 Verification of the pressure drop constraint in side I | Eq. (31) | The cold fluid pressure drop is calculated, |

Step 3.1.3 Verification of the velocity constraint for the fluid in side II. | Analogous to Step 3.1.1. | |

Step 3.1.4 Verification of the pressure drop constraint in side II. | Eq. (31) | Analogous to Step 3.1.2. |

Stage 3.2 Take | Analogous to stage 3.1. | |

Stage 3.3: Combination of the configuration parameters. | For | |

Step 4: Calculate the effectiveness in pure countercurrent flow, | If | |

Step 5: Verification of the thermal effectiveness constraint. The local optimal set of configurations (OS) is determined. | The selected configurations in Step 4 are simulated in a crescent order of the number of channels to find the possible local optimum set (OS). The remaining configurations do not need to be simulated. Both modeling can be used. | |

Step 6: Find global optimum. | By comparing all local optima, the global optimum is found. |

However, a drawback is the highly complex implementation of the simulation algorithm (see Table 5), in contrast to the second model, which is very simple.

Steps | Equations and tables | Comments |

Step 1. Tri-diagonal matrix coefficients are computed. | As | |

Step 2. Tri-diagonal matrix construction. | Eq. (6) | |

Step 3. Eigenvalues and eigenvectors are computed. | If one is using Matlab, one can use build-in functions. | |

Step 4. Linear system generation. | The boundary condition equations in algorithmic form of inlet fluids and change of pass are used. | |

Step 4.1. Generation of eigenvalues and eigenvectors matrix, | The resulting matrices are the sum of the matrices of both sides of the PHE. | |

Step 5. Determination of | ||

Step 6. Determination of the output dimensionless temperatures. | Table 6 | The boundary condition of output fluid in algorithmic form is used. |

Step 7. Computation of thermal effectiveness. | It may be obtained considering any side of the PHE, because the energy conservation is obeyed only if |

### 5.1. Simulation algorithm for the model using differential equations

For the development of this algorithm, the boundary conditions equations are used in the algorithm form described previously [6] (see Table 6). The simulation algorithm is applied separately for each value of

## 6. Case study

A case study was used to test the developed algorithm and compare the two mathematical models. Data were taken from examples presented in [18]. A cold water stream exchanges heat with a hot water stream of process. As the closed-form equations only consider configurations with a maximum of 4 passes for each fluid, a case was chosen in which the reduced set only had configurations with less than 4 passes for each stream. Table 8 presents the data used. Only one type of plate was considered.

The RS was obtained by applying the optimization algorithm up to Step 3. The optimal set was found by applying Step 5. As only one type of plate was considered, the local optimum was also the global optimum. The same optimal set was found with both approaches: two heat exchanger configurations with 144 channels and a 3-2 asymmetric pass arrangement (as presented in Table 9).

Plate characteristics | |

Process-water | Cooling-water |

Tin,cold= 20.0 °C | |

Constraints | |

80 | |

10 | |

It can be seen from the Tables that the simulations using the two models resulted in values that were very close. It is important to point out that the closed-form equations are only applicable when the end effects can be neglected (in the present case, when the number of thermal plates was greater than 40).

E (differential equation model) | E (closed-form model) | |||||||||||

1 | 80 | 1 | 2 | 0 | 80.1 | 80.3 | 80.3 | 80.1 | 80.4 | 80.4 | 80.4 | 80.4 |

2 | 80 | 2 | 1 | 1 | 80.3 | 80.1 | 80.3 | 80.1 | 80.4 | 80.4 | 80.4 | 80.4 |

3 | 81 | 1 | 2 | 0 | 80.3 | 80.3 | 80.3 | 80.3 | 80.6 | 80.6 | 80.6 | 80.6 |

4 | 83 | 2 | 1 | 1 | 80.5 | 80.5 | 80.5 | 80.5 | 80.8 | 80.8 | 80.8 | 80.8 |

5 | 84 | 1 | 2 | 0 | 80.5 | 80.8 | 80.8 | 80.5 | 80.9 | 80.9 | 80.9 | 80.9 |

6 | 84 | 2 | 1 | 1 | 80.8 | 80.5 | 80.8 | 80.5 | 80.9 | 80.9 | 80.9 | 80.9 |

7 | 85 | 1 | 2 | 0 | 80.8 | 80.8 | 80.8 | 80.8 | 81.0 | 81.0 | 81.0 | 81.0 |

8 | 87 | 2 | 1 | 1 | 80.9 | 80.9 | 80.9 | 80.9 | 81.2 | 81.2 | 81.2 | 81.2 |

9 | 88 | 1 | 2 | 0 | 80.9 | 81.2 | 81.2 | 80.9 | 81.3 | 81.3 | 81.3 | 81.3 |

10 | 88 | 2 | 1 | 1 | 81.2 | 80.9 | 81.2 | 80.9 | 81.3 | 81.3 | 81.3 | 81.3 |

11 | 89 | 1 | 2 | 0 | 81.2 | 81.2 | 81.2 | 81.2 | 81.4 | 81.4 | 81.4 | 81.4 |

12 | 91 | 2 | 1 | 1 | 81.4 | 81.4 | 81.4 | 81.4 | 81.6 | 81.6 | 81.6 | 81.6 |

13 | 92 | 1 | 2 | 0 | 81.3 | 81.6 | 81.6 | 81.3 | 81.7 | 81.7 | 81.7 | 81.7 |

14 | 92 | 2 | 1 | 1 | 81.6 | 81.3 | 81.6 | 81.3 | 81.7 | 81.7 | 81.7 | 81.7 |

15 | 93 | 1 | 2 | 0 | 81.6 | 81.6 | 81.6 | 81.6 | 81.8 | 81.8 | 81.8 | 81.8 |

16 | 95 | 2 | 1 | 1 | 81.7 | 81.7 | 81.7 | 81.7 | 82.0 | 82.0 | 82.0 | 82.0 |

17 | 96 | 1 | 2 | 0 | 81.7 | 81.9 | 81.9 | 81.7 | 82.1 | 82.1 | 82.1 | 82.1 |

18 | 96 | 2 | 1 | 1 | 81.9 | 81.7 | 81.9 | 81.7 | 82.1 | 82.1 | 82.1 | 82.1 |

19 | 97 | 1 | 2 | 0 | 81.9 | 81.9 | 81.9 | 81.9 | 82.2 | 82.2 | 82.2 | 82.2 |

20 | 144 | 2 | 3 | 0 | 71.8 | 71.7 | 92.8 | 92.9 | 71.8 | 71.8 | 93.0 | 93.0 |

21 | 144 | 3 | 2 | 1 | 71.7 | 71.8 | 92.8 | 92.9 | 71.8 | 71.8 | 93.0 | 93.0 |

22 | 149 | 3 | 2 | 1 | 71.7 | 71.7 | 93.0 | 93.0 | 71.7 | 71.7 | 93.2 | 93.2 |

## 7. Conclusions

In this chapter it was presented the development of two models for the design and optimization of plate heat exchangers. Both mathematical models were used to accomplish the heat exchanger design simulations. These methods use differential equations and closed-form equations based on the notion that a multi-pass PHE can be reduced to an arrangement consisting of assemblies of single-pass PHEs.

As a case study, an example obtained from the literature was used. The optimal sets were the same for both approaches, and agreement was achieved between the effectiveness values. The model using algebraic equations has the limitation of only being applicable to PHEs sufficiently large not to be affected by end channels and channels between adjacent passes. However, industrial PHEs generally possess more than 40 thermal plates. The major advantage of using this model is its general applicability to any configuration, without having to derive a specific closed-form equation for each configuration. However, its drawback is the highly complex implementation of the simulation algorithm, unlike the second approach, which is very simple.

## Nomenclature

Effective plate heat transfer area (m^{2}) | |

Plate effective area, m^{2} | |

Eigenvalues and eigenvectors matrix | |

Average thickness channel, m | |

Binary vector | |

Specific heat, J/kg·K | |

Heat capacity ratio | |

Equivalent diameter of the channel(m) | |

Port diameter of the plate(m) | |

Exchanger effectiveness, % | |

Thermal effectiveness in purely countercurrent flow (%) | |

Fanning factor | |

Gravitational acceleration (m/s^{2}) | |

Channel mass velocity (kg/m^{2}s) | |

Port mass velocity (kg/m^{2}s) | |

IS | Initial set of configurations |

Thermal conductivity(W/mK) | |

Plate thermal conductivity, W/m·K | |

Plate length, m | |

Horizontal port distance (m) | |

Vertical port distance (m) | |

mass velocity per channel (kg/s) | |

Mass flow rate, kg/s | |

Tri-diagonal matrix | |

Number of channels per pass | |

Number of channels | |

Number of plates | |

Number of transfer units | |

Nusselt number | |

OS | Optimal set of configurations |

Number of passes | |

Temperature effectiveness | |

Prandtl number | |

Heat transfer rate(J/s) | |

Heat capacity rate ratio | |

Fouling factor, K/W | |

Reynolds number | |

RS | Reduced set of configurations |

Binary parameter for flow direction | |

Plate thickness, m | |

Overall heat transfer coefficient, W/m^{2}·K | |

Fluid velocity inside channels, m/s | |

Plate width, m | |

Binary parameter for hot fluid location | |

Binary parameter for type of channel flow | |

Eigenvector of the tri-diagonal matrix |

Greek symbols

Heat transfer coefficient | |

Chevron corrugation inclination angle, degrees | |

Pressure drop, Pa | |

Log-mean temperature difference (K) | |

Effective mean temperature difference (K) | |

Normalized plate length | |

Dimensionless fluid temperature | |

Eingevalue of the tri-diagonal matrix | |

viscosity (Pa.s) | |

Viscosity at wall temperature (Pa.s) | |

density (kg/m^{3}) | |

Enlargement factor of the plate area | |

Parameter for feed connections position |

Subscripts

Cold fluid | |

Countercurrent | |

Hot fluid | |

Generic element | |

Inlet | |

Generic element | |

Outlet |

Superscripts

Odd channels of the heat exchanger | |

Even channels of the heat exchanger | |

Maximum | |

Minimum |