• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.4
• /
• pp.345-352
• /
• 2008

The motion of a projectile initiated by the release of highly pressurized air is simulated presuming the flow field as a two dimensional one. The effects of water compressibility on projectile movements are investigated, comparing results based on the Fluent VOF model where water is treated as an incompressible medium with those from the presently developed VOF scheme. The present model considers compressibility of both air and water. The Fluent results show that the body moves farther and at higher speeds than the present ones. As time proceeds, the relative difference of speed and displacement between the two results drops substantially, after acoustic waves in water traverse and return the full length of the tube several times. To estimate instantaneous accelerations, however, requires implementation of the water compressibility effect as discrepancies between them do not decrease even after several pressure wave cycles.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.12 no.8
• /
• pp.3374-3383
• /
• 2011

In this study, numerical analyses for slamming impact phenomena have been carried out using a 2-dimensional wedge shaped structure having finite deadrise angles. Fluid is assumed incompressible and entry speed of the structure is kept constant. Geo-reconstruct(or PLIC-VOF) scheme is used for the tracking of the deforming free surface. Numerical analyses are carried out for the deadrise angles of $10^{\circ}$, $20^{\circ}$ and $30^{\circ}$. For each deadrise angle, variations are made for the grid size on the wedge bottom and for the entry speed. The magnitude and the location of impact pressure and the total drag force, which is the summation of pressure distributed at the bottom of the structure, are analyzed. Results of the analyses are compared with the results of the Dobrovol'skaya similarity solutions, the asymptotic solution based on the Wagner method and the solution of Boundary Element Method(BEM).

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.575-589
• /
• 2008

In the present work, we investigate the hydrodynamic behavior of a turbulent bore, such as tsunami bore and tidal bore, generated by the removal of a gate with water impounded on one side. The bore generation system is similar to that used in a general dam-break problem. In order to the numerical simulation of the formation and propagation of a bore, we consider the incompressible flows of two immiscible fluids, liquid and gas, governed by the Navier-Stokes equations. The interface tracking between two fluids is achieved by the volume-of-fluid (VOF) technique and the M-type cubic interpolated propagation (MCIP) scheme is used to solve the Navier-Stokes equations. The MCIP method is a low diffusive and stable scheme and is generally extended the original one-dimensional CIP to higher dimensions, using a fractional step technique. Further, large eddy simulation (LES) closure scheme, a cost-effective approach to turbulence simulation, is used to predict the evolution of quantities associated with turbulence. In order to verify the applicability of the developed numerical model to the bore simulation, laboratory experiments are performed in a wave tank. Comparisons are made between the numerical results by the present model and the experimental data and good agreement is achieved.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.219-235
• /
• 2014

With the development of the technology of underwater moving bodies, the need for developing the knowledge of surface effect interaction of free surface and underwater moving bodies is increased. Hence, the two-phase flow is a subject which is interesting for many researchers all around the world. In this paper, the non-linear free surface deformations which occur during the water-exit of a circular cylinder due to its buoyancy are solved using finite volume discretization based code, and using Volume of Fluid (VOF) scheme for solving two phase flow. Dynamic mesh model is used to simulate dynamic motion of the cylinder. In addition, the effect of cylinder mass in presence of an external force is studied. Moreover, the oblique exit and entry of a circular cylinder with two exit angles is simulated. At last, water-exit of a circular cylinder in six degrees of freedom is simulated in 3D using parallel processing. The simulation errors of present work (using VOF method) for maximum velocity and height of a circular cylinder are less than the corresponding errors of level set method reported by previous researchers. Oblique exit shows interesting results; formation of waves caused by exit of the cylinder, wave motion in horizontal direction and the air trapped between the waves are observable. In 3D simulation the visualization of water motion on the top surface of the cylinder and the free surface breaking on the front and back faces of the 3D cylinder at the exit phase are observed which cannot be seen in 2D simulation. Comparing the results, 3D simulation shows better agreement with experimental data, specially in the maximum height position of the cylinder.

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

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.21 no.2
• /
• pp.182-190
• /
• 2009

