• 대한기계학회논문집B
• /
• 제27권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.

• 한국전산유체공학회지
• /
• 제14권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
• /
• 제12권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.

• 대한기계학회논문집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.

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

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

• 한국전산유체공학회지
• /
• 제15권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.

• 한국해양공학회지
• /
• 제14권1호
• /
• pp.11-16
• /
• 2000

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

• 한국수자원학회논문집
• /
• 제46권2호
• /
• pp.139-153
• /
• 2013

고르지 않은 바닥을 지나는 천수 흐름을 해석하기 위해 천수방정식의 흐름률 경사항과 바닥 경사 생성항에 대해 HLLL 기법과 DFB(Divergence Form for Bed slope source term) 기법을 각각 적용하여 유한체적 모형을 구성하였다. 또한, PSC(Partially Submerged Cell)의 고려를 위해 VFR(Volume/Free-surface Relationship)도 이용하였다. MUSCL에서 WSDGM(Weighted Surface-Depth Gradient Method)을 보다 단순하게 고쳐도 원래의 방법과 정확도가 동등함을 1차원 정상 흐름에 대해 확인하였다. 1차원 PSC에 대한 VFR를 통해 흐름률 경사항과 바닥 경사 생성항의 선평형성이 정확하게 충족됨을 입증하였다. 2차원 PSC에서 DFB 기법으로는 지배방정식의 선평형성이 충족되지 않은 문제를 삼각형 격자에 대한 VFR를 이용하여 해소하였다. 삼각형 턱과 둥근 융기를 지나는 2차원 댐 붕괴 흐름에 대한 모의에서 실험실 실험 결과와 잘 부합됨을 확인하였다. 또한, 부분 댐 붕괴 흐름에 대한 모형의 적용에서 경사면은 물론 불규칙 바닥에서도 요철의 잠김이 성공적으로 모의되었다. 따라서 고르지 않은 실제 하천 지형에 대한 이 모형의 적용성이 기대된다.