• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.7
• /
• pp.723-733
• /
• 2011

A numerical method of the CSF(Continuum Surface Force) model is presented for the calculation of the surface tension force and implemented in an in-house solution code(PowerCFD). The present method(code) employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with volume capturing method(CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The application of the present method to a 2-D liquid drop problem is illustrated by an equilibrium and nonequilibrium oscillating drop calculation. It is found that the present method simulates efficiently and accurately surface tension-dominant multiphase flows.

• Journal of the Society of Naval Architects of Korea
• /
• v.53 no.5
• /
• pp.400-409
• /
• 2016

In recent, shipping companies have made an enormous effort to improve the operation of vessel in various approaches, due to recession of shipping market and increasing competition among shipping companies. One of important parameters for improving the efficiency of vessel is the resistance performance that consist of friction and residual resistance. Especially, it is recognized that the friction resistance tends to be affected by conditions of vessel’s surface and occupies approximately 70~90% of the total resistance for slow speed ships. In general, the surface of vessel is covered with various type of paint to reduce fouling and corrosion. As time goes by, however, it is so hull roughness would be increased by fouling over the wetted surface that anti-fouling paints, such as CDP(Controlled Depletion Paint), Tin-Free SPC(Self Polishing Co-polymer) or Foul Release, are applied evenly on the hull surface. Nevertheless, these anti-fouling paints could not prevent fouling absolutely. A fundamental study on evaluating ship resistance performance variation due to hull roughness has been performed using a commercial software, Star-CCM+, which solves the continuity and Navier-Stokes equations for incompressible and viscous flow. The results of present simulation for plate are compared with some experimental data available and the effect of surface roughness to ship resistance performance is discussed.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.3
• /
• pp.91-100
• /
• 1988

Assuming that tidal level is constantly changed with an amplitude of 10 meters and a cycle of 12 hours, the slope stability for a typical dyke is analysed. The variation of pore water pressure within the dyke during the tidal change is obtained using a computer program, FLUMP, which is incorporated with saturated-unsaturated and transient flow. The results show that the variation of free water surface and distribution of pore water pressure within the dyke during the tidal fluctuations can be clearly predicted with the computer program. When a tide is lowered to the minimum level, the predicted pressure head is higher than the level of the free water surface in some parts of the dyke; that is, excess pore water pressure is generated in a zone affected by the tidal change. Also an unsaturated zone which shows negative pore water pressure is temporally created when a tide is lowered. The shear strength of the zone can be predicted based on the proposal suggested by Fredlund et al. It is emphasized that the excess pore water pressure generated during tidal fluctuations and strength parameters for the unsaturated zone should be considered in analyzing the slope stability of dykes. When those are considered, the critical slip surface seems to be located below the free water surface obtained when a tide is at the lowest.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.5
• /
• pp.109-117
• /
• 2002

An off-surface flow visualization experiments have been performed to investigate the flow field over a delta wing with the leading edge extension(LEX). The model is a flat wing with $65^{\circ}$ sweepback angle. The free stream velocity is 6.2 m/s, which corresponds to Reynolds number of $4.4\times10^5$ based on the wing root chord. The angle of attack and sideslip angle range from $16^{\circ}\sim28^{\circ}$ and $0^{\circ}\sim-15^{\circ}$, respectively. The visualization technique of using the micro water-droplet and the laser beam sheet enabled to observe the vortical flow structures, which can not be obtained by 5-hole probe measurements. With sideslip angle, the interaction and breakdown of the LEX and wing vortices was promoted in the windward side, whereas, it was suppressed in the leeward side.

• Journal of Korea Water Resources Association
• /
• v.32 no.3
• /
• pp.237-246
• /
• 1999

A finite element model is studied to simulate unsteady free surface flow based on dynamic wave equation and collocation method. The collocation method is used in conjunction with Hermite polynomials, and resulting matrix equations are solved by skyline method. The model is verified by applying to hydraulic jump, nonlinear disturbance propagation and dam-break flow in a horizontal frictionless channel. The computed results are compared with those by Bubnov-Galerkin and Petrov-Galerkin methods. It is also applied to the North Han River to simulate the floodwave propagation. The computed results have good agreements with those of DWOPER model in terms of discharge hydrographs. The suggested model has proven to be one of the promising scheme for simulating the gradually and rapidly varied unsteady flow in open channels.

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

Simulation of flow past a complex marine structure requires a fine resolution in the vicinity of the structure, whereas a coarse resolution is enough far away from it. Therefore, a lot of grid cells may be wasted, when a simple Cartesian grid system is used for an Immersed Boundary Method (IBM). To alleviate this problems while maintaining the Cartesian frame work, we adopted an Adaptive Mesh Refinement (AMR) scheme where the grid system dynamically and locally refines as needed. In this study, We implemented a moving IBM and an AMR technique in our basic 3D incompressible Navier-Stokes solver. A Volume Of Fluid (VOF) method was used to effectively treat the free surface, and a recently developed Lagrangian Dynamic Subgrid-scale Model (LDSM) was incorporated in the code for accurate turbulence modeling. To capture vortex induced vibration accurately, the equation for the structure movement and the governing equations for fluid flow were solved at the same time implicitly. Also, We have developed an interface by using AutoLISP, which can properly distribute marker particles for IBM, compute the geometrical information of the object, and transfer it to the solver for the main simulation. To verify our numerical methodology, our results were compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. Using the verified code, we investigated the following cases. (1) simulating flow around a floating sphere. (2) simulating flow past a marine structure.

• Journal of Korea Water Resources Association
• /
• v.38 no.11
• /
• pp.917-925
• /
• 2005

LSIV (Large Scale Image Velcocimetry) Is one of the image-based velocity measurement techniques. Since it measures surface velocities, it gives simple and inexpensive way to measure velocity, compared to other methods. Because of these advantages, there have been many studies to apply LSIV to the river discharge measurement in the field. Measuring the discharge by using LSIV requires a method which converts a surface velocity to a mean velocity In the present study, experiments and analysis of vortical velocity profile of open-channel flow in various conditions were performed to develop methods which estimate a mean velocity from a surface velocity. The result of this experiment reveals that velocity-dip phenomena occur at free-surface layer in open channel flow and Froude number affects more than bed roughness does. Two methods for estimating the mean velocity were proposed. One is to correct the wake law's profiles by using the difference of surface velocity from the mean velocity, and the other is to use the ratio of mean and surface velocities. The result of applying these methods in an experiment shows that they are quite accurate having an error of approximately $6\%$ only.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.1
• /
• pp.81-92
• /
• 1992

This paper deals with the procedure for developing a computer program which can predict the pressure fluctuation on the ship hull by solving the boundary value problem on the hull subject to the influence of the unsteady propeller and cavity motions. The program is applied to the solution of flow around a sphere under the influence of point sources simulating the propeller cavity, and then is compared with the analytic solution based on Butler's sphere theorem. The effect of free surface condition, either pressure-free or rigid-wall, upon the pressure distribution is studied. The computer code is also applied to a RO-RO ship, leading to the conclusion that the package may be useful for the analysis of excitation forces on the ship hull induced by the propeller in the design process.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.8 no.3
• /
• pp.338-350
• /
• 1996

A numerical analysis is carried out to study the two-dimensional steady state natural convection from vertical wires immersed in cold pure water. The surface of the wire is $0^{\circ}C$ unifrom temperature. Results of the analysis are presented for free stream temperature from $0^{\circ}C$ to $25^{\circ}C$ and the aspect ratio N from $5.26{\times}10^{-3}$ to $1.0{\times}10^{-3}$. The effects of the density extremum and aspect ratio on the flow pattern and the heat transfer characteristics are discussed As the aspect ratio N becomes larger, in the range of $1.0^{\circ}C{\leq}T_{\infty}{\leq}4.4^{\circ}C$ and $6{^{\circ}C}{\leq}T_{\infty}{\leq}17^{\circ}C$, the effect of Pr number on the heat transfer is shown to be more significant than the aspect ratio. Investigating into the effect of the density extremum on the heat transfer from wires, the new heat transfer correlations are suggested with the relation of average Nu mumber vs. modified Ra number. Here, the coefficient values C of correlations are presented as the function of density extremum parameter $R^*$. The effects of the density extremum parameter are also discussed.

• Proceedings of the Korean Vacuum Society Conference
• /
• 1999.07a
• /
• pp.109-109
• /
• 1999

Recently, low energy ion beam irradiation has been adopted for surface modification. Low energy ion beam irradiation has many advantages in polymer engineering such as weak damage, good adhesion, noticeably-enhanced wettability(less than 15 degree), good reproducibility and so on. In this experiment, chemical reactions between free radicals and environment gas species have been investigated using angle-resolved XPS and TRIM code. In the case of low ion beam energy (around 1 keV), energy loss in polymer is mainly originated from atomic collisions instead of electronic interference. Atomic collisions could generated displaced atoms and free radicals. Cold cathode-ion source equipped with 5cm convex grid was used in an O2 environment. Base and working pressure were 5$\times$10-6 and 2.3$\times$10-4 Torr. Flow rates of argon and oxygen were fixed at 1.2 and 8 sccm. target materials are polyethylene polyvinyidenefluoride and polytetrafluoroethylene.