• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.5
• /
• pp.628-634
• /
• 2003

We present numerical and experimental results of periodic flows inside a rectangular container under a background rotation. The periodic flows are generated by changing the speed of rotation periodically so that a time-periodic body forces produce the unsteady flows. In numerical computation, a parallel-computation technique with MPI is implemented. Flow visualization and PIV measurement are also performed to obtain velocity fields at the free surface. Through a series of numerical and experimental works, we aim to clarify, if any, the fundamental reasons \ulcornerf discrepancy between the two-dimensional computation and the experimental measurement, which was detected in the previous study for the same flow model. Specifically, we check if the various assumptions prerequisite for the validity of the classical Ekman pumping law are satisfied for periodic flows under a background rotation.

• Journal of computational fluids engineering
• /
• v.14 no.2
• /
• pp.68-76
• /
• 2009

The sloshing tank causes the instability of the fluid flows and the fluctuation of the impact pressure by the liquid on the tank. These flow characteristics inside the sloshing tank can generate the uncomfortable sloshing noise. In the present study, a numerical analysis for the reduction of a fuel tank sloshing noise was performed. To simulate the flow characteristics in a sloshing tank with partially filled liquid, a VOF method was used for interfacial flows by applying a momentum source term for the sloshing motion in a non-inertial reference frame. This numerical method was verified by comparing its results with the available experimental data. For the reduction of the sloshing noise, the horizontal and vertical baffles and porous media inside a sloshing tank were considered and numerically analyzed in the present study. For various installations of these baffles and porous media, the characteristics of the liquid behavior in the sloshing tank were obtained along with the impact pressure on the wall and the height of the free surface along the wall. These basic results can be used for the design of the actual vehicular fuel tank with the reduced sloshing noise.

• Journal of Ship and Ocean Technology
• /
• v.12 no.1
• /
• pp.16-27
• /
• 2008

The interaction between advancing ships and the waves generated by them plays important roles in wave resistances and ship motions. Wave breaking phenomena near the ship bow at different speeds are investigated both numerically and experimentally. Numerical simulations of free surface profiles near the fore bodies of ships are performed and visualized to grasp the general trend or the mechanism of wave breaking phenomena from moderate waves rather than concentrating on local chaotic irregularities as ship speeds increase. Navier-Stokes equations are differentiated based on the finite difference method. The Marker and Cell (MAC) Method and Marker-Density Method are employed, and they are compared for the description of free surface conditions associated with the governing equations. Extra effort has been directed toward the realization of extremely complex free surface conditions at wave breaking. For this purpose, the air-water interface is treated with marker density, which is used for two layer flows of fluids with different properties. Adaptation schemes and refinement of the numerical grid system are also used at local complex flows to improve the accuracy of the solutions. In addition to numerical simulations, various model tests are performed in a ship model towing tank. The results are compared with numerical calculations for verification and for realizing better, more efficient research performance. It is expected that the present research results regarding wave breaking and the geometry of the fore body of ship will facilitate better hull form design productivity at the preliminary ship design stage, especially in the case of small and fast ship design. Also, the obtained knowledge on the impact due to the interaction of breaking waves and an advancing hull surface is expected to be applicable to investigation of the ship bow slamming problem as a specific application.

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

• International Journal of Fluid Machinery and Systems
• /
• v.4 no.1
• /
• pp.33-46
• /
• 2011

Concave surface boundary-layer flows are subjected to centrifugal instability which results in the formation of streamwise counter-rotating vortices. Such boundary layer flows have been experimentally investigated on concave surfaces of 1 m and 2 m radius of curvature. In the experiments, to obtain uniform vortex wavelengths, thin perturbation wires placed upstream and perpendicular to the concave surface leading edge, were used to pre-set the wavelengths. Velocity contours were obtained from hot-wire anemometer velocity measurements. The most amplified vortex wavelengths can be pre-set by the spanwise spacing of the thin wires and the free-stream velocity. The velocity contours on the cross-sectional planes at several streamwise locations show the growth and breakdown of the vortices. Three different vortex growth regions can be identified. The occurrence of a secondary instability mode is also shown as mushroom-like structures as a consequence of the non-linear growth of the streamwise vortices. Wall shear stress measurements on concave surface of 1 m radius of curvature reveal that the spanwise-averaged wall shear stress increases well beyond the flat plate boundary layer values. By pre-setting much larger or much smaller vortex wavelength than the most amplified one, the splitting or merging of the streamwise vortices will respectively occur.

