Development of Dynamic LES-SPH Method for Thermal Ablation Analysis with Turbulent Natural Convection

Abstract

The primary objective in nuclear severe accident mitigation strategies is to retain and efficiently cool corium to prevent the release of fission products. However, the decay heat generated from corium can lead to thermal ablation of the pressure boundary, compromising structural integrity and potentially leading to the release of fission product into the environment. Prominent examples include RPV ablation in IVR-ERVC strategy, concrete structure breaching in MCCI and thermal erosion of sacrificial material caused by corium jet impingement above core-catcher. The extent of thermal ablation is determined by the flow characteristics of corium and heat flux at boundary. Therefore, a comprehensive understanding of thermal ablation and the behavior of corium is crucial from a severe accident mitigation perspective.

The analysis of turbulent natural convection behavior of corium and thermal ablation has been based on empirical approaches using simulant experiments or conventional computational fluid dynamics (CFD). Researches has been conducted to refine and validate of turbulence models for analyzing the behavior of corium flow, and recently, comprehensive analysis incorporating grid-based phase change analysis has been carried out. Numerical assumptions have been introduced, to simulate phenomena such as transient interface tracking between corium layers and domain deformation due to phase change.

With advancements in hardware and software for large-scale parallel computing, particle-based CFD methods, such as Smoothed Particle Hydrodynamics (SPH), have been applied to nuclear safety analysis. The particle-based CFD methods analyze flow using moving particles with physical quantities making it suitable for analyzing multi-phase, multi-fluid, free surface and phase change effectively.

Thermal ablation and corium thermal flow involve solid-liquid phase change and multi-fluid heat transfer analysis, making the particle-based CFD an effective framework by utilizing advanced phase change model. However, there has been no cases of analysis using sophisticated turbulence models based on particle-based CFD, and therefore, no research on thermal ablation accompanied by turbulent thermal flow.

Therefore, in this study, a turbulence model applicable to transient laminar, transitional, and turbulent flow conditions was developed based on Smoothed Particle Hydrodynamics (SPH), which is a representative particle-based fluid analysis method. By combining the common features of SPH, which uses spatial weight functions, and Large Eddy Simulation (LES), which utilizes spatial filter functions for modeling turbulent flow, a dynamic LES-SPH framework was developed without additional filtering processes. Numerical corrections were applied to SPH discretization to ensure the LES models numerical accuracy when analyzing low Prandtl number fluids, such as corium and liquid metals. Turbulent Prandtl number modeling and the dynamic Vreman model were introduced for this purpose. The capabilities of the developed dynamic LES-SPH model to analyze thermal flow under transitional and turbulent conditions were validated through comparisons with various numerical and experimental studies.

Furthermore, an advanced phase change model based on the Enthalpy-Porosity Model (EPM) was developed to analyze thermal ablation using the SPH method. By incorporating the assumption of porosity effects in the momentum equation, the EPM enabled the analysis of phase change without the need for numerical modification related to viscosity and other properties. The developed SPH phase change model was validated through comparisons with various numerical and experimental studies, confirming its capability to analyze phase interfaces and heat transfer.

To demonstrate the utility of the developed LES-SPH with phase change model, an analysis was conducted on the behavior of corium pool and the thermal ablation of the pressure vessel under in-vessel retention conditions. The thermal ablation of the pressure vessel at different times and locations was compared to CFD and previous SPH approaches. The analysis results confirmed that the dynamic LES-SPH and phase change models developed in this study effectively analyze turbulent corium flow accompanied by thermal ablation.

The particle-based LES-SPH phase change framework developed in this study provides an analysis methodology for phenomena that were difficult to analyze using conventional CFD in the context of nuclear severe accidents. It is expected to complement conventional CFD methods and provide a more comprehensive understanding. Furthermore, rigorous validation and analysis of the developed LES-SPH and phase change models will be conducted to ensure their reliability. This will facilitate the improvement and proposal of correlations in areas without experimental data related to nuclear severe accidents, contributing to the field of nuclear safety.

원자력 중대사고 완화전략에서는 노심용융물의 노내/노외에서의 효율적인 냉각을 통해 사고 진전 완화 및 핵분열물질의 유출방지를 목표로 한다. 그러나 노심용융물 붕괴열로 인해 압력경계의 열적 침식이 발생하며, 구조건전성로 인해 핵분열 생성물이 노외로 유출될 수 있다. 대표적인 사례로는 노내억류 전략에서 압력용기 용발, 노외억류 전략에서 MCCI 및 코어캐처 희생물질 열적 침식이 있다. 이때 압력경계의 열적 침식 정도는 노심용융물 열유동과 경계에서의 열유속에 의해 결정된다. 따라서 중대사고 완화 관점에서 열적 침식 현상과 상경계에서 노심용융물 거동에 대한 종합적인 이해가 필요하다. 그러나 극한 조건을 동반한 현상의 특징으로 인해, 실험적 해석에는 제한이 있고, 난류 자연대류를 동반한 노심용융물 거동, 고체-액체 상변화, 성층화된 노심용융물간 열유동 등의 거동 자체의 복잡성으로 인해 예측 및 평가에 불확실성이 존재한다.

