Plasma and Fusion Research

Volume 17, 2405073 (2022)

Regular Articles

Computational Investigation of 2-D Temperature Distribution in Static Liquid Metals Exposed to Steady State Plasmas
Department of Physics, Faculty of Science, Kasetsart University, Chatuchak, Bangkok 10900, Thailand
(Received 10 January 2022 / Accepted 24 April 2022 / Published 6 June 2022)


Understanding the temperature distribution in a liquid metal under plasma bombardment is required for characterizing and controlling their own impurity releases and heat transfer. To achieve this, a 2-D heat conduction inside a static liquid bombarded by a plasma is numerically solved in this study. Thin layers consist of Al, Li, In, Sn and Ga. Plasma constituents are D and Ar. The liquid temperature is initialized by melting temperature. The upper surface is heated by ions and electrons but cooled by evaporation and thermal radiation. The lower surface is in contact with the boundary conditions: the fixed melting temperature for characterizing the temperature spread influenced by different thermal diffusivities and ion masses; and the floating temperature governed by convective cooling and thermal radiation for characterizing the heat transfer across the liquid provided by the implementation of a coolant. It appears that Al conducts heat well so the temperature distribution is smoothed out, probably a choice for excessive heat flushing during abnormal events. Sn and In may be good undesirable impurity collectors because of low evaporation with less coolant concern, but not for Li and Ga.


plasma facing component, plasma surface interaction, liquid metal, Maxwell-Boltzmann distribution, floating potential, heat conduction, thermal radiation, convective cooling, Hertz-Knudsen-Langmuir (HKL) equation, Forward Time Central Space (FTCS) difference

DOI: 10.1585/pfr.17.2405073


  • [1] R.E. Nygren and F.L. Taberés, Nucl. Mater. Energy 9, 6 (2016).
  • [2] S.V. Mirnov and V.A. Evtikhin, Fusion Eng. Des. 81, 113 (2006).
  • [3] G.Z. Zuo et al., Fusion Eng. Des. 89, 2845 (2014).
  • [4] R.B. Gomes et al., Fusion Eng. Des. 83, 102 (2008).
  • [5] F.A. Hernández and P. Pereslavtsev, Fusion Eng. Des. 137, 243 (2018).
  • [6] E. Pickering et al., Entropy 23, 98 (2021).
  • [7] P.C. Stangeby, The Plasma Boundary of Magnetic Fusion Devices (IOP Publishing, Bristol, 2000).
  • [8] N. Somboonkittichai, J. Phys.: Conf. Ser. 2145, 012022 (2021).
  • [9] J.P. Holman, Heat Transfer (McGraw-Hill, Inc., New York, 1981).
  • [10] A. Chamber, Modern Vacuum Physics (CRC Press, Florida, 2005).
  • [11] Y. Jaluria and K.E. Torrance, Computational Heat Transfer (Taylor & Francis, New York, 2003).
  • [12] W.H. Press et al., Numerical Recipes: The Art of Scientific Computing (Cambridge University Press, New York, 2007).
  • [13] W.M. Hayes et al., CRC Handbook of Chemistry and Physics (CRC Press, Florida, 2014).
  • [14] N. Somboonkittichai, Manuscript in preparation (2022).
  • [15] Python Software Foundation, Python 3.8.10 documentation (Available at
  • [16] C.R. Harris et al., Nature 585, 357 (2020).
  • [17] P. Virtanen et al., Nature Methods 17(3), 261 (2020).
  • [18] J.D. Hunter, Comput. Sci. Eng. 9(3), 90 (2007).