Next Article in Journal
Effect of Centrifugal Force on Power Output of a Spin-Coated Poly(Vinylidene Fluoride-Trifluoroethylene)-Based Piezoelectric Nanogenerator
Next Article in Special Issue
Heat Transfer and Heat Recovery Systems
Previous Article in Journal
Experimental Investigation on the Deformation Characteristics and Mechanical Behaviors of Tectonic Coal under Complex Unloading Confining Pressure
Previous Article in Special Issue
Heat Integration for Phenols and Ammonia Recovery Process of Coal Gasification Wastewater Considering Optimization of Process Parameters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Blanket Cooling of a Fusion Reactor

Competence Center for Thermal Energy Storage (CCTES), Lucerne University of Applied Sciences and Arts, Technikumstrasse 21, 6048 Horw, Switzerland
*
Author to whom correspondence should be addressed.
Energies 2023, 16(4), 1890; https://doi.org/10.3390/en16041890
Submission received: 13 January 2023 / Revised: 1 February 2023 / Accepted: 10 February 2023 / Published: 14 February 2023
(This article belongs to the Special Issue Volume II: Heat Transfer and Heat Recovery Systems)

Abstract

:
Nuclear fusion is the gateway to a whole new paradigm of energy and is a strong candidate for the decarbonization of electricity generation on a global scale. With recent developments in high-temperature super-conducting magnets, the race is on to develop sub-systems which will support a commercially viable fusion reactor for use as a thermal power plant. The fusion of lighter elements creates an enormous amount of heat which must be transferred away from the reactor core. These intense conditions require novel approaches to efficiently transfer very high heat loads into useable thermal energy without compromising the structural integrity of the reactor core and the surrounding components. This report outlines the concept of a fundamental approach to solve the heat transfer problem as proposed by Commonwealth Fusion System’s design for a fusion reactor. A literature review was conducted for other applications that could serve as inspirations, as well as material properties and machining methods for the proposed power exhaust system. A dive into the theoretical thermodynamic and fluid dynamic characteristics of plate heat exchangers and finned surfaces was conducted from a fundamental perspective. A laminar flow regime was studied for the purpose of setting the floor for energy needed to pump coolant while simultaneously representing the least favorable heat transfer regime between a solid surface and a fluid. The results served as a basis for dimensioning and executing numerical simulations as a means for a first look into a solution of this heat transfer problem. The results were compared with the theoretical conclusions and judged based on constraints of the system. Recommendations were made for the continued development of a corresponding system.

1. Introduction

Energy from the nuclear fusion of two lighter elements represents the next paradigm in energy generation. It represents the most common type of observable reaction in the universe, the enormous amount of energy released by stars, and the incredible potential it could provide to humanity [1,2,3]. These intense reactions create extremely high specific heat fluxes. The Competency Center Thermal Energy Storage (CCTES) is conducting a study based on principle ideas created by Commonwealth Fusion Systems (CFS) to develop a concept for a part of the power exhaust system (PEXS) necessary for the removal of heat from a fusion reactor core. In a former CFS publication, the basic engineering principles for a blanket cooling system within an affordable compact robust (ARC) reactor were described. However, a design proposal for heat exchange surfaces and heat transfer fluid has not been published yet [4]. A reactor core consists of fusion fuels reacting in a magnetically confined vacuum vessel, requiring cooling. Vacuum vessel heat flux density reaches a peak of 12 MW/m2, which needs to be removed via FLiBe (Fluoride–Lithium–Beryllium) molten (liquid) salt to maintain material integrity. The intention of this work is to suggest and examine one potential solution approach: a novel design proposal for removing high-end vacuum vessel heat fluxes using a FLiBe molten salt heat transfer fluid. The design proposal is based on data from CFS and former research work at CCTES. The design proposal includes the incorporation of manufacturing possibilities and a general design to keep the inner (thin) shell mechanically stable, to allow for high heat transfer coefficients by forced convection (or even free convection), and for mass transfer exchange with the molten salt reservoir (for the transport of the heat to the bulk FLiBe tank). A general layout of the design is developed, general assumptions of the required flow rates of FLiBe and a possible way to evenly distribute this over all 300 m2 of the blanket surface are shown, the blanket surface is structured with a grid consisting of fins, and a detailed numerical calculation of the local heat transfer under such conditions is executed. Starting with a parameter study of the influence channel dimensions and fin width, the general feasibility in terms of heat transfer and pressure drop for FLiBe and laminar flow is examined.

2. Materials, Theory, and Methods

2.1. Materials

