• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.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.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.435-449
• /
• 2014

The effects of the gravity field and the free surface on the cavity shape and the drag are investigated through a numerical analysis for the steady supercavitating flow past a simple two-dimensional body underneath the free surface. The continuity and the RANS equations are numerically solved for an incompressible fluid using a $k-{\epsilon}$ turbulence model and a mixture fluid model has been applied for calculating the multiphase flow of air, water and vapor using the method of volume of fluid and the Schnerr-Sauer cavitation model. Numerical solutions have been obtained for the supercavitating flow about a two-dimensional $30^{\circ}$ wedge in wide range of depths of submergence and inflow velocities. The results are presented for the cavity shape, especially the length and the width, and the drag of the wedge in comparison with those of the case for the infinite fluid flow neglecting the gravity and the free surface. The influences of the gravity field and the free surface on the aforementioned quantities are discussed. The length and the width of the supercavity are reduced and the centerline of the cavity rises toward the free surface due to the effects of the gravity field and the free surface. The drag coefficient of the wedge, however, is about the same except for shallow depths of submergence. As the supercavitating wedge is approaching very close to the free surface, it is found the length and the width of a cavity are shorten even though the cavitation number is reduced. Also the present result suggests that, under the influence of the gravity field and the free surface, the length of the supercavity for a certain cavitation number varies and moreover is proportional to the inverse of the submergence depth Froude number.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.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.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2003.06a
• /
• pp.979-982
• /
• 2003

In this study, because of change in electromagnetic force, deformation of the free surface motion of a magnetic fluid is changed. Deformation of the free surface motion of a magnetic fluid for the change in electromagnetic force is discussed and carried out theoretically and experimentally on the basis of Rosensweig Ferrohydrodynamic Bernoulli Equation. Objective of this study explicates free surface motion by electromagnetic force and planes to designed controller. To control free surface of magnetic fluid, it embody designed two-dimensional free surface form of magnetic fluid. By using this characteristics, they applied to oscillator for surface control, flow control, boundary layer control. Strength of magnetic field and height of free surface of magnetic fluid measure as a hall-effect sensor. As performing height control of magnetic fluid, the result will be presented possibility of free surface deformation control.

• Journal of Korea Water Resources Association
• /
• v.44 no.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.

• 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 computational fluids engineering
• /
• v.19 no.4
• /
• pp.20-28
• /
• 2014

An efficient solution algorithm for simulating free surface problem is presented. Navier-Stokes equations for variable density incompressible flow are employed as the governing equation on Cartesian meshes. In order to describe the free surface motion efficiently, VOF(Volume Of Fluid) method utilizing THINC(Tangent of Hyperbola for Interface Capturing) scheme is employed. The most time-consuming part of the current free surface flow simulations is the solution step of the linear system, derived by the pressure Poisson equation. To solve a pressure Poisson equation efficiently, the PCG(Preconditioned Conjugate Gradient) method is utilized. This study showed that the proper application of the preconditioner is the key for the efficient solution of the free surface flow when its pressure Poisson equation is solved by the CG method. To demonstrate the efficiency of the current approach, we compared the convergence histories of different algorithms for solving the pressure Poisson equation.

• Nuclear Engineering and Technology
• /
• v.41 no.8
• /
• pp.1045-1052
• /
• 2009

A calculation model is developed to predict the transient free surface flow on the containment floor following a loss-of-coolant accident (LOCA) of pressurized water reactors (PWR) for the use of debris transport evaluation. The model solves the two-dimensional Shallow Water Equation (SWE) using a finite volume method (FVM) with unstructured triangular meshes. The numerical scheme is based on a fully explicit predictor-corrector method to achieve a fast-running capability and numerical accuracy. The Harten-Lax-van Leer (HLL) scheme is used to reserve a shock-capturing capability in determining the convective flux term at the cell interface where the dry-to-wet changing proceeds. An experiment simulating a sudden break of a water reservoir with L-shape open channel is calculated for validation of the present model. It is shown that the present model agrees well with the experiment data, thus it can be justified for the free surface flow with accuracy. From the calculation of flow field over the simplified containment floor of APR1400, the important phenomena of free surface flow including propagations and interactions of waves generated by local water level distribution and reflection with a solid wall are found and the transient flow rates entering the Holdup Volume Tank (HVT) are obtained within a practical computational resource.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.5 s.236
• /
• pp.609-616
• /
• 2005

While the SIMPLE algorithm is most widely used for the simulations of flow phenomena that take place in the industrial equipments or the manufacturing processes, it is less adopted for the simulations of the free surface flow. Though the SIMPLE algorithm is free from the limitation of time step, the free surface behavior imposes the restriction on the time step. As a result, the explicit schemes are faster than the implicit scheme in terms of computation time when the same time step is applied to, since the implicit scheme includes the numerical method to solve the simultaneous equations in its procedure. If the computation time of SIMPLE algorithm can be reduced when it is applied to the unsteady free surface flow problems, the calculation can be carried out in the more stable way and, in the design process, the process variables can be controlled based on the more accurate data base. In this study, a modified SIMPLE algorithm is presented fur the free surface flow. The broken water column problem is adopted for the validation of the modified algorithm (MoSIMPLE) and for comparison to the conventional SIMPLE algorithm.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.1
• /
• pp.104-114
• /
• 1993

In this paper, effects of free surface on viscous flow is investigated. Continuity equation coupled with Navier-Stokes equations are solved numerically by using an artificial compressibility method[1, 2]. The body-fitted generalized curvilinear coordinate system is employed to deal with arbitrary body shape. The IAF scheme with finite difference method is used to solve the equations, and a diagonal algorithm is applied to time-varying Jacobian matrices for the computational economics. Free surface shape is obtained by applying zero pressure condition to still water surface at each time step. A numerical test is made for larminar flow around a circular cylinder vertically piercing the free surface. Computed flow patterns are largely affected by the existance of free surface in low Reynolds number flows treated in this paper. Free surface causes viscous pressure drag to vary much in depth direction in accordance with the variations of flow pattern.