Thermal Energy Storage by PCM using fins (Paper CFD Validation), ANSYS Fluent
The present problem simulates heat transfer in a triplex tube heat exchanger containing phase change material (PCM) using ANSYS Fluent software. This simulation is based on the information of a reference article “Internal and external fin heat transfer enhancement technique for latent heat thermal energy storage in triplex tube heat exchangers” and its results are compared and validated with the results in the article.
This product includes a Mesh file and a comprehensive Training Movie.
There are some free products to check our service quality.
To order your ANSYS Fluent project (CFD simulation and training), contact our experts via [email protected], online support, or WhatsApp.
The present problem simulates heat transfer in a triplex tube heat exchanger containing phase change material (PCM) using ANSYS Fluent software. This simulation is based on the information of a reference article “Internal and external fin heat transfer enhancement technique for latent heat thermal energy storage in triplex tube heat exchangers” and its results are compared and validated with the results in the article. The triplex tube heat exchanger consists of two coaxial tubes that divide the heat exchanger cross-section into three sections or three tubes.
In the designed heat exchanger in this simulation, the inner and outer tubes are made of aluminum and the space between these two tubes is the place where the phase change material flows. To reduce the computational cost, the modeling is done in two dimensions and only on a section of this heat exchanger. Also, two different models of cross section of this heat exchanger have been designed, the first model has no fin (case A in the article) and the second model has four fins in its interior (case B in the article). Also, the PCM used in the simulation is RT82, whose properties are defined based on the data in Table 4 of the article.
Thus, this material has a density equal to 770 kg.m-3 and a specific heat capacity equal to 2000 j.kg-1.K-1 and a thermal conductivity equal to 0.2 Wm-1.K-1 and a viscosity equal to 0.03499 kg .m-1.s-1. Also, according to the data in this table, the maximum temperature at which the solid phase temperature is (solidus temperature) is equal to 350.15 K and the minimum temperature at which the liquid phase is dominant (liquidus temperature) is 358.15 K and the Pure solvent melting heat is defined as 176,000 j.kg-1.
In this heat exchanger, the inner wall of the inner tube and the outer wall of the outer tube are considered as thermal insulation; While the outer wall of the inner tube and the inner wall of the outer tube, which are in direct contact with the PCM material, have a constant temperature of 363.15 K.
Geometry & Mesh
The present model is designed in two dimensions using Design Modeler software. The model is a cross-section of a three-tube heat exchanger drawn in two different modes. A circular cross section is designed to be drawn in both models as two pipes coaxial with each other; Thus, four circles with different diameters are created. The inner and outer radii of the inner tube are 25.4 mm and 26.6 mm, respectively, and the inner and outer radii of the outer tube are 75 mm and 76 mm, respectively. In the first model, no fins are designed, but in the second model, two fins are designed to face each other on the inner tube’s outer wall, and the other two fins are designed to meet each other on the inner wall of the outer tube.
We carry out the meshing of the model using ANSYS Meshing software, and the mesh type is structured. The element number is 23908 for the first case study, and 24492 for the second one with fins. The following figure shows the mesh.
Thermal Energy Storage CFD Simulation
We consider several assumptions to simulate the present model:
- We perform a pressure-based solver.
- The simulation is unsteady. Because the purpose of the present work is to study the process of Solidification and Melting over time.
- The gravity effect on the fluid is equal to -9.81 m.s-2 along the Y-axis.
The following table represents a summary of the defining steps of the problem and its solution:
|Solidifiacation & Melting||On|
|Boundary conditions (PCM)
|Outer Wall of Inner Tube & Inner Wall of Outer Tube||Wall|
|wall motion||stationary wall|
|Outer Wall of Outer Tube & Inner Wall of Inner Tube||Wall|
|wall motion||stationary wall|
|heat flux||0 W.m-2|
|momentum||second order upwind|
|energy||second order upwind|
|gauge pressure||0 pascal|
|x-velocity & y-velocity||0 m.s-1|
The present validation is based on the diagram in Figure 14 of the mentioned article. This diagram is related to the changes in the melting fraction over time. This diagram includes various geometric designs based on the arrangement of heat exchanger fins. In the present simulation, two case studies are examined, which are related to case A and case B from Figure 2 of the paper. Based on the resulting images, it can be said that almost all of the PCM material turns into a liquid phase after a specific time.
At the end of the solution process, two-dimensional contours related to the temperature and liquid mass fraction are obtained for both different models (the fin model and the without fin model). It is observed that as the temperature of the phase change material increases, the amount of liquid volume fraction increases.
There are a Mesh file and a comprehensive Training Movie that presents how to solve the problem and extract all desired results.