• Journal of Korea Water Resources Association
• /
• v.45 no.4
• /
• pp.349-360
• /
• 2012

A horizontally curvilinear non-hydrostatic free surface model that was applicable to three-dimensional viscous flows was developed. The proposed model employed a top-layer equation to close kinematic free-surface boundary condition, and an isotropic k-${\varepsilon}$ model to close turbulence viscosity in the Reynolds averaged Navier-Stokes equation. The model solved the governing equations with a fractional step method, which solved intermediate velocities in the advection-diffusion step, and corrects these provisional velocities by accounting for source terms including pressure gradient and gravity acceleration. Numerical applications were implemented to the wind-driven currents in a two-dimensional closed basin, the flow in a steep-sided trench, and the flow in a strongly-curved channel accounting for secondary current by the centrifugal force. Through the numerical simulations, the model showed its capability that were in good agreement with experimental data with respect to free surface elevation, velocity, and turbulence characteristics.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.3 no.1
• /
• pp.13-24
• /
• 1983

A mathematical model for the two-dimensional unsteady flow was developed by introducing Abbott's implicit finite difference operator and double sweep algorithm, which could be applied to simulate the respose of a harbor against the intrusion of long waves through the entrance connected to open sea. In order to improve its accuracy corresponding to the field phenomena, bottom resistance, Coriolis force, wind effect terms were included and wave direction and radiating effect was considered. The result of seiche test was always stable and the amplitude was accurate. Some phase shift was occured, but it could be reduced by using small values of Courant number and many points per a wave length as well. A comparision with the Ippen and Goda's theoritical and hydraulic experimental works was fulfilled.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.3B
• /
• pp.175-183
• /
• 2012

In this study, Galerkin scheme and SU/PG scheme of Petrov-Galerkin family were applied to the shallow water equations and a finite element model for shallow water flow was developed. Numerical simulations were conducted in several flumes with convection-dominated flow condition. Flow simulation of channel with slender structure in the water course revealed that Galerkin and SU/PG schemes showed similar results under very low Fr number and Re number condition. However, when the Fr number increased up to 1.58, Galerkin scheme did not converge while SU/PG scheme produced stable solutions after 5 iterations by Newton-Raphson method. For the transcritical flow simulation in diverging channel, the present model predicted the hydraulic jump accurately in terms of the jump location, the depth slope, and the flow depth after jump, and the numerical results showed good agreements with the hydraulic experiments carried out by Khalifa(1980). In the oblique hydraulic jump simulation, in which convection-dominated supercritical flow (Fr=2.74) evolves, Galerkin scheme blew up just after the first iteration of the initial time step. However, SU/PG scheme captured the boundary of oblique hydraulic jump accurately without numerical oscillation. The maximum errors quantified with exact solutions were less than 0.2% in water depth and velocity calculations, and thereby SU/PG scheme predicted the oblique hydraulic jump phenomena more accurately compared with the previous studies (Levin et al., 2006; Ricchiuto et al., 2007).

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.169-169
• /
• 2015

유사 입자의 크기는 유사의 특성 및 그에 따른 거동 변화에 중요한 영향을 미친다. Stokes 침강 속도 모형에서 유사의 침강 속도에 가장 많은 영향을 주는 인자는 유사의 크기인 것 또한 확인된다. 유사 입자의 크기가 약 $60{\mu}m$보다 작은 유사들은 알갱이 사이의 점착력을 무시할 수 없다. 이로 인해 유사들은 응집 현상을 겪으며 입자 본래의 크기보다 크기가 큰 플럭을 형성하는 점착성 유사로 분류된다. 응집 현상이란, 흐름 내 점착성을 띠는 일차입자(Primary Particle)가 응집과 파괴를 반복하며 플럭을 형성하는 현상을 뜻한다. 입자 간의 충돌을 통해 응집이 진행되며 난류 전단으로 인해 형성된 플럭의 파괴가 발생한다. 많은 연구에서 점착성 유사의 충돌을 야기하는 가장 지배적인 원리는 난류라 알려져 있다. 이러한 응집 현상으로 인하여 플럭의 크기와 밀도는 지속적으로 변화를 겪으며 비점착성 유사와 다른 특징들을 보인다. 흐름에 존재하는 유사의 이동은 이송-확산 방정식을 통해 표현된다. 이송-확산 방정식은 시간 변화에 따른 농도의 변화를 입자의 침강과 난류 및 유사 자체의 특징에 의한 확산으로 해석한다. 침강속도로 대변되는 이송과 달리, 확산은 난류흐름 내에서 유사가 확산되는 정도를 정량화하기 위한 인자가 요구된다. 난류에 의한 유사의 확산은 유사 자체 특성에 따른 물질 확산에 비하여 매우 큰 값을 가지며, 이를 확산 계수로 개념화 한다. 확산계수는 와점성계수와 Schmidt 수(${\sigma}_c$)의 비로 정의된다. ${\sigma}_c$는 난류의 점성과 난류로 인한 부유과정에 의해 유사가 확산되는 정도를 나타낸다. 이에 따라 ${\sigma}_c$의 변화가 유사의 부유 및 침강거동에 많은 영향을 미칠 것이라 판단되나, 국내외에서 수행된 연구 동향에서는 ${\sigma}_c$를 0.5부터 1.0 사이의 상수를 적용하여 수행되었다. 이에 본 연구에서는 ${\sigma}_c$의 크기에 따라 달라지는 유사의 부유 및 침강 변화에 의한 총 부유량을 살펴보고자 한다. 유사의 점착성을 고려할 수 있는 1DV 수치 모형을 이용하여 비점착성 유사와 점착성 유사를 대상으로 수치연구를 수행하며, 유사의 크기 및 ${\sigma}_c$의 변화에 따른 총 부유량 경향을 살펴본다. 그 결과, 점착성 유사는 ${\sigma}_c$의 증가에 따라서 유사의 총 부유량이 증가하는 현상이 나타난 반면 비점착성 유사는 ${\sigma}_c$의 증가에 따라 유사의 총 부유량이 감소하는 경향이 나타났다. 그러나 크기가 아주 작은 비점착성 유사를 대상으로 수치 연구를 수행한 결과, ??에 따른 총 부유량의 경향은 유사의 점착성에서 기인하는 것이 아닌 입자의 크기로부터 야기되는 특성이라는 결론이 도출되었다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.883-888
• /
• 2007

