• 대한기계학회논문집B
• /
• 제24권12호
• /
• pp.1570-1579
• /
• 2000

Finite element analysis of fluid flow with moving free surface has been carried out in two and tree dimensions. The new VOF-based numerical algorithm that has been proposed by the present authors was applied to several 2-D and 3-D free surface flow problems. The proposed free surface tracking scheme is based on two numerical tools that have been newly introduced by the present authots; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed numerical algorithm has been applied to 2-D and 3-D cavity filling and sloshing problems, which demonstrated versatility and effectiveness of the new free surface tracking scheme as well as the overall solution procedure. The proposed numerical algorithm resolved successfully the interacting free surface with each other. The simulated results demonstrated the applicability of proposed numerical algorithm to the practical problems of large free surface motion. Also, it has been demonstrated that the proposed free surface tracking scheme can be easily implemented in any irregular non-uniform grid systems and can be extended to the 3-D free surface flow problem without additional efforts.

• 대한기계학회논문집B
• /
• 제24권12호
• /
• pp.1555-1569
• /
• 2000

Numerical simulation of fluid flow with moving free surface has been carried out. For the free surface flow, a VOF(Volume of Fluid)-based algorithm utilizing a fixed grid system has been investigated. In order to reduce numerical smearing at the free surface represented on a fixed grid system, a new free surface tracking algorithm based on the donor-acceptor scheme has been presented. Novel features of the proposed algorithm are characterized as two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed algorithm can be easily implemented in any irregular non-uniform grid systems that are usual in finite element method (FEM). Moreover, the proposed algorithm can be extended and applied to the 3-D free surface flow problem without additional efforts. For computation of unsteady incompressible flow, a finite element approximation based on the explicit fractional step method has been adopted. In addition, the SUPG(streamline upwind/Petrov-Galerkin) method has been implemented to deal with convection dominated flows. Combination of the proposed free surface tracking scheme and explicit fractional step formulation resulted in an efficient solution algorithm. Validity of the present solution algorithm was demonstrated from its application to the broken dam and the solitary wave propagation problems.

• 대한조선학회논문집
• /
• 제54권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.

• 한국수자원학회논문집
• /
• 제44권6호
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• pp.65-67
• /
• 2011

The moment-of-fluid (MOF) method is a new volume-tracking method that accurately treats evolving material interfaces. The MOF method uses moment data, namely the material volume fraction, as well as the centroid, for a more accurate representation of the material configuration, interfaces and concomitant volume advection. In this paper, unstructured mesh extension of the MOF method is to be presented. The MOF method is coupled with a stabilized finite element incompressible Navier-Stokes solver for two materials. The effectiveness of the MOF method is demonstrated with a free-surface dam-break problem.

• Journal of Mechanical Science and Technology
• /
• 제19권7호
• /
• pp.1517-1525
• /
• 2005

Liquid sloshing displays the highly nonlinear free surface fluctuation when either the external excitation is of large amplitude or its frequency approaches natural sloshing frequencies. Naturally, the accurate tracking of time-varying free surface configuration becomes a key task for the reliable prediction of the sloshing time-history response. However, the numerical instability and dissipation may occur in the nonlinear sloshing analysis, particularly in the long-time beating simulation, when two simulation parameters, the relative time-increment parameter a and the fluid mesh pattern, are not elaborately chosen. This paper intends to examine the effects of these two parameters on the potential-based nonlinear finite element method introduced for the large amplitude sloshing flow.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.228-230
• /
• 2009

The moment-of-fluid (MOF) method is a new volume-tracking method that accurately treats evolving material interfaces. Based on the moment data (volume and centroid) for each material, the material interfaces are reconstructed with second-order spatial accuracy in a strictly conservative manner. The MOF method is coupled with a stabilized finite element incompressible Navier-Stokes solver for two fluids, namely water and air. The effectiveness of the MOF method is demonstrated with a free-surface dam-break problem.

• 제어로봇시스템학회논문지
• /
• 제20권8호
• /
• pp.789-794
• /
• 2014

This paper presents some real performance of two typical bilateral teleoperation benchmark tasks: peg-in-hole and surface tracking tasks. The tasks are performed by an energy-bounding algorithm in the master control and position-based impedance algorithm in the slave control. Performance is analyzed for the position-force tracking capabilities from free space motion to surface contacting motion. In addition, preliminary user performance is evaluated by measuring the completion time and maximum/average contact forces. The quality of the measured performance is also compared with that of other existing approaches.

• 한국해안·해양공학회논문집
• /
• 제32권1호
• /
• pp.1-10
• /
• 2020

자유수면을 포함하는 파동장과 같이 단상의 경계가 시간발전에 따라 지속적으로 변화하는 경우나 액상과 기상이 혼합되는 문제에 있어서는 다상유동(multiphase flow) 문제를 적용하는 예가 증가하고 있다. 특히, 파동장과 같은 자유수면의 문제를 취급하는데 있어서는 혼합되지 않는 액상과 기상의 비압축성 뉴턴유체를 고려한 혼상류 모델이 적용되는 경우가 많다. 일반적으로 혼상류 모델은 각상의 경계면에 대한 시간기반 거동추적이 필수적이며, 궁극적으로는 계산의 정도를 좌우한다. 본 연구는 다양한 CFD 수치해석코드에 적용되고 있는 대표적인 VOF-type의 경계면 추적기법들의 이류성능을 평가하였다. 특히, 기존의 전통적인 VOF-type의 경계면 추적기법 및 이류계산에서 발생하는 수치확산을 최소화하기 위해 수치유속(numerical flux)을 제어하는 FCT 법의 효용성을 평가하고, 더불어 CIP 법을 활용한 자유수면 추적성능의 가능성을 고찰하였다. 그 결과, 본 연구에서 적용한 제한된 조건하에서는 수치확산 방지를 위해 수치확산방지 유속을 도입한 FCT-VOF 법이 가장 높은 경계면의 추적성능을 보였다. 본 연구에서 도출되는 결과는 다양한 수치해석코드에 적용되는 자유수면의 추적기법을 선택함에 있어서 중요한 기초자료로 활용될 것으로 기대된다.

• 한국해양공학회지
• /
• 제4권1호
• /
• pp.55-61
• /
• 1990

For the tracking of three dimensional free surface motions, a method referred to as the Volume of Fluid(VOF) algorithm is extended. In order to calculate the slope of three dimensional free surface which is the most important for the advection algorithm that decides the amount of fluid from cell to cell and for the application of free surface boundary condition, a simple method utilizing two dimensional slope informations is introduced. The extended algroithm is tested by demonstrating the simulation of a propagating sinusoidal wave through the channel whose width changes abruptly.