• Nuclear Engineering and Technology
• /
• 제54권1호
• /
• pp.162-176
• /
• 2022

In Lead-based reactor (LBR) severe accident, the meltdown and migration inside the reactor core will lead to fuel fragment concentration, which may further cause re-criticality and even core disintegration. Accurately predicting the migration and solidification behavior of melt in LBR severe accidents is of prime importance for safety analysis of LBR. In this study, the Moving Particle Semi-implicit (MPS) method is validated and used to simulate the migration and solidification behavior. Two main surface tension models are validated and compared. Meanwhile, the MPS method is validated by the L-plate solidification test. Based on the improved MPS method, the migration and solidification behavior of melt in LBR severe accident was studied furthermore. In the Pb-Bi coolant, the melt flows upward due to density difference. The migration and solidification behavior are greatly affected by the surface tension and viscous resistance varying with enthalpy. The whole movement process can be divided into three stages depending on the change in velocity. The heat transfer of core melt is determined jointly by two heat transfer modes: flow heat transfer and solid conductivity. Generally, the research results indicate that the MPS method has unique advantage in studying the migration and solidification behavior in LBR severe accident.

• Nuclear Engineering and Technology
• /
• 제52권6호
• /
• pp.1172-1179
• /
• 2020

As a Lagrangian particle method, Moving Particle Semi-implicit (MPS) method has great capability to capture interface/surface. In recent years, the multiphase flow simulation using MPS method has become one of the important directions of its developments. In this study, some key methods for multiphase flow have been introduced. The interface tension model in multiphase flow is modified to maintain the smooth of the interface and suitable for the three-phase flow. The mass transfer at immiscible liquid interface entrained by single bubble which could occur in Molten Core-Concrete Interaction (MCCI) has been investigated using this particle method. With the increase of bubble size, the height of entrainment column also increases, but the time of film rupture is slightly different. With the increase of density ratio between the two liquids, the height of entrained column decreases significantly due to the decreasing buoyancy of the denser liquid in the lighter liquid. In addition, the larger the interface tension coefficient is, the more rapidly the entrained denser liquid falls. This study validates that the MPS method has shown great performance for multiphase flow simulation. Besides, the influence of physical parameters on the mass transfer at immiscible interface has also been investigated in this study.

• 한국해양공학회지
• /
• 제31권6호
• /
• pp.388-396
• /
• 2017

A Lagrangian approach based computational fluid dynamics (CFD) was used to simulate large and/or sharp deformations and fragmentations of interfaces, including free surfaces, through tracing each particle with physical quantities. According to the concept of the particle-based CFD method, it is possible to apply it to both fluid particles and solid particles such as sand, gravel, and rock. However, the presence of more than two different phases in the same domain can make it complicated to calculate the interaction between different phases. In order to solve multiphase problems, particle interaction models for multiphase problems, including surface tension, buoyancy-correction, and interface boundary condition models, were newly adopted into the moving particle semi-implicit (MPS) method. The newly developed MPS method was used to simulate a typical validation problem involving dam breaking. Because the soil and other particles, excluding the water, may have different viscosities, various viscosity coefficients were applied in the simulations for validation. The newly developed and validated MPS method was used to simulate the mobile beds induced by broken dam flows. The effects of the viscosity on soil particles were also investigated.

• Nuclear Engineering and Technology
• /
• 제53권6호
• /
• pp.1810-1820
• /
• 2021

Study of single bubble behavior under rolling motions can prove useful for fundamental understanding of flow field inside the modern small modular nuclear reactors. The objective of the present study is to simulate the influence of rolling conditions on single rising bubble in a liquid using multiphase Moving Particle Semi-implicit (MPS) method. Rolling force term was added to 2D Navier-Stokes equations and a computer program was written using C language employing OpenACC to port the code to GPU. Computational results obtained were found to be in good agreement with the results available in literature. The impact of rolling parameters on trajectory and velocity of the rising bubble has been studied. It has been found that bubble rise velocity increases with rolling amplitude due to modification of flow field around the bubble. It has also been concluded that the oscillations of free surface, caused by rolling, influence the bubble trajectory. Furthermore, it has been discovered that smaller vessel width reduces the impact of rolling motions on the rising bubble. The effect of liquid viscosity on bubble rising under rolling was also investigated and it was found that effects of rolling became more pronounced with the increase of liquid viscosity.

• 한국전산유체공학회지
• /
• 제16권2호
• /
• pp.1-10
• /
• 2011

