• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.92-92
• /
• 2019

최근 기존의 격자방식의 해석 방법을 벗어나 해석 영역에 대한 분할을 별도로 고려치 않는 수치기법의 실무적 적용사례가 증가 하고 있으며, 이러한 방식중 SPH(Smoothed Particle Hydrodynamics) 방식이 근자에 수자원 분야에서도 도입되어 관수로 및 개수로 해석 또는 복합해석 등에 활용된 바 있다. 최초 도입된 무격자방식의 모형들은 주로 복잡한 형상을 지니는 유체기계 등에 활용성이 높았던 바, 큰 규모의 해석 도메인을 가지는 수자원 분야에서의 SPH의 실무적용 평가와 효율성의 확보를 위해서, 본 연구는 국내 댐을 시범 대상으로 하여 SPH 수치해석 툴을 적용하고자 하였다. 분석 대상댐은 국내 P댐으로서 관리수위의 변동은 크지 않으나, 댐 직상류의 만곡이 심하고 다수의 대규모 취수구를 가진 곳으로 상시 발전방류 및 수시 댐 수문방류에 의해 유체의 흐름이 2,3차원의 복잡성을 띄고있기 때문에, 3차원 전산유체역학 Tool의 적용이 적절한 것으로 판단하였다. 해석을 위해 하류경계조건을 댐축과 문비로 설정하였고 상류 1km까지를 해석의 도메인으로 설정하였다. 소요시간을 줄이기 위해 여러 번의 모의를 거쳐 입자의 평균 입경은 0.6m로 제안하고 시격은 1초 미만(평균 0.5초)로 설정하였다. 수문 및 발전방류는 해당댐의 1~2년 빈도 수준에 해당하는 $5,000m^3/s$ 이하의 유량을 기준으로 하여 모의를 수행하였다. 모의의 안정성을 확보한 이후에는 해당 댐지역의 하류영향을 고려한 문비제어를 반영한 다양한 방식의 수문운영 및 취수지점의 순간 수위 영향을 검토하였다. 그 결과로 본 모의에서는 특정한 수문의 운영 조건에서는 댐수위 계측지점과 인근 취수지점 간에도 0.2m 수준의 순간 수위차가 발생할 수 있음을 보였으며, 이는 경우에 따라 취수시설의 일시적 장애요소로 발생할 수 있음을 의미한다. 따라서, 현재의 취수구조물과 문비운영 특성에 따라 발생가능한 취수장애를 줄일 수 있는 운영조건의 탐색을 위해서 수치모의를 추가로 하였으며, 이 때 댐축 상류의 유속분포에 대한 추가 검토도 수행하였다. 다만, 댐에서 방류시 하류조건에 대한 검토는 추후 보강되어야 할 것으로 판단된다.

• Journal of Korea Water Resources Association
• /
• v.46 no.6
• /
• pp.597-611
• /
• 2013

This paper presents the methodology for construction of time-area curve via the width function and thereby rational estimation of time of concentration and storage coefficient of Clark model within the framework of method of moments. To this end time-area curve is built by rescaling the grid-based width function under the assumption of pure translation and then the analytical expressions for two parameters of Clark model are proposed in terms of method of moments. The methodology in this study based on the analytical expressions mentioned before is compared with both (1) the traditional optimization method of Clark model provided by HEC-1 in which the symmetric time-area curve is used and the difference between observed and simulated hydrographs is minimized (2) and the same optimization method but replacing time-area curve with rescaled width function in respect of peak discharge and time to peak of simulated direct runoff hydrographs and their efficiency coefficient relative to the observed ones. The following points are worth of emphasizing: (1) The optimization method by HEC-1 with rescaled width function among others results in the parameters well reflecting the observed runoff hydrograph with respect to peak discharge coordinates and coefficient of efficiency; (2) For the better application of Clark model it is recommended to use the time-area curve capable of accounting for irregular drainage structure of a river basin such as rescaled width function instead of symmetric time-area curve by HEC-1; (3) Moment-based methodology with rescaled width function developed in this study also gives rise to satisfactory simulation results in terms of peak discharge coordinates and coefficient of efficiency. Especially the mean velocities estimated from this method, characterizing the translation effect of time-area curve, are well consistent with the field surveying results for the points of interest in this study; (4) It is confirmed that the moment-based methodology could be an effective tool for quantitative assessment of translation and storage effects of natural river basin; (5) The runoff hydrographs simulated by the moment-based methodology tend to be more right skewed relative to the observed ones and have lower peaks. It is inferred that this is due to consideration of only one mean velocity in the parameter estimation. Further research is required to combine the hydrodynamic heterogeneity between hillslope and channel network into the construction of time-area curve.

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

