LithiumMolybdenum Infused Trenches Li MIT David Ruzic Wenyu

Lithium/Molybdenum Infused Trenches (Li. MIT) David Ruzic, Wenyu Xu, Vijay Surla PFC meeting UCLA, August 2010 Nuclear, Plasma, and Radiological Engineering Center for Plasma-Material Interactions Contact: druzic@illinois. edu PFC/FNST UCLA August 3 -6, 2010

Outline Thermoelectric Infused 2 MHD theory trenches divertor plate design Calculations § Flow needed to remove heat § Flow predicted from TEMHD § 3 -D computer simulation Ejection Future issue and force balance work at Illinois Conclusion PFC/FNST UCLA August 3 -6, 2010

3 Thermoelectric effect Lithium j Molybdenum j Seebeck Coefficient of Li Black: our measurements Lithium Like a thermocouple, a voltage is created at a junction of two metals dependent on the temperature. [1] A current will flow based on that voltage difference: where σ is the conductivity, and ΔS is the difference in Seebeck coefficients. [2] There is a large difference in S between Li and most other metals and it increases with T. [3] PFC/FNST UCLA August 3 -6, 2010

4 Does it work? Yes! SLi. DE Results Sheet electron beam (like a divertor heat flux stripe) hits a well diagnosed tray of molten Li in the presence of a magnetic field Theoretical vs experimental velocities : Movie is of a case where the “Jaworski” number is near 1 and TEMHD and TCMHD (Marangoni effect) are balanced, so flow osscilates between swirling and splitting. [4] M. A. Jaworski, et al. Phys. Rev. Lett. 104, 094503 (2010) PFC/FNST UCLA August 3 -6, 2010

5 Li. MIT Design Inner divertor shelf Divertor Strikepoint Stripe Left is a cross-section of NSTX showing the “shelf-like” inner divertor plates. Right is the Li. MIT concept: molybdenum tiles with radial trenches containing lithium. The trenches run in the radial (polodial) direction such that they lie primarily perpendicular to the torroidal magnetic field. PFC/FNST UCLA August 3 -6, 2010

Lithium Flow in the Trenches is Self-Pumping The top surface of the Li is hotter than the surface that is deeper. Therefore there is a VERTICAL temperature gradient Outlet Heat flux Hot F 6 J B Li flow Cooling channels Inlet Passive Li replenishment Concept for heat removal using TEMHD. The Li flows in the slots of the Mo plate powered by the vertical temperature gradient. This vertical temperature gradient generates vertical current, which when “crossed” by the torroidal magnetic field, will create a radial force on the Li driving it along the slot. This flow will transfer the heat from the strike point to other portions of the divertor plate. The bulk of the Mo plate could be actively cooled for a long-pulsed device or passively cooled for something like NSTX. Under the plate the Li flows back naturally. PFC/FNST UCLA August 3 -6, 2010

Thermoelectric Driven Flow Calculation 7 Heat flux: divertor strike point width, 1 cm Z direction t w h Grad T TE current Flow B L Look just at one channel. All channels are equivalent The width of Li (w) is 1 mm and the width of the Mo (t) trench is 1 mm. The depth of Li (h) is 5 mm. The length of the Li (L) is 100 mm. The magnetic field is 1 T. The heat flux on top surface is 20 MW/m 2 and the footprint of the heat flux is 2 mm*10 mm. This represents the divertor heat flux maximum value for the calculations here, but the concept should work for even higher values. . PFC/FNST UCLA August 3 -6, 2010

Heat load on Li. MIT: What Velocity is Needed? 8 The assumed heat flux means the 1 mm Li slot needs to handle 400 W heat load. (400 W = 20 MW/m 2 X 2 mm X 10 mm) To put this in the proper perspective, the temperature rise over the 5 second shot of an uncooled one-half inch thick (1. 2 cm) Mo plate the same size as our imagined Li/Mo trench/finger can be calculated using as follows: E = 400 W*5 s = 20, 000 J. V = (0. 2 cm)(10 cm) = 2. 4 cm 3 Cv=2. 57 J cm-3 K-1 for Mo. Therefore the temperature rise for a passive plate is 324 K. To see if a flowing surface can remove that much power with a lower temperature rise, we can use the following equation [5] : Here the first term is the heat conduction term and the second is the heat convection which transfers the heat along the trench. ΔT+ is the temperature difference across the depth of Li and ΔTll is the temperature difference between the inlet and outlet of the Li trench. When the heat flux is assumed to be 20 MW/m 2 on 1 cm, and both temperature differences are assumed to be 200 K, the velocity needed is only u=9. 9 cm/s. With this flow the temperature rise would only be 200 K. PFC/FNST UCLA August 3 -6, 2010

What Velocity is Generated by TEMHD ? The Seebeck coefficient for Mo is about 13μV/K (at 400 C). The value for Li reported in the literature [3] is 23μV/K (at 400 C). However, we have measured this number ourselves and obtained much higher values closer in line with those of Bidwell. This higher number 43μV/K (at 400 C) will be used first, for a difference of 30μV/K, but then results with one-third of that value will also be included. The current and velocity can be calculated with these two equations [1] Here P is the difference of Seebeck coefficient (thermoelectric power) between two material. Ha is the Hartmann number. C is a parameter coming from the conservation of current. The velocity with a higher thermoelectric power is 118 cm/s and the velocity is 39 cm/s with the smaller thermoelectric power, still enough to remove the heat. Both of these numbers are much greater than the 9. 9 cm/s on the last slide. 9 PFC/FNST UCLA August 3 -6, 2010

