• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.64-64
• /
• 2016

하천은 인간에게 용수의 이용 및 하천호안의 휴식처로써의 이용을 통해 직접적인 영향을 주고, 하천구조물의 심미적 영향, 랜드마크로써의 역할을 통해 간접적인 영향을 준다. 또한, 하천은 하천생태계에 서식하는 동 식물에게 영향을 준다. 그러나 하천유사로 인해 통수능이 감소하고, 하천구조물 주변에 침식을 야기할 뿐만 아니라, 댐과 저수지에 유사의 퇴적으로 저수용량의 감소시킨다. 그러므로 이를 예측하는 것은 경제적, 환경적으로 중요하다. 하상변동의 모의를 위해 기존의 2차원 모형은 만곡흐름에서 유동의 helical flow를 고려하지 않아 예측이 부정확하였다. 본 연구에서는 천수방정식을 이용한 하상변동 수치모의에 helical flow의 영향을 고려하였다. 하천과 같은 천수영역에서의 흐름 및 하상변동을 해석하기 위하여 수심평균 된 Navier-Stokes equations인 천수방정식을 이용하였다. 지배방정식은 곡선 좌표계에서 유한체적법으로 차분하였고, 비엇갈림격자를 사용하였다. 지배방정식의 닫힘 문제를 해결하기위해 0-방정식 난류모형을 사용하였고, "time marching" 기법의 적용을 위해 계산단계분할 방법을 이용하였다. 비엇갈림격자의 사용으로 인해 검사체적의 면에서의 유속이 필요하여 pressure-velocity coupling을 사용하여 유속의 진동을 줄였다. 또한, 만곡부의 helical flow를 모의하기위해 helical flow intensity model을 도입하였다. 앞에서 계산한 흐름을 바탕으로 유사량 산정공식과 Exner 방정식을 이용하여 하상변동을 모의하였다. 흐름의 검증, helical flow의 영향에 대한 확인, 하상변동의 적용을 위해 선행연구의 실험이 사용되었다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2022.05a
• /
• pp.245-245
• /
• 2022

