Plate heat exchanger can obtain higher thermal performance because of its advantages in high heat transfer coefficient, small scale, and can realize pure counter current flow. It has been widely applied in HVAC industry. In this paper, the numerical research of plate heat exchanger in solar heating system has been proposed. Aimed at the type of herringbone corrugated plate which has better thermal performance and been widely used, the three dimensional model is established by Gambit software. Using FLUENT software for numerical calculation, by studying the effect of corrugated inclination angle, corrugated depth, corrugated spacing and inlet velocity of heat exchanger on internal temperature, pressure, velocity distributions of domains, the relationship between the above parameters and the Nusselt Number and the pressure drop was obtained by simulation data. Heat transfer coefficient and pressure drop correlations used to measure the overall performance of the heat exchanger. The result shows that the optimal structure parameters is corrugated angle 60°, corrugated depth 4 mm and corrugated spacing 16 mm.
The availability of resources and energy is a serious problem of the whole world. Because of the environment, resources and energy crisis in recent years, the requirement of saving energy and reducing consumption is increasingly high, and using efficient heat exchange equipment to reduce the energy consumption has become the focus of modern industry attention and study. Compared with the traditional shell and tube heat exchanger, the development of plate heat exchanger has the obvious advantage [
Plate heat exchanger is a kind of efficient and compact heat exchange equipment and its application involves almost all industrial fields. In recent years, brazing type plate heat exchanger has been known by people widely by its superiority of good compactness, light weight, good heat transfer performance and low initial cost. Therefore, it is very necessary to find an effective and feasible method for the study of flow and heat transfer between heat exchanger plates. It is generally recognized that herringbone corrugated plate have higher heat transfer efficiency, big resistance and bearing capacity, because the complex change in plate flow channel section is easy to induce turbulence and the fluid flow consumes more energy in this kind of variable flow [
Each plate shall be four angle hole, there are two angle hole provide flow channel on the each side of the plate. When assembling plate, angle hole is arranged into bypass duct of two fluids in proper order. In a group of channel, fluid flow between the two plate after entering angle hole and outflow from the other hole. Each kind of fluid has a single port because of sheet gasket alternant of angle hole. The most common material of plate heat exchanger is stainless steel and titanium, sometimes using nickel alloy steel [
How to use the remaining energy reasonably and improve energy utilization has become a research direction for experts and scholars from all over the world. Authors of [
In this paper, numerical simulation is conducted for the most widely used herringbone corrugated plate heat exchanger. The temperature field and pressure field of the fluid in the channel between the plates were simulated and the main factors influencing the heat transfer and flow were analyzed. At the same time, by analyzing the performance of the plate heat exchanger and analyzing the heat transfer performance and resistance characteristics, the relationship between plate heat exchanger performance and geometrical parameters of corrugated plate are researched to provide the theoretical basis for the optimization and development of the plate heat exchanger plate.
In the process of numerical simulation, we must first establish the mathematical model which reflects the nature of the problem, that is, reflect the differential equation and the definite condition between the various problems, and then find the high efficiency and high precision calculation method.
Combined with the characteristics of flow and heat transfer in the actual plate heat exchanger, the mathematical model based on the numerical simulation is mainly used as follows:
1) The working medium is incompressible Newtonian fluid;
2) Gravity and floating capacity due to density differences are ignored;
3) Due to the lower flow velocity of the fluid in the heat exchanger, ignore the thermal effect of the viscous dissipation when the fluid flows;
4) Assuming that the fluid inlet speed is known, the outlet pressure is constant. The inlet pressure and outlet speed are free boundary conditions.
The control equations for the fluid motion and heat transfer in the calculation region are as follows:
Conti Nuity equation:
∂ u i ∂ x i = 0 (1)
Mass conservation equation:
∂ u ∂ x + ∂ v ∂ y + ∂ w ∂ z = 0 (2)
where u, v, w is the component of flow velocity.
The momentum conservation equation in the i-direction Cartesian coordinate system:
u ∂ U i ∂ x + v ∂ U i ∂ y + w ∂ U i ∂ z = − 1 ρ ∂ p ∂ x 2 + v ( ∂ 2 U i ∂ x 2 + ∂ 2 U i ∂ y 2 + ∂ 2 U i ∂ z 2 ) (3)
where ρ is the fluid density, p is the pressure, v is the kinematic viscosity, and Ui is the velocity component in the i direction.
Energy conservation equation:
u ∂ t ∂ x + v ∂ t ∂ y + w ∂ t ∂ z = a ∂ 2 t ∂ x 2 + ∂ 2 t ∂ y 2 + ∂ 2 t ∂ z 2 (4)
where a is the thermal diffusion coefficient.
There are two main categories of turbulence models: Reynolds stress model and vortex model. In the vortex model, the Reynolds stress term is not directly processed, but the turbulence viscosity is introduced, and then the turbulence stress is expressed as a function of the turbulence viscosity, which reduces the difficulty of solving. The vortex model includes zero-equation model, one-equation model, and two-equation model. The standard k-ε model is a typical two-equation model, which is formed by introducing the turbulence dissipation rate ε on the basis of the k-equation, which is the most widely used turbulence model. In this paper, the RNG model in the k-ε turbulence model is used as the renormalization group model, can deal better with high strain rate and the flow with larger flow line bending degree, due to correct the turbulence viscosity, take into account the rotational flow in the average flow, add an item to ε, reflect the mainstream time-dependent change rate Eij. Because the RNG model is only for the turbulence with full development, so wall function method is adopted for low Reynolds Number movement near the wall, the dense and a large Number of grid layout on the wall saves memory and time, which is widely used in engineering turbulence calculation.
The object of this study is EATB55 brazing plate heat exchanger, the size of which is 539 mm long, width of which is 123 mm, and the plate thickness of which is 2.34 mm. The model created using gambit is shown in
In order to explore the main factors that affect the performance of the plate heat exchanger, the corrugated angel, the corrugated depth and the corrugated spacing are defined as shown in
ripple depth is 3 mm to 6 mm, and the ripple pitch is 10 mm to 20 mm.
The geometrical parameters of the plate for the stimulation of the single corrugated plate and the selection of the velocity of the plate are shown in
Boundary condition setting:
1) Fluid inlet adopts speed inlet boundary conditions: Assuming inlet speed is 0.6 m/s and inlet temperature is set to 300 K;
2) The outlet of the fluid adopts pressure outlet boundary conditions, set the outlet pressure as 101,325 Pa;
3) The upper and lower wall is set to constant wall, temperature 340 K;
4) The remaining walls are non-slip adiabatic walls, and the fluid medium is water.
The pressure distribution at different angle, when the spacing λ = 16 mm, depth h = 4 mm, v = 0.6 m/s is shown in
The temperature distribution at different angles is shown in
The velocity distribution at different angles is shown in
The local vector velocity at different angles is shown in
The corrugated angle is smaller, the local vector velocity around the corner is higher, the pressure is higher, the temperature is higher as well, the highest temperature appears at the corner. When β = 30˚ the details of local vector velocity is shown in
and the whole plate velocity distribution uniform, the plate will not cause damage.
When angle β = 60˚, depth h = 4 mm, v = 0.6 m/s, the pressure distribution at different spacing is shown in
The temperature distribution at different spacing is shown in
The velocity distribution at different spacing is shown in
The local vector velocity at different spacing is shown in
When angle β = 60˚, spacing λ = 16 mm, v = 0.6 m/s, pressure distribution at different depths is shown in
The temperature distribution at different depths is shown in
The velocity distribution at different depths is shown in
The local vector velocity at different depths is shown in
When angle β = 60˚, spacing λ = 16 mm, h = 4 mm, the pressure distribution at different inlet velocity is shown in
The temperature distribution at different inlet velocity is shown in
The velocity distribution at different inlet velocity is shown in
The local vector velocity at different inlet velocity is shown in
The overall performance of the heat exchanger is determined by the heat transfer performance and hydraulic characteristics. Where Nu can be used to measure the heat transfer performance, ΔP can be used to measure the hydraulic characteristics. In this paper, the relationship between the heat transfer coefficient and the pressure drop is used to measure the overall performance of the heat exchanger, said η:
η = K / Δ P 1 / 3 (5)
Change the corrugation angle, the relationship of β with ΔP, Nu, the average heat transfer coefficient and η is shown as follow
With the increase of β, the flow of fluid is smoother at each inflection point between the plates, resulting in the decrease of heat transfer performance. Nu and average heat transfer coefficient decrease with the same trend, and the decrease is more obvious at 60˚. The resistance through the plate is reduced, which causes the reduction of pressure drop, and it declines rapidly when the temperature is 60˚. When the inclination is less than 60˚, the cross section of the heat exchanger is crossed, and then the transition is a meandering flow, and the flow velocity decreases at the inflection point. The overall performance of the heat exchanger is increased first and then decreased, and the maximum value of which appears at 60˚.
Change the corrugation spacing, the relationship of λ with ΔP, Nu, the average heat transfer coefficient and η is shown in
With the increase of λ, the Number of corrugations decreases, the contact of
the fluid between the plates decreases, the vortex region between the plates increases, the secondary flow caused by the contact decreases, resulting in reduced heat transfer performance, Nu Number and the average thermal coefficient decreases steadily with the same trend, and the pressure of the fluid flowing through the plate decreases, and the amplitude increases at 16 mm, and then tends to be stable. The overall performance of the heat exchanger increased first and then became smooth, and the maximum value of which appears at 16 mm.
Change corrugation depth, the relationship of h with ΔP, Nu Number, average heat transfer coefficient and η is as follows
With the increase of h, the distance and time of the fluid flow between the plates increases, resulting in an increase in heat transfer performance. The Nu Number and the average heat transfer coefficient increase first and then decrease at the same trend, and reach the maximum value at 4 mm, while the fluid need more flow resistance to flow through the plate, there is a larger increase when in the 5 mm, and then it becomes stable. The velocity field of different h is similar, and the flow dead zone is reduced. The overall performance of the heat exchanger first increases and then becomes stable, the maximum value of which appears at 4 mm Department.
Change the velocity of the fluid, the relationship of v with ΔP, Nu Number, average heat transfer coefficient and η is as follows
With the increase of velocity, the turbulence degree of the fluid in the plates increases, which leads to the enhancement of the heat transfer performance. The Nu Number and the average heat transfer coefficient increase with the same trend, and the fluid requires more flow resistance to flow through the plate. The overall performance of the heat exchanger fluctuates little and has little effect.
In this paper, CFD method is used to predict the performance of plate heat exchanger in solar heating system. The influence of structural parameters like ripple inclination, ripple depth, ripple spacing and fluid velocity on the temperature field, pressure field and velocity field of the heat exchanger were studied by using FLUENT6.3 software, further analyzing its influence on the heat transfer and resistance of the chevron plate heat exchanger. The results show that the heat transfer effect at the contact point of corrugated plate heat exchanger is the best, the turbulence is the highest, the heat transfer is strengthened, but the pressure loss after passing the contact point of the fluid reaches the maximum, which leads to the increase of the pressure drop. With the increase of corrugation angle, there are two different flow patterns appear between the plates, and the overall performance first increases and then decreases, and it is optimal at 60˚. With the increase of ripple spacing, the heat transfer coefficient
and pressure drop are reduced at the same time, but the amplitude is different, and the overall performance reaches the maximum at 16 mm. The increase of the ripple depth leads to heat transfer coefficient increasing first and decreasing then, and the decrease of the pressure drop, and the overall performance is the best at 4 mm. The increase of the fluid flow rate leads to the increase of the heat transfer coefficient and the pressure drop, and the overall performance is gradually increased. In this paper, the optimal structure of the plate heat exchanger is ripple angle of 60˚, ripple depth of 4 mm, ripple spacing of 16 mm. The results of the above analysis are in good agreement with the conclusion of the references, indicating that the model used in this simulation can well reflect the flow and heat transfer characteristics of the chevron plate heat exchanger.
The author would like to acknowledge the following projects for financial support. Wind and solar energy utilization technology of Inner Mongolia University of Technology key laboratory open fund project authorized by China’s ministry of education: Wind Speed Prediction Research Based on Geographic Information System. The National Natural Science Foundation of China: Performance Analysis of Functioned Brackish Water for Lighting-Induced Evaporation Type Solar Desalination System. And Science & technology plan projects in Inner Mongolia (ZD2012MS0703).
Jia, Y., Wang, C., Chang, Z.H., Li, W.X. and Zhang, C. (2017) Heat Transfer Modelling of Plate Heat Exchanger in Solar Heating System. Open Journal of Fluid Dynamics, 7, 426-447. https://doi.org/10.4236/ojfd.2017.73029