일반적인 하천의 흐름방향으로 발생하는 주흐름(primary flow)에 중첩하여 주흐름 방향의 수직단면에 이차류(secondary flow)가 발생하게 되며 이러한 이차류의 발달은 투입된 오염물질의 횡혼합을 증대시킨다. 오염물질의 혼합은 이송(advection)과 확산(diffusion) 또는 분산(dispersion)의 과정으로 설명되며 본 연구에서는 수로전체의 혼합과정을 설명하기 위해서 이송 확산 방정식을 적용하였다. 본 연구에서는 실험수조를 $150^{\circ}$의 중심각을 갖는 S자 형태의 만곡수로를 제작하여 유량조건은 15, 30, $60\;{\ell}l/sec$의 세 가지 경우로, 수심은 15, 20, 30, 40 cm의 경우로 총 12 케이스의 실험을 수행하였다. 유속장의 측정은 Sontek사의 3차원 micro-ADV(Acoustic Doppler Velocimeter)를 이용하였다. 오염물질 확산실험은 소금물 용액에 주변수와의 밀도차를 없애기 위해서 메탄올 용액을 첨가하여 추적자로 이용하여 농도장의 분석을 일본 KENEK사의 전기전도도계(conductivity meter)와 Gartner사의 DAS(data acquisition system)를 이용하여 횡방향 유속장의 분포와 오염운의 거동을 비교하여 다음과 같은 결론을 얻게 되었다. 주 흐름은 직선구간에서는 중앙에서 최대 유속을 나타내며, 좌우대칭적인 유속분포의 모습을 보이고, 만곡부에서는 수로안쪽을 따라 최대유속이 발생하였다. 수로의 직선구간에서는 최대유속이 발생하는 즉, 중앙에서의 오염물질의 분산이 가장 활발하게 이뤄졌으며 농도의 퍼짐형상인 오염운 역시 만곡부에서는 수로만곡부의 안쪽을 따라 확산 이동함을 알 수 있었다. 만곡부 외측에서는 오염물질의 정체현상이 일시적으로 발생하며, 유속구조의 횡방향 비대칭구조로 인한 종 횡방향의 분리현상이 발생하고, 오염운의 중첩현상이 종방향으로 연속되게 나타난다. 향후 수심방향 거동을 포함한 3차원적 분석이 요구되며 이 연구결과는 2차원적 수치해석의 적용 및 분석 자료로써 이용이 가능하다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.6B
• /
• pp.507-514
• /
• 2011

A numerical model has been developed by employing a finite element method to simulate the depth-averaged 2-D dispersion of the heat pollutant, which is an important pollutant material in natural streams. Among the finite element methods, the Streamline Upwind/Petrov Galerkin (SUPG) method was applied. Also both linear and quadratic elements can be applied so that irregular river boundaries can be easily represented. To show the movement of heat pollutants, the reaction term describing heat transfer was represented as an equation in which sink/source term is proportional to the difference between the equilibrium temperature and water surface temperature. The equation was expressed so that the water surface temperature changes according to the temperature transfer coefficient and the equilibrium temperature. For the calibration of the model developed, analytic and numerical results from a case of rectangular channel with full width continuous injection have been compared in a steady state. The comparisons showed that the numerical results were in good agreement with analytical solutions. The application site was selected from the downstream of Paldang dam to Jamsil submerged weir, and overall length of this site is about 22.5 km. The change of water temperature caused by the discharge from the Guri sewage treatment plant has been simulated, and results were similar to the observed data. Overall it is concluded that the developed model can represent the water temperature changes due to heat transport accurately. But the verification using observed data will further enhance the validity of the model.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.174-174
• /
• 2021