최근 국내 기후변화에 따른 국지성 집중호우로 인한 시간당 강우량의 증가로 도로부 유출량의 증가와 배수관거에서의 내수배제 불량에 따른 도심지 내수침수 피해가 증가함에 따라 이를 해결하기 위한 우수유출저감시설이 설치되고 있다. 그러나 대단위의 지하 저류시설의 지속적인 설치는 과밀화된 도심지에서 설치 지하공간의 구조적인 한계 및 적정 설치 위치의 미확보 등의 다양한 문제가 발생하여 저류시설의 침수저감 효과에 대한 추가적이고 새로운 저류시설에 대한 연구가 필요한 실정이다. 이에 내수 침수 저감 및 배수 능력 향상을 위한 도로 배수시설과 연계된 도로 측구부 저류시스템 구축이 필요하다. 이를 위해 역류 방지 및 노면수 저류 빗물받이에 적용되는 부력식 역류차단장치를 개발하였으며, 역류차단장치의 최적 형상 개발을 위해서 기존 빗물받이 연결관과의 통수능 비교 및 분석이 필요한 실정이다. 따라서 본 연구에서는 기존 빗물받이 연결관 및 연결관 내에 역류차단장치가 적용된 역류차단 빗물받이의 흐름분석을 위해 Fluent 모형을 이용하여 3차원 수치모의를 수행하였다. 수치모의 구성으로는 전체 형상을 40×50cm의 빗물받이 유입부와 50×50cm의 빗물받이로 결정하고 격자는 빗물받이 내부의 복잡한 3차원 흐름을 모의하기 위해 1.2~2mm 크기로 생성하였다. 다상유동 해석을 위해 VOF(Volume of Fluid)방법을 적용하였고, 수치해석 방법으로는 비정상류, 난류 모형으로는 SST k-𝜔모형을 적용하였다. 해석조건으로는 김정수(2021) 등이 제시한 4차선 기준 설계빈도별(5~30년) 빗물받이 유입유량을 산정하여 빗물받이 유입조건으로 선정하였으며, 빗물받이와 연결관에서의 통수능력 분석 조건으로는 빗물받이에 기존 연결관이 부착된 조건과 연결관 내에 역류차단장치가 설치되어 역류차단장치가 개방된 조건에서의 통수능을 비교하였으며, 역류상황을 가정한 연결관에서의 통수능을 비교하기 위하여 역류차단장치의 개폐정도를 15도(통수단면 33%감소) 닫힌 상태 및 30도(통수단면 67% 감소) 닫힌 상태 조건을 대상으로 빗물받이와 연결관에서의 흐름을 모의하였다. 수치모의 결과 역류차단장치의 계폐조건에 상관없이 5년 빈도유입량 조건에서는 완전 배수가 되었으며, 개폐조건 15도에서는 10년 빈도의 유입량에서는 완전 배수가 되었으나 20년 빈도 이상의 유입량 조건에서 빗물받이 유입부로의 역류가 발생하였으며, 개폐조건 30도에서는 5년 빈도 이상 유입량 조건에서 빗물받이 유입부로 역류가 발생하는 것으로 나타났다. 특히, 30년 빈도 이상의 유입량부터는 빗물받이 연결관 내에 역류차단장치 개페조건과 관계없이 빗물받이 유입부로의 역류로 인한 도로 침수가 발생하기 때문에 유휴공간인 도로 측구부를 저류공간으로 활용할 수 있는 도로 측구부 저류시스템의 구축은 필수적이라고 판단되며, 유량 조건에 따른 빗물받이 내부 와 흐름과 유출부에서의 유속 변화 특성을 확인하였다. 그러므로 측구 저류조 개발 형상과 연결한 3차원 흐름의 구현 및 분석에 Fluent 모형의 적용이 가능하다고 판단된다.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.1
• /
• pp.178-189
• /
• 1989

An experimental study has been performed to study the effect of the base slant angle of a 1/10 scale two-dimensional vehicle-like body on its wake flow including the recirculating region, where the simplified shape of the body has been originated from a profile of a domestic passenger car. In the case of a Reynolds number based on the length of the model R=7.96*10$^{5}$ , the surface pressure coefficient, the mean velocity and the turbulent stresses have been measured, while the flow visualization technique using wool tuft has been adopted as well. When the base slant angle of the model is 15.deg., the free stream flowing parallel to the slant is observed to be separated from the lower edge of the slant, thus forming the smallest recirculating region. When the base slant angles are 30.deg. and 45.deg., the free streams are separated from the upper edge of the slant and the sizes of the recirculating zones are observed to be almost the same as when the base slant angle is 0.deg. From these observations, it is conjectured that between the base slant angles of 15.deg. and 30.deg. there exists a critical angle at which the size of the recirculating region becomes minimum and as the slant angle becomes larger than this critical angle the separation line moves along the slant towards the rear edge of the roof. Through the flow visualization technique, the existence of the two counter-rotating bubbles in the recirculating region has been clearly observed and verified.

• Journal of computational fluids engineering
• /
• v.15 no.3
• /
• pp.39-45
• /
• 2010

The turbulent flows in a tunnel mock-up($10L{\times}0.5W{\times}0.25H$ m3 : scale reduction 1/20) with rectangular cross section were investigated. The instantaneous velocity fields of Re = 49,029, 89,571 were measured by the 2-D PIV system which is consisted of double pulsed Nd:Yag laser and the tracer particles in the straight-duct mock-up where the flows were fully developed. The mean velocity profiles were taken from the ensemble averages of 1,000 instantaneous velocity fields. Simultaneously, numerical simulations(RANS) were performed to compare with experimental data using STREAM code. Non-linear eddy viscosity model (NLEVM : Abe-Jang-Leschziner Eddy Viscosity Model) was employed to resolve the turbulent flows in the duct. The calculated mean velocity profiles were well compared with PIV results. In the log-law profiles, the experimental data were in good agreement with numerical simulations all the way to the wake region except the viscous sub-layer (near wall region).

