• Journal of Korea Water Resources Association
• /
• v.53 no.spc1
• /
• pp.755-764
• /
• 2020

In South Korea, 850 standard watersheds and 7,807 KRF catchment areas have been used as basic maps for water resources policy establishment, however it becomes necessary to set up new standard maps with a more appropriate scale for the integrated managements of surface water-groundwater as well as water quantity-quality in the era of integrated water management. Since groundwater has a slow flow velocity and also has 3-D flow properties compared to surface water, the sub-catchment size is more effective than the regional watershed for the evaluation of surface water-groundwater interaction. The KRF catchment area, which has averagely a smaller area than the standard watershed, is similar to the sub-catchment area that generally includes the first-order or second-order tributaries. Some KRF catchment areas, which are based on the surface reach, are too small or large in a wide plain or high mountain area. Therefore, it is necessary to revise the existing KRF area if being used as a unit area for integrated management of surface-water and groundwater. A unit watershed with a KRF area of about 5 to 15 ㎢ can be effective as a basic unit for water management of local government considering a tributary composition and the location of groundwater wells, and as well it can be used as a basic tool for water demand-supply evaluation, hydrological observation system establishment, judgment of groundwater permission through a total quantity management system, pollution assessment, and prioritizing water policy, and etc.

• Journal of Korea Water Resources Association
• /
• v.47 no.11
• /
• pp.1061-1066
• /
• 2014

A stable second-order finite volume method was proposed to predict sediment transport under rapidly varied flow conditions such as transcritical flow. For the use under unsteady flow conditions, a sediment transport model was coupled with shallow water equations. HLLC approximate Riemann solver based on a monotone upstream-centered schemes for conservation laws (MUSCL) reconstruction was used for the computation of the flux terms. From the comparisons of dam break flow experiments on erodible beds in one- and two-dimensional channels, good agreements were obtained when proper parameters were provided. Lastly, dam surface erosion problem by overtopped water was simulated. Overall, the numerical solutions showed reasonable results, which demonstrated that the proposed numerical scheme could provide stable and physical results in the cases of subcritical and supercritical flow conditions.

• Journal of Environmental Science International
• /
• v.11 no.1
• /
• pp.41-56
• /
• 2002

In order to estimate the quantitative roles of the tide induced residual currents, density currents and wind driven currents on the residual currents in Kwangynag Bay Korea, numerical experiments with a diagnostic multi-level model were carried out. Density currents were calculated from the temperature and salinity data observed in January, May, August and November 1998. Anti-clockwise circulations are formed at the western inner part of the bay, the location from the east of Myodo Is. to the south of POSCO Co. and the place between Yeosu peninsula and Namhae Is. from the calculation results of the tide induced residual currents. Velocities of the density currents are less than 3cm/s at the western inner part of the bay and about 5cm/s at the southern entrance of the bay. Density currents get strong in summer and weak in autumn. Wind driven currents at the surface layer flow in the directions of the given winds which are the daily mean winds when the temperature and salinity observations are carried out. In the middle and lower layers, however the wind driven currents flow in the opposite direction to the surface currents as supplementary currents. The surface wind driven currents are greater than the tide induced residual currents or the density currents. The calculated residual currents including the tide induced residual currents, density currents and wind driven currents agree with the results of the current observations approximately. In the Bay, the wind driven currents affect on the residual currents greatestly and tide induced residual currents and density currents do in the second place and the third place.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.311-316
• /
• 2001

The thin-layer Navier-Stokes equations are solved for the hypersonic flow over blunt cone configurations with applications to laminar as well as turbulent flows. The equations are expressed in the forms of flux-vector splitting and explicit algorithm. The upwind schemes of Steger-Warming and van Leer are investigated in their ability to accurately predict the heating loads along the surface of the body. A comparison with the second order extensions of these schemes is made and a hybrid scheme incorporating a combination of central differencing and flux-vector-splitting is presented. This scheme is also investigated in its ability to accurately predict heat transfer distributions.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.41-48
• /
• 2013

Ice accretion on the solid surface is an importance factor in assessing the performance of aircraft and wind turbine blade. Changes in the external shape due to ice accretion can greatly deteriorate the aerodynamic performance. In this study, a three-dimensional upwind-type second-order positivity-preserving finite volume CFD scheme based on the unstructured mesh topology is developed to simulate two-phase flow in atmospheric icing condition. The code is then validated by comparing with NASA IRT experimental data on the sphere. The present results of the collection efficiency are found to be in close agreement with experimental data and show improvement near the stagnation region.

• Journal of Hydrospace Technology
• /
• v.2 no.2
• /
• pp.1-13
• /
• 1996

