Abstract

Design of an additively manufactured molten salt (MS) to supercritical carbon dioxide (sCO2) primary heat exchanger (PHE) for solar thermal power generation is presented. The PHE is designed to handle temperatures up to 720 °C on the MS side and an internal pressure of 200 bar on the sCO2 side. In the core, MS flows through a three-dimensional periodic lattice network, while sCO2 flows within pin arrays. The design includes integrated sCO2 headers located within the MS flow, allowing for a counterflow design of the PHE. The sCO2 headers are configured to enable uniform flow distribution into each sCO2 plate while withstanding an internal pressure of 200 bar and minimizing obstruction to the flow of MS around it. The structural integrity of the design is verified on additively manufactured (AM) 316 stainless steel sub-scale specimens. An experimentally validated, correlation-based sectional PHE core thermofluidic model is developed to study the impact of flow and geometrical parameters on the PHE performance, with varied parameters including the mass flowrate, surface roughness, and PHE dimensions. A process-based cost model is used to determine the impact of parameter variation on build cost. The model results show that a heat exchanger with a power density of 18.6 MW/m3 (including sCO2 header volume) and effectiveness of 0.88 can be achieved at a heat capacity rate ratio of 0.8. The impact of design and AM machine parameters on the cost of the PHE are assessed.

1 Introduction

Solar thermal power generation with thermal storage offers the possibility of using a renewable resource for baseload or peak-load electrical power generation. Over the past decade, there has been an impetus to use the supercritical carbon dioxide (sCO2) Brayton cycle for power generation. To couple the power cycle with the solar field and storage, three pathways are being explored by several research groups based on different receiver designs. These include the gas, liquid, and a solid particle medium pathway. The liquid pathway is mature, with roots in the Solar II demonstration project, and has been successfully commercialized over three decades [14]. However, in place of the nitrate salt mixture, to achieve higher turbine inlet temperatures of 700 °C, chloride salts are being considered for sCO2 solar thermal cycles with the liquid pathway. A two-tank system is typically employed, with the molten salt (MS) pumped from the cold tank through the receiver and stored in a hot tank. An MS-sCO2 primary heat exchanger is needed to transfer the heat into the power cycle. The inlet temperature of MS into this primary heat exchanger (PHE) from the hot tank would be 720 °C, and the sCO2 inlet temperature would be 500 °C at a pressure of 200 bar. The pressure drop across the sCO2 side is to be kept below 2 % of the line pressure (i.e., < 4 bar) from a cycle efficiency consideration. The PHE needs to be designed for high effectiveness as well as for a low exit temperature of MS to reduce parasitic losses in the receiver sub-system. In addition, the PHE must be able to resist corrosion of the MS. There is a dearth of literature on designs of PHEs that meet these challenging operating conditions and constraints. Printed circuit heat exchangers (PCHE) represent a common architecture for high-effectiveness and compact heat exchangers that lend to high conductance.

The PCHEs in literature have typically been designed and demonstrated for the 500–600 °C range on the high side. While the designs are applicable to higher temperatures, fabrication difficulties pertaining to etching and diffusion bonding need to be addressed. Aakre and Anderson [5] experimentally tested a 316L stainless steel diffusion-bonded heat exchanger using solar salt (0.6 NaNO3—0.4KNO3) on the hot side and 160 bar sCO2 on the cold side. The temperature on the salt side remained below 550 °C, because of salt degradation around 600 °C. Yao et al. [6] designed a 316 stainless steel MS-to-sCO2 pillow plate heat exchanger consisting of a pattern of concave and convex regions to enhance heat transfer through boundary layer thickness reduction. The heat exchanger was designed for MS entering at a temperature of 520 °C and exchanging energy with sCO2 at an inlet temperature of 335 °C. Shi et al. [7] developed microlaminated heat exchanger for this application with rectangular parallel channels on the cold side and airfoil-shaped fins on the hot side. While the design was for chloride MS to sCO2, they fabricated the device using stainless steel and demonstrated it at lower temperatures in Wang et al. [8] using a ternary salt (Hitec) and a synthetic oil. While the design shows improvement over traditional PCHE designs with zig-zag channels, it should be noted that the challenges for such a PHE are in using the microlamination method of chemical etching and diffusion bonding for nickel superalloys such as IN740H and H282 that are needed for this application owing to their superior creep strength. Caccia et al. [9] developed a design for a cermet PCHE for this application. The Zirconium–Tungsten composite material considered in this study allowed for heat transfer at temperatures ranging from 600 to 800 °C, an operating pressure of 200 bar on the sCO2 side and 1 bar for the molten salt, and an allowable stress of 125 MPa giving a factor of safety of 3 with respect to the material flexure strength of 370 MPa. They reported power densities of around 8.2 MW/m3 for the core, assuming straight parallel channels having a semi-circular cross section. Montes et al. [10] investigated the optimization of the PCHE at the junction between a solar field and supercritical CO2 power cycle under three different configurations recompression, intercooling, and partial-cooling. Their techno-economic analysis targeted an MS-to-sCO2 printed circuit heat exchanger made of Haynes 242, and it was concluded that the investment cost for the solar plant could be significantly reduced by optimizing the heat exchanger. The authors were able to achieve an effectiveness of 0.83 and thermal power of 103.4 MW-th for the optimized heat exchanger for the partial-cooling cycle configuration. In this configuration, they also observed a reduction in costs from 40.456 million dollars for the baseline to 7.924 million dollars for the optimized PCHE.

In recent years, metal and ceramics additive manufacturing (AM) has been investigated as a viable technique for the fabrication of heat exchangers combining high effectiveness and compactness, especially in applications with extreme conditions and high strength requirements. Singh et al. [11] designed a ceramic heat exchanger for concentrating solar power (CSP) applications. The use of ceramics for this application was justified by their superior physical properties and corrosion/oxidation resistance that make them desirable candidates for surviving applications involving corrosive fluids at high temperatures (up to 750 °C). The modeled geometry was optimized not only for heat transfer performance but also for stress constraints imposed by high pressures and temperatures. One notable advantage of the AM process is the ability to design integrated headers—a feature that eases the scalability of the design while reducing cost and increasing the power density of heat exchangers. In a follow-on work [12], the group printed silicon carbide heat exchangers’ samples, characterized their properties, and experimentally validated the thermal performance using streams of cold and hot air at temperatures varying from 42 to 250 °C. Gerstler and Erno [13] presented a multi-furcating heat exchanger printed by the method of direct metal laser melting with aluminum, Titanium 6–4, cobalt chrome, and Inconel 718. The AM prototype was tested and compared to a conventional heat exchanger for aviation applications, and the authors showed that the AM heat exchanger was able to achieve the same thermal performance as the conventional one with 50 and 66 percent reductions in volume and weight, respectively. This AM heat exchanger also has the added advantage of increased reliability due to the absence of braze joints. Zhang et al. [14] experimentally validated the ability to use a compact AM heat exchanger for a high-temperature application. The prototype was an Inconel 718 manifold heat exchanger that used nitrogen at 600 °C and air at 38 °C. A good agreement between numerical and experimental results was observed, along with a 25% heat transfer density enhancement compared to plate-fin heat exchanger. Methods of heat transfer enhancement through additive manufacturing have been investigated for various applications. El Achkar et al. [15] studied the performance of an additively manufactured plate-fin heat exchanger for phase change applications. Ning et al. [16] compared the performance of three additively manufactured compact heat exchangers for industrial applications, with different levels of design complexity. The study concluded that the highest heat transfer improvement was observed with the most complex design whose printing was made possible by additive manufacturing. This thermal improvement comes, however, with the penalty of pressure losses up to twice those of the conventional design. Kelly et al. [17] demonstrated the use of binder jet additive manufacturing for the fabrication of a triply periodic minimal surface (TPMS) heat exchanger with an ultra-high-temperature ceramic. The team estimated an achievable 50 MW/m3 power density for their MS-sCO2 heat exchanger. Tano et al. [18] described the design of a compact MS-sCO2 PHE for this application, fabricated using Haynes 282 and metal AM (Fig. 1(a)). The design consisted of multiple pairs of mm-scale pin arrays on the sCO2 (cold) side and parallel channels on the MS side. The sCO2 headers were located outside the MS path, resulting in a partial cross-flow region for sCO2 flow at the entrance and exit of the core (see Fig. 2(b)), and a counterflow region in the rest of the core. The feasibility of printing the features in the core, as well as surface roughness characterization, was shown. A correlation-based model was used to present a parametric study of the impact of flow and geometrical parameters on the performance and cost of the PHE. Results indicated that designs with an effectiveness of 90% and power densities of 10.8 MW/m3, including headers, can be obtained for an exit hot side temperature of 550 °C. Because of the low flowrates on the MS side (Re < 50), the PHE design was heat transfer limited on the MS side, requiring core lengths of 80 cm to obtain an effectiveness of 0.90 and a high power density of 11.4 MW/m3. Scalability was possible in the lateral direction by increasing the number of pairs of hot and cold flow paths; however, scalability was not possible in the vertical direction of the build due to the location of the sCO2 headers. The design by Tano et al. [18] will be henceforth referred to as the Gen 1 design.

Fig. 1
Overview of (a) the AM PHE design proposed by Tano et al. [18] (Gen 1 design), with rectangular channels on the Molten Salt side (MS) and (b) the Gen 2 design with a lattice of slanted pins on the molten salt side, (c) schematic of the pin fin array on the sCO2 side, (d) cut-away view of the Gen 1 design with rectangular channels on the MS side, (e) lattice of slanted pins on the MS side in the Gen 2 design, schematic of the cross section of (f) Gen 1 and (g) Gen 2 designs
Fig. 1
Overview of (a) the AM PHE design proposed by Tano et al. [18] (Gen 1 design), with rectangular channels on the Molten Salt side (MS) and (b) the Gen 2 design with a lattice of slanted pins on the molten salt side, (c) schematic of the pin fin array on the sCO2 side, (d) cut-away view of the Gen 1 design with rectangular channels on the MS side, (e) lattice of slanted pins on the MS side in the Gen 2 design, schematic of the cross section of (f) Gen 1 and (g) Gen 2 designs
Close modal
Fig. 2
Gen 1 PHE design features: (a) MS flow in rectangular channels, (b) sCO2 flow in pin array plate with cross-flow region at the entrance, (c) pin array schematic; Gen 2 design features (d) MS flow in 3D lattice, (e) velocity vectors illustrating flow path in the MS lattice, and (f) cross-view of sCO2 flow passages
Fig. 2
Gen 1 PHE design features: (a) MS flow in rectangular channels, (b) sCO2 flow in pin array plate with cross-flow region at the entrance, (c) pin array schematic; Gen 2 design features (d) MS flow in 3D lattice, (e) velocity vectors illustrating flow path in the MS lattice, and (f) cross-view of sCO2 flow passages
Close modal