지표면에서 소규모의 요철로 이루어진 미세지형 (Micro Topography)은 자연 지형의 일반적인 특징이다. 지표면에서의 흐름해석 시 이러한 미세지형에 의한 효과를 정확히 고려하기 위해서는 고해상도의 지형격자를 사용하여야 하나, 대부분의 지표유출모델은 미세지형이 존재하지 않는 저해상도의 거친 표면으로 간주함으로써 바닥 마찰 효과를 증가시켜 미세지형에 의한 영향을 대략적으로만 반영한다. 본 연구에서는 보다 정확한 강우-유출 및 유사-유출 해석을 목적으로 고해상도의 지형격자를 사용하여 미세지형을 포함한 건조 및 반건조 지대에서의 흐름을 수치적으로 모의하였다. 미세지형의 형태는 마루와 골 사이가 직선으로 이루어진 파동의 형태로 이상화 하였으며, 파동의 진폭과 파장을 조절하여 다양한 형태의 미세지형을 고려하였다. 수치모형은 흐름의 움직임에 대한 Saint-Venant 방정식과 침식 및 유사이송에 대한 Hairsine-Rose 방정식을 함께 계산하는 통합모형인 tRIBS-FEaST를 사용하였다. 수치 모의 결과에서 나타난 미세지형의 핵심 영향은 최대 유량 및 유사량 도달시간의 지연, 유사 입자 크기별 유사-유출량 증감, 그리고 하천 유출(stream flow)의 생성이었다. 또한, 미세지형의 형태에 따라 미세지형과 강우-유출 및 유사-유출 사이에 비례 혹은 반비례 관계가 성립함을 보였다. 수치 모의 결과를 종합적으로 검토하여 미세지형이 강우-유출과 유사-유출에 미치는 영향에 대해 논의하였으며, 기존 Manning 거칠기 계수를 통한 해석 방법의 적부를 판단하였다.

• Journal of the Korean Geotechnical Society
• /
• v.29 no.12
• /
• pp.45-56
• /
• 2013

Abrams equation could be effectively applied to predict strength of cement-admixed clay and clay-water content to cement content ratio is a fundamental parameter for governing strength. This paper analyses unconfined compression strength varying with $w_c/C$ and curing time using laboratory test results. An attempt is made to identify strength of composite soil of cement and clay according to variation of Abrams coefficients and curing time. The value B, which was considered to be constant value in past researches, needs to be considered as parameter variable with curing time. From Abrams equation a correlation was formed for unconfined compression strength with mixing conditions by $w_c/C$ and curing time as dependent variable. Regression results in this paper could be used to predict strength of cement-admixed clay at various mixing conditions.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.43-43
• /
• 2015

천수방정식을 사용하는 초기 수치모형은 프로드수($F_4$)가 변화하는 흐름 즉, 상류방향과 하류방향으로 전파하는 홍수파를 동시에 해석하기 위해 중앙 차분기법이 필요한 상류(sub-critical flow)와 흐름방향에 따른 상류이송(upwinding)기법이 필요한 사류(super-critical flow)가 나타나는 흐름해석에서 어려움이 있었다. 하지만, 근사 Riemann 해법의 등장으로 흐름방향에 관계없이 특성선을 따라 정확한 상향가중기법의 적용이 가능하게 되어, 천수방정식을 지배방정식으로 하는 수치모형이 더욱 실용적으로 적용될 수 있도록 하였다. 따라서, 현재 근사 Riemann 해법은 Godunov 형 유한체적 기법, 불연속 Galerkin 혹은 Petrov-Galerkin 유한요소기법 그리고 Boussinesq 기법에도 적용되고 있으며, 특히 Godunov 형 유한체적기법과 결합한 근사 Riemann 해법은 댐 붕괴, 하천 범람 그리고 도시 및 해안지역 침수에 이르기까지 여러 가지 문제에 폭넓게 적용되고 있다. 지금까지 홍수 모델링에 적용된 Godunov형 유한체적모형은 정형 사각격자나 비정형 삼각격자 중에서 한가지의 격자 종류만을 적용한 연구가 주로 수행되었으며, 유한요소모형과 같이 이 두 가지 격자를 동시에 적용한 연구는 거의 이루어지지 않고 있다. 일반적으로, 삼각격자는 사각격자와 는 달리 연구유역의 경계나 지형이 복잡한 경우에도 큰 노력없이 격자의 생성이 가능하나, 격자와 노드의 수가 사각격자보다 많아 계산시간이 많이 소요되는 단점이 있다. 반면, 사각격자는 하천과 같이 선형으로 변하는 지형에 대해서는 표현하기가 용이하며 계산시간의 효율성도 뛰어나다. 본 연구에서는 하천, 도시 그리고 해안지역에서의 효율적이고 정확한 홍수 모델링을 위해 삼각 및 사각격자 그리고 이 두 격자를 동시에 고려한 하이브리드 격자의 적용이 가능한 Godunov형 2차원 유한체적 모형을 개발하였다. 그리고 개발모형을 정확해가 있는 댐 붕괴 문제, 실측치가 존재하는 실험하도 및 실제하도에 삼각, 사각 그리고 혼합격자를 생성하여 모의를 수행하고, 각 적용 격자에 따른 정확성과 효율성 및 장점과 단점을 연구하였다.

• 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.