• 한국해양공학회지
• /
• 제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.

• Computers and Concrete
• /
• 제24권6호
• /
• pp.571-578
• /
• 2019

This study aimed to investigate the permeable nature of pervious concretes (PC) through the moving particle simulation (MPS) method. In the simulation, the complex structure of a pervious concrete was virtually demonstrated as a lattice model (LM) of spherical beads, where the test of permeability was conducted. Results of the simulation were compared with the experimental ones for validation. As a result, MPS results showed the permeability index of the LM as almost twice as big as the actual PCs. A proposed virtual model was created to prevent the stuck of water flow in the MPS simulation of PC or LM. Successful simulation results were demonstrated with the model.

• 한국해양공학회지
• /
• 제24권5호
• /
• pp.31-38
• /
• 2010

The Floating storage and re-gasification unit (FSRU), which has large cargo storage tanks, is a floating liquefied natural gas (LNG) import terminal. The sloshing motion in tanks that are partially filled with LNG can cause impact pressure on the containment system and affect the global motion of the FSRU. Therefore, the accurate prediction of sloshing motion has been a significant issue in the offshore gas production industry. In this paper, a particle method based on the moving particle semi-implicit (MPS) method proposed by Koshizuka and Oka (1996) has been modified to predict sloshing motion accurately in a rectangular tank with the filling ratio of water. The simulation results, including the violent sloshing of the fluid, were validated by comparison with the original MPS method.

• Nuclear Engineering and Technology
• /
• 제55권12호
• /
• pp.4395-4407
• /
• 2023

Bubble columns are widely encountered in several industries, especially in the field of nuclear safety. The Eulerian-Eulerian and the Eulerian-Lagrangian methods are commonly used to investigate bubble columns. Eulerian approaches require additional tasks such as strict volume conservation at the interface and a predefined well-structured grid. In contrast, the Lagrangian approach can be easily implemented. Hence, we introduce a fully Lagrangian approach for the simulation of bubble columns using the discrete bubble model (DBM) and moving particle semi-implicit (MPS) methods. Additionally, we propose a rigorous method to estimate the volume fraction accurately, and verified it through experimental data and analytical results. The MPS method was compared with the experimental data of Dambreak. The DBM was verified by analyzing the terminal velocity of a single bubble for each bubble size. It agreed with the analytical results for each of the four drag correlations. Additionally, the improved method for calculating the volume fraction showed agreement with the Ergun equation for the pressure drop in a packed bed. The implemented MPS-DBM was used to simulate the bubble column, and the results were compared with the experimental results. We demonstrated that the MPS-DBM was in quantitative agreement with the experimental data.

• Ocean Systems Engineering
• /
• 제7권3호
• /
• pp.299-318
• /
• 2017

The paper aims at illustrating several key issues and ongoing efforts for development of a reliable fully-Lagrangian particle-based solver for simulation of hydroelastic slamming. Fluid model is founded on the solution of Navier-Stokes along with continuity equations via an enhanced version of a projection-based particle method, namely, Moving Particle Semi-implicit (MPS) method. The fluid model is carefully coupled with a structure model on the basis of conservation of linear and angular momenta for an elastic solid. The developed coupled FSI (Fluid-Structure Interaction) solver is applied to simulations of high velocity impact of an elastic aluminum wedge and hydroelastic slammings of marine panels. Validations are made both qualitatively and quantitatively in terms of reproduced pressure as well as structure deformation. Several remaining challenges as well as important key issues are highlighted. At last, a recently developed multi-scale MPS method is incorporated in the developed FSI solver towards enhancement of its adaptivity.

• Ocean Systems Engineering
• /
• 제8권2호
• /
• pp.167-182
• /
• 2018

RTI (Rayleigh-Taylor instability) is investigated by a multi-liquid MPS (Moving Particle Semi-implicit) method for both viscous and inviscid flows for various density differences, initial-disturbance amplitudes, viscosities, and surface tensions. The MPS simulation can be continued up to the late stage of high nonlinearity with complicated patterns and its initial developments agree well with the linear theoretical results. According to the relevant linear theory, the difference between inviscid and viscous fluids is the rising velocity at which upward-mushroom-like RTI flow with vortex formation is generated. However, with the developed MPS program, significant differences in both growing patters and developing speeds are observed. Also, more dispersion can be observed in the inviscid case. With larger Atwood (AT) number, stronger RTI flows are developed earlier, as expected, with higher potential-energy differences. With larger initial disturbances, quite different patterns of RTI-development are observed compared to the small-initial-disturbance case. If AT number is small, the surface tension tends to delay and suppress the RTI development when it is sufficiently large. Interestingly, at high AT number, the RTI-suppressions by increased surface tension become less effective.

• 한국해양공학회지
• /
• 제28권2호
• /
• pp.93-101
• /
• 2014

In general, particle simulation methods such as the MPS(Moving Particle Simulation) or SPH(Smoothed Particle Hydrodynamics) methods have some serious drawbacks for pressure solutions. The pressure field shows spurious high fluctuations both temporally and spatially. It is well known that pressure fluctuation primarily occurs because of the numerical approximation of the partial differential operators. The MPS and SPH methods employ a pre-defined kernel function in the approximation of the gradient and Laplacian operators. Because this kernel function is constructed artificially, an accurate solution cannot be guaranteed, especially when the distribution of particles is irregular. In this paper, we propose a particle simulation method based on the moving least-square technique for solving the partial differential operators using a Taylor-series expansion. The developed method was applied to the hydro-static pressure and dam-broken problems to validate it.

• International Journal of Ocean System Engineering
• /
• 제1권1호
• /
• pp.37-45
• /
• 2011

The nonlinear free-surface motions interacting with a floating body were investigated using the Moving Particle Semi-implicit (MPS) method proposed by Koshizuka and Oka [6] for incompressible flow. In the numerical method, more realistic Lagrangian moving particles were used for solving the flow field instead of the Eulerian approach with a grid system. Therefore, the convection terms and time derivatives in the Navier-Stokes equation can be calculated more directly, without any numerical diffusion, instabilities, or topological failure. The MPS method was applied to a numerical simulation of predicting the efficiency of floating-type breakwater interacting with waves.

• 한국해양공학회지
• /
• 제27권4호
• /
• pp.90-97
• /
• 2013

Oil booms are one of the most widely used types of equipment for the protection of coastal areas against oil spills. In some situations, however, there are several types of oil leaks from the oil boom. Important factors regarding these phenomena include the surrounding ocean environment, such as waves, the density and viscosity of oil, the length of the oil boom skirt, etc. To estimate the performance of the oil boom, it is necessary to predict the behavior of the spilled oil and oil boom. In the present study, the prediction of oil boom performance in waves was carried out using the Pusan-National-University-modified Moving Particle Semi-implicit (PNU-MPS) method, which is an improved version of the original MPS proposed by Koshizuka and Oka (1996). The governing equations, which consist of continuity and Navier-Stokes equations, are solved by Lagrangian moving particles, and all terms expressed by differential operators in the governing equations are replaced by the particle interaction models based on a kernel function. The simulation results were validated through a comparison with the results of Violeau et al. (2007)..

• 한국해양공학회지
• /
• 제21권6호
• /
• pp.42-46
• /
• 2007

The impact of a single wave generated by a dam break with a tall structure is modeled with a three-dimensional version of the Moving particle semi-implicit (MPS) method. The particle method is more feasible and effective than methods based on grid connection problems involving violent free surface motions. In the present study, the Tsunami impact load and the change of longitudinal velocity component around the structure, which are obtained from the numerical simulation, are compared to those from experiments.