2 Objectives

The objective of this work is to present a design that is more scalable and compact than the Gen 1 design by innovations in the header design and the MS channel flow path. The overall PHE design is first presented highlighting its unique attributes, followed by a description of the details of the header design. One of the key advantages of AM is that it permits an integrated header-core design that enables compact designs and the elimination of failure at the joint. Integrity of the header–core interface is demonstrated by pressure tests on a sub-scale PHE. A simplified thermofluidic model, developed in prior work [18], is used to assess the performance of the PHE core using correlations developed from experimental and computational fluid dynamics (CFD) simulations. A cost model, also developed in prior work [18], is used to assess the variation in cost based on these parameter variations. The PHE performance is discussed over a range of flow and geometric parameter variations. The impact of AM machine features and design on cost is assessed.

3 Design

A comparative view of the two designs is presented in Fig. 1 with the Gen 1 and 2 PHEs in Figs. 1(a) and 1(b), respectively. Both models are designed with nominal dimensions to fit the base plate of an EOS M290 (EOS GmbH) laser powder bed fusion printer. Details of the sCO2 and MS flows in the core are shown in Fig. 2. The modifications relative to the Gen 1 design shown in Figs. 1 and 2 address the low heat transfer rate on the MS side, which hinders the compactness of the PHE, and the scalability of the design in the vertical build direction. As shown in Figs 1(d) and 1(e), the rectangular channels on the MS side in the Gen 1 design were replaced by a 3D lattice of pins that provide additional surface area for heat transfer while increasing the residence time and mixing of the fluid in the channel. This design modification is intended to decrease the thermal resistance on the hot side of the PHE; however, the more intricate flow path in the lattice structure (see Fig. 2(e)) is anticipated to increase pressure drop when compared to the parallel channel path in the Gen 1 design shown in Fig. 1(d).

The second design change pertained to the headers, which are different in shape but also location. As illustrated in Fig. 2(d), the new sCO2 headers are aerodynamically shaped and integrated into the MS flow passages. By comparing Fig. 2(b) with Fig. 2(f), the less efficient cross-flow zones observed in Gen 1 are minimized in the Gen 2 PHE. Furthermore, the integration of the headers in the flow path allows for the scalability of the unit. The Gen 1 design (Fig. 1(a)) can only be scaled up by adding more plates in the depth (horizontal scalability) but with the Gen 2 architecture, it is now possible to achieve vertical scalability by building the stacks on top of the other. The overall length of a single Gen 2 stack, which conforms to the dimensions of the EOS M290 build plate, is 24 cm, including the headers. The sCO2 pin array plates are 18 cm in total length and include 2.3-cm-long diverging and converging sections near the headers followed by a 13.4-cm-long rectangular section that is 5 cm tall. Twenty-six such pin array plates are assembled in the width to make up the stack with a nominal rating of 20.7 kW under the current design conditions. The MS side channel width is 1.8 mm, with a lattice of slanted pin fins, tilted 45 deg in both the flow direction and the transverse direction to form the three-dimensional periodic flow network (see Fig. 2(e)).

Figure 3 shows a summary of the design methodology which includes structural, thermofluidic, and cost considerations. A unit cell of the core was used for mechanical simulations, while a correlation-based model was used for fluid flow and heat transfer. Headers for the sCO2 side were designed to ensure uniform flow distribution across each pin array and between the pin arrays, while also satisfying structural integrity. A process-based cost model was used to assess the impact of variables on cost.

Fig. 3
Summary of the design methodology
Fig. 3
Summary of the design methodology
Close modal

In order for the header design to be viable, the sCO2 headers would have to be made as compact as possible and of aerodynamic shape to reduce pressure drop on the MS side. The challenge with the compactness of the sCO2 header are threefold: (1) there needs to be sufficient wall thickness to withstand the 200 bar internal pressure, (2) the flow distribution among multiple sCO2 pin arrays needs to be uniform, and (3) the flow must be distributed uniformly across the lateral height within each pin array. Iterative simulations were performed between structural and fluid flow to determine the geometry of the header that would meet these requirements. The allowable stress in the design was restricted to a value of about 150 MPa, based on creep considerations using Haynes 282 [19] and the Larson Miller parameter approach [20] with a maximum temperature of 720 °C and a lifetime of 100,000 h. After conducting a static structural simulation of the header using linear periodic boundary conditions and applying the expected 200 bar pressure on the internal walls of the model, the final stress distribution showed a maximum equivalent stress of 168 MPa with most of the header subjected to stresses lower than 100 MPa.

Fluid flow simulations were performed in Ansys Fluent to ensure a uniform flow distribution from the sCO2 header into each pin array. The modeled geometry is shown in Fig. 4(a). To manage the computational complexity of incorporating multiple pin array flow paths of the PHE core into the header flow distribution analysis, a porous media approach was employed, wherein the pressure drop in the core was modeled using inertial and viscous terms [21] based on CFD simulations of the pin array. The red boxed area in Fig. 4(a) identifies the region of the model defined as porous media. CFD simulations were run for different flowrates on a representative stripe of the pin array plate. A correlation between velocity and pressure drop was obtained. From this analysis, it was possible to derive the inertial and viscous parameters establishing an equivalence between the fluid flow in the pin array plate and in a section of porous media [18]. Due to the presence of bends and recirculation regions in the wake of the flow distribution features at the entrance of each pin array (see Fig. 2(f)), the k-omega shear stress transport (SST) turbulence model was used. Figure 4(b) shows the flow uniformity in terms of mass flowrate in each plate (primary ordinate) and the percent difference between the local and average mass flowrates (secondary ordinate) for total incoming flowrates of 50, 75, and 100 g/s. For reference, the 50 g/s incoming flowrate corresponds to a Reynolds number of 86,160 based on the inlet diameter and a Reynolds number based on the pin diameter of 1314 for the pin array. For all the flowrates considered, the mass flowrate into all pin arrays of the core was within 5% except for the one closest to the inlet, for which it was about 6 percent lower than average.

Fig. 4
(a) Fluid domain for the study of flow distribution within the header, (b) flow distribution in the PHE plates for varied flowrates, (c) flow domain showing location of velocity sampling, and (d) flow distribution across the height of the core section in terms of total velocity vector magnitude (array dimensions such as pin diameter and spacing are constant at the baseline values in Table 5)
Fig. 4
(a) Fluid domain for the study of flow distribution within the header, (b) flow distribution in the PHE plates for varied flowrates, (c) flow domain showing location of velocity sampling, and (d) flow distribution across the height of the core section in terms of total velocity vector magnitude (array dimensions such as pin diameter and spacing are constant at the baseline values in Table 5)
Close modal

Upon verification of uniform flow into each sCO2 pin array, the design of each pin array for flow uniformity across its width was performed. The final iteration of the design that ensured structural integrity and uniform flow distribution is shown in Fig. 4(c). Elongated ribs near the entrance of the pin array were added to ensure uniform flow distribution across the width. The streamwise velocity magnitude across the width of the pin array is plotted in Fig. 4(d) at three different axial locations corresponding to x/L = 0.1, 0.25 and 0.5. The percent difference relative to the mean velocity at the respective location is shown in the secondary ordinate in Fig. 4(d). In the diverging section near the inlet, at location x/L = 0.1, some flow non-uniformity is observed, with regions of higher velocity in the center of the pin array. However, once the flow enters the straight rectangular section of the pin array plate, the uniformity in flow distribution improves significantly at all locations across the plate height such that the highest value percent difference drops to below 3% (see x/L = 0.25 and 0.5 locations), indicating uniform flow distribution exists for a majority of the core.

4 Header-Core Mechanical Integrity Verification

One of the stated advantages of AM is the ability to design and fabricate integral header and core architectures for heat exchangers. In traditionally manufactured heat exchangers such as PCHE, the core is welded or brazed to the headers, and this could lead to potential failure locations at the header-core location, in addition to inefficient heat exchange in cross-flow sections. To experimentally verify the mechanical integrity of the integrated header-core architecture of the AM PHE, sub-scale PHEs specimens, with three pairs of hot and cold plates and headers, were designed and fabricated for accelerated testing at temperature. Additively manufactured 316L SS was used instead of H282 in order to accelerate time to failure and obtain results in a timely manner. Since the Larson–Miller analysis is an accurate description of the material behavior of metallic systems, the choice of material was dictated by testing expediency and powder availability. The design was intentionally weakened in the core to determine whether failure would occur at this weakened region or the header-core intersection. Both location and time to failure were of interest. If the failure occurred at the weakened section in the core rather than in the header-core section, it would corroborate the finite element analysis model and confirm that the PHE can be designed for the appropriate life considering the high-stress regions. If the time to failure at operating temperature matched the predicted lifetime based on a Larson Miller approach, confidence in the printing parameters and AM approach would be obtained.

4.1 Design and Modeling for Mechanical Integrity Testing.

Two different models of the sub-scale PHE were designed with varying degrees of weakened core regions. The first sub-scale PHE model for structural analysis, along with the imposed boundary conditions, is shown in Fig. 5(a) and the second design is shown in Fig. 5(d). For the first model, an entire row of pins was removed from the core, while for the second design, every other pin in the same row was removed. For both designs, a symmetry boundary was imposed on the top half pin plate and a 200 bar internal pressure was imposed on the sCO2 pin arrays while atmospheric pressure was imposed on the MS side.