The flow in the meandering channel is characterized by the spiral motion of secondary currents that typically cause the erosion along the outer bank. Hydraulic structures, such as spur dike and groyne, are commonly installed on the channel bottom near the outer bank to mitigate the strength of secondary currents. This study is to investigate the effects of submerged vanes installed in a $90^{\circ}$ meandering channel on the development of secondary currents through three-dimensional numerical modeling using the hybrid RANS/LES method for turbulence and the volume of fluid method, based on OpenFOAM open source toolbox, for capturing the free surface at the Froude number of 0.43. We employ the second-order-accurate finite volume methods in the space and time for the numerical modeling and compare numerical results with experimental measurements for evaluating the numerical predictions. Numerical results show that the present simulations well reproduce the experimental measurements, in terms of the time-averaged streamwise velocity and secondary velocity vector fields in the bend with submerged vanes. The computed flow fields reveal that the streamwise velocity near the bed along the outer bank at the end section of bend dramatically decrease by one third of mean velocity after the installation of vanes, which support that submerged vanes mitigate the strength of primary secondary flow and are helpful for the channel stability along the outer bank. The flow between the top of vanes and the free surface accelerates and the maximum velocity of free surface flow near the flow impingement along the outer bank increases about 20% due to the installation of submerged vanes. Numerical solutions show the formations of the horseshoe vortices at the front of vanes and the lee wakes behind the vanes, which are responsible for strong local scour around vanes. Additional study on the shapes and arrangement of vanes is required for mitigate the local scour.

• Journal of Korea Water Resources Association
• /
• v.34 no.1
• /
• pp.19-30
• /
• 2001

A laterally averaged 2-dimensional hydrodynamic model is newly developed. The coordinate system is first transformed to minimize the effects of irregularity of bottom and surface. The advection terms of the governing equations are then discretized by an upwind scheme. By employing an explicit scheme for longitudinal direction and an implicit scheme for vertical direction, the model is free from restriction of temporal step size caused by a relatively small grid ratio. To demonstrate the applicability of the model, calculated time histories of free surface displacements and distributions of velocity and salinity are compared with the field measurements of the Keum River Estuary before construction of the estuary dam. A reasonable agreement is observed between them.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.12
• /
• pp.1383-1389
• /
• 2011

The current study focuses on the consistent analysis of heat transfer in multichannel volumetric solar receivers used for concentrating solar power. Changes in the properties of the absorbing material and channel dimensions are considered in an optical model based on the Monte Carlo ray-tracing method and in a one-dimensional heat transfer model that includes conduction, convection, and radiation. The optical model results show that most of the solar radiation energy is absorbed within a very small channel length of around 15 mm because of the large length-to-radius ratio. Classification of radiation losses reveals that at low absorptivity, increased reflection losses cause reduction of the receiver efficiency, notwithstanding the decrease in the emission loss. As the average temperature increases because of the large channel radius or small mass flow rate, both emission and reflection losses increase but the effect of emission losses prevails.

• Journal of the Korean GEO-environmental Society
• /
• v.11 no.1
• /
• pp.5-11
• /
• 2010

The objective of this study is to evaluate the effect of operational parameters such as rotation speed and vibrating milling time for the fabrication of photocatalytic $TiO_2$ thin film using aerosol deposition methods. $TiO_2$ powders produced in the range of 1,000-3,000 rpm of rotation speed of centrifugal separator are ineffective on the fabrication of $TiO_2$ thin film by aerosol deposition due to the problem of nozzle powder jam. $TiO_2$ powders controlled by vibrating milling had about 420 nm of average diameter after 2 hr of vibrating milling time. The result of XRD analysis indicated that $TiO_2$ powders had a anatase phase. Vibrating milling methods was considered to be an effective pre-treatment process for $TiO_2$ powder control. Consequently $TiO_2$ photocatalytic thin film with dispersion of anatase crystallites controled by vibrating milling was successfully fabricated by aerosol deposition. The permeation flux of $TiO_2$ photocatalytic thin film with the immobilized $TiO_2$ powder was higher than that of suspended $TiO_2$ powder. Therefore, $TiO_2$ photocatalytic thin film promises to be one of the effective methods for enhancing filtration performance for the treatment of organic pollutants.