• Journal of the Society of Naval Architects of Korea
• /
• v.45 no.6
• /
• pp.631-636
• /
• 2008

In this paper, numerical analysis based on the RANS equation and the Realizable ${\kappa}-{\varepsilon}$ turbulence model is carried out for flows around an axisymmetric body at three Reynolds numbers($1.22{\times}10^7$, $1.0{\times}10^8$, $1.5{\times}10^8$) and the numerical results are compared with experiments data. Computed velocity distributions agree well with experiments as the Reynolds number increases. Pressure distributions agree well with the results of the potential flow except the tail region but differ from experiments for the parallel middle body as well as tail region. Pressure gradients show a good agreement with those of potential flow and experiment except the tail region. Friction coefficients show that the numerical results generally are lower than the experimental results estimated from the measured velocity. The difference of friction coefficients between the calculation and the experiment increases with growing of a boundary layer.

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

여름철의 강우시 저수지에 주로 유입되는 고탁수층은 저수지내의 수체에 대하여 밀도차와 유속차를 가지며, 난류혼합을 일으키면서 주변수를 유입하며 하류방향으로 이송 확산되어 간다. 임하댐과 같은 대형 저수지는 수심이 30 m가 넘는 관계로 여름철에 2개의 수온약층이 존재한다. 따라서 댐체에 도달한 고탁수층은 저수지내의 밀도성층으로 인하여 주로 중층에 분포하게 되며, 가을철에 발생하는 수평확산과 전도현상으로 인해 저수지 전역에 분포하게 되어 탁질입자의 분포에 따라 탁수현상의 장기화를 유발하기도 한다. 이와 같은 탁수문제의 저수지내 대책으로는 홍수기에 고탁수층을 우선적으로 취수하는 방법이 있다. 고탁수층은 이와같이 수온약층과 밀접한 관계를 가지며, 취수탑 주변의 선택취수 결과에 영향을 받는다. 형성된 탁수층을 효율적으로 선택배제하기 위해서는 취수시설에 접근하는 성층흐름을 정확하게 이해해야 한다. 본 연구에서는 3차원 CFD를 사용하여 밀도성층에 따른 취수탑 주변의 접근흐름을 수치해석하여 방류수심, 방류유속 및 밀도성층구조의 선택배제에 대한 영향을 분석하였다. 임하댐 취수설비에 적용한 결과에 의하면 상층, 중층, 하층에서 취수탑 문비를 개도하였을 경우, 취수탑 주변에서만 유속이 증가하였을 뿐 저수지 내에서는 유속이 크게 증가되지 않은 것으로 나타났다. 하지만 문비 개도구간의 변화에 따라 수심별 유속분포는 크게 변화하는 것을 확인할 수 있다. 따라서 문비를 조절하는 것은 저수지 전체의 유동을 깨뜨리지 않으면서 취수탑 인근에서 선택취수가 가능하다는 것을 확인할 수 있었다. 또한 본 연구에서는 정밀한 현장조사를 수행하여 저수지와 방류수의 탁도변화를 모니터링하였으며, 취수탑 개구부 주변의 연직 유속분포를 측정하였다. 3차원 수치모의 결과와 현장에서 관측한 유속장을 비교함으로써 본 연구에서 제시한 실제 탁수배제능력을 검증하였다.를 구축하였다는데 의의가 있다.로와 접하는 건물의 경우 모서리부 광고 효과가 지배적이며 대부분 곡선돌출형이 사용되고 있었다. 그러므로 모서리 저층부를 필로티로 계획하여 보행흐름을 원활하게 하고 대신 입면을 투명하게 하여 간접광고(내부전시) 효과를 유도하는 것이 좋다. 특히 원형모서리는 건물 특화 성격이 강하므로 불가피할 경우 소형 액센트 광고 위치를 미리 벽면으로 할애하는 것이 경관 및 입면계획에 유리한 것으로 분석되었다. 불확실도 해석모형 등의 새로운 기능을 추가하여 제시하였다. 모든 입출력자료는 프로젝트 단위별로 운영되어 data의 관리가 손쉽도록 하였으며 결과를 DB에 저장하여 다른 모형에서도 적용할 수 있도록 하였다. 그리고 HyGIS-HMS 및 HyGIS-RAS 모형에서 강우-유출-하도 수리해석-범람해석 등이 일괄되게 하나의 시스템 내에서 구현될 수 있도록 하였다. 따라서 HyGIS와 통합된 수리, 수문모형은 국내 하천 및 유역에 적합한 시스템으로서 향후 HydroInformatics 구현을 염두에 둔 특화된 국내 수자원 분야 소프트웨어의 개발에 기본 토대를 제공할 것으로 판단된다.았다. 또한 저자들의 임상병리학적 연구결과가 다른 문헌에서 보고된 소아 신증후군의 연구결과와 큰 차이를 보이지 않음을 알 수 있었다. 자극에 차이가 있지 않나 추측되며 이에 관한 추후 연구가 요망된다. 총대장통과시간의 단축은 결장 분절 모두에서 줄어들어 나타났으나 좌측결장 통과시간의 감소 및 이로 인한 이 부위의 통과시간 비율의 저하가 가장 주요하였다. 이러한 결과는 차가운 생수 섭취가 주로 결장 근위부를 자극하는 효과를 발휘하는 것이 아닌가 해석된다. 이와 같은 연구결과를 통해 생다시마를 주원료로 개발된 생다시마차와 생다시마 음료가 만성 기능성 변비 증세를 개선하는 효능이 잠재적으로 있음을 확인하였다. 그러나 생약제재의 변비약 수준으로

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.1
• /
• pp.81-96
• /
• 2007