Three main materials are needed to construct a heat exchanger or PEXS. Tungsten (W, Wolfram) is the proposed material for the plasma-facing surface responsible for transferring heat to the molten FLiBe salt in the ARC reactor design. It has a high melting temperature of 3400 °C, low vapor pressure, and high tensile strength [5,6]. Tungsten is the material of choice in several divertor designs, including an ITER (International Thermonuclear Experimental Reactor), for its high thermal conductivity (kW ~175 W/(m·K)), high temperature operating range (~1500 K), and resistance to plasma erosion [4,7,8]. It is proposed to use a plasma-facing layer that is as thin as possible to minimize thermal resistance while maintaining structural integrity during the operating life cycle.
Above the Tungsten layer is the fluid responsible for capturing and removing heat from the plasma-facing surface. The fluid proposed by the CFS paper [4] is FLiBe (LiF—BeF2; 66 mol%, 34 mol%). FLiBe has several engineering advantages for use as a coolant in a magnetically confined plasma fusion reactor. It exhibits a favorable operating temperature range between 800 and 875 K, performs well as a radiation shield, and has favorable flow characteristics [9,10,11]. For the proposed FLiBe operating temperature (800~875 K), the dynamic viscosity of FLiBe is approximately 10 times higher than water at room temperature (293 K). Table 1 outlines the recommended values for FLiBe over a liquid temperature range between melting (732 K) and boiling (1703 K) points.
The design proposed in this paper was a super-structure of Inconel 718 tasked with distributing FLiBe to the tungsten layer and providing structural support to the “cool” side of the PEXS. It operates at a lower temperature than plasma-facing surfaces and provides a structural level of support to the inner tungsten layer. Inconel was chosen for its high melting temperature (~1480–1610 K), very high yield and tensile strength, and excellent weldability with itself and tungsten [14,15,16]. Figure 1 below illustrates the conceptual design of the channeled tungsten plasma-facing surface (orange in the back) attached to a second set of Inconel 718 channels (dark gray). FLiBe capillary channels are shown over the top of the entire construction.

2.2. Theory

A full understanding of heat transfer and the effects of material extended away from a surface is requisite to properly evaluating designs which best conform to boundary conditions. Table 2 below shows the boundary conditions of the PEXS which must not be exceeded. This set of boundary conditions is considered nominal operating conditions, as outlined by CFS.

2.2.1. Heat Transfer

Heat from the inner core to the tungsten shell is transferred mainly by radiation. The following theory assumes the maximum tungsten temperature will be maintained by an equilibrium of radiation from the core to tungsten and convective heat transfer from tungsten to FLiBe. The boundary conditions were conceptualized in terms of heat transfer using the general theory for plate heat exchangers. Figure 2 below qualitatively shows the evolution of the FLiBe temperature as it passes through the tungsten structure. It is governed by the following equation:
Q ˙ = m ˙ C p T o u t T i n = h A Δ T l m
Q ˙ is the heat flux in watts, m ˙ is the FLiBe mass flow in kg/s, C p is the specific heat capacity J/(kg·K), h is the overall heat transfer coefficient W/(m2·K), and A is the area of the plasma-facing surface. Δ T l m is the logarithmic mean temperature difference (LMTD) between the heat source and heat sink. It is defined in Equation (2) below and accounts for a non-uniform temperature difference along the path of the fluid. T p l a s m a is the temperature of the tungsten shell on the inner plasma-facing side.
Δ T l m = T p l a s m a T i n T p l a s m a T o u t ln T p l a s m a T i n T p l a s m a T o u t   K
Equation (1) can be rearranged to give us the minimum required heat transfer coefficient h m i n and minimum area-specific FLiBe mass flow m ˙ m i n / A , as shown in Equations (3) and (4), respectively.
h m i n = Q ˙ A   1 Δ T l m   W / m 2 K
m ˙ m i n A = Q ˙ A 1 C p T o u t T i n kg / m 2 s
Using Equations (2)–(4), and the boundary conditions in Table 2, a Δ T l m ~660 K, h m i n ~18 kW/(m2·K), and a m ˙ m i n / A of ~68 kg/(m2·s) were found. These calculated quantities represent the heat flow characteristics to remain within the boundary conditions set forth by the ARC reactor design. Figure 3 below shows the proposed ARC reactor operating point in terms of fluid temperature change, LMTD, and h m i n . This helps to maximize the thermal driving force, keep plasma facing surfaces cooler than the maximum temperatures allowed, and minimize asymmetric thermal stresses across the tungsten elements.
The minimum required overall heat transfer coefficient ( h m i n )   is a function of tungsten thickness δ W and the minimum average heat transfer coefficient for the fluid surface α m i n . Equation (5) below defines α m i n in terms of these new quantities and is plotted against tungsten layer thickness in Figure 4.
1 α m i n = 1 h m i n δ W k W     m 2 K W

2.2.2. Area Ratio