• Analytical Science and Technology
• /
• v.28 no.1
• /
• pp.33-39
• /
• 2015

CdS-QD particles are a nano-sized semiconducting crystal that emits light. Their optical properties show great potential in many areas of applications such as disease-diagnostic reagents, optical technologies, media industries and solar cells. The wavelength of emitting light depends on the particle size and thus the quality control of CdS-QD particle requires accurate determination of the size distribution. In this study, CdS-QD particles were synthesized by a simple ${\gamma}$-ray irradiation method. As a particle stabilizer polyvinyl pyrrolidone (PVP) were added. In order to determine the size and size distribution of the CdS-QD particles, sedimentation field-flow fractionation (SdFFF) was employed. Effects of carious parameters including the the flow rate, external field strength, and field programming conditions were investigated to optimize SdFFF for analysis of CdS-QD particles. The Transmission electron microscopy (TEM) analysis show the primary single particle size was ~4 nm, TEM images indicate that the primarty particles were aggregated to form secondary particles having the mean size of about 159 nm. As the concentration of the stabilizer increases, the particle size tends to decrease. Mean size determined by SdFFF, TEM, and dynamic light scattering (DLS) were 126, 159, and 152 nm, respectively. Results showed SdFFF may become a useful tool for determination of the size and its distribution of various types of inorganic particles.

• Journal of the Korean Ceramic Society
• /
• v.39 no.10
• /
• pp.994-1000
• /
• 2002

Oxygen permeability and the mechanical properties of mixed ionic-electronic conductive $La_{0.6}Sr_{0.4}Co_{0.2}Fe_{0.8}O_{3-{\delta}}$ perovskite-type membrane, fabricated by solid state reaction, were investigated with regard to microstructure. The microstructure of the membrane was controlled by changing the sintering temperature and holding time. The average grain size and relative density were evaluated as a function of sintering conditions. As the fraction of grain boundary decreased, oxygen permeability showed a tendency to increase. Especially the maximum oxygen flux of 0.37 ml/$cm^2$${\cdot}$min was measured for the specimen sintered at 1300${\circ}C$ for 10 h, which has high density and relatively large grain size. Fracture strength was dependent on the relative density of sintered body, while fracture toughness increased with average grain size.

• 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 Korean Society of Marine Environment & Safety
• /
• v.23 no.2
• /
• pp.194-199
• /
• 2017

When sea tide and wave velocity change, the behavior of silt protectors underwater changes, and a hydraulic force exceeding the anchor wave force is applied. In this study, the movement mechanism of a silt protector has been analyzed using the mass-spring method. The initial position of the silt protector was in the Jindo area near Gwangpo Port (742-1, Gyupori, Chongdo-myeon, Jindo-gun, Jeonnam, Korea). The tension required to exceed the holding power of the anchor was 0.05 m/s at 318 sec., 0.15 m/s at 77 sec., 0.25 m/s at 43 sec., and 0.3 m/s at 37 sec.. As the anchor started to move from the sea floor and the tide speed increased to 0.01 m/s, anchor movement start time shortened by an average of 11.2 sec.. Compared with when tide was the only affecting factor, the silt protector and anchor were found to have moved 19.7 % at 0.1 m/s, 7.6 % at 0.15 m/s, 5.8 % at 0.2 m/s, 4.3 % at 0.25 m/s and 2.8 % at 0.3 m/s, showing an increase. When wave effect was added to the tide, anchor movement started when the flow rate was slow 7.6 % of the time. With a high flow velocity, anchor movement started without any significant difference less than 4.3 % of the time. When tide speed exceeded 0.13 m/s and the direction of the waves matched, the silt protector was not able to perform due to collisions with surrounding sea structures. When installing a silt protector, the fluid flow situation and the silt protector situation must be carefully analyzed using the mass-spring method to apply the result found in this study.