• Journal of Advanced Marine Engineering and Technology
• /
• v.32 no.8
• /
• pp.1215-1220
• /
• 2008

To obtain the shape of the free surface more accurately, computations are carried out by a finite volume method using unstructured meshes and an interface capturing method. Free-surface flow, which is very important in the fields of ship and marine engineering, is numerically simulated for flows of both water and air. Control volumes are used with an arbitrary number of faces and allows a local mesh refinement. The integration is of second order, with a midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation. The solution method of pressure-correction type solves sequentially equations of momentum, continuity, conservation, and two-equations turbulence model. Comparison are quantitatively made between the computation and experiment in order to confirm the solution method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.32 no.1
• /
• pp.23-29
• /
• 2008

Two high resolution compressive schemes, CICSAM(Ubbink, 1997) and HRIC(Muzaferija & Peric, 1999), in interface capturing method are reviewed briefly with respect to the extended forms suitable for unstructured meshes. And then those are applied to three typical test cases of translation test, shearing flow test and collapsing water problem with an obstacle. It is accomplished by implementing the high resolution schemes in the in-house CFD code(PowerCFD) for computing 3-D flow with an unstructured cell-centered method, which is based on the finite-volume technique and fully conservative. The calculated results show that CICSAM is better than HRIC with respect to accuracy and robustness, although either scheme can be used as a good choice for free surface or two-phase flow simulation.

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

• Journal of Advanced Research in Ocean Engineering
• /
• v.3 no.2
• /
• pp.75-82
• /
• 2017

With the development of computational fluid dynamics (CFD), studies on shipbuilding and maritime issues including free-surface wave flow have been conducted. Although the volume of fluid (VOF) and level-set methods are widely used to study the free-surface wave flow, disadvantages exist. In particular, it takes a long time to obtain solutions. In this study, a free-surface capturing code is developed for ship and offshore structures. The developed code focuses on accuracy and computation time. Open source CFD libraries, termed OpenFOAM, are used to develop the code. The results obtained using the developed code are compared with those obtained using interFoam. The results show that the developed code could be used to capture the free-surface wave flow without numerical diffusion; moreover, the accuracy of the developed code is largely the same as that of interFoam.

• Journal of Ship and Ocean Technology
• /
• v.6 no.2
• /
• pp.37-50
• /
• 2002

In the present numerical simulation of viscous free surface flow around a ship, two-fluids in-compressible Reynolds-averaged Navier-Stokes equations with the standard $\textsc{k}-\varepsilon$turbulence model are discretized on a regular grid by using a finite volume method. A local level-set method is introduced for capturing the free surface movement and the influence of the viscous layer and dynamic boundary condition of the free surface are implicitly considered. Partial differential equations in the level-set method are discretized with second order ENO scheme and explicit Euler scheme in the space and time integration, respectively. The computational results for the Series-60 model with $C_B=0.6$ show a good agreement with the experimental data, but more validation studies for commercial complicated hull forms are necessary.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2005.05a
• /
• pp.43-46
• /
• 2005

Micro molding technology is a promising mass production technology for polymer based microstructures. Mass production technologies such as the micro injection/compression molding, hot embossing, and micro reaction molding are already in use. In the present study, we have developed a numerical analysis system to simulate three-dimensional non-isothermal cavity filling for hot embossing, with a special emphasis on the free surface capturing. Precise free surface capturing has been successfully accomplished with the level set method, which is solved by means of the Runge-Kutta discontinuous Galerkin (RKDG) method. The RKDG method turns out to be excellent from the viewpoint of both numerical stability and accuracy of volume conservation. The Stokes equations are solved by the stabilized finite element method using the equal order tri-linear interpolation function. To prevent possible numerical oscillation in temperature Held we employ the streamline upwind Petrov-Galerkin (SUPG) method. With the developed code we investigated the detailed change of free surface shape in time during the mold filling. In the filling simulation of a simple rectangular cavity with repeating protruded parts, we find out that filling patterns are significantly influenced by the geometric characteristics such as the thickness of base plate and the aspect ratio and pitch of repeating microstructures. The numerical analysis system enables us to understand the basic flow and material deformation taking place during the cavity filling stage in microstructure fabrications.

• 한국전산유체공학회:학술대회논문집
• /
• 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 computational fluids engineering
• /
• v.13 no.1
• /
• pp.7-13
• /
• 2008

Several high resolution schemes such as OSHER, MUSCL, SMART, GAMMA, WACEB and CUBISTA are comparatively studied with respect to the accurate capturing of fluid interfaces throughout the application to two typical test cases of a translation test and a collapsing water column problem with a return wave. It is accomplished by implementing the high resolution schemes in the in-house CFD code(PowerCFD) for computing 3-D flow with an unstructured cell-centered method and an interface capturing method, which is based on the finite-volume technique and fully conservative. The calculated results show that SMART scheme gives the best performance with respect to accuracy and robustness.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.7
• /
• pp.723-733
• /
• 2011

A numerical method of the CSF(Continuum Surface Force) model is presented 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. The application of the present method to a 2-D liquid drop problem is illustrated by an equilibrium and nonequilibrium oscillating drop calculation. It is found that the present method simulates efficiently and accurately surface tension-dominant multiphase flows.