For RANS equation model using VOF scheme Lin and Liu (1999) developed internal wave-maker method to generate target wave trains by using designed mass source functions of the continuity equation. By using this method studies on various numerical wave experiments has been achieved without the problem caused by wave reflection due to an external wave-maker. In this study, the wave-maker method to generate target wave trains by using a momentum source function was proposed. The computational results obtained by applying the mass and momentum source functions into FLUENT were compared with each other. To see its applicability, the hydraulic experiment of Luth et al. (1994) were numerically simulated and their measurements are compared with the computations, and the vertical variations of computed results were shown and investigated.

• Journal of Ocean Engineering and Technology
• /
• v.26 no.6
• /
• pp.7-13
• /
• 2012

This paper presents a numerical study of the wave loads acting on offshore structures using a Cartesian-grid-based flow simulation method. Finite volume discretization with a volume-of-fluid (VOF) method is adopted to solve two-phase Navier-Stokes equations. Among the many variations of the VOF method, the CICSAM scheme is applied. The body boundary conditions are satisfied using a porosity function, and wave generation is carried out by using transient (wave or damping) zone approaches. In order to validate the present numerical method, three different basic offshore structures, including a sphere, Pinkster barge, and Wigley model, are numerically investigated. First, diffraction and radiation problems are solved using the present numerical method. The wave exciting and drift forces from the diffraction problems are compared with potential-based solutions. The added mass and wave damping forces from the radiation problems are also compared with the potential results. Next, the wave-induced motion responses of the structures are calculated and compared with the existing experimental data. The comparison results are fairly good, showing the validity of the present numerical method.

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

In a numerical simulation of open channel turbulent flows using RANS (Reynolds averaged Navier-Stokes) equations model equipped with VOF (Volume of Fluid) scheme, the determination of wall roughness for wall function was studied. The roughness constant, based on the law-of-the-wall for flow on rough walls, obtained by experimental works for pipe flows is employed in general wall functions. However, this constant of wall function is the function of Froude number in open channel flows. Thus, the wall roughness should be determined by taking into account the effect of Froude number. In addition, the wall roughness should be corresponding to Manning's roughness coefficient widely used for open channels. In this study, the relation between wall roughness height as an input condition and Manning's roughness coefficient was investigated, and an equation for effective wall roughness height considering the characteristics of numerical models was proposed as a function of Manning's roughness coefficient.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.1-10
• /
• 1998

서로 섞이지 않는 두 비압축성 유체의 유동을 해석하기 위하여 VOF 방법에 기초한 수치 기법을 개발하였다. 유체간의 계면형상의 거동은 유동장내의 유체의 점유체적비의 변화에 의해 묘사되는데 이를 지배하는 이동방정식을 풀기 위한 새로운 대류항 계산법을 고안하였다. 대류항은 유체계면의 방향에 따라 풍상법과 역풍상법의 적절한 조합을 취하여 계산하는데 여기에 대각방향의 상류효과를 포함시켜 시간에 대한 2차 정확도를 갖도록 하였다. 또한 이 방법을 유량보정수송(FCT)법과 결합시켜 해의 단조성을 보장하였다. 몇 가지 단순 문제에 대한 시험 결과 이 기법이 수치오차에 의한 계면형상의 변형과 파손을 감소시킴을 확인하였다.

• Proceedings of the KSME Conference
• /
• 2007.05a
• /
• pp.532-537
• /
• 2007

Nano imprint lithography(NIL) is a cost-efficient, high-throughput processing technique to transfer nano-scale patterns onto thin polymer films. Polymers used as the resist include UV cured resins as well as thermoplastics such as polymethyl-methacrylate(PMMA). In this study, an analytic investigation was performed for the NIL process of transferring nano scale patterns onto polymeric films. Process optimization calls for a thorough understanding of resist flow during the process. We carried out 2D and 3D numerical analyses of resist flow during NIL process. The simulation incorporated continuum-hypothesis and the effects of surface tension were taken into account. For a more effective prediction of free surface, fixed grid scheme with the volume of fluid (VOF) method were used. The simulation results were verified with experimental results qualitatively. And the parametric study was performed for various process conditions.