This paper presents the result of a computational study on the wake characteristics of two tanker models, i.e. HSVA and Mystery hull forms. The focus of the study is on the distributions of axial, radial and tangential velocities of the two hull forms in way of the propeller, especially over the propeller disk. The effect of bilge vortices on the velocity distribution is also concerned. For the computation of stern and wake flows of the two hull farms, the incompressible Reynolds-Averaged Wavier-Stokes(RANS) equations are numerically solved by the second order finite difference method, which employs a four stage Runge-Kutta scheme with a residual averaging technique and the Baldwin-Lomax model. The calculated pressure distributions on the hull surface and the axial, radial and tangential velocity distributions over the propeller disk are presented for the two hull forms. Finally, the result of wake analysis for the computed wake distribution over the propeller disk is given in comparison with those for the experimental wake distribution fur the both hull forms.

• Journal of Korean Society of Water and Wastewater
• /
• v.18 no.5
• /
• pp.656-665
• /
• 2004

An Ozone reaction model combined with CFD(Computational Fluid Dynamics) technique was developed in this research, in the simulation of ozonation, hydrodynamic behavior as well as reaction model is important because ozone is supplied to treated water as gas ozone. In order to evaluate hydrodynamic behavior in an ozone contactor, CFD technique was applied. CFD technique elucidated hydrodynamic behavior in the selected ozone contactor, which consisted of three main chambers. Three back-mixing zones were found in the contactor. The higher velocities of water were observed in the second and third compartments than that in the first compartment. The flow of the opposite direction to the main flow was observed near the water surface. Based on the results of CFD simulation, the ozone contactor was divided into small compartments. Mass balance equations were established were established in each compartment with reaction terms. This reaction model was intended to predict dissolved ozone concentration, especially. We concluded that the model could predict favorably the mass balance of ozone, namely absorption efficiency of gaseous ozone, dissolved ozone concentration and ozone consumption. After establishing the model, we discussed the effect of concentration of gaseous ozone at inlet, temperature and organic compounds on dissolved ozone concentration.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.9 no.4
• /
• pp.403-413
• /
• 1985

Hot-wire measurements of second and third-order mean products of velocity fluctuations have been made in the separated, reattached, and redeveloping boundary layer behind a vertical fence. Mean velocity, wall static pressure distributions have also been measured in the whole flow field. Upstream of the reattachment point, the separated shear layer developes as a free mixing layer, but the gradient of the maximum slope thickness, turbulent intensities and the Reynolds shear stress are higher than that of the mixing layer due to initial streamline curvature and the effects of highly turbulent recirculating flow region. In the reattachment region, Reynolds shear stress and triple products near the surface is far more rapid than the decrease of the shear stress; that is the presence of the solid wall has a marked effect on the apparent gradient diffusivity of intensity or shear stress and throws doubts upon the usefulness of the simple gradient diffusivity model in this region.

• Journal of Korean Society of Water and Wastewater
• /
• v.33 no.2
• /
• pp.131-140
• /
• 2019

This study used a packed column reactor and a horizontal flow mesh reactor to examine the removal of copper ions from aqueous solutions using pine bark, a natural adsorbent prepared from Korean red pine (Pinus densiflora). Both equilibrium and nonequilibrium adsorption experiments were conducted on copper ion concentrations of 10mg/L, and the removals of copper ions at equilibrium were close to 95%. Adsorption of copper ions could be well described by both the Langmuir and Freundlich adsorption isotherms. The bark was treated with nitric acid to enhance efficiency of copper removal, and sorption capacity was improved by about 48% at equilibrium; mechanisms such as ion exchange and chelation may have been involved in the sorption process. A pseudo second-order kinetic model described the kinetic behavior of the copper ion adsorption onto the bark. Regeneration with nitric acid resulted in extended use of spent bark in the packed column. The horizontal flow mesh reactor allowed approximately 80% removal efficiency, demonstrating its operational flexibility and the potential for its practical use as a bark filter reactor.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.38 no.3
• /
• pp.92-100
• /
• 1996

In order to control the water quality in rivers or lakes, it is needed to evaluate accurate amount of pollutant loadings from watersheds. The daily pollutant loadings were simulated using the pollutant loading calculation model which was composed of mathematical equations superimposed on the TANK model. The calibration of runoff and pollutant loading parameters were carried out with observed data, using a trial-and-error method. In addition, the proposed model was applied to evaluate its applicability for the representative watershed, the Bokha river watershed, Icheon city, Korea. The parameters of SS and T-P showed large values in the first tank while T-N showed large in the second tank. As a result of simulating the daily pollutant loadings by the pollutant loading calculation model, all of SS, T-N and T-P loadings were increased or decreased according to the amount of runoff discharge. Especially, it was apparent that SS and T-P loadings were significantly influenced by the runoff variation when it was rain. These results could partly explain that SS and T-P would occur mainly from the surface runoff while T-N would occur from both surface and subsurface flow.