Further Analysis: 3 D – FLUENT Calculation Mo plate 10 Li channel 20 MW/m 2 heat flux Back flow of Li What will really happen to the temperature of the Li when the strike point heat flux hits Li. MIT? Is the 200 K gradient assumption reasonable? A simple 3 D heat transfer model is run with FLUENT has convection flow and 3 D heat transfer included. Boundary Conditions: The Li flows through the slot between two Mo plate and flow back below the Li slot. The inlet velocity is set to be 100 cm/s. The initial temperature is 500 K. The top heat flux is 20 MW/m 2 with a 2 mm*10 mm rectangle shape. The bottom temperature is set to be 500 K as a constant. Will we remain under 700 K as the one-D calculation indicated? PFC/FNST UCLA August 3 -6, 2010

11 Answer: Almost, but that is OK Left figure is the temperature distribution of the top surface of Li slot and right figure is the temperature change along the center line of the top surface. The temperature difference between the inlet and the outlet is less than 100 K – great! However the temperature immediately at the strike point reaches a maximum of 750 K (480 C) over 1 mm. There will be evaporation there – which gives additional cooling not included in the model. Some evaporation (radiative cooling) is acceptable and even desired. Also, a higher gradient in T gives a larger thermoelectric current. The temperature of most of the Li surface is under 650 K (380 C) which proves the ability of the Li. MIT concept to handle a high heat flux through TEMHD. PFC/FNST UCLA August 3 -6, 2010

12 Does the Radial Temp. Gradient Cause Ejection? Heat flux Li Grad T Lorentz Force B TE current Mo One concern about using free surface Li is the ejection problem. The temperature gradient along the Li flowing direction will generate a thermoelectric current along the same direction and the Lorentz force may could eject the Li into the plasma. Similar to the capillary porous system (CPS) [6] which effectively has very narrow channels, Li. MIT’s trench design has very narrow slots (1 mm) to utilize the capillary force to hold the Li in place. PFC/FNST UCLA August 3 -6, 2010

Does it Work? Capillary Force Balance 13 The thermoelectric current parallel to the Li flowing direction is [1]. Here the temperature difference is also assumed to be 200 K and the temperature gradient d. T/dy=2000 K/m. Under these conditions j. TEMHD=1. 78*105 A/m 2. So total current along the Li trench is 0. 89 A and the force from the TEMHD is 0. 089 N. Here P is chosen to be 30μV/K (the high value) The capillary force is 2ΣL and Σ=0. 3 N/m [7] at 600 K. So the capillary force is about 0. 06 N. Close, but slightly less than the ejection force. However this is the worse case. Mitigating factors: § The capillary force is not effected by thermoelectic power, P. It P is actually smaller, the forces will balance. Also, if the trenches were 0. 66 mm wide instead of 1 mm wide, then the two forces would balance. § The temperature difference in the radial direction could be less than 200 K. § The leading edges of the Mo fingers could be Mo-sprayed, greatly increasing the surface area and therefore the capillary force. This would easily hold the Li in. § The part of the Li with the highest ejection force is also the hottest Li and has likely already evaporated anyway! PFC/FNST UCLA August 3 -6, 2010

14 Future Work at Illinois on SLi. DE Heat Stripe for SLi. DE’s Electron Beam Magnet coil E-beam line source B Li Tray The SLi. DE experiment at Illinois is being reconfigured to test this concept. We expect to be able to show radial flows of Li along radial trenches in a Mo plate and measure the flow velocity compared to calculations. An electron beam is used to provide the heat flux while the magnet can generate about 800 Gauss magnetic field parallel to the tray surface. The temperature rise of the Li will also be monitored and compared to theory. Return channels (lithium hydraulic engineering) will also be tested to find a design compatible with tokamak operation. PFC/FNST UCLA August 3 -6, 2010

15 Conclusion A Mo trench structure with flowing Li is proposed as a potential method to absorb the high heat flux at the divertor region. Thermoelectric effect is utilized to drive the Li flowing along the radius trench direction. The heat transfer ability is estimated based on the 1 -D heat transfer model and simulated in the 3 -D domain. Both give positive results showing the ability of flowing Li to mitigate the peak heat flux and to transfer the heat without an unacceptable temperature increase. The ejection problem is analyzed and should be able to be suppressed by the capillary force. A simple trench structure system is under construction at UIUC to validate this design. This type of divertor structure should “stay on the radar” for future fusion devices. There is significant interest for using it on HT-7 in China. PFC/FNST UCLA August 3 -6, 2010

16 Reference [1] M. A. Jaworski, Ph. D. thesis, University of Illinois, Urbana, IL (2009) [2] J. A. Shercliff, Thermoelectric magnetohydrodynamics, J. Fluid Mech. 91, 231 (1979) [3] P. Ioannides, et al. Journal of Physics E 8, 315 (1975) [4] M. A. Jaworski, et al. Phys. Rev. Lett. 104, 094503 (2010) [5] M. A. Jaworski, et al. Journal of Nuclear Materials 390– 391, 1055– 1058(2009) [6] V. A. Evtikhin, et al. , Fusion Engineering and Design 49– 50 (2000) 195– 199. [7] M. A. Abdou, et al. On the exploration of innovative concepts for fusion chamber technology: APEX interim report. 282 Technical Report UCLA-ENG-99206, University of California, Los Angeles, November 1999. PFC/FNST UCLA August 3 -6, 2010