Fig. 5
Setup and boundary conditions for (a) Design 1 and (d) Design 2, stress distribution for (b) Design 1 and (e) Design 2, high-stress region for (c) Design 1 and (f) Design 2
Fig. 5
Setup and boundary conditions for (a) Design 1 and (d) Design 2, stress distribution for (b) Design 1 and (e) Design 2, high-stress region for (c) Design 1 and (f) Design 2
Close modal
In order to determine the time to failure, the creep-rupture data of printed SS316 generated by Li et al. [20] were used to generate a relationship between rupture stress, time to failure, and temperature using the Larson–Miller Parameter (LMP) approach
(1)
where T is in Kelvin and the time t in seconds.

The results of the structural simulations in ANSYS Mechanical for the first design are shown in Figs. 5(b) and 5(c) and life estimations at different locations of Design 1 and 2 are summarized in Table 1. For the first design, the peak and average stresses in the unit are approximately 510 and 48 MPa, respectively, and it is seen that regions of higher stress in the 170 MPa range (and the peak of 220 MPa) exist at the entrance of the pin array in between the elongated flow distribution structures (see Fig. 5(b)). The stress in the region of the core with the removed row of pins is also similarly high; the regions of high stress are located in the surrounding pins facing the open row of pins. The high-stress regions in the pins surrounding the open core region are shown in Fig. 5(c). Based on the peak stress of 510 MPa in the design, the expected time for failure for a pressure of 200 bar and a temperature of 743 °C (highest operational temperature recorded during the experiment) was estimated to be less than 0.1 s, while the estimate based on the global average stress yield a rupture time of 8471 min (141 h). Since the high stresses observed in the simulation will relax in practice during the experiment, it was expected that the correct metric to predict the time to rupture of the unit should be a value bracketed by the 48 and 510 MPa stress estimates. To further this analysis, the average surface stress in the missing row region of each plate was calculated to estimate the expected life.

Table 1

Estimated rupture time at different locations of Designs 1 and 2

Design 1Design 2
#Stress casesFailure stress (MPa)Life (minutes)Failure stress (MPa)Life (minutes)
1Maximum5100.000153760.026265
2Global average4884714210684
3Average in the weakened region (3 plates)772760546717
4Average surface stress—top plate150164113686
5Average surface stress—middle plate17758115635
6Average surface stress—bottom plate16785112713
Design 1Design 2
#Stress casesFailure stress (MPa)Life (minutes)Failure stress (MPa)Life (minutes)
1Maximum5100.000153760.026265
2Global average4884714210684
3Average in the weakened region (3 plates)772760546717
4Average surface stress—top plate150164113686
5Average surface stress—middle plate17758115635
6Average surface stress—bottom plate16785112713

The average stress in row #3 of Table 1 was calculated as the average equivalent stress in the weakened region. It includes both the stress at the surface of the pins and in the solid. As another measure of the local stress, the surface stress at the surface of the pins in the top, middle, and bottom pin arrays in the weakened regions are reported in rows 4–6 of Table 1. The average stress in this region was 77 MPa with local surface-averaged values of 150, 177, and 167 MPa for the stress-concentration region in the top, middle, and bottom plates, respectively. Given the location of the absolute maximum stress and the local peak in each plate, the expectation of failure for this design was in the core region with the missing row of pins. Among the three local results considered, the average surface stress in the middle and bottom plates gave the most realistic estimate of time to rupture of 58 and 85 min, which compare well with the experimental results described in the following section.

The stress distribution in the second design is shown in Fig. 5(e). As with the first design, it is seen that there are regions of high stress near the entrance to the pin array in between the elongated flow distribution structures and the peak stress in this region is similar to that in the first design. The high-stress regions in the core row with every other pin removed are shown in Fig. 5(f). As with the first design, the high-stress zones are in the pins surrounding the regions where the pins have been removed. The peak stress in this region is 376 MPa with an average of 115 MPa in the middle plate. Compared to the first design, the second design exhibits less stress variation from one plate to an adjacent one. Based on the average surface stresses in this region, the time to failure for a pressure of 200 bar and temperature of 743 °C was estimated to be 10.6 h.

The two designs were printed using SS316 in a Trumpf TruPrint 3000 laser powder bed fusion printer. The majority of the volume was printed with laser power, scanning velocity, hatch spacing, and layer thickness of 220 W, 700.0 mm/s, 150 μm, and 30 μm, respectively, with additional contour, up-skin, and down-skin laser scans. A porous, net-like structure was printed below the specimens to support the bottom surface of the specimens during printing.

4.2 Facility Description.

A picture of the pressure and temperature test facility is shown in Fig. 6. The facility consisted of a 500,000 BTU/h (165 kW) Maxon® burner (Honeywell Process Solutions, Phoenix, AZ) that was connected to a test chamber using 21-in. circular stainless steel ducting lined with high-temperature rigid insulation. The circular duct wall was equipped with an active cooling scheme consisting of water-cooled copper coils integrated around the ducting, a water pump, and a small air-cooled radiator to dissipate the heat from the water. The test chamber was 3 ft × 3 ft × 2 ft (width × height × depth) and was made with ¼” wall thickness steel plate with high-temperature insulation lining the inner walls. The chamber was mounted on a strut channel frame with locking castors. The back and top walls were adjustable to permit the placement of the test article to be tested and for exhausting the combustion gases, respectively. The test article to be pressure tested was placed within the chamber, and it was connected to a high-pressure nitrogen source with high-pressure-rated tubing integrated with an electronically controlled pressure regulator and an electronically controlled three-way solenoid valve assembly. The solenoid valve assembly could be used to isolate the test article from the nitrogen cylinder upon application of pressure. Furthermore, it could be used to relieve pressure and re-pressurize the test article during cyclic pressure testing. A proportional-integral-derivative (PID) controller on the burner was wired to the control station located in an adjacent office space container where the experiments were monitored remotely. The burner controller enabled the adjustment of the fuel–air ratio and hence the test chamber set temperature.

Fig. 6
(a) Experimental facility and (b) AM stainless steel unit before and after testing
Fig. 6
(a) Experimental facility and (b) AM stainless steel unit before and after testing
Close modal

During the pressure test, since the units were not stress-relieved, the test specimens were left attached to the build plate and sectioned apart as individual pieces for ease of testing. The experiments started by pressurizing the units at ambient temperature with a hold time at each set pressure. During the hold time, the unit was isolated from the pressure source to monitor the pressure for any decay indicative of a leak. All the above units passed the pressure test at ambient with 30 min hold time at the max pressure of 200 bar.

The elevated temperature testing started with the weakest design (Design 1), and a target temperature and pressure of 720 °C and 200 bar, respectively. The Design 1 unit was placed in the burner chamber, and five k-type thermocouples were inserted at different locations of the unit (as shown in Fig. 6). The time series of this unit body temperature (the average of all five thermocouples) and the unit internal pressure throughout the first round of testing are plotted in Figs. 7(a) and 7(b) shows a picture of the failed unit. First, the unit was heated up to the target temperature of 720 °C at atmospheric pressure. The unit was then pressurized in steps of 50 bar and 15 min hold time until the target value of 200 bar was achieved. The test proceeded until any potential failure sign could be noticed.

Fig. 7
(a) Time series of weakened unit (Design 1) body temperature and internal pressure throughout the test and (b) weakened sub-scale unit showing leak in the middle section after testing at high temperature and pressure. The internal features dimensions are kept at the baseline values in Table 5.
Fig. 7
(a) Time series of weakened unit (Design 1) body temperature and internal pressure throughout the test and (b) weakened sub-scale unit showing leak in the middle section after testing at high temperature and pressure. The internal features dimensions are kept at the baseline values in Table 5.
Close modal

The first unit withstood a cumulative testing time of around 150 min before failing, with the last 35 min being at the highest pressure of 200 bar. Failure was characterized by an abrupt decay in the pressure observed while the unit was still at the target temperature, which is an indication of a leak. The test was stopped, and after the cool-down process, the unit was removed from the chamber for inspection. No visible sign of defect was observed on the external walls, but after pressurizing the unit with compressed air (6.9 bar) and applying soapy water, the approximate location of the leak was revealed to be adjacent to the location of the missing rows at the center of the part and between the sCO2 plates (see Fig. 7(b)). The test was conducted a second time with a second identically printed sub-scale unit to ensure the repeatability of the experiment. The unit failed after about 57 min at the highest temperature and pressure. The leak test revealed that this unit failed at the core region and away from the header core. The second sub-scale design with alternate pins missing in two rows (Fig. 5(d)) was also pressure tested at temperature. Over a span of two days, the unit was under 200 bar pressure at 730 °C for over 11 h without any loss of pressure, at which point the test was concluded.

When compared with the simulation results, these experimental results reveal the following: (i) The prediction of failure location for the first design is consistent between the simulation and experiments and (ii) the predicted lifetime with the LMP approach using the maximum stress is an underestimate of the actual failure time while using the average stress is an overestimate. A local average of the surface stresses near the weakened location from the mechanical simulations seems to be the most accurate metric to predict the time to rupture. The experimental results also established confidence in the printing parameters and process.

5 Primary Heat Exchanger Core Model

Upon verification of the design of the header for flow uniformity and stresses in the prior sections, thermofluidic modeling and process-based cost modeling of the PHE core were performed. The correlations in the thermofluidic model were validated through pressure drop and heat transfer experiments on additively manufactured pin arrays for the cold side. Simulations were used to determine the flow and heat transfer characteristics on the hot side. The reliability of these correlations and simulations will also be validated in future work by conducting experiments on the full-scale PHE.

5.1 Thermofluidic Model.

