• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.409-413
• /
• 2005

Even though the relative importance of length scale of flow system allow us to simplify three dimensional flow problem to one or two dimensional representation, many systems still require three dimensional analysis. The objective of this study is to develop an efficient and accurate finite element model for analyzing and predicting three dimensional flow features in natural rivers and to offend to model spreading of pollutants and transport of sediments in the future. Firstly, three dimensional Reynolds averaged Navier-Stokes equations with the hydrostatic pressure assumption in generalized curvilinear coordinates were combined with the kinematic free-surface condition. Secondly. to simulate realistic high Reynolds number flow, the model employed the Streamline Upwind/Petrov-Galerkin(SU/PG) scheme as a weighting function for the finite element method in conjunction with an appropriate turbulence model(Smagorinsky scheme for the horizontal plain and Mellor-Yamada scheme for the vertical direction). Several tests is performed for the purpose of validation and verification of the developed model. A simple rectangular channel, 5-shaped and U-shaped channel are used for tests and comparisons are made with RMA-10 model. Runs for each case is converged stably without a oscillation and calculated water-surface deformation, longitudinal and transversal velocities, and velocity vector fields are in good agreement with the results of RMA-10 model.

• 한국가시화정보학회:학술대회논문집
• /
• 2004.12a
• /
• pp.169-175
• /
• 2004

The characteristics of flow around a rotating propeller were investigated using PIV technique. For each of four different blade phases of $0^{\circ},\;18^{\circ},\;36^{\circ}\;and\;54^{\circ}$four hundred instantaneous velocity fields were ensemble averaged to investigate the spatial evolution of the flow around a propeller. The phase-averaged mean velocity fields show that the viscous wake formed by the boundary layers developed on the blade surfaces and the slipstream contraction in the near-wake region. The out-of-plane velocity component and strain rate had large values at the locations of the tip and trailing vortices. The boundary layer developed along the ship hull bottom surface of the ship stern provides a strong turbulent shear layer, affecting the vortex structure in the propeller near-wake. As the flow develops in the downstream direction, the trailing vortices formed behind the propeller hub move upward slightly due to the presence of the hull wake and free surface. The turbulence intensity has large values around the tip and trailing vortices. As the wake moves downstream, the strength of the vorticity diminishes and the turbulence intensity increases due to turbulent diffusion and active mixing between the tip vortices and adjacent wake flow.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.19 no.2
• /
• pp.13-18
• /
• 1982

The method of Rankine source distribution is emerging as a powerful yet simple alternative for the solution of complicated free surface problems. But it has been uncertain whether the radiation condition could be satisfied exactly by distributing the simple sources on the free surface only. In this paper, it is proved rigorously that the Rankine sources, whose intensities are varying sinusoidally along the axis satisfying the free surface boundary condition, generate the radiation waves both in the infinite and finite-depth flows. A formula is derived to give the upper and lower bounds of the errors in the induced velocity computation that will be introduced by truncating the extent of source distribution on the free surface. Since the truncation is inevitable in the numerical analysis, this formula may be used as a criterion to limit the position of the field points, where velocity computation is made, away from the truncation boundary. A typical analysis shows that the maximum error will be 3.4 percent of the exact induced velocity when the field point is on the free surface two wave lengths away from the truncation boundary.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.99-102
• /
• 2002

The present study describes a numerical analysis for simulation of the sloshing of flows with free-surface which contained in a rectangular tank The SOLA-VOF (Volume of fluid) method uses a fixed mesh for calculating the motion of flow and the free-surface. This Eulerian approach enables the VOF method to use only a small amount of computer memory for simulating sloshing problems with complicated free-surface contours. The VOF function, representing the volume fraction of a cell occupied by the fluid, is calculated for each cells, which gives the locating of the free-surface filling any some fraction of cells with fluid. Using SOLA-VOF method, the study describes visualization about simulation of the sloshing of flows and damping effect by baffle. Translation and pitching motion of the forms have been investigated The time-dependent changes of free-surface height are used for visualization subject to several conditions such as fluid height horizontal acceleration, sinusoidal motion, and viscosity. The free-surface heights were used for comparing wall-force, which is caused by sloshing of flows. Baffle was Installed to reduce the force on the wall by sloshing of flows. Damping effects was extensively expressed under the conditions such as baffle shape and position.

• 한국해양학회지
• /
• v.17 no.2
• /
• pp.69-73
• /
• 1982

