• International Journal of Ocean System Engineering
• /
• v.1 no.3
• /
• pp.136-143
• /
• 2011

In this paper, three-dimensional free-surface flows are simulated by using two different numerical methods, the constrained interpolation profile (CIP)-based and finite volume (FV)-based methods. In the CIP-based method, the governing equations are solved on stationary staggered Cartesian grids by a finite difference method, and an immersed boundary technique is applied to deal with wave-body interactions. In the FV-based method, the governing equations are solved by applying collocated finite volume discretization, and body-fitted meshes are used. A free-surface boundary is considered as the interface of the multi-phase flow with air and water, and a volumeof-fluid (VOF) approach is applied to trace the free surface. Among many variations of the VOF-type method, the tangent of hyperbola for interface capturing (THINC) and the compressive interface capturing scheme for arbitrary meshes (CICSAM) techniques are used in the CIP-based method and FV-based method, respectively. Numerical simulations have been carried out for dam-breaking and wave-body interaction problems. The computational results of the two methods are compared with experimental data and their differences are observed.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.10
• /
• pp.754-760
• /
• 2008

Computation of moving interface by the level set method typically requires the reinitialization of level set function. An inaccurate estimation of level set function $\phi$ results in incorrect free-surface capturing and thus errors such as mass gain/loss. Therefore, an accurate and robust reinitialization process is essential to the simulation of free-surface flows. In the present paper, we pursue further development of the reinitialization process, which evaluates level set function directly using a normal vector on the interface without solving there-distancing equation of hyperbolic type. The Taylor-Galerkin approximation and P1P1 splitting/SUPG (Streamline Upwind Petrov-Galerkin) FEM are adopted to discretize advection equation of the level set function and the incompressible Navier-Stokes equation, respectively. Advection equation and re-initialization process of free surface capturing are validated with benchmark problems, i.e., a broken dam flow and timereversed single vortex flow. The simulation results are in good agreement with the existing results.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6B
• /
• pp.691-696
• /
• 2008

In this paper, a three-dimensional non-hydrostatic pressure model for free surface flows using a normalized vertical coordinate system is presented. To strongly couple the free surface and non-hydrostatic pressure in the momentum equations, a double predictor-corrector method is employed. This research is especially focused on implementing the dynamic boundary condition (a zero pressure condition) at the free surface. This boundary condition can be specified accurately with a small modification to existing models. Numerical results with and without this modification clearly show that a precise implementation of the dynamic boundary condition is paramountly important.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.369-372
• /
• 2002

Under the assumption of potential flow, free-surface flows around a 2-dimensional hydrofoil are calculated by high-order spoctral/boundary-integral method. This method is one of the most efficient numerical methods by which the nonlinear interactions between hydrofoil and free-surface can be simulated in time-domain. Comparisons of the calculated free-surface profiles with other experimental results show relatively good agreements. As another example, free-surface flow generated by the heaving and translating hydrofoil is calculated and discussed.

• Journal of Ship and Ocean Technology
• /
• v.5 no.4
• /
• pp.29-43
• /
• 2001

In this paper, free surface flows around an advancing two-dimensional rectangular cylinder piercing the free surface are studied using numerical and experimental methods. Especially, wave breaking phenomenon around the cylinder is treated in detail. A series of numerical simulations and experiments were performed for the purpose of comparison. For the numerical simulations, a finite difference method was adopted with a rectangular grid system, and the variation of the free surface was computed by the marker density method. The computational results are compared with the experiments. It is confirmed that the present numerical method is useful for the numerical simulation of nonlinear free surface waves around a piercing body.

• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.1
• /
• pp.26-33
• /
• 2017

In this study, the surface marker method, one of the particle tracking methods, used to track the free surface is extended to cover the more general cases easily including the collision and separation of the free surface. In surface particle method to redistribute particles effectively using the grid system, the free surface is composed of the sum of quadrilaterals having four curves where fixed markers are placed at ends of each curve. Fixed markers are used to know how curves are connected to each other. The position of fixed markers can move as the free surface deforms but all fixed markers cannot be deleted during all time of simulation to keep informations of curve connection. In the case of the collision or separtion of the free surface where several curves can be intersected disorderly, severe difficulties can occur to define newly states of curve connection. In this study, the adaptable surface parTicle method without fixed markers is introduced. Intersection markers instead of the fixed markers are used to define quadrilaterals. The position of the intersection markers is defined to be the intersection point between the free surface and the edge of the grid and it can be added or deleted during the time of simulation to allow more flexibilities. To verify numerical schemes, two flow cases are simulated and the numerical results are compared with other's one and shown to be valid.

• Journal of the Korean Society of Visualization
• /
• v.2 no.2
• /
• pp.38-44
• /
• 2004

Theoretical and numerical studies are given to the sink flows within a rotating circular tank driven by the fluid withdrawal from a bottom circular hole. It was found that, when the upper free surface was set with no-slip boundary conditions, the Ekman boundary-layer develops not only above the bottom surface but under the top surface. The sink fluid is coming from the two Ekman layers, and the mass transfer from the bulk, inviscid region is dependent on the rotational speed. It is also remarkable to see that all the fluid gathered along the axis flows in a form of rapidly rotating fluid column haying almost the same diameter as the bottom hole.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.2-5
• /
• 2011

This paper presents a numerical study on multiphase flows induced by wall adhesion The CSF(Continuum Surface Force} model is used for the calculation of the surface tension force and implemented in an in-house solution code(PowerCFD). The present method(code) employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with volume capturing method(CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing As an application of the present method, the effects of wall adhesion are numerically simulated with the CSF model for a shallow pool of water located at the bottom of a cylindrical tank. Two different cases are computed, one in which the water wets the wall and one in which the water does not wet the wall. It is found that the present method simulates efficiently and accurately surface tension-dominant multiphase flows induced by wall adhesion.

• Journal of Korean Port Research
• /
• v.13 no.1
• /
• pp.175-182
• /
• 1999

The ocean wave is hydrodynamically investigated to get more reliable solution. To improve the computational accuracy more fine grids are used with relatively less computer storage on the free surface. One element of the free surface is discretized into more fine grids because the free-surface waves are much affected by the grid size in the finite difference scheme. Here the multi-grid method is applied to confirm the efficiency for the S103 ship model by solving the Navier-Stokes equation for the turbulent flows. According to the computational result approximately 30% can be improved in the free surface generation, Finally the limiting streamlines show numerical result is similar to the experiment by twin tuft.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.11-22
• /
• 1999

A numerical prediction method has been proposed to predict non-linear free surface oscillation in a three-dimensional container. The fluid motions are numerically predicted with Navier-Stokes equations discretized in a Lagrangian scheme with sufficient numerical accuracy. The profile of a free surface is precisely represented with three-dimensional body-fitted coordinates (BFC), which are regenerated in each computational step on the basis of the arbitrary Lagrangian-Eulerian (ALE) formulation. In order to confirm the reliability of the computational method, it was firstly applied to three-dimensional flows within complicated-shaped rigid boundaries, such as curved pipes and ducts. Than it was applied to benchmark computations related to free surface oscillations. Following these basic verifications, non-linear sloshings in a cylindrical tank and transitions from sloshing to swirling motions were numerically predicted. Throughout these computations, the applicability of the present computational method has been confirmed and some of the predicted free surface motions were visualized as sequential images and animations to understand their dynamic futures.