Geometric features added to a surface increase the surface area and direct fluid flow in an advantageous manner. These geometric features take the form of extended elements, “fins”, which extend away from the plate. The additional area is an advantage as the required α m i n decreases as the total area of the finned surface A f i n increases relative to the flat plate area A p l a s m a . The area ratio A r is defined as A f i n / A p l a s m a and can be used to adjust the required α m i n as a guide for design. Figure 5 below shows the general concept of increasing the area of a surface with a rectangular fin relative to a non-finned surface.
Figure 6 below shows the change in α m i n with an increasing area ratio A r . Theoretical thermal performance was analyzed to calculate the effects of different thicknesses of tungsten for the vacuum vessel. It was found that by reducing the thickness of the tungsten to 2 mm and increasing the area ratio significantly reduced α m i n . Figure 6 below shows the effect of reducing the thickness of the tungsten layer and increasing the area ratio using Equation (5).
The higher the area ratio, the better the chance a design has at successfully meeting the task of transferring the required amount of heat. However, there are limiting factors which constrain the geometry of the fluid surface to a range which is achievable. The assessment of the proposed design used a tungsten plate of 2 mm in thickness and an area ratio >2. This represented the minimum vacuum vessel wall thickness and a minimum in thermal resistance to test the thermal performance of the design concept. A thermodynamic design point of an α m i n ~12,000 [W/(m2·K)] was considered in the scope of this study as it is typical of forced convection for non-metallic, single-phase liquids [17].

2.2.3. Fin and Surface Efficiency

The use of fins on a surface constitutes the use of more material above the plate itself. The added material incurs thermal resistance to heat flow and thus factors in when determining optimal fin geometry. The use of very long fins is not required to reach the optimal heat transfer rate [18]. Fin efficiency η f i n is a term which describes the ability to transfer the heat of a flat surface or a tube with fins relative to the same surface without fins. It is the ratio of the average temperature difference between the end of geometry extending away from the base and the fluid to the average temperature difference between the base of the surface and the fluid [17]. Equation (6) below summarizes this definition.
η f i n = T ¯ f i n T ¯ f l u i d T ¯ b a s e T ¯ f l u i d ; 0 η f i n 1
The efficiency of the fins can be further expanded into an overall surface efficiency η o . The overall performance of a finned plate is a function of the fin efficiency η f i n , the number of fins N , the area of each fin A f , the total area of the finned plate A t , and η f i n . Equation (7) below is the overall efficiency of a finned plate flat plate.
η o = 1 N A f A t 1 η f i n
The surface efficiency of rectangular fins is plotted below in Figure 7 for an area of 1 m2 against A r , fin thickness t h , fin length L , and a channel width of 1 mm.
The primary purpose of optimizing fin geometry is to maximize thermal performance using the least amount of material, while the secondary purpose is to promote the mixing of cool fluid down to the surface of the finned structure. Rectangular fins were chosen as the best solution, as their shape allows for the optimal addition of surface area with vertical walls. The five other finned surfaces considered (rectangular pins, triangular/parabolic fins, and pins) exhibited poorer overall fin efficiencies with the same area ratios.
The fluid surface geometry is of critical importance to design a structure capable of handling the proposed high heat flux density of the ARC divertor leg. An increased area ratio was achieved by extending material away from the plate surface in the form of rectangular fins, as shown in Figure 5 and Figure 7. A logical starting point for analysis was using a narrow fin thickness (1 mm), narrow fluid channels (1 mm), and tall fins. Table 3 below shows the dimensions and their corresponding overall fin efficiencies, area ratios, and hydraulic diameters d h .

2.2.4. Pumping Power

About 300 m2 is needed to envelop an ARC vacuum vessel [4]. Using the calculated area specific mass flow rate m ˙ m i n / A , ~20,400 kg/s is required to provide the required amount of FLiBe to cool the vacuum vessel shell under the proposed conditions. The minimization of energy use by auxiliary systems to circulate a coolant is a priority. The pumping power required has the potential to demand a significant proportion of generated electricity and must be kept to a minimum. Equation (8) below shows the hydrodynamic relationship between pressure loss, fluid mass flow, and electrical power needed to move a fluid with pump efficiency.
P e l = m ˙ Δ p l o s s ρ η p u m p
Table 4 below displays the estimated pumping power required for the entire proposed PEXS for pressure losses between 1 and 16 bar a pump efficiency η p u m p of 0.85 and an overall coolant mass flow rate of 20,400 kg/s for a ~200 MWel ARC fusion reactor [19].

2.3. Methods

The boundary conditions relevant for this study were the heat densities to be transferred to the FLiBe heat transfer fluid. The design had to be built to minimize fluid pressure losses and thus minimize energy consumption for cooling. Furthermore, uniform temperature gradients within the PEXS material were crucial to maintaining symmetric thermal stresses and material longevity.
A literature review was performed to find any high-power heat transfer applications either in commercial use or tested in a lab. The results of the literature review served as a basis of inspiration for the initial designs tested in a numerical setting using COMSOL 6.0 [20]. Laminar flow conditions were analyzed as they represent flow conditions which exhibit the lowest pressure losses and worst (in terms of low) heat transfer conditions. This provided a safe-sided starting point in evaluating the feasibility of certain designs with regard to the boundary conditions.