The correlation-based thermofluidic model presented by Tano et al. [18] for the core design shown in Fig. 2(a) was modified to reflect the changes made to the design in the PHE core. The objective of this simplified numerical model was to evaluate the thermofluidic performance of the heat exchanger over a wide design space efficiently. The inlet temperature and pressure of the fluid streams were kept fixed at 720 °C and pressure of 1 bar for the molten salt, and 500 °C and 200 bar for the sCO2. As shown in Fig. 8, the PHE was divided into three main sections- diverging entrance section, a converging exit section, and the main counterflow region. The entrance and exit regions were not truly in cross-flow with MS; however, 5% of the length at the entrance and exit was modeled as cross-flow to account for a slight reduction in effectiveness in these sections relative to a pure counterflow design. The main counterflow section was further discretized to incorporate the effect of temperature-dependent properties.

Fig. 8
(a) Schematic of the Gen 2 PHE core and (b) differential energy balance for a control volume with half control volume on each side
Fig. 8
(a) Schematic of the Gen 2 PHE core and (b) differential energy balance for a control volume with half control volume on each side
Close modal
For the sake of brevity, the model details presented by Tano et al. [18] are omitted in the current work, and only salient modifications are discussed. The heat exchanger equations are derived for a half control volume on the hot and cold sides (see rectangular section shown in Fig. 8(b)), and the overall conductance and heat rates are multiplied by 2Nplates to determine the scaled PHE parameters. For the entrance and exit regions, the HX effectiveness of a cross-flow HX with both sides mixed was used [22]
(2)
The heat capacity rate ratio Cr in Eq. (2) was calculated as
(3)
where Cmin in Eq. (3) is the minimum of the hot and cold heat capacity rates, Ch and Cc, respectively.
The number of transfer units for the cross-flow section, NTUcr, in Eq. (2) was obtained as
(4)
in which the overall heat transfer conductance in the cross-flow section was estimated as
(5)
where Aw is the unfinned area. The factor of 2 in the above equation is required because the analysis is formulated for a control volume consisting of half hot and half cold channels with solid walls separating them.
Due to the complex geometry of the lattice structure on the MS side, the thermal resistance on this side was written in terms of thermal resistance rather than heat transfer coefficient. An analogous modification is implemented in the counterflow section and with the assumption of incompressible flow, the specific enthalpy i in Fig. 8 can be expressed in terms of the specific heat at constant pressure such that the temperature gradient on the hot side is given by
(6)

A similar equation can be written for the cold side. The thermal resistance to heat flux on the hot side Rh″, and the heat transfer coefficient on the cold side, obtained from CFD and experiments, respectively, are discussed in Sec. 5.2. The local temperatures on the hot and cold sides were obtained through an iterative process. First, a guess was made on the outlet sCO2 temperature, and the average temperatures on the hot and cold sides are obtained for the cross-flow 2 region. Next, a differential equation solver in matlab was used to obtain the temperature distribution in each control volume of the counterflow region, and finally, a procedure analogous to that used in the cross-flow 2 domain was employed to calculate the average temperatures in the cross-flow 1 section. At this stage, the inlet sCO2 temperature could be calculated, and the result was compared to the known inlet temperature. The iterative process was terminated when the two values were within the specified tolerance of 0.01 °C.

5.2 Pressure Drop and Heat Transfer Correlations

5.2.1 Cold (sCO2) Side.

Experimental testing of pin arrays with supercritical carbon dioxide: The heat transfer correlation used in the thermofluidic model on the sCO2 side was experimentally verified using the setup illustrated in Fig. 9. The design of AM pin array heat sink used for sCO2 heat transfer characterization is shown in Fig. 9(b). The pin array was laid out on a 5 cm × 10 cm footprint area connected to plena with fluidic ports to measure temperature and pressure on either side.

Fig. 9
(a) AM heat sink attached to the build plate, (b) AM cold-side pin array heat sink design, (c) heat sink and fluxmeter assembly, and (d) heat transfer test experimental setup line schematic
Fig. 9
(a) AM heat sink attached to the build plate, (b) AM cold-side pin array heat sink design, (c) heat sink and fluxmeter assembly, and (d) heat transfer test experimental setup line schematic
Close modal

The line schematic of the CO2 loop is depicted in Fig. 9(d). The main components of the test assembly included a tank containing liquid CO2, a circulating pump designed specifically for high-pressure liquid CO2 application and a preheater that converted liquid CO2 to a supercritical state and sensibly increased the temperature to a desired inlet temperature before it entered the heat sink which rested on top of a heat flux meter. This flux meter provided heat to the test section through a set of cartridge heaters inserted at the bottom and controlled by an AC variable transformer (see Fig. 9(c)). When the heated CO2 exited the heat sink, it was sent to a coil heat exchanger immersed in a glycol bath and cooled by a chiller. Upon becoming a liquid, the flow was directed through a Coriolis flowmeter and eventually back to the liquid CO2 tank located upstream of the CO2 pump. The fluxmeter and heat sink assembly were thoroughly insulated to minimize heat losses during the experiments.

The heat flux meter was made out of aluminum (6061 with material certificate) and had three sections separated by a 1.6 mm gap. Three thermocouples were embedded along the height of the heat flux meter in each section distanced 16 mm apart to measure the temperature gradient and calculate the heat flow toward the base of the pin array heat sink that was generated by the cartridge heaters. The AM pin array heat sink was designed with three thermocouple holes along the flow direction 2 mm below the base of the pin array (Fig. 9(b)) for surface temperature measurement. The fluid bulk temperature was monitored by thermocouples inserted into the side plena. The temperature data combined with the heat flux information were used to determine the heat transfer coefficient and compare it against the correlations used in the model
(7)
where q is the total amount of heat quantified by the heat flux meter, ΔT is the difference between the pin array base temperature (average of three measurements along the length) and the bulk mean fluid temperature, Asurf is the pin array footprint area, Npin is the total number of pins, and Dpin, Hpin, and ηpin are the pin fin diameter, pin fin height, and pin fin efficiency, respectively. During data reduction, the measurements from thermocouples inserted into the heat sink bottom wall were corrected to account for the heat conduction between the probe location in the wall and the surface in contact with the fluid. The corrected wall temperatures were then averaged to obtain a single value representative of the entire array. The bulk mean temperature was obtained by taking an average of the inlet and outlet fluid temperatures in the plena.
5.2.1.1 Results and comparison against pin array correlations.

The dimensions of two AM pin array heat sinks tested in this study are listed in Table 2. The pin fin rows, as shown in Fig. 2(c), were arranged in a staggered arrangement with respect to the flow direction. The non-dimensional heat transfer coefficient, Nusselt number, as a function of ReAmin for the flowrates considered in the experiment as well as a comparison with three correlations in the literature—Prasher et al. [23], vortex and non-vortex shedding correlations by Rasouli et al. [24] are shown in Fig. 10. The Prasher et al. [23] correlation was developed for micro pin fin heat sink with staggered circular and square pin fins, and the Rasouli et al. [24] correlations were established using PF-5060 and LN2 flows through diamond-shaped pin fin heat sinks (Table 3). The heat sink HS2 has dimensions designed so that vortex shedding was expected while HS1 has an aspect ratio that prevents the onset of vortex shedding at any location. As seen in Fig. 10, for the range of tested Reynolds numbers, the Nu number showed an increasing trend as predicted by the correlations, and for any given ReAmin, the HS2 (designed to develop vortex shedding) resulted in higher Nu compared with HS1.

Fig. 10
Nusselt number varying with Reynolds number for (a) HS1 and (b) HS2
Fig. 10
Nusselt number varying with Reynolds number for (a) HS1 and (b) HS2
Close modal
Table 2

AM Pin array heat sink geometry dimensions

ParametersHS1HS2
Dpin (mm)1.651.1
α (−)0.731.64
βT (−)1.492.43
βL (−)1.292.11
ParametersHS1HS2
Dpin (mm)1.651.1
α (−)0.731.64
βT (−)1.492.43
βL (−)1.292.11
Table 3

Literature correlations for single phase heat transfer and pressure drop in micro pin array heat sinks

ReferenceRemarks/Correlation
Prasher et al.a [23]
  • Micro pin fin heat sink

  • Staggered circular and square pin fin

  • Water

2.48 < α < 2.8
2.4 < βT = βL < 3.6
NuDh=0.132(βL1)0.256ReDh0.84(PrtestfluidPrwater)0.36 if ReDh < 100
NuDh=0.281(βL1)0.63ReDh0.73(PrtestfluidPrwater)0.36 if ReDh > 100
1.3 < α < 2.8
2.0 < βT = βL < 3.6
2.4 < βL < 4.0

f = 4 × 0.295(α)1.249(βL − 1)−0.7(βT − 1)−0.36Re−0.1 if ReDh > 100
Rasouli et al. [24]
  • Micro pin fin heat sink

  • Staggered diamond pin fin

  • LN2 and PF-5060

0.7 < α< 3.2
1.7 < βT < 3.0
0.8 < βL < 1.5

NuAmin,noVS=0.039(βT1)0.19ReAmin0.837Pr0.557

NuAmin,VS=0.253(α)0.39(βL1)0.062ReAmin0.747 for α > 1, βT > 2, and βL > 1
f=1373.6(βL1)0.16ReAmin1.3 if ReAmin < 100
f=9.2(α)0.43(βL1)0.07(βT1)0.07ReAmin0.15 if ReAmin > 100
Moores and Joshi [25]
  • Meso pin fin heat sink

  • Staggered circular pin fin

0.5 < α < 1.1
1.3 < βT < 1.36
1.13 < βL < 1.18
1000 < ReDh < 10000

f = 4 * 3.2(α)−0.138Re−0.42
Short et al.b [26]
  • Macro pin fin heat sink

  • Staggered circular pin fin

  • Air

1.9 < α < 7.2
2.0 < βT < 6.4
1.8 < βL < 3.2
ReDh > 1000
f=4×0.221(α)0.056(βL)1.4(βT)0.54ReDh0.08×(LDhNrow)
ReferenceRemarks/Correlation
Prasher et al.a [23]
  • Micro pin fin heat sink

  • Staggered circular and square pin fin

  • Water