A new computational program, which is based on the CIP/CCUP(Constraint Interpolation Profile/CIP Combined Unified Procedure) method, has been developed to numerically analyse sloshing phenomena dealt as multiphase-flow problems. For the convection terms of Navier-Stokes equations, the RCIP(Rational function CIP) method was adopted and the THINC-WLIC(Tangent of Hyperbola for Interface Capturing-Weighted Line Interface Calculation) method was used to capture the air/water interface. To validate the present numerical method, two-dimensional dam-breaking and sloshing problems in a rectangular tank were solved by the developed method in a stationary Cartesian grid system. In the case of sloshing problems, simulations by using a improved MPS(Moving Particle Simulation) method, which is named as PNU-MPS(Pusan National University-MPS), were also carried out. The computational results are compared with those of experiments and most of the comparisons are reasonably good.

• 한국해양환경ㆍ에너지학회지
• /
• 제20권1호
• /
• pp.37-44
• /
• 2017

하나의 시스템 내에 2개 이상의 상이 다른 유체가 존재할 시에는 다상유동에 의한 복장성이 존재하며, 이는 해석의 어려움이 따른다. 두 개 이상의 상이 다른 다상유동은 유동 및 경계면에 영향을 끼치지 때문에, 불안정성과 같은 비선형 유동이 나타나게 된다. 여러 종류의 불안정성 중 레일리히-테일러 불안정성은 대표적인 예로 알려져 있다. 본 연구에서는 밀도차가 레일리히-테일러 불안정성에 미치는 영향을 조사하기 위해 다양한 Atwood 수를 선정하였으며, 초기 경계면 형상 역시 다양한 형태를 설정하고 시뮬레이션 하였다. 본 연구에서 사용된 입자법인 MPS(Moving particle simulation)은 이러한 다상유동에서 널리 쓰이지는 않았으나, 다상유동을 위한 입자간 상호 연성 모델인 자가-부력 항, 표면 장력 항과 경계면 경계 조건 항을 추가로 사용하여 수치해석이 가능하게 하였다. 본 연구에서 새로이 개발된 다상유동형 입자법을 이용하여 고려된 경우들에 대해 수치해석을 수행하였으며, 각각의 결과들을 비교 분석하였다. 또한 레일리히-테일러 불안정성에 기인한 유동의 속도를 측정하여 포텐셜 기반의 이론값과의 비교를 통해 경향성이 일치함을 알 수 있었다. 이론값과의 크기의 차는 포텐셜 기반의 이론값에서는 고려가 힘든 비선형성에 기인한다고 사료된다.

• 청정기술
• /
• 제25권3호
• /
• pp.206-222
• /
• 2019

최근에는 다양한 다상오염물 거동 흐름 모델들이 개발되었고 일부는 상용화되기도 하였으나, 대부분이 압력기저접근방식을 갖고 개발된 프로그램들이므로 다양한 수치적 어려움을 내재하고 있다. 이러한 수치적 어려움을 극복하기 위해서는 분율흐름접근방식을 따르는 기존 다상흐름거동 수치모델로 개발된 MultiPhaSe flow (MPS) 모델에 계면활성제에 의한 용해 현상을 모사할 수 있는 오염물 거동 모듈을 결합해서 MultiPhaSe flow and TranSport (MPSTS) 프로그램을 본 연구에서 개발하였다. 개발된 모델은 Clement의 해석 해를 사용하여 검증하였다. 여기서 MPSTS프로그램은 입자추적법과 결합한 라그랑지안-율러리안 기법을 이용해서 상간물질전달 효과와 다상내 오염물 거동 기능을 결합한 계면활성제 활용 복원과정을 모사할 수 있는 프로그램이다. 본 연구에서는 개발된 모델을 이용해서 소수성 액체(non aqueous phase liquid, NAPL)로 오염된 지역의 계면활성제에 의한 오염 정화 시 층상구조를 가지는 수리지질학적 불 균질성이 복원효율에 미치는 영향을 수치 모의 하였다. 수치모의 결과, 하부 층의 수리전도도가 상부 층의 수리전도도보다 10배, 20배, 50배로 큰 경우에 대해서 하부에서 물속에 용해된 디젤의 농도가 높게 나타난다. 왜냐하면 계면활성제가 하부 층을 따라서 좀 더 빨리 움직여서 하부 층에서 잔류 소수성 액체를 좀 더 많이 용해시켰기 때문이다.