Author Contributions
Conceptualization, J.K. and M.F.; methodology, J.K. and M.F.; software, M.F.; validation, J.K.; formal analysis, M.F.; investigation, J.K.; resources, A.P.S.; data curation, M.F. and A.S.; writing—original draft preparation, M.F. and J.K.; writing—review and editing, M.F.; visualization, M.F.; supervision, A.P.S.; project administration, A.P.S.; funding acquisition, A.P.S. All authors have read and agreed to the published version of the manuscript.
Figure 1.
Cross section schematic of double pipe heat exchanger with and without PCM thermal capacitor.
Figure 1.
Cross section schematic of double pipe heat exchanger with and without PCM thermal capacitor.
Figure 2.
Validation of numerical model with experiments [
15].
Figure 2.
Validation of numerical model with experiments [
15].
Figure 3.
Case A and Case B configuration of hot and cold streams within the triplex heat exchanger.
Figure 3.
Case A and Case B configuration of hot and cold streams within the triplex heat exchanger.
Figure 4.
Cold outlet temperature fluctuation range and average temperature at various Reynolds numbers for different PCMs in Case A (hot inner) and Case B (cold inner) configuration.
Figure 4.
Cold outlet temperature fluctuation range and average temperature at various Reynolds numbers for different PCMs in Case A (hot inner) and Case B (cold inner) configuration.
Figure 5.
Cold outlets for base case and thermal capacitor cases with four PCMs at Re = 1000, cold inlet temperature is present to show the absolute temperature changes
Figure 5.
Cold outlets for base case and thermal capacitor cases with four PCMs at Re = 1000, cold inlet temperature is present to show the absolute temperature changes
Figure 6.
Cold outlet temperature and liquid fraction for eicosane ( = 309.35 K, = 310.35 K) at Reynolds from 200 to 2000. Reynolds 200–1000 on the left (a,c) and 1200–2000 on the right (b,d). Liquid fractions are shown in (c,d) show low usage of the PCMs full phase change capacity with (b,d) showing a trend of lower average temperature of the flow with higher PCM use.
Figure 6.
Cold outlet temperature and liquid fraction for eicosane ( = 309.35 K, = 310.35 K) at Reynolds from 200 to 2000. Reynolds 200–1000 on the left (a,c) and 1200–2000 on the right (b,d). Liquid fractions are shown in (c,d) show low usage of the PCMs full phase change capacity with (b,d) showing a trend of lower average temperature of the flow with higher PCM use.
Figure 7.
Cold outlet temperature and liquid fraction for myristic acid ( = 330.65 K, = 331.65 K) at Reynolds from 200 to 2000. Reynolds 200–1000 on the left (a,c) and 1200–2000 on the right (b,d). The PCM is not active as seen in (c,d), as the material in the PCM layer never reaches the much higher phase change temperature, thus (a,b) show results for sensible heat storage only.
Figure 7.
Cold outlet temperature and liquid fraction for myristic acid ( = 330.65 K, = 331.65 K) at Reynolds from 200 to 2000. Reynolds 200–1000 on the left (a,c) and 1200–2000 on the right (b,d). The PCM is not active as seen in (c,d), as the material in the PCM layer never reaches the much higher phase change temperature, thus (a,b) show results for sensible heat storage only.
Figure 8.
Cold outlet temperature and liquid fraction for octadecane ( = 300.65 K, = 301.65 K) at Reynolds from 200 to 2000. Reynolds 200–1000 on the left (a,c) and 1200–2000 on the right (b,d). The octadecane remains in a partial phased changed state during the full oscillation cycle. This is reflected in the shape of the cold outlet variations which have a stretched out shape on both the positive and negative variation of the cycle unlike with the eicosane which only shows the structure on the positive component.
Figure 8.
Cold outlet temperature and liquid fraction for octadecane ( = 300.65 K, = 301.65 K) at Reynolds from 200 to 2000. Reynolds 200–1000 on the left (a,c) and 1200–2000 on the right (b,d). The octadecane remains in a partial phased changed state during the full oscillation cycle. This is reflected in the shape of the cold outlet variations which have a stretched out shape on both the positive and negative variation of the cycle unlike with the eicosane which only shows the structure on the positive component.
Figure 9.
Cold outlet temperature and liquid fraction for paraffin wax ( = 319.15 K, = 320.15 K) at Reynolds from 200 to 2000. Reynolds 200–1000 on the left (a,c) and 1200–2000 on the right (b,d). The paraffin wax is shown to only phase change a small amount at all Reynolds numbers, the small amounts do not noticeably affect the outlet temperature and the results are very similar to myristic. Despite the phase change temperature being close to the average temperature of the oscillating hot stream, the material is not very active demonstrating the difficulty in selecting an appropriate phase change temperature.
Figure 9.
Cold outlet temperature and liquid fraction for paraffin wax ( = 319.15 K, = 320.15 K) at Reynolds from 200 to 2000. Reynolds 200–1000 on the left (a,c) and 1200–2000 on the right (b,d). The paraffin wax is shown to only phase change a small amount at all Reynolds numbers, the small amounts do not noticeably affect the outlet temperature and the results are very similar to myristic. Despite the phase change temperature being close to the average temperature of the oscillating hot stream, the material is not very active demonstrating the difficulty in selecting an appropriate phase change temperature.
Figure 10.
Average temperature and oscillation range for the cold stream outlet as a function of Reynolds number and phase change material. As Reynolds increases a general trend in decrease of the average temperature is seen due to decreased residence time within the pipe. Oscillations also become larger as heat transfer rates increase with no change in the amount of thermal mass being engaged. In the case of eicosane, more material phase changes at higher Reynolds which causes the oscillations to dampen.
Figure 10.
Average temperature and oscillation range for the cold stream outlet as a function of Reynolds number and phase change material. As Reynolds increases a general trend in decrease of the average temperature is seen due to decreased residence time within the pipe. Oscillations also become larger as heat transfer rates increase with no change in the amount of thermal mass being engaged. In the case of eicosane, more material phase changes at higher Reynolds which causes the oscillations to dampen.
Figure 11.
Results of eicosane at Re = 1000 with amplitudes of 10, 20, 30 K and periods of 5, 10, 20 min. Variations due to amplitude differences are seen to be clipped at the maximum at the same point even with higher temperature variations. Higher period oscillations allow for higher max temperature with more pronounced clipping.
Figure 11.
Results of eicosane at Re = 1000 with amplitudes of 10, 20, 30 K and periods of 5, 10, 20 min. Variations due to amplitude differences are seen to be clipped at the maximum at the same point even with higher temperature variations. Higher period oscillations allow for higher max temperature with more pronounced clipping.
Figure 12.
Results of myristic acid at Re = 1000 with amplitudes of 10, 20, 30 K and periods of 5, 10, 20 min. No phase change occurs even with the higher oscillation range and thus the wave simply is scaled. Higher and lower periods result in stretched out and compressed oscillations with higher maximums due to longer diffusion times.
Figure 12.
Results of myristic acid at Re = 1000 with amplitudes of 10, 20, 30 K and periods of 5, 10, 20 min. No phase change occurs even with the higher oscillation range and thus the wave simply is scaled. Higher and lower periods result in stretched out and compressed oscillations with higher maximums due to longer diffusion times.
Figure 13.
Results of N-octadecane at Re = 1000 with amplitudes of 10, 20, 30 K and periods of 5, 10, 20 min. Octadecane remains at all times in a semi-phase changed state, when the oscillation range is increased the top end of the temperature variations remain similar due to the constant temperature phase change. The high period case showed very constant PCM activity, whereas the low period case showed very rapid swings in PCM use. In all cases the absolute change in oscillation remains small with oscillations less than 1 K.
Figure 13.
Results of N-octadecane at Re = 1000 with amplitudes of 10, 20, 30 K and periods of 5, 10, 20 min. Octadecane remains at all times in a semi-phase changed state, when the oscillation range is increased the top end of the temperature variations remain similar due to the constant temperature phase change. The high period case showed very constant PCM activity, whereas the low period case showed very rapid swings in PCM use. In all cases the absolute change in oscillation remains small with oscillations less than 1 K.
Figure 14.
Results of paraffin wax at Re = 1000 with amplitudes of 10, 20, 30 K and periods of 5, 10, 20 min. When the temperature oscillations increase the PCM becomes slightly active which results in slightly lower peak temperature. In the case of variations in period, longer periods did not result in phase change of the material thus more heat diffused through the material due to the longer period at a high temperature in a fashion similar to the myristic acid.
Figure 14.
Results of paraffin wax at Re = 1000 with amplitudes of 10, 20, 30 K and periods of 5, 10, 20 min. When the temperature oscillations increase the PCM becomes slightly active which results in slightly lower peak temperature. In the case of variations in period, longer periods did not result in phase change of the material thus more heat diffused through the material due to the longer period at a high temperature in a fashion similar to the myristic acid.
Table 1.
Geometric parameters of heat exchanger.
Table 1.
Geometric parameters of heat exchanger.
Parameters | Value | Unit | Description |
---|
Di | 3.00 × | m | diameter of inner pipe |
Dm | 4.00 × | m | diameter of phase change pipe |
Do | 5.66 × | m | diameter of outer pipe |
Z | 2 | m | length of heat exchange domain |
Table 2.
Fitting constants for water material characteristic functions.
Table 2.
Fitting constants for water material characteristic functions.
Parameters | Value | Unit |
---|
| −3.57 × | kg |
| 1.88 | kg |
| 753.2 | kg |
| 2.591 × | kg |
| 238.3 | K |
| 143.2 | K |
| −8.354 × | W |
| 6.53 × | W |
| −0.5981 | W |
| 4200 | |
Table 3.
Fitting constants for PCM material characteristic functions.
Table 3.
Fitting constants for PCM material characteristic functions.
| 2160 (n-octadecane) | |
| 2210 (n-eicosane) | |
| 1590 (myristic acid) | |
| 2890 (wax) | |
| 300.15, 300.65 (n-octadecane) | K |
| 309.15, 309.65 (n-eicosane) | |
| 331.15, 331.65 (myristic acid) | |
| 319.15, 320.65 (wax) | |
| 0.358 (n-octadecane) | |
| 0.358 (n-eicosane) | |
| 0.365 (myristic acid) | |
| 0.321 (wax) | |
| 0.148 (n-octadecane) | |
| 0.148 (n-eicosane) | |
| 0.159 (myristic acid) | |
| 0.120 (wax) | |
| 0.001 | |
| −4.649 (n-octadecane) | - |
| −4.300 (n-eicosane) | |
| −4.659 (myristic acid) | |
| −4.250 (wax) | |
| 1.790 × | K |
L | 2.43 × (n-octadecane) | |
| 2.47 × (n-eicosane) | |
| 1.99 × (myristic acid) | |
| 1.73 × (wax) | |
| 774 (n-octadecane) | kg |
| 778 (n-eicosane) | |
| 860 (myristic acid) | |
| 750 (wax) | |
| 9.00 × | |
Table 4.
Base case boundary conditions and parameters.
Table 4.
Base case boundary conditions and parameters.
Parameters | Value | Unit |
---|
| 313.15, 20 | K |
| 293.15 | K |
(gauge) | 0 | Pa |
| 600 | s |
Table 5.
Inlet velocity boundary condition for hot and cold streams.
Table 5.
Inlet velocity boundary condition for hot and cold streams.
Reynolds Number | Inner (m/s) | Annulus (m/s) |
---|
200 | 0.0044 | 0.00246 |
400 | 0.0088 | 0.00493 |
600 | 0.0132 | 0.00740 |
800 | 0.0176 | 0.00987 |
1000 | 0.0219 | 0.01234 |
1200 | 0.0263 | 0.01481 |
1400 | 0.0307 | 0.01727 |
1600 | 0.0351 | 0.01974 |
1800 | 0.0395 | 0.02221 |
2000 | 0.0439 | 0.02468 |
Table 6.
Comparison of results with and without natural convection enabled.
Table 6.
Comparison of results with and without natural convection enabled.
| 2D | 3D (Convective) | % Error |
---|
Annulus Inlet Temp (K) | 293.15 |
Inner Inlet Temp (K) | 313.15 |
Bulk Annulus Temperature (K) | 293.49 | 294.01 | 0.18% |
Bulk Inner Temperature (K) | 312.52 | 311.74 | 0.25% |
Bulk PCM Temperature (K) | 300.25 | 300.30 | 0.02% |
PCM Liquid Fraction | 0.3077 | 0.3128 | 1.63% |
Flux Through PCM (W) | 121.11 | 119.69 | 1.19% |
Table 7.
Range and Average temperature for the cold outlet of the PCM capacitors at different Reynolds numbers.
Table 7.
Range and Average temperature for the cold outlet of the PCM capacitors at different Reynolds numbers.
Cold Stream Outlet Average Temperature (K) |
---|
Reynolds | Eicosane | Myristic | Octadecane | Wax |
200 | 300.46 | 300.52 | 299.48 | 299.19 |
400 | 297.80 | 297.85 | 297.08 | 296.82 |
600 | 296.60 | 296.66 | 296.03 | 295.82 |
800 | 295.91 | 295.98 | 295.44 | 295.27 |
1000 | 295.46 | 295.54 | 295.07 | 294.91 |
1200 | 295.15 | 295.24 | 294.81 | 294.67 |
1400 | 294.92 | 295.01 | 294.61 | 294.49 |
1600 | 294.74 | 294.84 | 294.46 | 294.36 |
1800 | 294.60 | 294.70 | 294.35 | 294.25 |
2000 | 294.48 | 294.58 | 294.25 | 294.16 |
Cold Stream Outlet Temperature Range (K) |
Reynolds | Eicosane | Myristic | Octadecane | Wax |
200 | 0.11 | 0.14 | 0.07 | 0.09 |
400 | 0.45 | 0.51 | 0.04 | 0.30 |
600 | 1.05 | 1.14 | 0.09 | 0.68 |
800 | 1.40 | 1.69 | 0.14 | 1.04 |
1000 | 1.50 | 1.99 | 0.17 | 1.24 |
1200 | 1.49 | 2.12 | 0.18 | 1.33 |
1400 | 1.45 | 2.15 | 0.19 | 1.36 |
1600 | 1.39 | 2.14 | 0.19 | 1.35 |
1800 | 1.34 | 2.10 | 0.18 | 1.31 |
2000 | 1.29 | 2.05 | 0.18 | 1.27 |