2.48 < α < 2.8
2.4 < βT = βL < 3.6
NuDh=0.132(βL1)0.256ReDh0.84(PrtestfluidPrwater)0.36 if ReDh < 100
NuDh=0.281(βL1)0.63ReDh0.73(PrtestfluidPrwater)0.36 if ReDh > 100
1.3 < α < 2.8
2.0 < βT = βL < 3.6
2.4 < βL < 4.0

f = 4 × 0.295(α)1.249(βL − 1)−0.7(βT − 1)−0.36Re−0.1 if ReDh > 100
Rasouli et al. [24]
  • Micro pin fin heat sink

  • Staggered diamond pin fin

  • LN2 and PF-5060

0.7 < α< 3.2
1.7 < βT < 3.0
0.8 < βL < 1.5

NuAmin,noVS=0.039(βT1)0.19ReAmin0.837Pr0.557

NuAmin,VS=0.253(α)0.39(βL1)0.062ReAmin0.747 for α > 1, βT > 2, and βL > 1
f=1373.6(βL1)0.16ReAmin1.3 if ReAmin < 100
f=9.2(α)0.43(βL1)0.07(βT1)0.07ReAmin0.15 if ReAmin > 100
Moores and Joshi [25]
  • Meso pin fin heat sink

  • Staggered circular pin fin

0.5 < α < 1.1
1.3 < βT < 1.36
1.13 < βL < 1.18
1000 < ReDh < 10000

f = 4 * 3.2(α)−0.138Re−0.42
Short et al.b [26]
  • Macro pin fin heat sink

  • Staggered circular pin fin

  • Air

1.9 < α < 7.2
2.0 < βT < 6.4
1.8 < βL < 3.2
ReDh > 1000
f=4×0.221(α)0.056(βL)1.4(βT)0.54ReDh0.08×(LDhNrow)
a

The authors’ correlation does not have a Pr term. However, they used a Pr correction based on the study by Zukauskas [27], to compare their experimental data with literature. Based on the reported test conditions, water Pr was calculated at 52.5 °C.

b

The factor of L/DhNrow was applied because of the difference in the definition of f used by Short et al. [26] with this study.

Comparing experimental results with predicted values by the correlations, it can be seen that the Prasher et al. [23] correlation better predicts heat transfer in cases where vortex shedding is expected, while the Rasouli et al. [24] non-VS correlation is more accurate for non-vortex shedding scenarios.

The results of pressure drop experiments in non-dimensional form (friction factor) are presented in Fig. 11 along with predictions from various pin array correlations from the literature. The literature correlations are listed in Table 3. The experimental friction factor in HS1 and 2, in agreement with correlation predictions, has a decreasing trend with an increase in ReAmin. Compared with predicted values, it can be concluded that none of the correlations considered has an acceptable predictive capability, except the one presented by Rasouli et al. [24]. This correlation was developed based on results from eight-pin array heat sinks with diamond pin shape, machined on aluminum substrate with 2 cm × 2 cm footprint area. Seven of the eight samples had an aspect ratio between 1.5 and 3.2. The aspect ratio of HS2 falls in this range, and the friction factor is hence predicted with better accuracy (Fig. 11(b)) compared to HS1, whose aspect ratio is outside the range. The Rasouli et al. [24] correlation was therefore selected to predict the friction factor of the sCO2 in the thermofluidic model.

Fig. 11
Comparison of experimental friction factor against values predicted by available correlations in literature for (a) HS1 and (b) HS2
Fig. 11
Comparison of experimental friction factor against values predicted by available correlations in literature for (a) HS1 and (b) HS2
Close modal
As a summary, the correlations used for this section are the Prasher et al. [23] correlation for geometries in which the pin array was tall and sparse (α > 1, βT > 2)
(8a)
(8b)
for α < 1 or βL < 1 or βT < 2.
And the Rasouli et al. [24] correlation, chosen for geometries in which the pin array was short and dense (α < 1, βT < 1)
(9)

Rasouli et al. [24] noted that for the tall and sparse pin arrays, enhancement in heat transfer was observed due to vortex shedding downstream of the pins while this effect was suppressed for denser or shorter pin arrays.

The friction factor was predicted using the Rasouli et al. [24] correlations

(10a)
if ReAmin < 100
(10b)
if ReAmin > 100.

5.2.2 Molten Salt Side.

Due to the complex nature of the 3D lattice architecture on the hot side, CFD simulations were performed to obtain the thermal resistance for various Reynolds numbers. Figure 12 shows the simulation setup that consists of two computational domains—the fluid domain with MS and the solid domain of Haynes 282. A representative stripe of the array of slanted pin fins was modeled in lieu of the complete array to reduce the computational time. A constant heat flux of 40,000 W/m2 was imposed on the top and bottom faces based on the estimated thermal rating of the full-scale PHE. The side faces were defined as symmetry boundaries. The flowrate at the inlet was varied to obtain simulation results as a function of the Reynolds number. The molten salt, whose temperature-dependent properties were obtained from Wang et al. [28], is a composition of 20 mol% NaCl, 40 mol% KCl, and 40 mol% MgCl2. The simulations were performed using the k-omega SST turbulence model, and convergence was deemed achieved when all residuals were lower than 10−6 except the turbulent kinetic energy, which was of the order of 10−4.

Fig. 12
(a) CFD simulation setup, side view of (b) the molten salt and (c) H282 domains, with pin diameter of 1.2 mm and channel height of 1.0 mm
Fig. 12
(a) CFD simulation setup, side view of (b) the molten salt and (c) H282 domains, with pin diameter of 1.2 mm and channel height of 1.0 mm
Close modal
The results of friction factor and thermal resistance are presented in Figs. 13(a) and 13(b), respectively. Two channel heights, 1.0 and 1.8 mm, were considered in the analysis. The smaller height would enhance the heat transfer coefficient, while the taller channel height would ease the powder removal process and reduce pressure drop. A comparison between the two options would establish whether the heat transfer gain can justify the selection of a 1.0 mm channel. The thermal resistance R″ is calculated at several locations along the length of the stripe as
(11)
where Twavg is the average temperature of the top and bottom walls at the location considered and Tf is the average temperature of the fluid region defined by these top and bottom walls.
Fig. 13
(a) Friction factor and (b) thermal resistance varying with Reynolds number for the 1.8 and 1.0 mm hot (MS) channel height. The pin diameter, channel height, and wall thickness are constant at the baseline value specified in Table 5.
Fig. 13
(a) Friction factor and (b) thermal resistance varying with Reynolds number for the 1.8 and 1.0 mm hot (MS) channel height. The pin diameter, channel height, and wall thickness are constant at the baseline value specified in Table 5.
Close modal

For the two channel heights considered, a second-order polynomial was used to relate the thermal resistance to the Reynolds number, and a piecewise model was used for the friction factor as a function of the Reynolds number. The resulting correlations based on CFD simulations for two channel heights are listed in Table 4. Figure 13 shows that for both channel heights; the friction factor decreases sharply at low ReAmin and at a slower rate after a ReAmin of 12 and 17 for 1.8 and 1.0 mm channels, respectively. The friction factor increases by 43% from 35.7 for a 1.8 mm gap to 51.2 for 1.0 mm at ReAmin of 20. The percent increase drops to 32% at ReAmin of 35. The thermal resistance of the 1.8 mm channel is however higher, such that going from the 1.8 to 1.0 mm channel, the value of Rh decreases by 27% on average over all ReAmin. The resistance of the 1.8 mm channel decreases more rapidly with ReAmin than the other configuration such that the percent change from 1.8 mm to 1.0 mm channel is only 16% at a Reynolds number of 35.

Table 4

Friction factor and thermal resistance correlations for MS

1.8 mm MS lattice corefh = 146.72 − 7.76ReAmin if ReAmin < 12
fh=86.63.47ReAmin+0.05ReAmin2 otherwise
R2 = 1
R2 = 0.994
Rh=3.85×1045.90×106ReAmin+3.06×108ReAmin2 in K-m2/WR2 = 0.999
1.0 mm MS lattice corefh = 166.7 − 6.49ReAmin if ReAmin < 17
fh=99.73.07ReAmin+0.033ReAmin2 otherwise
R2 = 1
R2 = 0.994
Rh=2.71×1043.05×106ReAmin+1.45×108ReAmin2 in K-m2/WR2 = 0.999
1.8 mm MS lattice corefh = 146.72 − 7.76ReAmin if ReAmin < 12
fh=86.63.47ReAmin+0.05ReAmin2 otherwise
R2 = 1
R2 = 0.994
Rh=3.85×1045.90×106ReAmin+3.06×108ReAmin2 in K-m2/WR2 = 0.999
1.0 mm MS lattice corefh = 166.7 − 6.49ReAmin if ReAmin < 17
fh=99.73.07ReAmin+0.033ReAmin2 otherwise
R2 = 1
R2 = 0.994
Rh=2.71×1043.05×106ReAmin+1.45×108ReAmin2 in K-m2/WR2 = 0.999

5.3 Cost Model Description.

A process-based cost model (PBCM) was developed to estimate the heat exchanger cost (in $/kW-thermal or $/kW-th) for every variation of geometric and flow parameters considered [29]. The estimation process was performed by evaluating the constituent cost at every step of the production line including expenses associated with the material, equipment, labor, consumables, facility, utilities, and overhead. The model determines the overall cost after consideration of parameters for the PHE design, the facility operation, each step of the production and the individual expenses they generate. The PBCM consisted of the four following steps. First, a laser powder bed fusion is employed to print the heat exchanger layer by layer by melting metal powder under the action of a laser and allowing the melt pool to solidify according to the desired design. Next, a heat treatment step was implemented to relieve the stress generated in the printed element during the repeated melting-cooling sequences. The third step consisted of removing the base plate and supporting elements by means of a computer-controlled bandsaw. Finally, the heat exchanger was internally cleaned and smoothed by flowing an abrasive medium through the internal passageways. The cost estimate presented in Sec. 6 uses a AM printing parameter set of 700 mm/s velocity, 330 W power, and 0.11 mm hatch spacing. This conservative parameter set was selected because it consistently achieved high part density, low surface roughness, and good dimensional tolerance [18]. A discussion of increasing laser speed on the cost is discussed in the last section. A detailed description of the cost model is found in the study by Ziev et al. [29].

