용융염 전해정련공정은 사용후핵연료로부터 전기화학적인 방법을 통해 음극에서 우라늄을 회수 하는 공정이다. 이 때 우라 늄은 약 30wt%의 염을 포함하고 있어 순수한 우라늄을 얻기 위해서는 염을 제거하는 Cathode Process (CP)가 필수적이다. CP는 대량의 우라늄을 처리해야 하므로 파이로공정의 난관중의 하나로 인식되고 있으며, 우라늄의 순도가 최종적으로 결정 되는 단계이므로 매우 중요한 공정이다. 현재, 이에 대한 연구는 주로 실험적 방법에 근거 하고 있어 염 제거 공정 중 온도, 압력, 염 가스의 거동을 관찰하기 어렵다. 따라서 본 연구에서는, 공정의 운전 조건에 대해 적합한 수학적 모델을 이용하여 전산모사 해석을 진행하였다. 본 연구는 증류부에서 염 가스의 증류 량, 확산계수에 의해 계산된 장치 내 염 가스의 이동 그 리고 응축부에서의 응결속도를 중점적으로 연구하였다. 장치내의 각각의 염 가스 거동을 정의하기 위해 Hertz-Langmuir 관 계식, Chapman-Enskog Theory, ANSYS-CFX의 상용 코드를 사용하였다. 그리고 HSC Chemistry에서 염의 물성 값을 이용 하여 모델을 구성하였다. 본 연구의 전산모사 해석을 통해 얻은 연구 결과를 이용하여 염 가스의 거동과 장치의 최적 운전 조건을 예측하였다. 따라서 본 해석 결과는 CP의 물리적 현상을 깊게 이해하는데 쓰일 뿐 아니라, 공학규모의 CP 장치를 상 용규모로 확장하는데 이용 할 수 있다.

Molten salt electrorefining process achieves uranium deposits at cathode using an electrochemical processing of spent nuclear fuel. In order to recover pure uranium from cathode deposit containing about 30wt% salt, the adhered salt should be removed by cathode process (CP). The CP has been regarded as one of the bottle-neck of the pyroprocess as the large amount of uranium is treated in this step and the operation parameters are crucial to determine the final purity of the product. Currently, related research activities are mainly based on experiments consequently it is hard to observe processing variables such as temperature, pressure and salt gas behavior during the operation of the cathode process. Hence, in this study operation procedure of cathode process is numerically described by using appropriate mathematical model. The key parameters of this research are the amount of evaporation at the distillation part, diffusion coefficient of gas phase salt in cathode processor and phase change rate at condensation part. Each of these conditions were composed by Hertz-Langmuir equation, Chapman-Enskog theory, and interphase mass flow application in ANSYS-CFX. And physical properties of salt were taken from the data base in HSC Chemistry. In this study, calculation results on the salt gas behavior and optimal operating condition are discussed. The numerical analysis results could be used to closely understand the physical phenomenon during CP and for further scale up to commercial level.

• 유병욱(충남대학교, 대전광역시 유성구 대학로 99) | Bung Uk Yoo (Chungnam National University, 99 Daehak-ro, Yuseong-gu, Daejeon, Korea)
• 박성빈(한국원자력연구원, 대전광역시 대덕대로 989번길 111) | Sung Bin Park (Korea Atomic Energy Research Institute, 111 Daedeok-daero 989, Yuseong-gu, Daejeon, Korea)
• 권상운(한국원자력연구원, 대전광역시 대덕대로 989번길 111) | Sang Woon Kwon (Korea Atomic Energy Research Institute, 111 Daedeok-daero 989, Yuseong-gu, Daejeon, Korea)
• 김정국(한국원자력연구원, 대전광역시 대덕대로 989번길 111) | Jeong Guck Kim (Korea Atomic Energy Research Institute, 111 Daedeok-daero 989, Yuseong-gu, Daejeon, Korea)
• 이한수(충남대학교, 대전광역시 유성구 대학로 99) | Han Soo Lee (Chungnam National University, 99 Daehak-ro, Yuseong-gu, Daejeon, Korea)
• 김인태(한국원자력연구원, 대전광역시 대덕대로 989번길 111) | In Tae Kim (Korea Atomic Energy Research Institute, 111 Daedeok-daero 989, Yuseong-gu, Daejeon, Korea)
• 이종현(충남대학교, 대전광역시 유성구 대학로 99) | Jong Hyeon Lee (Chungnam National University, 99 Daehak-ro, Yuseong-gu, Daejeon, Korea) Corresponding author