3. Results and Discussion

The results are presented in the following manner: the state of the art, design concept, and simulated output. Each section provides the design development chronologically, and therefore, the result of each section is built upon the previous result and presented as such. A short discussion of each step is included with the results.

3.1. State of the Art

A literature review was conducted to discover commercial and/or research applications which could serve as a starting point for a design concept. Microchips and integrated circuits have evolved to reach 2.25 MW/m² by employing one of the four types of heat transfer: radiation and free convection, forced air-cooling, forced liquid-cooling, and liquid evaporation [21]. Another study examined the thermal exhaust requirements for VLSI (very-large-scale-integrated) circuits [22]. Using a micro-channel structure, laminar flow regime, and employing water as the heat transfer fluid, a heat flux of up to 7.9 MW/m2 was achieved with a maximum temperature change of 71 K between substrate and inlet fluid temperature [21]. In this paper, it is further suggested a heat flux density of 10 MW/m2 is possible by scaling a liquid-cooled heat exchanger design to microscopic dimensions. Table 5 below shows the experimental results for several dimensions, fluid pressures p , flow rates V ˙ , and heat flux densities q . Note specific mass flow rates for water as a coolant, m ˙ w , were calculated from available data and are of the same order of magnitude as the specific mass flow rates calculated using FLiBe.
A review of heat sinks related to nuclear fission reactors was investigated for a possible corollary for a solution to this heat transfer problem [23]. Nuclear reactors designed to generate electricity exhibit high rates of heat transfer from fuel rods to coolant fluid. Various coolants are used depending on the reactor design. An overview is given in Table 6 below to outline potentially suitable solutions.
Coolants are put under forced convection flow regimes to maximize heat transfer and safety during operation. In addition, some designs allow for a phase change of the coolant to increase heat transfer rates.
Mirrors used for the concentration of high-power laser light experience significant heating. Light beam quality experiences a sharp decline when the mirror begins to distort due to heating [24,25]. To mitigate these negative thermal effects, active cooling of the mirror with water is applied. Table 7 below outlines the assumptions made for a simulation regarding the heat transfer area facing the fluid, and Figure 8 shows the proposed design. It was found that a microchannel heat sink was able to decrease the thermal deformation of a mirror. Similar to the VLSI circuit heat sink, microchannels were found to have a high depth to width ratio.
Another high-power heat sink application was designed to limit the temperature rise used in HVDC (high-voltage direct-current) thyristors in electrical grids. The heat flux associated with the use of this equipment could experience 2 MW/m2 under the condition of a very low acceptable temperature increase. The design of the heat transfer surface in this application is also of interest to us. Figure 9 below shows a double-layer cross-flow matrix of channels that is used to move heat away from the contact surface and limit the temperature rise in a thyristor [26].
This design proposes a second example of using small channels near the surface for the purpose of efficiently moving large amounts of heat away from the surface needing to be cooled.

3.2. Design Approach

The results of the literature review gave a deeper insight into the general design characteristic of high-power heat exchangers tasked with cooling a surface or device. Three of the four designs studied employ the use of small channels to guide a coolant near a surface to remove heat. In combination with the required geometry of the proposed ARC reactor design, the following concept was developed followed by an evaluation with numerical simulation.

3.2.1. Proposed ARC PEXS Design

The proposed design provided from the literature was used as a starting point for the ARC PEXS [4]. Theoretical principles and designs from other high-power heat exchange applications were used to develop a concept for a PEX capable of transferring heat away from the vacuum vessel per the stated boundary conditions. Figure 10 below shows a simplified schematic (a) and illustration (b) of how heat is removed from the fusion reactor.
FLiBe is circulated over the surface of the fusion reactor (vacuum vessel) to remove heat from the core and transfer heat to a secondary loop (a). Coolant is pumped from the point labeled T i n and distributed evenly over the surface of the vacuum vessel (dark blue/black dashed line). It removes heat Q ˙ from the reactor core and maintains operating temperatures below 1500 K on the plasma-facing side of the tungsten shell. The heated coolant is discharged locally (red arrows) into the FLiBe tank, which encloses the vacuum vessel. The heated FLiBe is taken from the tank to generate steam for electricity generation (b). The ARC reactor components relevant to this study are the divertor leg (2), FLiBe blanket (5), and vacuum vessel (8), shown in Figure 11 below.
A network of channels supplying coolant to the surface of the vacuum vessel is essential to maintain its structural integrity. Multiple inlets would be necessary to maintain a low temperature change in the coolant while reducing the pressure drop between pump and outlet to the FLiBe blanket. Figure 12 below shows a concept for cooling the divertor leg via a primary supply channel (circular blue pipe) branching into secondary capillary channels (blue longitudinal conduits) before entering a series of finned structures (gray lattice layer) attached to the vacuum vessel.

3.2.2. Fluid surface

