• Ocean Systems Engineering
• /
• v.8 no.2
• /
• pp.183-199
• /
• 2018

The open source software OpenFOAM is utilised to simulate the water entry and hydrodynamic impact process of 2D wedges and ship hull sections. Incompressible multiphase flow solver interDyMFoam is employed to calculate the free fall of structure from air into water using dynamically deforming mesh technique. Both vertical and oblique entry of wedges of various dead-rise angles have been examined. A convergence study of dynamics as well as kinematics of the flow problem is carried out on successively refined meshes. Obtained results are presented and compared to the experimental measurements showing good agreement and reasonable mesh convergence of the solution.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.20 no.1
• /
• pp.37-44
• /
• 2017

Complexity of multiphase flows due to existence of more than two interface including free-surface in one system, cannot be simulated easily. Since more than two fluids affect to flows and disturb interface, non-linearities such as instabilities can be appeared. Among several instabilities on multiphase flows, one of representative is Rayleigh-taylor instability. In order to examine in importance of density disparity, several cases with numerous Atwood number are set. Moreover, investigation of influence on initial disturbance were also considered. Moving particle simulation (MPS) method, which was employed in this paper, was not widely used for multiphase problem. In this study, by adding new particle interaction models such as self-buoyance correction, surface tension, and boundary condition at interface models, MPS were developed having more strength of physics and robust. By applying newly developed multiphase MPS, considered cases are performed and compared each other. Additionally, though existence of disagreement of magnitude of rising velocity between theoretical values from linear potential theory and that of numerical simulation, agreement of tendency can be proved of similarity of result. the discordance of magnitude can be explained due to non-linear effects on numerical simulation which was not considered in theoretical result.

• Journal of computational fluids engineering
• /
• v.16 no.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.

• Coupled systems mechanics
• /
• v.8 no.4
• /
• pp.339-350
• /
• 2019

A numerical simulation of the incompressible multiphase hydraulic jump flow was performed to compare the interface prediction through the use of the three RANS turbulence models: $k-{\varepsilon}$, $RNGk-{\varepsilon}$ and SST $k-{\omega}$. A three dimensional no submerged hydraulic jump and a two dimensional submerged hydraulic jump were modeled. Both the geometry and the mesh were created using the open source Gmsh code. The project's geometry consists of a rectangular channel with length and height differences between the two dimensional and three dimensional simulations. Uniform hexahedral cells were used for the mesh. Three refining meshes were constructed to allow to verify simulation convergence. The Volume of Fluid (abbr. VOF) method was used for treatment of the air-water surface. The turbulence models were evaluated in three distinct set up configurations to provide a greater accuracy in the flow representation. In the two-dimensional analysis of a submerged hydraulic jump simulation, the turbulence model RNG RNG $k-{\varepsilon}$ provided a better interface adjust with the experimental results than the model $k-{\varepsilon}$ and SST $k-{\omega}$. In the three-dimensional simulation of a no-submerged hydraulic jump the k-# showed better results than the SST $k-{\omega}$ and RNG $k-{\varepsilon}$ capturing the height and length of the ledge with a better fit with the experimental results.

• The KSFM Journal of Fluid Machinery
• /
• v.14 no.6
• /
• pp.45-53
• /
• 2011

Multiphase lattice Boltzmann method (LBM) has been used to simulate the nucleate pool boiling directly. For the phase change model, the thermal model and the Stefan boundary condition were introduced to the isothermal LBM. The phase change model was validated by the bubble growth in a superheated liquid under no gravity. The bubble growth on and departure from a superheated wall has been simulated successfully. The preliminary results showed that the detail process of nucleate pool boiling was in good agreement with the experimental results.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.1
• /
• pp.495-509
• /
• 2019

The objective of this study is to analyze the compressibility effects of multiphase cavitating flow during the water-entry process. For this purpose, the water-entry of a projectile at transonic speed is investigated computationally. A temperature-adjusted Tait equation is used to describe the compressibility effects in water, and air and vapor are treated as ideal gases. First, the computational methodology is validated by comparing the simulation results with the experimental measurements of drag coefficient and the theoretical results of cavity shape. Second, based on the computational methodology, the hydrodynamic characteristics of flow are investigated. After analyzing the cavitating flow in compressible and incompressible fluids, the characteristics under compressible conditions are focused upon. The results show that the compressibility effects play a significant role in the development of cavitation and the pressure inside the cavity. More specifically, the drag coefficient and cavity size tend to be larger in the compressible case than those in the incompressible case. Furthermore, the influence of entry velocities on the hydrodynamic characteristics is investigated to provide an insight into the compressibility effects on cavitating flow. The results show that the drag coefficient and the impact pressure vary with the entry velocity, and the prediction formulas for drag coefficient and impact pressure are established respectively in the present study.