노심용융물의 난류 자연대류 거동 및 압력경계 열적 침식 해석은 다른 중대사고와 마찬가지로, 모사실험 기반 경험적 접근 방식이나 성층화된 노심용융물 층의 유동특성을 고려한 전통적 전산 유체해석 방법론에 기반하였다. 노심용융물의 거동해석을 위해 고도화된 난류모델 개발 및 검증 연구가 진행되고 있으며, 최근에는 격자기반 상변화 해석을 포함한 종합적인 해석연구가 수행되었다. 이때 성층화된 노심용융물 간 계면변화, 노심용융물의 자유표면 해석 및 상변화로 인한 해석영역의 변형 등의 영향을 모사하기 위해 수치적 가정들이 도입되었다.

한편, 최근 대규모 병렬계산을 위한 장비 및 소프트웨어의 발전에 힘입어, 입자법 기반 전산유체 해석 방법론이 원자력 안전 해석에 적용하는 경우가 늘어나고 있다. 입자법은 격자 기반이 아닌, 물리량을 가진 입자의 움직임을 통해 유동을 해석하여, 다유체/다상 경계를 해석할 수 있으며, 상변화 및 자유표면을 다루는데 강점이 있다.

열적 침식 및 노심용융물 거동 관련 물리현상은 고체-액체 상변화현상으로 인한 해석계면 변형은 물론, 다물질 혼합물인 노심용융물 내에서 밀도차에 의한 성층화로 인한 다유체 열유동 해석등을 포함한다. 이러한 특징 때문에 입자법 기반 상변화 모델을 활용한 전산유체 해석방법론을 통해 효율적인 해석체계를 구축할 수 있다. 그러나 입자법 기반 전산유체 해석 방법론에서는 고도화된 난류모델을 통한 해석이 진행된 사례가 없으며, 따라서 난류 열유동을 동반한 열적 침식에 대한 해석 연구가 진행된 사례가 없다.

따라서, 본 연구에서는 대표적인 입자 기반 유체 해석 방법론인 완화입자유체동역학 (Smoothed Particle Hydrodynamics, SPH) 기반으로, 과도적인 층류, 전이 및 난류 유동 조건에 범용적으로 적용할 수 있는 난류모델을 개발하였다. 공간 가중함수를 활용하는 SPH 차분화 방법과, 공간 필터함수를 활용해 난류 유동 모델링을 수행하는 대와류모델 (Large Eddy Simulation, LES)의 공통점을 이용해, 추가적인 필터링 과정이 필요 없는 LES-SPH 해석체계를 구축하였다. LES 모델의 수치 정확도 확보를 위해, SPH 차분화 수치 보정을 수행하였다. 또한 저 Pr 유체인 노심용융물 및 액체금속 열유동 해석을 위해, 난류 Pr수 모델링 및 동적 Vreman 모델을 도입하였다. 구축된 동적 LES-SPH 모델의 전이 및 난류 조건에서의 열유동 해석 능력 검증은 다양한 수치, 실험 연구들과 비교를 통해 수행되었다.

또한 양해법 SPH 해석체계를 활용한 열적 침식 해석을 위해, 엔탈피-다공성(Enthalpy-Porosity Model, EPM) 모델을 기반으로 상변화 모델을 고도화하였다. 상변화 구간에서 다공성 가정을 활용해 운동량 방정식 구성함으로써, 점성 등의 물성치에 대한 수치적 수정 없이 상변화 해석을 수행할 수 있다. 개발된 SPH 상변화 모델은 다양한 수치, 실험 연구와의 비교검증을 통해 상경계 및 열유동 해석 능력을 검증하였다.

마지막으로 개발된 LES-SPH 기반의 상변화 모델의 유용성 입증하기 위해 원자로 중대사고 완화전략인 노내억류 조건에서 노심용융물의 거동 및 압력용기의 열적 용발에 대한 해석을 수행하였다. 시간 및 위치에 따른 압력용기의 열적 용발 정도를 격자법 전산유체해석 방법론 및 선행 SPH 방법과 비교 검증하였다. 분석 결과, 본 연구에서 개발한 동적 LES-SPH 및 상변화 모델이 열적 침식을 동반한 난류 다유체 유동을 정성적으로 잘 해석하는 것을 확인하였다.

본 연구에서 개발한 입자법 기반 LES-SPH 상변화 해석체계는 원자로 중대사고 관점에서 기존 전산유체해석 기법을 통해 해석하기 어려웠던 현상들에 대한 직접 해석 방법론을 제안하였다는 의의가 있으며, 추후 전통적 유체해석 기법과 상호보완적 역할을 수행할 것으로 기대한다. 또한 개발된 LES-SPH 모델 및 상변화 모델에 대한 보다 엄밀한 검증 및 분석을 통해 신뢰성을 확보할 예정이다. 이를 통해 원자력 중대사고에서 실험결과가 없는 영역에 대한 상관식 개선 및 제안할 수 있으며, 원자력 안전 분야에 기여할 것으로 기대한다.