Nonlinear shoaling characteristics over surf zone are numerically investigated based on spatially averaged NavierStokes equation. We also test the validity of gradient model for turbulent stresses due to wave breaking using the data acquainted during SUPERTANK LABORATORY DATA COLLECTION PROJECT(Krauss et al., 1992). It turns out that the characteristics length scale of breaking induced current is not negligible, which firmly stands against ever popular gradient model, ${\kappa}-{\varepsilon}$ model, but favors Large Eddy Simulation with finer grid. Based on these observations, we model the residual stress of spatially averaged NavierStokes equation after Lagrangian Dynamic Smagorinsky(Meneveau et al., 1996). We numerically integrate newly proposed wave equations using SPH with Gaussian kernel function. Severely deformed water surface profile, free falling water particle, queuing splash after landing of water particle on the free surface and wave finger due to structured vortex on rear side of wave crest(Narayanaswamy and Dalrymple, 2002) are successfully duplicated in the numerical simulation of wave propagation over uniform slope beach, which so far have been regarded very difficult features to mimic in the computational fluid mechanics.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.4
• /
• pp.244-252
• /
• 2018

This paper provides quantification of the effects of the turbulence model and grid refinement on the analysis of tip vortex flows by using the RANS(Reynolds averaged Navier-Stokes) method. Numerical simulations of the tip vortex flows of the NACA $66_2$-415 elliptic hydrofoil were conducted, and two turbulence models for RANS closure were tested, i.e., the Realizable $k-{\varepsilon}$ model and the Reynolds stress transport model. Numerical results were compared with available experimental data, and it was shown that the data for the Reynolds stress transport model that were computed on the finest grid system had better agreement in reproducing the development and propagation of the tip vortex. The Realizable $k-{\varepsilon}$ model overestimated the turbulence level in the vortex core and showed a diffusive behavior of the tip vortex. The tip vortex cavitation on the hydrofoil and its trajectory also showed good agreement between the current numerical results that were obtained using the Reynolds stress transport model and the results observed in the experiment.