• Nuclear Engineering and Technology
• /
• v.52 no.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.

• Journal of Ocean Engineering and Technology
• /
• v.21 no.4
• /
• pp.8-14
• /
• 2007

A particle method recognized as one of gridless methods has been developed to investigate incompressible viscous flaw. The method is more feasible and effective than conventional grid methods for solving the flaw field with complicated boundary shapes or multiple bodies. The method is consists of particle interaction models representing pressure gradient, diffusion, incompressibility and the boundary conditions. In the present study, the models in case of various simulation condition were checked with the analytic solution, and applied to the two-dimensional Poiseuille flow in order to validate the developed method.

• Journal of Ocean Engineering and Technology
• /
• v.22 no.2
• /
• pp.20-27
• /
• 2008

A moon-pool is a vertical well in a floating barge, drilling ship, or offshore support vessel. In this study, numerical simulation of two-dimensional moon-pool flaw coupled with a ship's motion in waves is carried out using a particle method, the so-called MPS method. The particle method, which is recognized as one of the gridless methods, was developed to investigate nonlinear free-surface motions interacting with structures. The method is more feasible and effective than convectional grid methods in order to solve a flaw field with complicated boundary shapes.

• Nuclear Engineering and Technology
• /
• v.53 no.3
• /
• pp.833-841
• /
• 2021

During the postulated severe accident of nuclear reactor, eutectic reaction leads to low-temperature melting of fuel cladding and early failure of core structure. In order to model eutectic melting with the moving particle semi-implicit (MPS) method, the eutectic reaction model is developed to simulate the eutectic reaction phenomenon. The coupling of mass diffusion and phase diagram is applied to calculate the eutectic reaction with the uniform temperature. A heat transfer formula is proposed based on the phase diagram to handle the heat release or absorption during the process of eutectic reaction, and it can combine with mass diffusion and phase diagram to describe the eutectic reaction with temperature variation. The heat transfer formula is verified by the one-dimensional melting simulations and the predicted interface position agrees well with the theoretical solution. In order to verify the eutectic reaction models, the eutectic reaction of uranium and iron in two semi-infinite domains is simulated, and the profile of solid thickness decrease over time follows the parabolic law. The modified MPS method is applied to calculate Transient Reactor Test Facility (TREAT) experiment, the penetration rate in the simulations are agreeable with the experiment results. In addition, a hypothetical case based on the TREAT experiment is also conducted to validate the eutectic reaction with temperature variation, the results present continuity with the simulations of TREAT experiment. Thus the improved method is proved to be capable of simulating the eutectic reaction in the severe accident.

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.1
• /
• pp.1-9
• /
• 2008

The free-surface motions interacting with structures are investigated numerically using the Moving Particle Semi-implicit (MPS) method proposed by Koshizuka et al. (1996) for solving incompressible flow. In the method, Lagrangian moving particles are used instead of Eulerian approach using grid system. Therefore the terms of time derivatives in Navier-Stokes equation can be directly calculated without any numerical diffusion or instabilities due to the fully Lagrangian treatment of fluid particles and topological failure never occur. The MPS method is applied to the numerical study on the fluid impact loads for wet-drop tests in a LNG tank, and the results are compared with experimental ones.

• Journal of the Korea Convergence Society
• /
• v.10 no.2
• /
• pp.201-206
• /
• 2019

The flooding and overtopping due to unexpected large ocean wave may occur serious problems to environments and structures. Generally fixed wave breakers and several structures were installed to prevent such damages, however, they may affect to environments and charted path of ships badly. In this regard, new type of initial wave breaker was investigated in both of experimentally and numerically. For the experiments, conceptual devices were built by authors with stereo vision system. The moving particle semi-implicit method was adopted for simulation. It is revealed that the initial wave breaker reduce the damages from ocean waves by energy dissipated earlier. Furthermore, the effects of position of the initial wave breaker was also considered.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.24 no.1
• /
• pp.119-125
• /
• 2018

The particle based computational fluid dynamics (CFD) method, which follow Lagrangian approach for fluid dynamics, fluid particle behavior by tracking all particle calculation physical quantities of each particle. According to basic concept of particle based CFD method, it is difficult to satisfy continuum theory and measure influences from neighboring particle. Article number density and weight function were used to solve aforementioned issue. Difficulties continuum mean simulate non-continuum particles such as solid including granular and sand. In this regard, the particle based CFD method modified solid particle problems by replacing viscous and surface tension forces friction and drag forces. In this paper, particle interaction model for solid particle friction model implemented to simulate solid particle problems. The broken dam problem, which is common to verify particle based CFD method, used fluid or solid particles. The angle of repose was observed in the simulation results the solid particle not fluid particle.

• Nuclear Engineering and Technology
• /
• v.53 no.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.

• Ecology and Resilient Infrastructure
• /
• v.7 no.4
• /
• pp.336-344
• /
• 2020

The prediction and evaluation of driftwood accumulation around river-crossing structures are essential because driftwood accumulation increases during flood disasters. In this study, the driftwood accumulation and behavior around bridge piers were evaluated via a numerical model that could be employed to analyze three-dimensional turbulent flow and driftwood motion. The moving particle semi-implicit-based model for driftwood motion was sensitive to the number of spheres. The numerical results showed that the approach velocity and the ratio of driftwood length to pier width were the key factors influencing driftwood accumulation, whereas the driftwood density had only a minor influence. Overall, it is expected that this study will contribute to the development of improved risk evaluation indexes for assessing driftwood accumulation around river-crossing structures.

• Nuclear Engineering and Technology
• /
• v.54 no.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.