6 Primary Heat Exchanger Core Techno-Economic Parametric Study

Based on the thermofluidic and cost models described in Sec. 5, a parametric study was conducted on the geometrical and flow parameters of the PHE to identify tradeoffs and design geometries that would lead to high power densities. The parameters were varied sequentially, holding all others at a baseline value, to isolate their impact on the performance. The varied parameters and their baseline values are summarized in Table 5. The baseline values for the hot and cold channel width (Wc and Wh) are the result of a compromise between the fin efficiency and the channel gap that will minimize manufacturing-induced flow obstruction. The pin diameter (Dpin), transverse pitch ratio (βT), and wall thickness (tw) were set from simulations on mechanical stress. The range of surface roughness was based on measurements of as-printed parts (upper limit) and the potential for reduction in roughness with surface treatments such as abrasive flow machining or other chemical treatments (lower limit). While all parameters were varied, only results for the parameters with the most impact on performance are presented in this section.

Table 5

Fixed and varied parameters in PHE optimization

ParametersBaseline valueRange
LHX (cm)18[16.00, 35.00]
Nplates26[21, 39]
Wmax (cm)12[10.00, 18.00]
H (cm)5[5.00, 6.00]
Wh (mm)1.801.00, 1.80
Wc (mm)1.80[0.75, 2.50]
βT (−)2.05[1.34, 2.05]
Dpin_c (mm)1.20[0.50, 1.50]
Dpin_h (mm)1.20
tw (mm)0.50[0.30, 1.50]
Rz (μm)60[30.00, 250.00]
m˙sCO2 (g/s)80 (ReAmin = 2480)
Cmin on this side at MS baseline flowrates
m˙MS (g/s)115 (ReAmin = 23)[97, 141] (ReAmin = [18, 30])
ParametersBaseline valueRange
LHX (cm)18[16.00, 35.00]
Nplates26[21, 39]
Wmax (cm)12[10.00, 18.00]
H (cm)5[5.00, 6.00]
Wh (mm)1.801.00, 1.80
Wc (mm)1.80[0.75, 2.50]
βT (−)2.05[1.34, 2.05]
Dpin_c (mm)1.20[0.50, 1.50]
Dpin_h (mm)1.20
tw (mm)0.50[0.30, 1.50]
Rz (μm)60[30.00, 250.00]
m˙sCO2 (g/s)80 (ReAmin = 2480)
Cmin on this side at MS baseline flowrates
m˙MS (g/s)115 (ReAmin = 23)[97, 141] (ReAmin = [18, 30])

6.1 Flowrate and Molten Salt Channel Variations.

Figure 14 shows the impact of varying the flowrate of the MS (hot) stream relative to the sCO2 (cold) stream, represented as heat capacity rate ratio Cr on the performance of the PHE with the baseline geometry indicated in Table 5. The impact of varying Cr on exit fluid temperatures and effectiveness is shown in Fig. 14(a), while the pressure drop and cost (in $/kW) are shown in Fig. 14(b). The sCO2 flowrate is kept at 80 g/s, and the arrows on Fig. 14 show the direction of increasing the MS flowrate. The PHE with a 1.8 mm channel gap width on the hot side is represented by the filled symbols, and the 1.0 mm gap width PHE is depicted with hollow symbols. At an equal overall width of the PHE, the 1.0 mm gap will result in more plates and an increased heat transfer surface area compared to the 1.8-mm gap configuration. The PHEs with the two different channel gap widths show a similar trends in terms of effectiveness and outlet temperatures (Fig. 14(a)), and thermal rating (Fig. 14(b)). The 1.0-mm channel PHE has an effectiveness that varies from 0.81 to 0.90 compared to that from 0.79 to 0.89 for the 1.8-mm channel case. The better thermal performance of the 1.0 mm channel PHE is observed in the form of higher thermal rating (Fig. 14(b)) and power density. For example, at a Cr ratio of 0.93, the thermal rating and power density are 20.2 kW and 18.6 MW/m3 compared to 19.9 kW and 15.3 MW/m3, respectively, for the 1.8 mm case. This trend can be explained by the higher fin efficiency on the hot side for the 1.0 mm channel and the increased local velocity with narrower MS channels. The abrupt change in effectiveness observed at the two highest Cr ratios is due to a switch in Cmin along the PHE. As the MS flowrate increases, however, Cmin is consistently on the sCO2 side for the subsequent cases. Figure 14(a) shows that an increase in MS flowrate is accompanied by an increase in Tco because the thermal resistance on the hot side decreases (as illustrated in Fig. 13). However, the increased velocity results in a decrease in residence time, which results in an increase in Tho. A tradeoff exists between a high effectiveness and a lower Tho, which is a desired feature to minimize parasitic losses in the system.

Fig. 14
PHE Performance prediction as a function of Cr —the filled and hollow symbols represent a hot channel width of 1.8 mm and 1.0 mm, respectively (for reference, the volumetric power density at 20 kW corresponds to 15.4 MW/m3). The arrows indicate the direction of increasing MS flowrate, and all fixed parameters are at their baseline value specified in Table 5.
Fig. 14
PHE Performance prediction as a function of Cr —the filled and hollow symbols represent a hot channel width of 1.8 mm and 1.0 mm, respectively (for reference, the volumetric power density at 20 kW corresponds to 15.4 MW/m3). The arrows indicate the direction of increasing MS flowrate, and all fixed parameters are at their baseline value specified in Table 5.
Close modal

The two PHE designs are comparable in terms of pressure drop on the cold side (Fig. 14(b)); the slight difference observed (6760 Pa for 1.8 mm hot side versus 6783 Pa for 1.0 mm) results from the counteracting effects of lower flow per plate in the 1.0 mm PHE and a higher temperature and thus viscosity on the sCO2 side. The pressure drop on the hot side, however, is several times greater for the narrower channel. For example, at a Cr of 0.93, the pressure drop through the 1.0 mm gap width is 7 kPa compared to 1.36 kPa for the 1.8 mm gap width. At a constant hot channel gap width, the volume of the heat exchanger remains unchanged; hence, the trends in cost (in $/kW-th) reflect the power rating variations. As the MS mass flowrate increases from 97 to 141 g/s, the power rating increases from 19.5 to 21.5 kW translates into a cost reduction from $460/kW-th to $419/kW-th for the 1.8-mm hot channel PHE. In the 1.0 mm case, a cost reduction from $434/kW-th to approximately $400/kW-th is achieved with a power rating enhancement from 19.9 to 21.7 kW.

6.2 Core Length Variation.

Figure 15 shows the PHE performance as a function of LHX, the total length of the PHE core. The exit temperatures of the streams and effectiveness trends with core length are shown in Fig. 15(a) while the pressure drops, thermal rating, and cost (in $/kW-th) are shown in Fig. 15(b). As the length increases, there is more surface area and residence time for heat exchange between the two streams as shown by the increase in temperature difference with Tho decreasing from 552 to 542 °C and Tco increasing from 704 to 716 °C. The effectiveness increases 12% from 0.84 to 0.94 with core length. Although the thermal rating increases from 20.4 to 21.7 kW as LHX changes from 16 to 35 cm, the volumetric power density decreases from 17.3 to 9.3 MW/m3 due to the increase in core volume from 1.18 to 2.32 dm3. The pressure drop follows a linear increase from 1360 to 2824 Pa on the hot side and 6327 to 11,400 Pa on the cold side. The non-monotonic trend in cost can be explained in terms of the different growth rates of the power rating and PHE volume with LHX. As the length increases, so does the heat exchanger volume and the expenses related to material cost, and the other processing parameters. The additional cost is justified by the associated increase in power rating, as evidenced by the initial decrease in cost per kW. However, a point is reached at LHX of 20 cm, where the power rating gain comes with a larger increase in material volume, thereby causing a rise in $/kW-th cost beyond LHX of 20 cm from $434/kW-th for LHX of 20 cm to $445/kW-th for LHX of 36 cm.

Fig. 15
PHE Performance varying with the total length of the heat exchanger from 16 to 35 cm: (a) ɛ, Tco, and Tho variation and (b) $/kW-th, ΔPc, ΔPh, and thermal rating variation with LHX. sCO2 and MS flowrates of 80 and 115 g/s, respectively. All other parameters are constant at their respective baseline value specified in Table 5.
Fig. 15
PHE Performance varying with the total length of the heat exchanger from 16 to 35 cm: (a) ɛ, Tco, and Tho variation and (b) $/kW-th, ΔPc, ΔPh, and thermal rating variation with LHX. sCO2 and MS flowrates of 80 and 115 g/s, respectively. All other parameters are constant at their respective baseline value specified in Table 5.
Close modal

6.3 Total Width of Primary Heat Exchanger.

The impact of the overall width (number of plates) on the PHE performance is shown in Fig. 16. The exit temperatures of the streams and effectiveness trends with core width are shown in Fig. 16(a) while the pressure drop, thermal rating, and cost (in $/kW) are shown in Fig. 16(b). The residence time increases with the number of plates in the PHE, since the total flowrate gets subdivided into more streams, and the PHE effectiveness increases from 0.848 to 0.88. The thermal rating improves from 20.5 to 21.1 kW, but the additional pairs of hot and cold streams translate into a volume increase, which causes a drop in power density from 19 to 10.9 MW/m3. The pressure drop decreases on both hot and cold sides due to a decrease in mass flowrate through each plate in the core with an increase in the overall width while keeping the total MS and sCO2 flowrates constant. A cost reduction from $439/kW-th to $427/kW-th is observed because the increase in material cost associated with adding more plates is lower than the increase in power rating.