The construction assumed each fin is attached to an upper cross-fin, as shown in Figure 1, at a 90° angle, contrary to the geometry in Figure 1. This was chosen to simplify manufacturing and assembly. Consequently, it presented an undesired thermal bridge to the upper channels. Thermal bridging was limited by attaching every other lower-channel fin with an upper-channel fin, resulting in a “∩-type” of lower-channel geometry. Figure 13 below shows a cross-section (a) and perspective (b) of the ∩-type lower channel.

3.3. Numerical Simulation

The construction shown in Figure 14 was chosen to take advantage of symmetry planes along the xz-planar vertical faces and the right-hand yz-planar vertical face. Figure 14 shows the fluid flow regime as simulated for all cases outlined in Table 3.
Inlets are defined by blue arrows and the outlets by red. Heat from the fusion reaction is applied on the bottom side and enters the fluid stream along the various surfaces of the channels. The illustration above is a rendering of the design concept with three channels in the simulated environment. Geometries with 1, 2, 3, 4, 5, 10, and 20 inlets for each fin length were analyzed for a total of 18 unique assessments. Figure 15 below shows a 20-inlet geometry.
Each combination of parameters was using a mass flow rate consistent with the area-specific mass flow rate calculated in Section 2.2.1.
The thermal performance was first assessed by analyzing the average outlet temperature of each configuration. Figure 16 below shows calculated average outlet temperatures for each fin length (1.5, 2.5, and 3.5 mm) shown in the legend.
Each successive addition of inlets showed a reduction in the average fluid temperature at the outlet. This was expected as the addition of inlets introduced more FLiBe to a channel with fixed cross-sectional areas. The maximum fluid velocity u m a x increased with the number of inlets while decreasing with fin length f L for geometries with the same number of inlets. The maximum domain temperature, T m a x , and average heat transfer coefficient on the fluid surface, α ¯ f , also decreased with the addition of inlets up to three to four inlets before increasing again. The maximum pressure-drop from inlet to outlet across each geometry, Δ p m a x , increased with increasing u m a x . Table 8 below shows the trends observed organized by fin length, area ratio, and average minimum heat transfer coefficient, α m i n , for each geometry by the number of inlets N i n .
The simulated results indicate a thermal performance satisfying the boundary conditions for geometries with one to five inlets based on T m a x . The minimum average heat transfer coefficient for a 2 mm thick tungsten wall and flat plate ( α m i n ) was not consistent with the maximum temperature results. Simulated geometries included thermal bridges and a cool side boundary condition not considered in the theoretical analysis. Errors associated with the model, the use of average domain values for calculations (e.g., average outlet temperature), and deviations from theory explained simulations where T m a x exceeded 1500 K while α ¯ f exceeded α m i n .

4. Conclusions

The results of this study have determined there is a viable path forward to develop a PEXS capable of fulfilling the thermodynamic constraints of the ARC reactor design. Theories regarding heat transfer, fin and surface efficiency, and fluid dynamics were used as design guides for the development of a concept analyzed using numerical simulation. The analysis of the concept performance indicated positive progress in developing a heat exchange surface capable of transferring a specific heat flux of 12 MW/m². Simulated heat transfer coefficients for the fluid surface approached values found using an energy balance for heat exchangers. The most feasible solution was a channeled block of tungsten with fins of 3.5 mm in height, fin and channel widths of 1 mm, and five sequential inlets. It represented the solution with the lowest pressure losses and thus only consumed 2.3% of the produced electrical power generated. However, the simulated results were not fully in agreement with the theory. Further investigation is necessary into the accuracy of the model by means of empirical results and tests in a laboratory setting.
Before moving forward with the proposed geometries, the calculated and simulated numerical results in this work need verification. The construction of an experimental test rig for measuring empirical performance is required to understand the accuracy of the calculated and simulated results. A design for an experiment must be developed using the materials (or their analogs) which accurately model the problem described in this work.

Author Contributions

Conceptualization, R.B. and L.F.; methodology, R.B. and L.F.; software, R.B.; formal analysis, R.B. and L.F.; writing—original draft preparation, R.B.; writing—review and editing, R.B. and L.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data may be accessed by a request placed with either of the authors of this manuscript.

Acknowledgments