Wind wave generation is generally accompanied by a strong wind drift currint near the surface. The current can not be considered irrotational. The classical eqration gorerning free surface wave is formulated on the basis of irrptationality. To deal with the situation proterly, an equation governing free surface fluctuation of small amplitude is derived when the mean flow can not be assumed irrotational. The equation encompass the [lassical one as a limiting case as expected.

• Journal of the Korean Society of Visualization
• /
• v.10 no.3
• /
• pp.21-24
• /
• 2012

Although the motion of the three-phase contact line on a solid substrate has been extensively studied thus far, the understanding of the dynamics of the contact line of liquid/liquid/gas phases is far from complete. Here we deposit a drop of isopropyl alcohol (IPA) on water and HFE-7100 whose free surfaces are exposed to air to observe the flow field around the contact line. By combining the shadowgraph and high-speed imaging techniques, we find that vortices are spontaneously generated at the contact line, which grow in size with time. The flow is attributed to the Marangoni stress that pulls a liquid of lower-surface tension toward a liquid surface having a higher surface tension. However, it is not still clear why the entrained lower-surface-tension liquid should whirl rapidly beneath the contact line. We also visualize the flow by the particle image velocimetry (PIV) to find out that the rotational velocity reaches the order of 1 mm/s near the free surface.

• Ocean Systems Engineering
• /
• v.6 no.3
• /
• pp.245-264
• /
• 2016

The iterative boundary element method (IBEM) developed originally before for cavitating two-dimensional (2-D) and three-dimensional (3-D) hydrofoils moving under free surface is modified and applied to the case of 2-D (two-dimensional) airfoils and 3-D (three-dimensional) wings over water. The calculation of the steady-state flow characteristics of an inviscid, incompressible fluid past 2-D airfoils and 3-D wings above free water surface is of practical importance for air-assisted marine vehicles such as some racing boats including catamarans with hydrofoils and WIG (Wing-In-Ground) effect crafts. In the present paper, the effects of free surface both on 2-D airfoils and 3-D wings moving steadily over free water surface are investigated in detail. The iterative numerical method (IBEM) based on the Green's theorem allows separating the airfoil or wing problems and the free surface problem. Both the 2-D airfoil surface (or 3-D wing surface) and the free surface are modeled with constant strength dipole and constant strength source panels. While the kinematic boundary condition is applied on the airfoil surface or on the wing surface, the linearized kinematic-dynamic combined condition is applied on the free surface. The source strengths on the free surface are expressed in terms of perturbation potential by applying the linearized free surface conditions. No radiation condition is enforced for downstream boundary in 2-D airfoil and 3-D wing cases and transverse boundaries in only 3-D wing case. The method is first applied to 2-D NACA0004 airfoil with angle of attack of four degrees to validate the method. The effects of height of 2-D airfoil from free surface and Froude number on lift and drag coefficients are investigated. The method is also applied to NACA0015 airfoil for another validation with experiments in case of ground effect. The lift coefficient with different clearance values are compared with those of experiments. The numerical method is then applied to NACA0012 airfoil with the angle of attack of five degrees and the effects of Froude number and clearance on the lift and drag coefficients are discussed. The method is lastly applied to a rectangular 3-D wing and the effects of Froude number on wing performance have been investigated. The numerical results for wing moving under free surface have also been compared with those of the same wing moving above free surface. It has been found that the free surface can affect the wing performance significantly.

• Journal of Advanced Marine Engineering and Technology
• /
• v.32 no.8
• /
• pp.1215-1220
• /
• 2008

To obtain the shape of the free surface more accurately, computations are carried out by a finite volume method using unstructured meshes and an interface capturing method. Free-surface flow, which is very important in the fields of ship and marine engineering, is numerically simulated for flows of both water and air. Control volumes are used with an arbitrary number of faces and allows a local mesh refinement. The integration is of second order, with a midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation. The solution method of pressure-correction type solves sequentially equations of momentum, continuity, conservation, and two-equations turbulence model. Comparison are quantitatively made between the computation and experiment in order to confirm the solution method.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.8 no.2
• /
• pp.153-168
• /
• 2016

Sloshing phenomenon is a complicated free surface flow problem that increases the dynamic pressure on the sidewalls and the bottom of the storage tanks. When the storage tanks are partially filled, it is essential to be able to evaluate the fluid dynamic loads on the tank's perimeter. In this paper, a numerical code was developed to determine the pressure distribution on the rectangular and trapezoidal storage tanks' perimeters due to liquid sloshing phenomenon. Assuming the fluid to be inviscid, the Laplace equation and the nonlinear free surface boundary conditions were solved using coupled boundary element - finite element method. The code performance for sloshing modeling was validated using Nakayama and Washizu's results. Finally, this code was used for partially filled rectangular and trapezoidal storage tanks and free surface displacement, pressure distribution and horizontal and vertical forces exerted on the tanks' perimeters due to liquid sloshing phenomenon were estimated and discussed.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.90-100
• /
• 2006

A digital wave tank (DWT) simulation technique has been developed by authors to investigate the interactions of fully nonlinear waves with 3D marine structures. A finite-difference/volume method and a modified marker-and-cell (MAC) algorithm have been used, which are based on the Navier-Stokes (NS) and continuity equations. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique or the Level-Set (LS) technique developed for one or two fluid layers. In this paper, some applications for various engineering problems with free-surface are introduced and discussed. It includes numerical simulation of marine environments by simulation equipments, fully nonlinear wave motions around offshore structures, nonlinear ship waves, ship motions in waves and marine flow simulation with free-surface. From the presented simulations, it seems that the developed DWT simulation technique can handle various engineering problems with free-surface and reliably predict hydrodynamic features due to the fully-nonlinear wave motions interacting with such marine structures.