Fig. 16
PHE performance varying with the total width of the heat exchanger from 10 to 18 cm: (a) ɛ, Tco, and Tho variation and (b) $/kW-th, ΔPc, ΔPh, and thermal rating variation with Wmax. sCO2 and MS flowrates of 80 and 115 g/s, respectively. All other parameters are constant at their respective baseline value specified in Table 5.
Fig. 16
PHE performance varying with the total width of the heat exchanger from 10 to 18 cm: (a) ɛ, Tco, and Tho variation and (b) $/kW-th, ΔPc, ΔPh, and thermal rating variation with Wmax. sCO2 and MS flowrates of 80 and 115 g/s, respectively. All other parameters are constant at their respective baseline value specified in Table 5.
Close modal

6.4 Roughness.

Figure 17 illustrates the influence of the peak-to-valley roughness Rz on the thermal performance of the heat exchanger. The Rz range was selected to span the values obtained from sample measurements on the sCO2 side, before and after abrasive flow machining, and was modeled so that Rz effectively results in a reduction in the height of the pin array on the cold side. Due to the complex nature of the flow path through the 3D lattice on the MS side, the impact of surface roughness was not modeled on this side. The reduction in the minimum cross-sectional flow area Amin due to the channel constriction by surface roughness appears to have a prominent effect on the thermal performance by increasing the heat transfer coefficient on the cold side and reducing the thermal resistance. As a result, the effectiveness increases from 0.85 to 0.88 and the rating from 20.6 to 21.1 kW. The pressure drop on the cold side, however, increases significantly from 6.4 kPa to 11.8 kPa due to the constriction in the open channel area.

Fig. 17
PHE performance varying with the cold side roughness from 30 to 250 µm: (a) ɛ, Tco, and Tho variation and (b) $/kW-th, ΔPc, ΔPh, and thermal rating variation with. sCO2 and MS flowrates of 80 and 115 g/s, respectively. All other parameters are constant at their respective baseline values specified in Table 5.
Fig. 17
PHE performance varying with the cold side roughness from 30 to 250 µm: (a) ɛ, Tco, and Tho variation and (b) $/kW-th, ΔPc, ΔPh, and thermal rating variation with. sCO2 and MS flowrates of 80 and 115 g/s, respectively. All other parameters are constant at their respective baseline values specified in Table 5.
Close modal

The effect of various parameters was considered in this study and based on the findings, it can be suggested that the most effective way to combine the advantages of high power density, low pressure drop and high effectiveness is by increasing the number of plates (larger Wmax). Increasing the length of the PHE does increase the thermal performance of the PHE and lower the cost, but only up to the point when the increase in material volume and associated manufacturing costs match the gain in thermal performance, after which there is only a marginal increase in thermal rating and the trend in cost (in $/kW-th) reverses. The effect of increasing the flow length in the PHE has the adverse effect of increasing the pressure drop whereas an increase in the number of plates, at a constant mass flowrate lowers the pressure losses. The effect of surface roughness is less straightforward to incorporate at the design stage because roughness is not a parameter that can easily be monitored and predicted before printing. Further study is required to formulate accurate design guidelines, but a reasonable operation window can be specified based on the average size of the powder and available surface smoothing methods.

6.5 Pathways to Additive Manufacturing Cost Reduction.

In this section, we explore the impact of changing PHE design and printing process parameters on cost. The cost in $/kW-th presented in the parametric study is based on the conservative set of parameters and a single-laser machine (print speed of 700 mm/s). However, high part density is obtained at higher print speeds and laser powers. A porosity map was generated from a build using H282 with various speeds and laser powers (see Fig. 18) with the motivation of operating within a “process window” for which low porosity is maintained regardless of small variations in local process conditions. Fourteen cube specimens printed with various speeds and laser powers were cross-sectioned, and around 20 to 50 optical micrographs were each taken from different parts of the cubes. The porosities which showed up as dark regions in the micrographs were segmented out, and the amount of bright region in the set of micrographs collected from a given cube was regarded as the part density for the given set process parameters. From the map, it is seen that very low porosities of 0.03%, obtained at the conservative printing speed of 700 mm/s can also be obtained at higher speeds of 959 mm/s and 1366 mm/s. Figure 19 shows a cost comparison between Gen 1 [18] and Gen 2 PHEs for two printing speeds (959 and 1366 mm/s) and three commercially available machines: the EOS M290 which is a single-laser machine currently used for the printing of the articles tested in this study, the EOS M400–4 machine with four lasers, and the SLMNGXII600 twelve-laser machine equipped with a system that reduces the setup/teardown and heat-up/cool-down time.

Fig. 18
Porosity of H282 parts printed via EOS M290 LPBF plotted on the process parameter map
Fig. 18
Porosity of H282 parts printed via EOS M290 LPBF plotted on the process parameter map
Close modal
Fig. 19
Cost in $/kW for (a) Gen 1 and (b) Gen 2 on three different machines and two printing speeds, 959 and 1366 mm/s, based on an annual production volume of 1500 units
Fig. 19
Cost in $/kW for (a) Gen 1 and (b) Gen 2 on three different machines and two printing speeds, 959 and 1366 mm/s, based on an annual production volume of 1500 units
Close modal

The results reveal that switching to the Gen 2 PHE design results in more than two-fold cost reduction for all the printers investigated. The cost reductions are primarily driven by reduced material used per kW and reduced printing costs per kW. The reduced printing costs per kW are related to the reduced material used because melting less material requires less time. Additionally, printing at higher speeds and with more advanced printers that have more lasers reduce printing costs by reducing the time required to melt the material. These reductions in cost show diminishing returns with respect to printing speed because the print time-related costs only account for a portion of the total cost. Costs associated with machine setup, pre-heating, cooldown and part removal time and post-processing costs are not affected by the reducing in printing time. It is seen that the Gen 2 design reduces the PHE cost by 2.8 times for the baseline print speed of 959 mm/s. A further 22 percent decrease in cost, to $185/kW-th, can be obtained by using the twelve-laser machine.

7 Conclusions

Design of a compact, scalable additively manufactured (AM) primary heat exchanger for solar thermal applications was presented. The design is compared against a recent design in literature from Tano et al. [18], referred to as Gen 1 design. The key design novelties include (a) a lattice structure in the MS flow, which increases the heat transfer rate on the side with the largest thermal resistance, and (b) incorporation of sCO2 headers within the MS path. The latter results in a near-counterflow arrangement for the PHE and a design that can be scaled in two dimensions.

The structural integrity of the header-core junction was assessed and validated by conducting experiments at a pressure of 200 bar and a temperature of 730–743 °C on AM 316 stainless steel printed sub-scale specimens. A comparison between the structural simulations and the experimental results revealed that the maximum stress used as design parameter is quite conservative. A more realistic value to consider is an average surface stress localized around the stress-concentration zone.

Thermofluidic and process-based cost models were used to determine the impact of variation of geometric and flow parameters on the performance and cost of the PHE. The results show that a factor of 2.8 times reduction in cost can be achieved by the Gen 2 design, with power densities as high as 18.6 MW/m3 (including headers) for an effectiveness of 0.88.

Pathways to cost reduction by increase in print speed and industrial-scale AM machines were presented. The largest cost reduction is achieved by the design changes, followed by increase in print speed from 700 mm/s to 959 mm/s. While the presented design was for MS to sCO2, the design can be readily adapted for solid particle medium to sCO2 heat exchangers, with particles flowing through the lattice channels.

Acknowledgment

This material is based upon work supported by the U.S. Department of Energy's Office of Energy Efficiency and Renewable Energy (EERE) under the Solar Energy Technologies Office Award Number DE-EE0008536. Simulations of the unit cell core were performed using computational resources sponsored by the Department of Energy's Office of Energy Efficiency and Renewable Energy and located at the National Renewable Energy Laboratory. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The data sets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

Nomenclature

f =

Darcy friction factor

h =

average heat transfer coefficient (W/m2-K)

t =

thickness (m)

C =

heat capacity rate (W/K)

H =

height (m)

L =

length (m)

m˙ =

mass flowrate (kg/s)

R″ =

thermal resistance (K-m2/W)

cp =

specific heat capacity (J/kg-K)

kw =

wall conductivity (W/m-K)

Acr =

footprint of cross-flow region, LcrH (m2)

Amin =

minimum flow area within the pin array (m2)

Apin =

pin half surface area, =0.5πDpinHpin (m2)

Asurf =

pin array footprint area (m2)

Aw =

wall surface area not including the pins (m2)

Cr =

heat capacity rate ratio

DAmin =

heat sink hydraulic diameter based on Amin and Pmin (m)

Dh =

hydraulic diameter based on pin fin size (m)

Dpin =

pin fin diameter (m)

Pmin =

wetted perimeter associated to Amin (m)

Nplates =

number of plates in the heat exchanger

Npin =

number of pin fins in the array

Nrow =

number of pin fins rows in the flow direction

Rz =

peak-to-valley roughness (µm)

SL =

longitudinal pitch (m)

ST =

transverse pitch (m)

Wmax =

overall width of heat exchanger core (m)

AM =

additive manufacturing

H282 =

Haynes 282

HX =

heat exchanger

IN740H =

Inconel 740H

MS =

molten salt

NuAmin =

average Nusselt number based on DAmin

NuDh =

average Nusselt number based on hydraulic diameter

NTU =

number of transfer units

PCHE =

printed circuit heat exchanger

PHE =

primary heat exchanger

Pr =

Prandtl number

ReAmin =

Reynolds number based on DAmin

ReDh =

Reynolds number based on hydraulic diameter

sCO2 =

supercritical carbon dioxide

UA =

overall heat transfer conductance (W/K)

Greek Symbols

α =

aspect ratio, =Hpin/Dpin

β =

pitch ratio, =SL/T/Dpin

ɛ =

heat exchanger effectiveness

  ηpin =

fin efficiency

Subscripts

c =

cold

co =

cold outlet

cr =

cross-flow

h =

hot

ho =

hot outlet

versus =

vortex shedding

w =

wall

References