The authors would like to thank all those involved in developing and supporting this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. EUROfusion. 2022. Available online: https://www.euro-fusion.org/fusion/history-of-fusion/#c2407 (accessed on 23 May 2022).
  2. Eddington, A. The Internal Constitution of the Stars. Nature 1920, 106, 14–20. [Google Scholar] [CrossRef]
  3. Oliphant, M.L.E.; Harteck, P.; Rutherford, O.M. Transmutation Effects Observed with Heavy Hydrogen. Proc. R. Soc. 1934, 144, 692–703. [Google Scholar] [CrossRef]
  4. Kuang, A.Q.; Cao, N.M.; Creely, A.J.; Dennett, C.A.; Hecla, J.; LaBomard, B.; Tinguely, R.A.; Tolman, E.A.; Hoffman, H.; Major, M.; et al. Conceptual design study for heat exhaust management in the ARC fusion pilot plant. Fusion Energy Des. 2018, 137, 221–242. [Google Scholar] [CrossRef]
  5. Tungsten, Midwest Tungsten Serv. 2022. Available online: https://www.tungsten.com/material-info/tungsten-w (accessed on 23 May 2022).
  6. Aldrich, R. Laser Fundamentals. Naval Surface Warfare Center, Dahlgren Division. Available online: https://man.fas.org/dod-101/navy/docs/laser/fundamentals.htm#LASERTYPES (accessed on 26 January 2022).
  7. Hu, L.; Wirth, B.D.; Maroudas, D. Thermal conductivity of tungsten: Effects of plasma-related structural defects from molecular-dynamics simulations. Appl. Phys. Lett. 2017, 111, 081902. [Google Scholar] [CrossRef]
  8. Suslova, A.; El-Atwani, O.; Sagapuram, D.; Harilal, S.S.; Hassanein, A. Recrystallization and grain growth induced by ELMs-like transient heat loads in deformed tungsten samples. Sci. Rep. 2014, 4, 1–11. [Google Scholar] [CrossRef] [PubMed]
  9. Haubenreich, P.N.; Engel, J.R. Experience with the Molten-Salt Reactor Experiment. Nucl. Appl. Technol. 1970, 8, 118–136. [Google Scholar] [CrossRef]
  10. Moir, R.W. HYLIFE-II Inertial Confinement Fusion Reactor Design. Fusion Technol. 1991, 19, 617–624. [Google Scholar] [CrossRef]
  11. Cadwallader, L.C.; Longhurst, G.R. Flibe Use in Fusion Reactors: An Initial Safety Assessment. Ida. Natl. Eng. Environ. Lab. 1999, 1, 1–25. [Google Scholar]
  12. Janz, G.J.; Gardner, G.L.; Krebs, U.; Tomkins, R.P.T. Molten Salts: Volume 4, Part 1, Fluorides and Mixtures Electrical Conductance, Density, Viscosity, and Surface Tension Data. J. Phys. Chem. Ref. Data 1974, 3, 31–36. [Google Scholar] [CrossRef]
  13. Romatoski, R.R.; Hu, L.W. Fluoride salt coolant properties for nuclear reactor applications: A review. Ann. Nucl. Energy 2017, 109, 635–647. [Google Scholar] [CrossRef]
  14. Inconel 718 Technical Data, High Temp. Metals. 2015. Available online: https://www.hightempmetals.com/techdata/hitempInconel718data.php (accessed on 27 January 2022).
  15. Cortés, R.; Barragán, E.R.; López, V.H.; Ambriz, R.R.; Jaramillo, D. Mechanical properties of Inconel 718 welds performed by gas tungsten arc welding. Int. J. Adv. Manuf. Technol. 2018, 94, 3949–3961. [Google Scholar] [CrossRef]
  16. Cole, N.C.; Gilliland, R.G.; Slaughter, G.M. Weldability of tungsten and its alloys. Weld. J. 1971, 50, 419–426. [Google Scholar]
  17. Kind, M.; Martin, H. Verein Deutscher Ingenieure (VDI)—Heat Atlas; Springer: Berlin/Heidelberg, Germany; Dordrecht, The Netherlands; London, UK; New York, NY, USA, 2010; ISBN 9783540778769. [Google Scholar]
  18. Bergman, T.L.; Lavine, A.S.; Incropera, F.P.; DeWitt, D.P. Fundamentals of Heat and Mass Transfer, 7th ed.; John Wiley & Sons: New York, NY, USA, 2011; ISBN 978-0470-50197-9. [Google Scholar]
  19. Sorbom, B.N.; Ball, J.; Palmer, T.R.; Mangiarotti, F.J.; Sierchio, J.M.; Bonoli, P.; Kasten, C.; Sutherland, D.A.; Barnard, H.S.; Haakonsen, C.B.; et al. ARC: A compact, high-field, fusion nuclear science facility and demonstration power plant with demountable magnets. Fusion Eng. Des. 2015, 100, 378–405. [Google Scholar] [CrossRef]
  20. COMSOL, COMSOL Multiphysics®. 2022. Available online: https://doc.comsol.com/6.0/docserver/#!/com.comsol.help.comsol/helpdesk/helpdesk.html (accessed on 20 July 2022).
  21. Scott, W.A. Cooling of Electronic Equipment; John Wiley & Sons: New York, NY, USA, 1974. [Google Scholar]
  22. Tuckerman, D.B.; Pease, R.F.W. High-Performance Heat Sinking for VLSI. IEEE Electron Device Lett. 1981, 2, 126–129. [Google Scholar] [CrossRef]
  23. Duderstadt, J.J.; Hamilton, L.J. Nuclear Reactor Analysis; John Wiley & Sons: New York, NY, USA, 1976; p. 675. [Google Scholar]
  24. Cao, H.; Chen, G. Optimization design of microchannel heat sink geometry for high power laser mirror. Appl. Therm. Eng. 2010, 30, 1644–1651. [Google Scholar] [CrossRef]
  25. Peng, Y.F.; Cheng, Z.H.; Zhang, Y.N.; Qiu, J.L. Laser-induced temperature distributions and thermal deformations in sapphire, silicon, and calcium fluoride substrates at 1.315 μm. Opt. Eng. 2001, 40, 2822. [Google Scholar] [CrossRef]
  26. Li, Q.; Fischer, L.; Qiao, G.; Mura, E.; Li, C.; Ding, Y. High performance cooling of a HVDC converter using a fatty acid ester-based phase change dispersion in a heat sink with double-layer oblique-crossed ribs. Int. J. Energy Res. 2020, 44, 5819–5840. [Google Scholar] [CrossRef] [Green Version]
  27. Chandler, D. A New Path to Solving a Longstanding Fusion Challenge, MIT News. 2018. Available online: https://news.mit.edu/2018/solving-excess-heat-fusion-power-plants-1009 (accessed on 13 February 2022).
