• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3B
• /
• pp.247-257
• /
• 2009

In the present paper, turbulent open-channel flows with alternate vegetated zones are numerically simulated using threedimensional model. The Reynolds-averaged Navier-Stokes Equations are solved with the ${\kappa}-{\varepsilon}$ model. The CFD code developed by Olsen(2004) is used for the present study. For model validation, the partly vegetated channel flows are simulated, and the computed depth-averaged mean velocity and Reynolds stress are compared with measured data in the literature. Comparisons reveal that the present model successfully predicts the mean flow and turbulent structures in vegetated open-channel. However, it is found that the ${\kappa}-{\varepsilon}$ model cannot accurately predict the momentum transfer at the interface between the vegetated zone and the non-vegetated zone. It is because the ${\kappa}-{\varepsilon}$ model is the isotropic turbulence model. Next, the open channel flows with alternate vegetated zones are simulated. The computed mean velocities are compared well with the previously reported measured data. Good agreement between the simulated results and the experimental data was found. Also, the turbulent flows are computed for different densities of vegetation. It is found that the vegetation curves the flow and the meandering flow pattern becomes more obvious with increasing vegetation density. When the vegetation density is 9.97%, the recirculation flows occur at the locations opposite to the vegetation zones. The impacts of vegetation on the flow velocity and the water surface elevation are also investigated.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.4
• /
• pp.206-216
• /
• 2017

The temporal distribution of the turbulence kinetic energy (TKE) and the vertical component of Reynolds stresses ($-{\bar{u^{\prime}w^{\prime}}}$) was measured during one wave period under high wave energy conditions. The wave data were obtained at Hujeong Beach in the east coast of Korea at January 14~18 of 2017 when an extratropical cyclone was developed in the East Sea. Among the whole thousands of waves measured during the period, hundreds of regular waves that had with similar pattern were selected for the analysis in order to give three representing mean wave patterns using the ensemble average technique. The turbulence properties were then estimated based on the selected wave data. It is interesting to find out that $-{\bar{u^{\prime}w^{\prime}}}$ has one clear peak near the time of flow reversal while TKE has two peaks at the corresponding times of maximum cross-shore velocity magnitudes. The distinguished pattern of Reynolds stress indicates that vertical fluxes of such properties as suspended sediments may be enhanced at the time when the horizontal flow direction is reversed to disturb the flows, supporting the turbulence convection process proposed by Nielsen (1992). The characteristic patterns of turbulence properties are examined using the CADMAS-SURF Reynolds-Averaged Navier-Stokes (RANS) model. Although the model can reasonably simulate the distribution of TKE pattern, it fails to produce the $-{\bar{u^{\prime}w^{\prime}}}$ peak at the time of flow reversal, which indicates that the application of RANS model is limited in the prediction of some turbulence properties such as Reynolds stresses.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.55-55
• /
• 2011

연안 및 해안공학의 발달과 더불어 부유식방파제의 기능적 효율성이 중요시 되고 있다. 흔히 사용되어오던 착저식방파제는 설치에 많은 시간과 경비가 소요되고 환경 및 생태계에 많은 변화를 줄 수 있으며, 설치 예정지의 수리학적 특성 등의 여건에 많은 제약을 받는 단점이 있다. 부유식방파제는 일본 등의 선진국을 중심으로 활용이 잦아지고 있는 방파제로서 수면 위에 설치되기 때문에 수중 생태계에 미치는 영향이 적은 친환경방파제이다. 또한 기존에 시공된 중력식방파제와는 달리 수심에 제한을 덜 받고, 공사기간이 짧기 때문에 경제적이다. 실제 시공사례로는 2007년 마산 원전항에 완공된 부유식방파제가 대표적이며, 지금까지도 부유식방파제에 대한 여러 연구자들의 관심이 증가하고 있는 추세이다. 방파제뿐만 아니라 우리나라처럼 국토의 면적이 작은 지역에서 증가하는 해상물동량을 소화하기 위해서 부유식방파제 등을 이용한 항만의 시공이 필요한 실정이다. 이러한 부유식방파제의 분석적인 측면에 있어서 수치해석은 파랑과 구조물의 상호작용을 해석하는 데 한계가 있으며, 부유식방파제 단면형상을 정확하게 재현할 수 없으므로, 수리모형실험을 통한 부유식방파제의 연구가 필요할 것으로 판단된다. 최근 기술의 발달로 인한 유동장 해명이 가능해 졌으며, PIV(Particle image velocimetry) 및 LDV시스템은 다양한 분야에서 응용되고 있다. 특히, LDV시스템은 측정하려는 한 지점에 대하여 레이저 빔을 단면(Cross-section)으로 만들고 입자의 산란광을 후방산란(Back scatter)으로 받아서 도플러 효과를 이용, 속도에 대한 주파수를 획득하며, 유속을 측정하는 장비로 매우 높은 정확도와 비접촉식 이라는 장점을 가지고 있다. 또한, PIV 시스템에 비하여 측정시간이 오래 걸리는 반면 데이터를 가공하지 않고 활용할 만큼 높은 정확성을 가지고 있다. 본 연구에서는 수리모형실험을 통하여 단독형, 2열형 및 3열형 부유식방파제의 형상, 흘수 및 거리를 변화시키며 유동장을 수집하였으며, 방파성능에 따른 와의 생성 및 소멸시점에서의 파랑변형과의 관계를 분석하였다. 방파제의 형상과 흘수를 달리하여 수리모형실험을 수행하였으며, 와류의 상관관계를 분석하였다. 또한, 연직 2차원 Navier-Stokes 방정식 모형을 이용하여 수치모형실험을 수행하였으며, 수치모형실험 결과와 수리모형실험 결과를 비교 분석하였다. 후방방파제에서 발생되는 파랑은 입사파의 주기가 길어질수록 상대적으로 커지는 현상을 보였으며, 흘수심이 깊어질수록 전방방파제 입사 면에서 자유 수면이 높게 관측되는 결과를 보였다. 또한, 비교적 장주기파랑에 해당하는 입사파랑의 경우 전달파고비 산정에 있어서 설계기준인 0.5를 대다수 초과하는 반면, 3열형 구조에서는 대부분이 0.5이하로 상당히 높은 방파성능 결과를 나타내었다.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.9 no.4
• /
• pp.878-885
• /
• 2008

The FPDs(Flat Panel Displays) such as LCD(Liquid Crystal Display) and PDP(Plasma Display Panel) and OLED(Organic Light Emitting Diode), recently, have been substituted for CRT(Cathode Ray Tube) displays because they have a convex surface, small volume, light weight and lower electric power consumption. The productivity of FPDs is greatly dependent on the area of thin glass panel with 0.6 - 0.8mm thickness because FPDs are manufactured by cutting a large-scaled thin glass panel with patterns to the required product dimensions. So FPD's industries are trying to increase the area of thin glass panel. For example, the thin glass panel size of the 8th generation is 2,200mm in width, 2,600mm in length and 0.7mm in thickness. The air flows both in the thin glass panel and in the porous PE-plate surface were modeled and analyzed, from which a working condition was estimated. The thin glass panel on the porous PE-plate surface with self-lubricating characteristics was investigated and compared with that on the square duct floating bar surface with many holes of 1mm diameter when the thin glass panel contacts the floating bar surface due to malfunction of electric power supply.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.4 no.4
• /
• pp.21-31
• /
• 2001

When an oil-spilling accident occurs at sea, it is of the primary importance to predict the amount of oil leakage for the swift response and decision-making. The simplest method of oil-leakage estimation is based on the hydrostatic pressure balance between oil inside the tank and seawater outside of leakage hole, that is the so-called Torricelli equilibrium relation. However, there exists discrepancy between the reality and the Torricelli relation, since the latter is obtained from the quasi-steady treatment of Bernoulli equation ignoring viscous friction. A preliminary experiment has been performed to find out the oil-leaking speed and shape. Soy-bean oil inside the inner tank was ejected into water of the outer tank through four different leakage holes to record the amount of oil leakage. Furthermore, a CFD (Computational Fluid Dynamics) method was utilized to simulate the experimental situation. The Wavier-Stokes equations were solved for two-density flow of oil and water. VOF method was employed to capture the shape of their interface. It is found that the oil-leaking speed varies due to the frictional resistance of the leakage hole passage dependent on its aspect ratio. The Torricelli factor relating the speed predicted by using the hydrostatic balance and the real leakage speed is assessed. For the present experimental setup, Torricelli factors were in the range of 35%~55% depending on the aspect ratio of leakage holes. On the other hand, CFD results predicted that Torricelli factor could be 52% regardless of the aspect ratio of the leakage holes, when the frictional resistance of leakage hole passage was neglected.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.5
• /
• pp.345-357
• /
• 2011

Recently, as part of diversifying energy sources and earth environmental issues, technology development of new renewable energy using wave energy is actively promoted and commercialized around Europe and Japan etc. In particular, OWC(Oscillating Water Column) wave power generation system using air flow induced by vertical movement of the water surface by waves in an air-chamber within caisson is known as the most efficient wave energy absorption device and therefore, is one of the wave power generation apparatus the closest to commercialization. This study examines air flow velocity, which operates turbine(Wells turbine) directly in oscillating water column type wave power generation structure from two-and three-dimensional numerical experiments and discusses optimal shape of oscillating water column type wave power generation structure by estimating the maximum flow rate of air according to change in shape. The three-dimensional numerical wave flume was applied in interpretation for this study which is the model for the immiscible two-phase flow based on the Navier-Stokes Equation. From this, it turned out that size of optimal shape appears differently according to the incident wave period and air flow is maximized at the period where minimum reflection ratio occurs.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.28 no.6
• /
• pp.332-349
• /
• 2016

OLAFOAM is the powerful CFD code and is an expanded version of $OpenFOAM^{(R)}$, for wave mechanics simulation. The $OpenFOAM^{(R)}$ does provide many solvers to correspond to each object of the numerical calculation in a variety of fields. OLAFOAM's governing equation bases on VARANS (Volume-Averaged Reynolds-Averaged Navier-Stokes) equation, and the finite volume method is applied to numerical techniques. The program is coded in C++ and run on the Linux operating system. First of all, in this study, OLAFOAM was validated for 1) wave transformation inside porous structure under bore and regular wave conditions, 2) wave transformation by submerged breakwater under regular wave condition, and 3) regular wave transformation and resultant vertical velocity distribution under current by comparison with existing laboratory measurements. Hereafter, this study, which is almost no examination carried out until now, analyzed closely variation characteristics of water surface level, wave height, frequency spectrum, breaking waves, averaged velocity and turbulent kinetic energy around porous submerged breakwater in the wave and current coexisting field for the case of permeable or impermeable rear beach. It was revealed that the wave height fluctuation according to current direction(following or opposing) was closely related to the turbulent kinetic energy, and others.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.1
• /
• pp.11-23
• /
• 1997

This paper contains a new approach to blade section design method for marine propellers. The hydrodynamic characteristics of 2-D section are highly influenced by its geometrical parameters i.e., thickness and camber distributions and leading edge radius etc. To consider fully turbulent flow field near 2-D section. the finite volume method with k-${\varepsilon}$ turbulent model which solve Reynolds time averaged Navier-Stokes(RANS) equation is applied. In this study, O-type grid system that can provide many calculation points on blade surface is used. The results were compared with those of the experiment of NACA0012 to confirm the accuracy of the developed codes. The goal of this study is the development of a blade section with high efficiency and low drag. To achieve this, we carried out the tests of lift, drag and cavitation characteristics in cavitation tunnel. The results of experiment were compared with numerical results in order to validate the proposed blades design method. By comparing the numerical results with the experiments, we found that the new blade section, KH28 allows superior performance in efficiency and cavitation avoidance characteristics. We further investigated the blade section design method and an application study of this section, KH28 to apply to the marine propeller. In order to improve the accuracy of numerical results on prediction of lift and drag, we conclude here that the 2-layer boundary model must be used.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.27 no.6
• /
• pp.699-710
• /
• 2009

If a country has her geoid model, it could be determine accurate orthometric height because the geoid model could provide continuous equi-gravity potential surface. And it is possible to improve the coordinates accuracy of national control points through geodetic network adjustment considering geoidal heights. This study aims to find the best way to develop geoid model in Philippines which have similar topographic conditions as like Malaysia and Indonesia in Eastsouth asia. So, in this study, it is surveyed the general theories of geoid determination and development cases of geoid model in Asia and it is computed that the geoidal heights and gravity anomalies by spherical harmonic analysis using EGM2008, the latest earth geopotential model. The results show that first, the development of gravimetric geoid model based on airborne gravimetry is needed and second, about 200 GPS surveying data at national benchmark is needed. It is concluded that it is the most reasonable way to develop the hybrid geoid model through fitting geometric geoid by GPS/leveling data to gravimetric geoid. Also, it is proposed that four band spherical Fast fourier transformation(FFT) method for evaluation of Stokes integration and remove and restore technique using EGM2008 and SRTM for calculation of gravimetric geoid model and least square collocation algorithm for calculation of hybrid geoid model.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.26 no.7
• /
• pp.931-941
• /
• 2020

There is a growing interest this paper for ocean sensing where autonomous vehicles can play an essential role in assisting engineers, researchers, and scientists with environmental monitoring and collecting oceanographic data. This study was conducted to develop a rigid sail for the autonomous sailing drone. Our study aims to numerically analyze the aerodynamic characteristics of curvy twin sail and compare it with wing sail. Because racing regulations limit the sail shape, only the two-dimensional geometry (2D) was open for an optimization. Therefore, the first objective was to identify the aerodynamic performance of such curvy twin sails. The secondary objective was to estimate the effect of the sail's spacing and shapes. A viscous Navier-Stokes flow solver was used for the numerical aerodynamic analysis. The 2D aerodynamic investigation is a preliminary evaluation. The results indicated that the curvy twin sail designs have improved lift, drag, and driving force coefficient compared to the wing sails. The spacing between the port and starboard sails of curvy twin sail was an important parameter. The spacing is 0.035 L, 0.07 L, and 0.14 L shows the lift coefficient reduction because of dramatically stall effect, while flow separation is improved with spacing is 0.21 L, 0.28 L, and 0.35 L. Significantly, the spacing 0.28 L shows the maximum high pressure at the lower area and the small low pressure area at leading edges. Therefore, the highest lift was generated.