1.
Mehos
,
M.
,
Turchi
,
C.
,
Vidal
,
J.
,
Wagner
,
M.
,
Ma
,
Z.
,
Ho
,
C.
,
Kolb
,
W.
,
Andraka
,
C.
, and
Kruizenga
,
A.
,
2017
, “
Concentrating Solar Power Gen3 Demonstration Roadmap
,”
National Renewable Energy Laboratory
, Report No. NREL/Tp-5500-67464.
2.
Achkari
,
O.
, and
El Fadar
,
A.
,
2020
, “
Latest Developments on TES and CSP Technologies—Energy and Environmental Issues, Applications and Research Trends
,”
Appl. Therm. Eng
,
167
, p.
114806
.
3.
Azouzoute
,
A.
,
Alami Merrouni
,
A.
, and
Touili
,
S.
,
2020
, “
Overview of the Integration of CSP as an Alternative Energy Source in the MENA Region
,”
Energy Strategy Rev.
,
29
, p.
100493
.
4.
Romero
,
M.
,
Buck
,
R.
, and
Pacheco
,
J. E.
,
2002
, “
An Update on Solar Central Receiver Systems, Projects, and Technologies
,”
ASME J. Sol. Energy Eng.
,
124
(
2
), pp.
98
108
.
5.
Aakre
,
S.
, and
Anderson
,
M.
,
2022
, “
Pressure Drop and Heat Transfer Characteristics of Nitrate Salt and Supercritical CO2 in a Diffusion-Bonded Heat Exchanger
,”
Int. J. Heat Mass Transfer
,
189
, p.
122691
.
6.
Yao
,
Y.
,
Ding
,
J.
,
Zhang
,
Y.
,
Wang
,
W.
, and
Lu
,
J.
,
2022
, “
Heat Transfer Performance of Pillow Plate Heat Exchanger With Molten Salt and Supercritical Carbon Dioxide
,”
Int. J. Heat Mass Transfer
,
183
(
Part C
), p.
122211
.
7.
Shi
,
H. Y.
,
Li
,
M.-J.
,
Wang
,
W.-Q.
,
Qiu
,
Y.
, and
Tao
,
W.-Q.
,
2020
, “
Heat Transfer and Friction of Molten Salt and Supercritical CO2 Flowing in an Airfoil Channel of a Printed Circuit Heat Exchanger
,”
Int. J. Heat Mass Transfer
,
150
, p.
119006
.
8.
Wang
,
W.-Q.
,
Qiu
,
Y.
,
He
,
Y.-L.
, and
Shi
,
H.-Y.
,
2019
, “
Experimental Study on the Heat Transfer Performance of a Molten-Salt Printed Circuit Heat Exchanger With Airfoil Fins for Concentrating Solar Power
,”
Int. J. Heat Mass Transfer
,
135
, pp.
837
846
.
9.
Caccia
,
M.
,
Tabandeh-Khorshid
,
M.
,
Itskos
,
G.
,
Strayer
,
A. R.
,
Caldwell
,
A. S.
,
Pidaparti
,
S.
,
Singnisai
,
S.
, et al
,
2018
, “
Ceramic–Metal Composites for Heat Exchangers in Concentrated Solar Power Plants
,”
Nature
,
562
, pp.
406
409
.
10.
Montes
,
M. J.
,
Linares
,
J. I.
,
Barbero
,
R.
, and
Rovira
,
A.
,
2020
, “
Proposal of a New Design of Source Heat Exchanger for the Technical Feasibility of Solar Thermal Plants Coupled to Supercritical Power Cycles
,”
Sol. Energy
,
211
, pp.
1027
1041
.
11.
Singh
,
D.
,
Yu
,
W.
,
France
,
D.
,
Allred
,
T.
,
Liu
,
I.-H.
,
Du
,
W.
,
Barua
,
B.
, and
Messner
,
M.
,
2020
, “
One Piece Ceramic Heat Exchanger for Concentrating Solar Power Electric Plants
,”
Renewable Energy
,
160
, pp.
1308
1315
.
12.
Du
,
W.
,
Yu
,
W.
,
France
,
D.
,
Singh
,
M.
, and
Singh
,
D.
,
2022
, “
Additive Manufacturing and Testing of a Ceramic Heat Exchanger for High-Temperature and High-Pressure Applications for Concentrating Solar Power
,”
Sol. Energy
,
236
, pp.
654
665
.
13.
Gerstler
,
W.
, and
Erno
,
D.
,
2017
, “
Introduction of an Additively Manufactured Multi-Furcating Heat Exchanger
,”
Proceedings of the 16th IEEE Intersociety Conference on Thermal and Thermomechanical Phenomena in Electronic Systems (ITherm)
,
Orlando, FL
,
May 30–June 2
,
IEEE
, pp.
624
633
.
14.
Zhang
,
X.
,
Tiwari
,
R.
,
Shooshtari
,
A.
, and
Ohadi
,
M.
,
2018
, “
An Additively Manufactured Metallic Manifold-Microchannel Heat Exchanger for High Temperature Applications
,”
Appl. Therm. Eng.
,
143
, pp.
899
908
.
15.
El Achkar
,
G.
,
Septet
,
C.
,
Le Metayer
,
O.
, and
Hugo
,
J.-M.
,
2022
, “
Experimental Thermohydraulic Characterisation of Flow Boiling and Condensation in Additive Manufactured Plate-Fin Heat Exchanger
,”
Int. J. Heat Mass Transfer
,
199
, p.
Article 123465
.
16.
Ning
,
J.
,
Wang
,
X.
,
Sun
,
Y.
,
Zheng
,
C.
,
Zhang
,
S.
,
Zhao
,
X.
,
Liu
,
C.
, and
Yan
,
W.
,
2022
, “
Experimental and Numerical Investigation of Additively Manufactured Novel Compact Plate-fin Heat Exchanger
,”
Int. J. Heat Mass Transfer
,
190
, p.
122818
.
17.
Kelly
,
J. P.
,
Finkenauer
,
L. R.
,
Roy
,
P.
,
Stolaroff
,
J. K.
,
Nguyen
,
D. T.
,
Ross
,
M. S.
,
Hoff
,
A. T.
, and
Haslam
,
J.
,
2022
, “
Binder Jet Additive Manufacturing of Ceramic Heat Exchangers for Concentrating Solar Power Applications With Thermal Energy Storage in Molten Chlorides
,”
Addit. Manuf.
,
56
, p.
102937
.
18.
Tano
,
I.-N.
,
Rasouli
,
E.
,
Ziev
,
T.
,
Wu
,
Z.
,
Lamprinakos
,
N.
,
Seo
,
J.
,
Schulze Balhorn
,
L.
,
Vaishnav
,
P.
,
Rollett
,
A.
, and
Narayanan
,
V.
,
2021
, “
An Additively-Manufactured Molten Salt-to-Supercritical Carbon Di-Oxide Primary Heat Exchanger For Solar Thermal Power Generation—Design And Techno-Economic Performance
,”
Sol. Energy
,
234
, pp.
152
169
.
20.
Li
,
M.
,
Zhang
,
X.
,
Chen
,
W.-Y.
, and
Heidet
,
F.
,
2020
, “
Progress Report on the Assessment of the Material Performance for TCR Applications
,” Report No. ANL/NSE-20/42.
21.
Ansys Fluent
,
Release 2023 R1, Help System, User Inputs for Porous Media, ANSYS, Inc
., http://www.ansys.com/academic/terms-and-conditions.
22.
Nellis
,
G.
, and
Klein
,
S.
,
2009
,
Heat Transfer
,
Cambridge University Press
,
New York
.
23.
Prasher
,
R. S.
,
Dirner
,
J.
,
Chang
,
J. Y.
,
Myers
,
A.
,
Chau
,
D.
,
He
,
D.
, and
Prstic
,
S.
,
2007
, “
Nusselt Number and Friction Factor of Staggered Arrays of Low Aspect Micropin–Fins Under Cross Flow for Water as Fluid
,”
ASME J. Heat Transfer
,
129
(
2
), pp.
141
153
.
24.
Rasouli
,
E.
,
Naderi
,
C.
, and
Narayanan
,
V.
,
2018
, “
Pitch and Aspect Ratios Effects on Single-Phase Heat Transfer Through Microscale Pin Fin Heat Sinks
,”
Int. J. Heat Mass Transfer
,
118
, pp.
416
428
.
25.
Moores
,
K. A.
, and
Joshi
,
Y. K.
,
2003
, “
Effect of Tip Clearance on the Thermal and Hydrodynamic Performance of a Shrouded Pin Fin Array
,”
ASME J. Heat Transfer
,
125
(
6
), pp.
999
1006
.
26.
Short
,
B. E.
,
Raad
,
P. E.
, and
Price
,
D. C.
,
2002
, “
Performance of Pin Fin Cast Aluminum Cold Walls, Part. 1: Friction Factor Correlations
,”
J. Thermophys. and Heat Transfer
,
16
(
3
), pp.
389
396
.
27.
Žukauskas
,
A.
,
1972
, “
Heat Transfer From Tubes in Crossflow
,”
Adv. Heat Transfer
,
8
(
1972
), pp.
93
160
.
28.
Wang
,
X.
,
Del Rincon
,
J.
,
Li
,
P.
,
Zhao
,
Y.
, and
Vidal
,
J.
,
2021
, “
Thermophysical Properties Experimentally Tested for NaCl-KCl-MgCl2 Eutectic Molten Salt as a Next-Generation High-Temperature Heat Transfer Fluids in Concentrated Solar Power Systems
,”
ASME J. Sol. Energy Eng.
,
143
(
4
), p.
041005
.
29.
Ziev
,
T.
,
Rasouli
,
E.
,
Tano
,
I.-N.
,
Wu
,
Z.
,
Yarasi
,
S. R.
,
Lamprinakos
,
N.
,
Seo
,
J.
,
Narayanan
,
V.
,
Rollett
,
A. D.
, and
Vaishnav
,
P.
,
2023
, “
Cost of Using Laser Powder Bed Fusion to Fabricate a Molten Salt-to-Supercritical Carbon Dioxide Heat Exchanger for Concentrating Solar Power
,”
3D Print. Addit. Manuf.
,
10
(
3
).