• Communications of the Korean Mathematical Society
• /
• v.16 no.1
• /
• pp.137-152
• /
• 2001

Steady two-dimensional flows due to an applied pressure distribution in water of finite depth are considered. Gravity is included in the dynamic boundary condition. Gravity is included in the dynamic boundary condition. The problem is solved numerically by using the boundary integral equation technique. It is shown that, for both supercritical and subcritical flows, solutions depend on three parameters: (i) the Froude number, (ii) the magnitude of applied pressure distribution, and (iii) the span length of pressure distribution. For supercritical flows, there exist up to two solutions corresponding to the same value of Froude number for positive pressures and a unique solution for negative pressures. For subcritical flows, there are solutions with waves behind the applied pressure distribution. As the Froude number decreases, these waves when the Froude numbers approach the critical values.

• Journal of computational fluids engineering
• /
• v.15 no.1
• /
• pp.71-80
• /
• 2010

The violent free-surface motions and the corresponding impact loads are numerically simulated by using the Moving Particle Semi-implicit (MPS) method, which was originally proposed by Koshizuka and Oka (1996) for incompressible flows. In the original MPS method, there were several shortcoming including non-optimal source term, gradient and collision models, and search of free-surface particles, which led to less-accurate fluid motions and non-physical pressure fluctuations. In the present study, how those defects can be remedied is illustrated by step-by-step improvements in respective processes of the revised MPS method. The improvement of each step is explained and numerically demonstrated. The numerical results are also compared with the experimental results of Martin and Moyce (1952) for dam-breaking problem. The current numerical results for violent free-surface motions and impact pressures are in good agreement with their experimental data.

• Journal of Ocean Engineering and Technology
• /
• v.14 no.1
• /
• pp.11-16
• /
• 2000

수중날개를 대상으로 임의로 변형하는 이층류의 수직계산을 수행하였다. 계산은 비구조격자와 접면포착법을 이용하여 자유표면의 형상을 결정하였다 본 연구에 의하면 낮은 물속 깊이와 높은 프루드수에서 쇄파현상이 발생하였다 수중날개에 대하여 쇄파형상의 효과를 연구하였고 실험값과 계산값을 상호비교함으로서 격자의 민감성을 수직적으로 확인하였다.

• Journal of Korea Water Resources Association
• /
• v.57 no.8
• /
• pp.493-507
• /
• 2024

RANS-based CFD analysis is widely applied in various engineering fields, including practical hydraulic engineering, due to its high computational efficiency. However, problems of non-physical behavior in the analysis of two phase flow, such as free surfaces, have long been raised. The two-equation turbulence models used in general RANS-based analysis were developed for single phase flow and simulate unrealistically high turbulence energy at the interface where there are abrupt changes in fluid density. To solve this issue, one of the methods recently developed is the buoyancy-modified turbulence model, which has been partially validated in coastal engineering, but has not been applied to open channel flows. In this study, the applicability of the buoyancy-modified turbulence model is evaluated using the VOF method in the open-source program OpenFoam. The results of the uniform flow showed that both the buoyancy-modified k-𝜖 model and the buoyancy-modified k-ω SST model effectively simulated the reduction of turbulence energy near the free surface. Specifically, the buoyancy-modified k-ω SST model accurately simulated the vertical velocity distribution. Additionally, the model is applied to dam-break flows to examine cases with significant surface variation and cavity formation. The simulation results show that the buoyancy-modified turbulence models produce varying results depending on the VOF method and shows non-physical behavior different from experimental results. While the buoyancy-modified turbulence model is applicable in cases with stable surface shapes, it still has limitations in general application when there are rapid changes in the free surface. It is concluded that appropriate adjustments to the turbulence model are necessary for flows with rapid surface changes or cavity formation.