Figure 1. General design concept with distribution channels (the left one semitransparent for better insight into the lower grid).
Figure 1. General design concept with distribution channels (the left one semitransparent for better insight into the lower grid).
Energies 16 01890 g001
Figure 2. Generalized heat exchanger temperature profile.
Figure 2. Generalized heat exchanger temperature profile.
Energies 16 01890 g002
Figure 3. Dependence of LMTD and hmin with ΔTfluid and m ˙ min/A.
Figure 3. Dependence of LMTD and hmin with ΔTfluid and m ˙ min/A.
Energies 16 01890 g003
Figure 4. Dependence of α m i n on the thickness of the inner layer of tungsten and Equation (5).
Figure 4. Dependence of α m i n on the thickness of the inner layer of tungsten and Equation (5).
Energies 16 01890 g004
Figure 5. Increased surface area relative to a flat surface by use of fins.
Figure 5. Increased surface area relative to a flat surface by use of fins.
Energies 16 01890 g005
Figure 6. α m i n as a function of area ratio Ar and tungsten shell thickness.
Figure 6. α m i n as a function of area ratio Ar and tungsten shell thickness.
Energies 16 01890 g006
Figure 7. Surface efficiency of rectangular fins on a flat plate.
Figure 7. Surface efficiency of rectangular fins on a flat plate.
Energies 16 01890 g007
Figure 8. Proposed heat sink for the cooling of a laser mirror.
Figure 8. Proposed heat sink for the cooling of a laser mirror.
Energies 16 01890 g008
Figure 9. HVDC thyristor heat sink [26].
Figure 9. HVDC thyristor heat sink [26].
Energies 16 01890 g009
Figure 10. (a) Simple schematic of PEX in operation; (b) detailed illustration of PEX cooling the vacuum vessel.
Figure 10. (a) Simple schematic of PEX in operation; (b) detailed illustration of PEX cooling the vacuum vessel.
Energies 16 01890 g010
Figure 11. Conceptual design of the ARC reactor [27].
Figure 11. Conceptual design of the ARC reactor [27].
Energies 16 01890 g011
Figure 12. Conceptual design for the PEXS on the divertor leg.
Figure 12. Conceptual design for the PEXS on the divertor leg.
Energies 16 01890 g012
Figure 13. (a) Cross-section of ∩-type lower channel; (b) perspective of ∩-type lower channel.
Figure 13. (a) Cross-section of ∩-type lower channel; (b) perspective of ∩-type lower channel.
Energies 16 01890 g013
Figure 14. Flow regime as set-up in COMSOL Multiphysics.
Figure 14. Flow regime as set-up in COMSOL Multiphysics.
Energies 16 01890 g014
Figure 15. Example of a 20-inlet geometry.
Figure 15. Example of a 20-inlet geometry.
Energies 16 01890 g015
Figure 16. Average outlet temperatures for each simulated geometry.
Figure 16. Average outlet temperatures for each simulated geometry.
Energies 16 01890 g016
Table 1. Material properties of FLiBe.
Table 1. Material properties of FLiBe.
ParameterCorrelationUnitsUncertainty
density [12] ρ = 2413–0.488 · T[K]kg/m32%
heat capacity [13] C p = 2386J/(kg·K)3%
dynamic viscosity [12] μ = 0.116 · e(3755/T[K])mPa·s20%
thermal conductivity [13] k = 1.1W/(m·K)10%
Table 2. Boundary conditions of vacuum vessel cooler [4].
Table 2. Boundary conditions of vacuum vessel cooler [4].
Boundary ConditionValueQuantity
max heat flux density~12 MW/m2 Q ˙ / A
Tmax tungsten1500 K T p l a s m a
ΔTmax FLiBe~75 K T o u t T i n
Table 3. Analyzed fluid surface geometries.
Table 3. Analyzed fluid surface geometries.
Fin Length [mm] η o [-] A r [-] d h [mm]
1.5~0.80~2.5~1.2
2.5~0.65~3.5~1.4
3.5~0.53~4.5~1.6
Table 4. Estimated electrical power consumption for pumping FLiBe with increasing PEXS pressure losses.
Table 4. Estimated electrical power consumption for pumping FLiBe with increasing PEXS pressure losses.
Δ p l o s s [Bar] P e l [MW] % of Output
11.80.9
23.61.8
47.13.6
814.27.1
1628.414.2
Table 5. Experimental VLSI heat sink with water filled microchannels [22].
Table 5. Experimental VLSI heat sink with water filled microchannels [22].
Exp. w c
[μm] 1
w w
[μm] 2
z
[μm] 3
p
B a r
V ˙     c m 3 s m ˙ w     k g m 2 s R t h , m a x   K W   4 q   M W m 2
156443201.0344.7470.1101.81
255452871.1726.5650.1132.77
350503022.1378.6860.0907.90
1 channel width; 2 wall width; 3 channel depth; 4 maximum thermal resistance.
Table 6. Experimental VLSI heat sink with water filled microchannels [23].
Table 6. Experimental VLSI heat sink with water filled microchannels [23].
Reactor Type 1PWR aBWR/6 bHTGR cLMFBR dGCFR eCANDU
PHW f
Manufacturer 2W aB&W bGE c-----
CoolantH2OH2OH2OH2OHeNaHeD2O
Avg. heat flux 30.6850.640.650.5030.2040.1050.0930.050
Max. heat flux 31.831.681.731.1150.5832.371.681.15
1a pressurized water reactor; 1b boiling water reactor; 1c high-temperature gas-cooled reactor; 1d liquid metal fast breeder reactor; 1e gas-cooled fast breeder reactor; 1f Canadian pressurized heavy water reactor. 2a Westinghouse; 2b Babcock and Wilcox; 2c General Electric. 3 (MW/m2).
Table 7. Assumed conditions for simulation of a laser disk mirror [24].
Table 7. Assumed conditions for simulation of a laser disk mirror [24].
Assumed ConditionsValueUnit
heat flux density2MW/m²
cooling region radius23.5mm
mirror thickness2mm
channel width1mm
channel depth2mm
fin width1.5mm
water flow rate500mL/min
specific mass flow rate4.8kg/(m²·s)
Table 8. Estimated electrical power consumption with increasing PEXS pressure losses.
Table 8. Estimated electrical power consumption with increasing PEXS pressure losses.
f L   [ mm ] / A r [-];
amin 3 W / m K
N i n [-] u m a x   1
m / s
  T m a x   1
