• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3B
• /
• pp.289-301
• /
• 2009

In this work, wave run-up heights and resultant wave forces on a vertical revetment due to tsunami (solitary wave) are investigated numerically using a numerical wave tank model called CADMAS-SURF (CDIT, 2001. Research and Development of Numerical Wave Channel (CADMAS-SURF). CDIT library, No. 12, Japan.), which is based on a 2-D Navier-Stokes solver, coupled to a volume of fluid (VOF) method. The third order approximate solution (Fenton, 1972. A ninth-order solution for the solitary wave. J. of Fluid Mech., Vol. 53, No.2, pp.257-271) is used to generate solitary waves and implemented in original CADMAS-SURF code. Numerical results of the wave profiles and forces are in good agreements with available experimental data. Using the numerical results, the regression curves determined from the least-square analysis are proposed, which can be used to determine the maximum wave run-up height and force on a vertical revetment due to tsunami. In addition, the capability of CADMAS-SURF is demonstrated for tsunami wave forces acting on an onshore structure using various configuration computations including the variations of the crown heights of the vertical wall and the position of the onshore structure. Based on the numerical results such as water level, velocity field and wave force, the direct effects of tsunami on an onshore structure are discussed.

• Journal of Ocean Engineering and Technology
• /
• v.35 no.1
• /
• pp.38-49
• /
• 2021

Slug flow is the most common multi-phase flow encountered in oil and gas industry. In this study, the hydrodynamic features of flow in pipes investigated numerically using computational fluid dynamic (CFD) simulations for the effect of slug flow on the vertical and bent pipeline. The compressible Reynold averaged Navier-Stokes (RANS) equation was used as the governing equation, with the volume of fluid (VOF) method to capture the outline of the bubble in a pipeline. The simulations were tested for the grid and time step convergence, and validated with the experimental and theoretical results for the main hydrodynamic characteristics of the Taylor bubble, i.e., bubble shape, terminal velocity of bubble, and the liquid film velocity. The slug flow was simulated with various air and water injection velocities in the pipeline. The simulations revealed the effect of slug flow as the pressure occurring in the wall of the pipeline. The peak pressure and pressure oscillations were observed, and those magnitudes and trends were compared with the change in air and water injection velocities. The mechanism of the peak pressures was studied in relation with the change in bubble length, and the maximum peak pressures were investigated for the different positions and velocities of the air and water in the pipeline. The pressure oscillations were investigated in comparison with the bubble length in the pipe and the oscillation was provided with the application of damping. The pressures were compared with the case of a bent pipe, and a 1.5 times higher pressures was observed due to the compression of the bubbles at the corner of the bent. These findings can be used as a basic data for further studies and designs on pipeline systems with multi-phase flow.

• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.3
• /
• pp.215-226
• /
• 2017

This paper presents numerical results of the performance of a marin propeller in cavitating and non-cavitating flow conditions. The geometry and experimental validation data of the propeller are provided in Potsdam Propeller Test Case(PPTC) in the framework of the second International Symposium on Marine Propulsors 2011(SMP'11) workshop. The PPTC includes open water tests, velocity field measurements and cavitation tests. The present numerical analysis was carried out by using the Reynolds averaged Navier-Stokes(RANS) method on a wall-resolved grid ensuring a y+=1, where the SST k-${\omega}$ model was mainly used for turbulence closure. The influence of the turbulence model was investigated in the prediction of the wake field under a non-cavitating flow condition. The propeller tip vortex flows in both cavitating and non-cavitating conditions were captured through adaptation of additional grids. For the cavitation flows at three operation points, Schnerr-Sauer's cavitation model was used with a Volume-Of Fluid(VOF) approach to capture the two-phase flows. The present numerical results for the propeller wake and cavitation predictions including the open water performance showed a qualitatively reasonable agreement with the model test results.

• Journal of Ocean Engineering and Technology
• /
• v.35 no.1
• /
• pp.50-58
• /
• 2021

The demand for eco-friendly energy is expected to increase due to the recently strengthened environmental regulations. In particular, the flow inside the pipe used in a cargo handling system (CHS) or fuel gas supply system (FGSS) of hydrogen transport ships and hydrogen-powered ships exhibits a very complex pattern of multiphase-thermal flow, including the boiling phenomenon and high accuracy analysis is required concerning safety. In this study, a feasibility study applying the boiling model was conducted to analyze the multiphase-thermal flow in the pipe considering the phase change. Two types of boiling models were employed and compared to implement the subcooled boiling phenomenon in nucleate boiling numerically. One was the "Rohsenow boiling model", which is the most commonly used one among the VOF (Volume-of-Fluid) boiling models under the Eulerian-Eulerian framework. The other was the "wall boiling model", which is suitable for nucleate boiling among the Eulerian multiphase models. Moreover, a comparative study was conducted by combining the nucleate site density and bubble departure diameter model that could influence the accuracy of the wall boiling model. A comparison of the Rohsenow boiling and the wall boiling models showed that the wall boiling model relatively well represented the process of bubble formation and development, even though more computation time was consumed. Among the combination of models used in the wall boiling model, the simulation results were affected significantly by the bubble departure diameter model, which had a very close relationship with the grid size. The present results are expected to provide useful information for identifying the characteristics of various parameters of the boiling model used in CFD simulations of multiphase-thermalflow, including phase change and selecting the appropriate parameters.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.2
• /
• pp.7-14
• /
• 2011

Numerical simulations were performed for the evaluation of wave and current loads on a fixed cylindrical substructure model for an ocean wind turbine using the ANSYS-CFX package. The numerical wave tank was actualized by specifying the velocity at the inlet and applying momentum loss as a wave damper at the end of the wave tank. The Volume-Of-Fluid (VOF) scheme was adopted to capture the air-water interface. An accuracy validation of the numerical wave tank with a truncated vertical circular cylinder was accomplished by comparing the CFD results with Morison's formula, experimental results, and potential flow solutions using the higher-order boundary element method (HOBEM). A parametric study was carried out by alternately varying the length and amplitude of the wave. As a meaningful engineering application, in the present study, three kinds of conditions were considered, i.e., cases with current, waves, and a combination of current and progressive waves, passing through a cylindrical substructure model. It was found that the CFD results showed reasonable agreement with the results of the HOBEM and Morison's formula when only progressive waves were considered. However, when a current was included, CFD gave a smaller load than Morison's formula.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.11
• /
• pp.871-884
• /
• 2015

Fluid transport is a key issue in the development of microfluidic systems. Recently, Myong (2014) has proposed a new concept for droplet transport without external power sources, and numerically validated the results for a hypothetical 2D shape, initially having a hemicylindrical droplet shape. Myong and Kwon (2015) have also examined the transport mechanism for an actual water droplet, initially having a 3D hemispherical shape, on a horizontal hydrophilic/hydrophobic surface, based on the numerical results of the time evolution of the droplet shape, as well as the total kinetic, gravitational, pressure and surface free energies inside the droplet. In this study, a 3D numerical analysis of an initially spherical droplet is carried out to establish a new concept for droplet transport. Further, the transport mechanism of an actual water droplet is examined in detail from the viewpoint of the capillarity force imbalance through the numerical results of droplet shape and various energies inside the droplet.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.36 no.6
• /
• pp.1037-1048
• /
• 2016

In this study, nonlinear wave interaction in the presence of a uniform current is studied using numerical model, named CADMAS-SURF which is based on the Navier-Stokes equations coupled with Volume of Fluid for tracking free surface deformation. The original CADMAS-SURF developed for interaction of wave with structure is modified/extended to simulate nonlinear fluid dynamic motions within wave-current coexisting field. The capability of Numerical Wave-Current Tank (NWCT) in this study is validated by comparing with available existing laboratory experiments for both wave-following and wave-opposing current. The numerical results for interaction between wave and current are shown to be in good agreement with experimental data. Then, this study focused on the dynamic motions of the water velocity, surface elevation and vorticity within combined wave-current field in demonstrating complex nonlinear physical phenomena due to interaction between wave and current. In addition, NWCT is applied to simulate a more complex wave-current-structure field for wave propagating over a submerged breakwater associated with current. Detailed discussion including characteristics of velocity and vorticity fields and the relation between free surface and vorticity are given.

• 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.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.6 s.237
• /
• pp.671-686
• /
• 2005

The transient nature and complex geometries of two-phase gas-liquid flows cause fundamental difficulties when measuring flow velocity using an electromagnetic flowmeter. Recently, a current-sensing flowmeter was introduced to obtain measurements with high temporal resolution (Ahn et al.). In this study, current-sensing flowmeter theory was applied to measure the fast velocity transients in slug flows. The velocity fields of axisymmetric gas-liquid slug flow in a vertical pipe were obtained using Volume-of-Fluid (VOF) method, and the virtual potential distributions for the electrodes of finite size were also computed using the finite volume method for simulating slug flow. The output signal prediction for slug flow was carried out from the velocity and virtual potential (or weight function) fields. The flowmeter was numerically calibrated to obtain the cross-sectional liquid mean velocity at an electrode plane from the predicted output signal. Two calibration parameters are proposed for this procedure: a flow pattern coefficient and a localization parameter. The flow pattern coefficient was defined by the ratio of the liquid resistance between the electrodes for two-phase flow with respect to that for single-phase flow, and the localization parameter was introduced to avoid errors in the flowmeter readings caused by liquid acceleration or deceleration around the electrodes. These parameters were also calculated from the computed velocity and virtual potential fields. The results can be used to obtain the liquid mean velocity from the slug flow signal measured by a current-sensing flowmeter.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1951-1956
• /
• 2004

The transient nature and complex flow geometries of two-phase gas-liquid flows cause fundamental difficulties when measuring flow velocity using an electromagnetic flowmeter. Recently, a current-sensing flowmeter was introduced to obtain measurements with high temporal resolution (Ahn et $al.^{(1)}$). In this study, current-sensing flowmeter theory was applied to measure the fast velocity transients in slug flows. To do this, the velocity fields of axisymmetric gas-liquid slug flow in a vertical pipe were obtained using Volume-of-Fluid (VOF) method and the virtual potential distributions for the electrodes of finite size were also computed using the finite volume method for the simulated slug flow. The output signal prediction for slug flow was carried out from the velocity and virtual potential (or weight function) fields. The flowmeter was numerically calibrated to obtain the cross-sectional liquid mean velocity at an electrode plane from the predicted output signal. Two calibration parameters are required for this procedure: a flow pattern coefficient and a localization parameter. The flow pattern coefficient was defined by the ratio of the liquid resistance between the electrodes for two-phase flow with respect to that for single-phase flow, and the localization parameter was introduced to avoid errors in the flowmeter readings caused by liquid acceleration or deceleration around the electrodes. These parameters were also calculated from the computed velocity and virtual potential fields. The results can be used to obtain the liquid mean velocity from the slug flow signal measured by a current-sensing flowmeter.