Development of GPU-Based SPH-DEM Code for Hydrodynamic Interactions between Multi-Phase Flow and Solid Debris
저자
발행사항
서울 : 서울대학교 대학원, 2020
학위논문사항
학위논문(박사)-- 서울대학교 대학원 : 에너지시스템공학부 2020. 8
발행연도
2020
작성언어
영어
주제어
DDC
622.33 판사항(22)
발행국(도시)
서울
기타서명
다상유동 및 고체 파편물의 수력학적 거동 해석을 위한 GPU 기반의 SPH-DEM 연계해석 코드 개발
형태사항
xi, 146 p. : 삽화, 표 ; 26 cm
일반주기명
참고문헌 수록
UCI식별코드
I804:11032-000000162918
소장기관
In the late phase severe accident of LWR, the massive corium releases out of the reactor pressure vessel (RPV) and falls to the coolant if the In-Vessel Retention (IVR) strategy fails. The melt jet can be fragmented into debris particles based on the assumption that the ex-vessel pool is sufficiently deep. It is known that there are various three-phase flow issues associated with the fragmented debris particles under the influence of phase change of cavity coolant. In such cases, the vapor phase forms a sharp and dynamic interface with the liquid phase while the transient relocation behavior of debris particles is the main concern. Thus, coupling Lagrangian-based multi-phase CFD techniques and Discrete Element Method (DEM) can be an effective approach in terms of numerical modeling of such behaviors. In this respect, an integrated numerical code for incompressible 3-phase flow has been developed in this study by two-way phase coupling of multi-phase Smoothed Particle Hydrodynamics (SPH) and DEM model.
Smoothed Particle Hydrodynamics (SPH) is one of the best-known meshless CFD methods in which the fluid system is represented as the finite number of Lagrangian particles. The SPH code developed in this study proposes a new density estimation model and improves the surface tension model for accurate simulation of incompressible two-phase flow behavior. The demonstration of its applicability has been performed through several V&V simulations including multi-phase dam-break and sloshing simulations.
Discrete Element Method (DEM) is a direct simulation method for a rigid body that can analyze the translation, rotation, and collision behavior of solid particles in detail. The soft-sphere collision model with Hertz-Mindlin contact force equations is adopted for developed DEM code in this study. To precisely estimate the wall boundary interactions of bed-formed debris particles, a versatile wall boundary model is newly proposed in this study that also covers the sliding and rolling behavior of solid particles. The inter-particle collision behavior and sliding & rolling motion of particles are well proven in several V&V cases.
The numerical code system for incompressible 3-phase flow is newly developed by two-way phase coupling of the above two models (SPH-DEM). The unresolved coupling approach between two methods was adopted for the analysis of the overall behavior of particulate solid debris. The fundamental validation of the phase coupled model was performed for both single-particle behavior and particulate granular flow such as dam-breaking motion of particle-fluid.
The SPH-DEM coupled code in this study has been parallelized based on Graphical Process Unit (GPU) in order to overcome the inherent efficiency problem of the Lagrangian-based numerical method. Parallel mapping and reduction are applied for solving discretized summation equations of each SPH particle, solving contact force equations for each DEM particle, and also for solving coupling equations between SPH and DEM particles. The efficiency of code parallelization was evaluated through the scalability analysis based on the benchmark calculation.
Finally, the simulation of the vapor-driven leveling behavior of spherical solids was performed as a case study to demonstrate the applicability of the developed code. The time-variant surface shape of solid particles was compared with the benchmark experiments both qualitatively and quantitatively. The effect of gas flow rate on the tendency of leveling behavior also has been analyzed.
The developed numerical system in this study is expected to be a good alternative for the simulation of such phenomena that were difficult to handle with traditional numerical methods since the numerical schemes used in the code have a high potential for simulation of complicated physics with highly deformable geometry. In addition, this validated code system can contribute to hydrodynamic modeling studies for severe accident technology by performing numerical experiments on conditions that hard to be conducted on a laboratory scale.
가압경수로 중대사고 후기 과정에서 핵연료 용융물 노내 억류 전략의 실패로 용융 핵연료가 용기 바깥쪽 침수 공동으로 재배치되면, 공동의 깊이나 제트의 속도 등의 조건에 따라 핵연료 파편화가 발생할 수 있다. 이러한 핵연료 파편물에서 발생하는 붕괴열에 의해 냉각재의 상변화가 이루어지면, 이상유동 자연대류와 고체 파편물이 공존하는 3상유동 시스템이 된다. 핵연료 용융물의 제트 파쇄, 하강 및 퇴적, 평탄화 등 일련의 과정에서 다양한 형태의 파편물 3상 냉각 거동이 발생할 수 있으며, 핵연료 파편물의 붕괴열이 충분히 제거되지 못하면 파편물의 응집 및 재용융이 발생할 수 있다. 중대사고 완화의 관점에서 이러한 3상 냉각 거동에 대한 상세한 이해가 필요하지만, 이상유동 자연대류, 이상유동 열전달, 고화, 용융, 비등, 유동-입자 상호작용 등을 포함하는 거동 자체의 복잡성으로 인해 관련 현상들의 예측 및 평가는 큰 불확실성을 내포하고 있다.
전통적인 핵연료 파편물을 포함한 3상 거동 해석은 다른 중대사고 해석과 마찬가지로 경험적 방법론에 의존한 보수적인 접근 방법이나 고체 입자와 이상 유체를 모두 연속체로 가정하는 다유체 모델을 기반으로 이루어졌다. 최근에는 이러한 방법론들이 가지는 본질적인 한계를 극복하고자 고체 입자들 사이의 충돌이나 회전을 별개로 다루는 이산요소법(Discrete Element Method, DEM)과 격자 기반의 오일러리안 전산유체해석(CFD) 기법을 연계하는 형태로 서로 간의 상호작용을 해석하는 연구들이 많이 이루어지고 있다.
한편, 최근에는 하드웨어 및 소프트웨어의 성능이 비약적으로 좋아지면서 격자(Grid)에 기반하지 않고 질점 하나하나의 움직임을 따라가면서 유동에 대한 지배방정식을 해석하는 라그랑지안 유체 해석기법의 응용이 늘어나고 있다. 라그랑지안 해석 기법에서는 다상유동 해석 시에 액체와 기체 두 상을 완전히 별개의 영역으로 해상하여 지배방정식을 풀기 때문에 계면 마찰력이나 항력, 양력 등에 대한 별도의 상관식 없이 제 1 원리 기반으로 유동을 해석할 수 있어, 이상유동에 대한 보다 근본적인 해석이 가능하다.
핵연료 파편물의 생성 및 3상 냉각 거동과 관련된 현상들은 대부분 기체 상의 생성이나 이상유동 자연대류에 영향을 받는 현상들로 액체 기체 사이의 인터페이스가 복잡하고 역동적인 경향이 있기 때문에, 액체-기체 다상유동에 효과적인 라그랑지안 기반의 유체해석 기법과 강체 이산요소법(DEM)을 연계하면 효과적인 3상 유동 해석 체계를 구축할 수 있다. 하지만, 핵연료 파편물을 포함한 3상 냉각 거동과 관련하여 라그랑지안 입자 기반 유체해석 기법을 활용한 연구는 세계적으로도 아직 수행된 바가 없다.
이러한 필요성에 따라, 본 연구에서는 대표적인 입자 기반의 유체해석 방법론 중 하나인 완화입자유체동역학(Smoothed Particle Hydrodynamics, SPH) 기법과 강체의 충돌, 병진, 회전 운동을 직접적으로 다루는 이산요소법(DEM)의 연계를 통해 고체 입자를 포함한 3상 유동 해석을 위한 라그랑지안 해석 체계를 구축하였다. 고체 입자와 유체 사이의 연계는 mm 이하의 스케일을 가지는 핵연료 파편물의 형상적 특성을 고려하여 두 상 사이의 겹침을 허용하여 운동량 교환을 모델링하는 비해상(unresolved) 방식으로 이루어졌다. 또한, SPH 유체 모델, DEM 강체 모델, SPH-DEM 연계 모델 각각에 대한 검증을 다양한 스케일에서 다양한 실험 연구들과의 비교를 통해 수행하였다.
한편, 완화입자유체동역학(SPH) 기법과 이산요소법(DEM) 모델은 라그랑지안 해석 기법의 특성상 컴퓨터 성능의 비약적인 발전에도 불구하고 오일러리안 해석 기법에 비해 계산 효율 및 시간에 대한 상대적인 제약이 존재한다. 특히 액체-기체의 이상유동 해석을 다룰 경우 기체 상의 밀도가 작기 때문에 라그랑지안 유체해석 기법에서는 더 작은 시간 간격이 요구된다. 이에 본 연구에서는 그래픽 처리 장치 (Graphics Processing Unit, GPU)를 활용하여 SPH 해석, DEM 해석, SPH-DEM 연계 해석이 모두 각 라그랑지안 질점에 대해 동시에 수행될 수 있도록 GPU 기반의 연계코드 병렬화 및 가속화를 수행하였다.
마지막으로 개발한 라그랑지안 기반의 3상유동 해석 체계의 유용성 입증을 위해 원자로 중대사고 후기 과정에서 발생할 수 있는 핵연료 파편물 층(debris bed)의 평탄화(self-leveling) 거동에 대한 검증 해석을 수행하였다. 시간에 따른 파편물 층 표면의 형상 변화를 타 기관에서 수행된 기체주입 실험과 비교하는 형태로 검증이 이루어졌다. 분석 결과, 본 연구에서 개발한 SPH-DEM 연계해석 코드가 고체 입자 상을 포함한 수력학적 3상 거동을 정량적으로, 정성적으로 잘 해석하는 것을 확인하였다.
본 연구에서 개발한 라그랑지안 기반의 SPH-DEM 3상유동 해석 체계는 원자로 중대사고의 해석적 관점에서 기존의 수치해석 기법들이 다루기 어려웠던 현상들에 대한 대안 또는 상호 보완의 역할을 할 수 있다. 또한, 본 연구에서 개발한 코드는 제 1원리 기반의 물리 법칙을 기반으로 유동 및 강체의 거동을 해석하기 때문에 실험으로 구현하기 어려운 조건이나 스케일에 대한 수치적 재현이 가능하고, 이를 바탕으로 기존의 스케일링 법칙을 검증하거나 실험 결과가 없는 영역에서 수치 데이터를 생성하여 기존의 상관식을 개선하는데 활용할 수 있다. 이러한 점에서 본 연구는 원자로 중대사고의 해석이나 안전성 평가와 관련하여 불확실성을 저감하는데 기여한다.
분석정보
서지정보 내보내기(Export)
닫기소장기관 정보
닫기권호소장정보
닫기오류접수
닫기오류 접수 확인
닫기음성서비스 신청
닫기음성서비스 신청 확인
닫기이용약관
닫기학술연구정보서비스 이용약관 (2017년 1월 1일 ~ 현재 적용)
학술연구정보서비스(이하 RISS)는 정보주체의 자유와 권리 보호를 위해 「개인정보 보호법」 및 관계 법령이 정한 바를 준수하여, 적법하게 개인정보를 처리하고 안전하게 관리하고 있습니다. 이에 「개인정보 보호법」 제30조에 따라 정보주체에게 개인정보 처리에 관한 절차 및 기준을 안내하고, 이와 관련한 고충을 신속하고 원활하게 처리할 수 있도록 하기 위하여 다음과 같이 개인정보 처리방침을 수립·공개합니다.
주요 개인정보 처리 표시(라벨링)
목 차
3년
또는 회원탈퇴시까지5년
(「전자상거래 등에서의 소비자보호에 관한3년
(「전자상거래 등에서의 소비자보호에 관한2년
이상(개인정보보호위원회 : 개인정보의 안전성 확보조치 기준)개인정보파일의 명칭 | 운영근거 / 처리목적 | 개인정보파일에 기록되는 개인정보의 항목 | 보유기간 | |
---|---|---|---|---|
학술연구정보서비스 이용자 가입정보 파일 | 한국교육학술정보원법 | 필수 | ID, 비밀번호, 성명, 생년월일, 신분(직업구분), 이메일, 소속분야, 웹진메일 수신동의 여부 | 3년 또는 탈퇴시 |
선택 | 소속기관명, 소속도서관명, 학과/부서명, 학번/직원번호, 휴대전화, 주소 |
구분 | 담당자 | 연락처 |
---|---|---|
KERIS 개인정보 보호책임자 | 정보보호본부 김태우 | - 이메일 : lsy@keris.or.kr - 전화번호 : 053-714-0439 - 팩스번호 : 053-714-0195 |
KERIS 개인정보 보호담당자 | 개인정보보호부 이상엽 | |
RISS 개인정보 보호책임자 | 대학학술본부 장금연 | - 이메일 : giltizen@keris.or.kr - 전화번호 : 053-714-0149 - 팩스번호 : 053-714-0194 |
RISS 개인정보 보호담당자 | 학술진흥부 길원진 |
자동로그아웃 안내
닫기인증오류 안내
닫기귀하께서는 휴면계정 전환 후 1년동안 회원정보 수집 및 이용에 대한
재동의를 하지 않으신 관계로 개인정보가 삭제되었습니다.
(참조 : RISS 이용약관 및 개인정보처리방침)
신규회원으로 가입하여 이용 부탁 드리며, 추가 문의는 고객센터로 연락 바랍니다.
- 기존 아이디 재사용 불가
휴면계정 안내
RISS는 [표준개인정보 보호지침]에 따라 2년을 주기로 개인정보 수집·이용에 관하여 (재)동의를 받고 있으며, (재)동의를 하지 않을 경우, 휴면계정으로 전환됩니다.
(※ 휴면계정은 원문이용 및 복사/대출 서비스를 이용할 수 없습니다.)
휴면계정으로 전환된 후 1년간 회원정보 수집·이용에 대한 재동의를 하지 않을 경우, RISS에서 자동탈퇴 및 개인정보가 삭제처리 됩니다.
고객센터 1599-3122
ARS번호+1번(회원가입 및 정보수정)