• Geomechanics and Engineering
• /
• v.13 no.1
• /
• pp.1-23
• /
• 2017

A coupled liquid-gas-solid three-phase model, linking two numerical codes (TOUGH2/EOS3 and $FLAC^{3D}$), was firstly established and validated by simulating an in-situ air flow test in Essen. Then the coupled model was employed to investigate responses of multiphase flow and soil skeleton deformation to compressed air or freshwater injection using the same simulation conditions in an aquifer of Tianjin, China. The simulation results show that with injecting pressurized fluids, the vertical effective stress in some area decreases owing to the pore pressure increasing, an expansion of soil skeleton appears, and land uplift occurs due to support actions from lower deformed soils. After fluids injection stops, soil deformation decreases overall due to injecting fluids dissipating. With the same applied pressure, changes in multiphase flow and geo-mechanical deformation caused by compressed air injection are relatively greater than those by freshwater injection. Furthermore, the expansion of soil skeleton induced by compressed air injection transfers upward and laterally continuously with time, while during and after freshwater injection, this expansion reaches rapidly a quasi-steady state. These differences induced by two fluids injection are mainly because air could spread upward and laterally easily for its lower density and phase state transition appears for compressed air injection.

• Nuclear Engineering and Technology
• /
• v.55 no.8
• /
• pp.3017-3029
• /
• 2023

In the late in-vessel phase of a nuclear reactor severe accident, the internal heat transfer and crust evolution during the debris bed melting process have important effects on the thermal load distribution along the vessel wall, and further affect the reactor pressure vessel (RPV) failure mode and the state of melt during leakage. This study coupled the phase change model and large eddy simulation to investigate the variations of the temperature, melt liquid fraction, crust and heat flux distributions during the debris bed melting process in the hypothetical severe accident of HPR1000. The results indicated that the heat flow towards the vessel wall and upper surface were similar at the beginning stage of debris melting, but the upward heat flow increased significantly as the development of the molten pool. The maximum heat flux towards the vessel wall reached 0.4 MW/m². The thickness of lower crust decreased as the debris melting. It was much thicker at the bottom region with the azimuthal angle below 20° and decreased rapidly at the azimuthal angle around 20-50°. The maximum and minimum thicknesses were 2 and 90 mm, respectively. By contrast, the distribution of upper crust was uniform and reached stable state much earlier than the lower crust, with the thickness of about 10 mm. Moreover, the sensitivity analysis of initial condition indicated that as the decrease of time interval from reactor scram to debris bed dried-out, the maximum debris temperature and melt fraction became larger, the lower crust thickness became thinner, but the upper crust had no significant change. The sensitivity analysis of in-vessel retention (IVR) strategies indicated that the passive and active external reactor vessel cooling (ERVC) had little effect on the internal heat transfer and crust evolution. In the case not considering the internal reactor vessel cooling (IRVC), the upper crust was not obvious.

• 한국전산유체공학회:학술대회논문집
• /
• 2002.10a
• /
• pp.144-149
• /
• 2002

This paper presents the result of numerical simulation of particles trace following melted zinc movement with nitrogen gas injection. The code of the computational fluid dynamics for numerical analysis was performed using FLUENT related to CFD. As application model, there was applied Eulerian multiphase model for simulation of melted zinc movement at first and then was used stochastic tracking technique for particles trace secondarily. Numerical simulation results are shown that particles move to the same direction as the movement of melted zinc.

• Journal of Ocean Engineering and Technology
• /
• v.28 no.4
• /
• pp.280-287
• /
• 2014

In this paper, the numerical simulation of the gas-liquid flow in a cylinder cyclone separator is performed to investigate the flow characteristics using a commercial software, FLUENT, which solves the Reynolds-averaged Navier-Stokes(RaNS) equations. First, a single-phase flow with water in a gas-liquid cylinder cyclone(GLCC) separator is simulated and compared with the experiments(Farchi, 1990) and numerical simulations(Erdal, 1997). Then, the characteristics of the multi-phase flow for water-air, mud-only, and mud-air cases are discussed in the view point of the feasibilities for a mud handling system.