[K]
α ¯ f   1
W / m K
Δ p m a x   1
P a
P e l   2,*
[MW]
P e l / P e l , A R C   4,*
[—]
1.5/2.5;
9160
10.134146510,7033914.47.2%
20.260143311,22213314.47.2%
30.393142911,26626913.36.7%
40.527143410,84244312.26.1%
50.666143410,85167511.85.9%
101.3501556854227749.64.8%
202.6641917847284948.04.0%
2.5/3.5;
6540
10.102146793242711.15.6%
20.199143510,1468210.65.3%
30.300143599271628.04.0%
40.402142510,1182857.83.9%
50.518144297404237.33.7%
101.0331621734015695.42.7%
202.0262131549548064.22.1%
3.5/4.5;
5090
10.08114869510229.34.7%
20.16014499364596.93.5%
30.246144891931145.62.8%
40.328146489221855.02.5%
50.413149085242734.62.3%
100.8551692659610323.61.8%
201.6752368454931702.81.4%
1 Simulated output. 2 Calculated with simulated output. 3 Calculated with theory Section 2.2.1. 4 Assuming power output of the ARC reactor is 200 MWel [19]. * Hydraulic power needed to supply FLiBe to the tungsten vacuum vessel only.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Beaufait, R.; Fischer, L. Blanket Cooling of a Fusion Reactor. Energies 2023, 16, 1890. https://doi.org/10.3390/en16041890

AMA Style

Beaufait R, Fischer L. Blanket Cooling of a Fusion Reactor. Energies. 2023; 16(4):1890. https://doi.org/10.3390/en16041890

Chicago/Turabian Style

Beaufait, Robert, and Ludger Fischer. 2023. "Blanket Cooling of a Fusion Reactor" Energies 16, no. 4: 1890. https://doi.org/10.3390/en16041890

APA Style

Beaufait, R., & Fischer, L. (2023). Blanket Cooling of a Fusion Reactor. Energies, 16(4), 1890. https://doi.org/10.3390/en16041890

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop