• Journal of Korea Water Resources Association
• /
• v.38 no.10 s.159
• /
• pp.871-883
• /
• 2005

This paper investigates the impacts of turbulent anisotropy on the mean flow and turbulence structures in vegetated open-channel flows. The Reynolds stress model, which is an anisotropic turbulence model, is used for the turbulence closure. Plain open-channel flows and vegetated flows with emergent and submerged plants are simulated. Computed profiles of the mean velocity and turbulence structures are compared with measured data available in the literature. Comparisons are also made with the predictions by the k-$\epsilon$ model and by the algebraic stress model. For plain open-channel flows and open-channel flows with emergent vegetation, the mean velocity and Reynolds stress profiles by isotropic and anisotropic turbulence models were hardly distinguished and they agreed well with measured data. This means that the mean flow and Reynolds stress is hardly affected by anisotropy of turbulence. However, anisotropy of turbulence due to the damping effect near the bottom and free surface is successfully simulated only by the Reynolds stress model. In open-channel flows with submerged vegetation, anisotropy of turbulence is strengthenednear the vegetation height. The Reynolds stress model predicts the mean velocity and turbulence intensity better than the algebraic stress model or the k-$\epsilon$ model. However, above the vegetation height, the k-$\epsilon$ model overestimates the mean velocity and underestimates turbulence intensity Sediment transport capacity of vegetated open-channel flows is also investigated by using the computed profiles. It is shown that the isotropic turbulence model underestimates seriously suspended load.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.8 no.4
• /
• pp.91-100
• /
• 2008

Numerical implementation with a Cartesian cut-cell method is conducted in this study. A Cartesian cut-cell method is an easy and efficient mesh generation methodology for complex geometries. In this method, a background Cartesian grid is employed for most of computational domain and a cut-cell grid is applied for the peculiar grids where the flow characteristics are changed such as solid boundary to enhance the accuracy, applicability and efficiency. Accurate representation of complex geometries can be obtained by using the cut-cell method. The cut-cell grids are constructed with irregular meshes which have various shape and size. Therefore, the finite volume method is applied to numerical discretization on a irregular domain. The HLLC approximate Riemann solver, a Godunov-type finite volume method, is employed to discretize the advection terms in the governing equations. The weighted average flux method applied on the Cartesian cut cell grid for stabilization of the numerical results. To validate the numerical model using the Cartesian cut-cell grids, the model is applied to the rectangular tank problem of which the exact solutions exist. As a comparison of numerical results with the analytical solutions, the numerical scheme well represents flow characteristics such as free surface elevation and velocities in x-and y-directions in a rectangular tank with the Cartesian and cut-cell grids.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.4
• /
• pp.491-498
• /
• 2018

The cost evaluation for voyage route planning in an ice-covered sea is one of the major topics among ship owners. Information of the ice properties, such as ice type, concentration of ice, ice thickness, strength of ice, and speed-power relation under ice conditions are important for determining the optimal route in ice and low operational cost perspective. To determine achievable speed at any designated pack ice condition, a model test of resistance, self-propulsion, and overload test in ice and ice-free water were carried out in a KRISO ice tank and towing tank. The available net thrust for ice and an estimation of the ice resistance under any pack ice condition were also performed by I-RES. The in-house code called 'I-RES', which is an ice resistance estimation tool that applies an empirical formula, was modified for the pack ice module in this study. Careful observations of underwater videos of the ice model test made it possible to understand the physical phenomena of underneath of the hull bottom surface and determine the coverage of buoyancy. The clearing resistance of ice can be calculated by subtracting the buoyance and open water resistance form the pre-sawn ice resistance. The model test results in pack ice were compared with the calculation results to obtain a correlation factor among the pack ice resistance, ice concentration, and ship speed. The resulting correlation factors were applied to the calculation results to determine the pack ice resistance under any pack ice condition. The pack ice resistance under the arbitrary pack ice condition could be estimated because software I-RES could control all the ice properties. The available net thrust in ice, which is the over thrust that overcomes the pack ice resistance, will change the speed of a ship according to the bollard pull test results and thruster characteristics (engine & propulsion combination). The attainable speed at a certain ice concentration of pack ice was determined using the interpolation method. This paper reports a procedure to determine the attainable speed in pack ice and the sample calculation using the Araon vessel was performed to confirm the entire process. A more detailed description of the determination of the attainable speed is described. The attainable speed in 1.0 m, 90% pack ice and 540 kPa strength was 13.3 knots.

• Korean Journal of Materials Research
• /
• v.11 no.5
• /
• pp.405-412
• /
• 2001

Corrosion tests were carried out in $360^{\circ}C$ water and $360^{\circ}C$ 70ppm LiOH solution to investigate the corrosion behavior of new zirconium alloys (Zr-0.4Nb-0.8Sn-xFeCrMn, Zr-0.2Nb-1.1Sn-xFeCrMn, Zr-1.0Nb-xFeCu). Microstructures of tested alloys were analyzed by optical microscope and TEM. The cross-sectional surface and crystalline structure of the oxide layer were analyzed by SEM and XRD. From the results of corrosion test, all the alloys showed higher corrosion rates in $360^{\circ}C$ 70ppm LiOH aqueous solution thats in $360^{\circ}C$ water. Especially, high Nb-containing alloy exhibited the acceleration of corrosion rate in LiOH solution. The low Nb- and Sn-added alloys showed better corrosion resistance than the Sn- free high Nb alloy. from the effect of final annealing on the corrosion, it was observed that the partially recrystallized alloys showed better corrosion resistance than fully recrystallized alloys. This would be related to the size and distribution of the second phase particles.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.2
• /
• pp.106-121
• /
• 2020

In order to make harbor outskirt facilities robust using the reliability-based design, probabilistic models of wave heights at varying stage of shoaling process optimized for Korean sea waves are prerequisite. In this rationale, we numerically simulate the nonlinear shoaling process of random waves over the beach with a sandbar at its foreshore. In doing so, comprehensive numerical models made of spatially filtered Navier-Stokes Eq., LES [Large Eddy Simulation], dynamic Smagorinsky turbulence closure were used. Considering the characteristics of swells observed at the east coast of Korean Peninsula, random waves were simulated using JONSWAP wave spectrum of various peak enhancement coefficients and random phase method. The coefficients of probabilistic models proposed in this study are estimated from the results of frequency analysis of wave crests and its associated trough detected by Wave by Wave Analysis of the time series of numerically simulated free surface displacements based on the threshold crossing method. Numerical results show that Modified Glukhovskiy wave height distribution, the most referred probabilistic models at finite water depth in the literature, over-predicts the occurring probability of relatively large and small wave heights, and under predicts the occurrence rate of waves of moderate heights. On the other hand, probabilistic models developed in this study show vary encouraging agreements. In addition, the discrepancy of the Modified Glukhovskiy distribution from the measured one are most visible over the surf zone, and as a result, the Modified Glukhovskiy distribution should be applied with caution for the reliability-based design of harbor outskirt facilities deployed near the surf-zone.

• Journal of Korea Water Resources Association
• /
• v.52 no.10
• /
• pp.697-706
• /
• 2019

A multiphase flow modeling approach equipped with a hybrid turbulence modeling method is applied to compute the gravity currents in a rectangular channel. The present multiphase solver considers the dense fluid, the less-dense ambient fluid and the air above free surface as three phases with separate flow equations for each phase. The turbulent effect is simulated by the IDDES (improved delayed detach eddy simulation), a hybrid RANS/LES, approach which resolves the turbulent flow away from the wall in the LES mode and models the near wall flow in RANS mode on moderately fine computational meshes. The numerical results show that the present model can successfully reproduce the gravity currents in terms of the propagation speed of the current heads and the emergence of large-scale Kelvin-Helmholtz type interfacial billows and their three dimensional break down into smaller turbulent structures, even on the relatively coarse mesh for wall-modeled RANS computation with low-Reynolds number turbulence model. The present solutions reveal that the modeling approach can capture the large-scale three dimensional behaviors of gravity current head accompanied by the lobe-and-cleft instability at affordable computational resources, which is comparable to the LES results obtained on much fine meshes. It demonstrates that the multiphase modeling method using the hybrid turbulence model can be a promising engineering solver for predicting the physical behaviors of gravity currents in natural environmental configurations.

• Korean Journal of Remote Sensing
• /
• v.34 no.6_2
• /
• pp.1273-1282
• /
• 2018

Recently, active research on the Arctic Ocean has been conducted due to the influence of global warming and new Arctic ship route. Although previous studies already calculated quantitative extent of sea ice using passive microwave radiometers, melting at the edge of sea ice and surface roughness were hardly considered due to low spatial resolution. Since Sentienl-1A/B data in Extended Wide (EW) mode are being distributed as free of charge and bulk data for Arctic sea can be generated during a short period, the entire Arctic sea ice data can be covered in high spatial resolution by mosaicking bulk data. However, Sentinel-1A/B data in EW mode, especially in HV polarization, needs significant radiometric correction for further classification. Thus, in this study, we developed algorithms that can correct thermal noise and scalloping effects, and confirmed that Arctic sea ice and open-water were well classified using the corrected dual-polarization SAR data.

• Journal of Korea Water Resources Association
• /
• v.55 no.2
• /
• pp.147-157
• /
• 2022

The backward-facing step (BFS) is a benchmark geometry for analyzing flow separation occurred at the edge and resulting development of shear layer and recirculation zone that are occupied by turbulent flow. It is important to accurately reproduce and analyze the mean flow and turbulence statistics of such flows to design physically stable and performance assurance structure. We carried out 3D RANS computations with widely used, two representative turbulence models, k-ω SST and RNG k-ε, to reproduce BFS flow at the Reynolds number of 23,000 and the Froude number of 0.22. The performance of RANS computations is evaluated by comparing numerical results with an experimental measurement. Both RANS computations with two turbulence models appear to reasonably well reproduce mean flow in the shear layer and recirculation zone, while RNG k-ε computation results in about 5% larger velocity between the outer edge of boundary layer and the free surface above the recirculation zone than k-ω SST computation and experiment. Both turbulence models underestimate the shear stress distribution experimentally observed just downstream of the sharp edge of BFS, while shear stresses computed in the boundary layer downstream of reattachment point are agree reasonably well with experimental measurement. RNG k-ε modeling reproduces better shear stress distribution along the bottom boundary layer, but overestimates shear shear stress in the approaching boundary layer and above the bottom boundary layer downstream of the BFS.

• The Journal of Engineering Geology
• /
• v.29 no.1
• /
• pp.51-68
• /
• 2019

This study performs to evaluate adaptability of seashore filtering type seawater-intake which adapts dua1 filter well alternative for direct seawater-intake. This study varies filter condition of seashore free surface aquifer which is composed of sand layer then installs real size dual filter well and single filter well to evaluate water permeability and proper pumping amount according to filter condition. According to result of step aquifer test, it is analysed that 110.3% synergy effect of water permeability coefficient is happened compare to single filter since dual filter well has better improvement. dual filter has higher water permeability coefficient compare to same pumping amount, this means dual filter has more improved water permeability than single filter. According to analysis result of continuous aquifer test, it is evaluated that dual filter well (SD1200) has higher water permeability than single filter well (SS800) by analysis of water permeability coefficient using monitoring well and gauging well, it is also analysed dual filter has 110.7% synergy effect of water permeability coefficient. As a evaluation result of pumping amount according to analysis of water level dropping rate, it is analysed that dual filter well increased 122.8% pumping amount compare to single filter well when water level dropping is 2.0 m. As a result of calculating proper pumping amount using water level dropping rate, it is analysed that dual filter well shows 136.0% higher pumping amount compare to single filter well. It is evaluated that proper pumping amount has 122.8~160% improvement compare to single filter, pumping amount improvement rate is 139.6% compare to averaged single filter. In other words, about 40% water intake efficiency can be improved by just installation of dual filter compare to normal well. Proper pumping amount of dual filter well using inflection point is 2843.3 L/min and it is evaluated that daily seawater intake amount is about $4,100m^3/day$ (${\fallingdotseq}4094.3m^3/day$) in one hole of dual filter well. Since it is possible to intake plenty of water in one hole, higher adaptability is anticipated. In case of intaking seawater using dual filter well, no worries regarding damages on facilities caused by natural disaster such as severe weather or typhoon, improvement of pollution is anticipated due to seashore sand layer acts like filter. Therefore, It can be alternative of environmental issue for existing seawater intake technique, can save maintenance expenses related to installation fee or damages and has excellent adaptability in economic aspect. The result of this study will be utilized as a basic data of site demonstration test for adaptation of riverside filtered water of upcoming dual filter well and this study is also anticipated to present standard of well design and construction related to riverside filter and seashore filter technique.

• Korean Journal of Soil Science and Fertilizer
• /
• v.18 no.3
• /
• pp.289-296
• /
• 1985

To clarify the soil characteristics and genesis of the clayey terrace soils distributed in Yeongcheon (inland region) and Yeongil (coastal region), 9 typifying pedons including a recent local deposits were studied. The physico-chemical properties found were summarized as follows; 1. The contents of available water of the terrace soils were ranged from 12.0 to 20.3%. They were increased with increaseing clay content up to 35%, and then decreased at clay content of 35 to 55%, but were maintained about 12 to 13% when the clay content exceeds 55%. 2. The amount of water stable aggregates in the surface soils were ranged from 55.0 to 81.1%, and were low in cultivated upland where water erosion were accelerated while high in the paddy where the double cropping were practiced. The index of aggregate in cultivated upland were higher in the surface soils but the paddy soils were higher in subsoils. 3. The colloidal reactions checked by the differences between pH in $0.01N-Na_2SO_4$ and $H_2O$ indicate exchangeable acidity to exchangeable neutrality, which means a sort of Red-yellow soil forming process had been prevailed during genesis of the terrace soils. 4. Cation exchange capacity of the clayey terrace soils were averaged about 15.8-20.2me/100g of fine earth and 35.6-52.6me/100g of clay. The rates of base saturation were 32.7 to 57.6% and the rates decreased with increasing the elevations of terraces. 5. The activities of free iron oxides were 0.061 to 0.739 and the values decreased with increasing elevation while the crystallinity index of iron oxides positively correlated with elevations the values ranged from 0.067 to 0.537.