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

트렌치 수로에서 동수압을 고려한 자유수면 흐름을 해석하고 난류 모형의 적용성을 제시하였다. 본 연구에 사용된 지배방정식으로 비정상 상태의 비압축성 유체에 대한 연속방정식과 비점성 Navier-Stokes 방정식을 사용하였다. 난류완결문제를 해결하기 위해서 $\kappa-\varepsilon$방정식을 사용하여서 난류 와점성을 구할 수 있다. 자유수면과 동수압을 고려하기 위하여 3단계로 나누어서 해석하였다. 제 1단계에서는 운동량방정식을 연직방향에 대해 음해적으로 차분하였다. 제2단계에서는 유속과 동수압 보정항으로 이루어진 식을 연속방정식에 대입하여 타원형 방정식인 동수압-포와송 방정식을 해석하여 얻어지는 유속은 질량보존법칙을 만족하게 된다. 마지막으로 자유수면과 최종유속을 보정 및 계산하였다. 본 연구에서 제시한 수치모형을 검증하기 위해서, 트렌치 수로에서 난류의 흐름에 변화를 고려하기 위해 수치모의를 하였다. 전반적으로 수치모의에 의한 결과와 실험 자료가 일치하는 경향을 보였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2020.06a
• /
• pp.149-149
• /
• 2020

하천에서 물 흐름이 보와 댐과 같은 수공구조물을 지날 때 일반적으로 흐름상태에 다양하고 급진적인 변화가 발생한다. 특히 흐름이 구조물을 지나면서 사류(supercritical flow)로 변하고 다시 상류(subcritical flow)로 복원되면서 일어나는 도수(hydraulic jump) 현상은 수위의 급변화, 흐름 에너지 소산, 변동성이 강한 압력 분포 등이 특징이다. 이러한 흐름 특성들은 보나 여수로와 같은 수공구조물 자체의 성능뿐만 아니라 이들 수공구조물의 하류에서 발생하는 국부세굴로 인해 구조물의 안정성에 부정적인 영향을 줄 수 있다. 따라서 수공구조물을 설계할 때는 이들 구조물을 통과하는 흐름의 비정상 난류 흐름 특성을 정확하게 해석하여 반영하여야 한다. 이 연구에서는 k-omega SST 난류 모형과 자유수면의 급격한 변동을 해석하기 위한 하이브리드-VOF(hybrid volume of fluid)기법을 이용하여 도수현상을 수치적으로 재현하고자 한다. 기존 CFD(computational fluid Dynamics) 모델링에서는 자유수면 변동의 영향을 고려하기 위해 VOF 기법을 많이 사용하였다. 하지면 전통적인 VOF 기법은 다상흐름(multiphase flow)을 오직 부피분율(volume fraction)의 함수로만 고려하며 모의함으로써 강한 수면변동뿐만 아니라 공기연행(air entrainment)를 동반하는 난류 흐름을 모의하는데는 한계가 있다. 이 연구에서 이용하는 Eulerian 기법인 하이브리드 VOF 기법은 물과 공기의 각 상에 대하여 흐름 특성들을 개별적으로 계산하기 때문에 공기연행을 포함한 급변류 흐름에서 전통적인 VOF 기법보다 적용성이 우수하다. 이와 같은 난류모형과 자유수면 포착기법을 이용하여 3차원 비정상 난류 흐름 수치모형을 구축하여 수공구조물 주변에서 발생하는 강한 공기연행과 난류 특성를 보이는 급변류를 수치적으로 재현한다. 이 연구는 계산된 수치해석 결과를 기존 수리실험 결과와 비교하여 수치모형의 적용성을 평가하고 도수 현상에서 발생하는 독특한 흐름 특성을 제시한다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2004.05b
• /
• pp.376-380
• /
• 2004

측벽이 존재하는 개수로 난류흐름에 대한 DNS 자료를 이용하여 난류의 비등방성을 해석하였다. 측벽의 2등분선(sidewall bisector)에서 난류강도의 분포를 통해 바닥과 자유수면에서의 비등방성을 분석하고, AIM을 도입하여 흐름장 전반에 걸쳐 세부적인 비등방성 해석을 수행하였다. 측벽의 2등분선에서의 난류강도의 분포를 통해 바닥과 자유수면 근처에서 난류강도가 강한 비등방성을 갖는 것을 볼 수 있었다. 또한 3개의 다른 영역에서 AIM의 도입을 통해 측벽이나 바닥근처에서는 난류의 비등방성이 구형 관수로 흐름과 유사한 것으로 나타났으나, 개수로 난류흐름의 주된 특성이 velocity-dip phenomena가 존재하는 영역에서는 구형 관수로 흐름과는 달리 천이영역이 존재한다는 것을 알 수 있었다.

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

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

사행수로에서의 원심력은 비선형적인 압력분포를 야기하여 이차류, 편수위 등과 같은 불규칙하고 복잡한 흐름을 발생시킨다. 일반적으로 이들 흐름은 난류이고 매우 3차원적이며 자유수면과의 상호작용이 중요한 역할을 할 수도 있다. 환경, 유사이동, 지형 변화와 관련된 환경 수리학적 관점에서 사행수로에서의 흐름을 이해하고 설계하기 위해서는 이러한 복잡한 3차원 흐름을 정확하게 계산하는 것이 매우 중요하다. 이 연구에서는 유한차분법에 근거한 3차원 흐름해석 모형을 이용하여 사행수로에서의 난류 흐름을 모의하고자 한다. 지배방정식은 3차원 비정상 RANS(Reynolds averaged Navier-Stokes) 방정식이며, 난류 해석을 위해서 공학적으로 널리 이용되고 있는 난류 모형 중 k-omega 모형을 이용한다. 수치모형은 시간과 공간에 대해서 2차 정확도의 이산화 기법을 적용한다. 자유수면의 변동은 이상(two-phase) VOF (volume of fluid) 기법을 이용하여 계산한다. 수치모형의 적용 대상은 기존 문헌에서 제시되어 있는 키노시타 사인곡선을 이용하여 만든 폭 60cm의 사행수로에서 후르드수 0.23 그리고 레이놀즈수 41,700의 조건에서 발생시킨 난류 흐름이다. 적용한 난류모형들을 이용하여 해석한 결과들을 유속벡터분포와 수위의 항으로 비교분석하여 사행수로에서 발생되는 이차류와 편수위 변화 재현에 대한 수치모형의 적용성을 평가하고 각 난류모형들의 특성을 제시한다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.434-438
• /
• 2005

본 연구에서는 개수로 흐름에서 전단속도의 분포를 조사하였다. 이를 위해 직사각형 개수로 흐름과 복단면 개수로 흐름에 대한 DNS 자료를 이용하였으며, 계산된 전단속도는 기존의 실험자료 및 수치해석결과와 비교$\cdot$분석되었다. 평균흐름장에서는 직사각형 개수로 흐름의 경우, 자유수면 부근에서 내측이차흐름이 발생하는 것을 볼 수 있었으며, 복단면 개수로 흐름에서는 주수로와 홍수터의 접합부 부근에서 쌍와가 생성되는 것으로 나타났다. 이차흐름의 최대크기는 체적평균된 주흐름방향 유속의 약 $5\%$로서, 직사각형 개수로 흐름에서 발생되는 것보다 더 강하다는 것을 알 수 있었다. 이러한 이차흐름의 영향으로 주흐름방향 평균유속의 등속선이 편향된 것으로 나타났다. 본 DNS에서 전단속도가 최대값을 갖는 지점은 직사각형 개수로 흐름의 경우 자유수면이며, 복단면 개수로 흐름의 경우 접합부인 것으로 나타났다. 이는 기존 LES 결과와 동일하지만, 실험연구와는 다른 결과이다. 접합부에서 이러한 차이를 보이는 이유는 실험에서 마찰속도를 산정하기 위해 사용하는 방법의 부적합성에 있는 것으로 나타났다.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.18 no.4
• /
• pp.352-359
• /
• 2012

An experiment was carried out to figure out the turbulence flow characteristics in the wake of the transom stern's 2-dimensional section by 2-frame grey level cross correlation PIV method at Re= $3.5{\times}10^3$, Re= $7.0{\times}10^3$. The angles of transom stern are $45^{\circ}$(Model "A"), $90^{\circ}$(Model "B") and $135^{\circ}$(Model "C") respectively. The depth of wetted surface is 40mm from free surface. Strong turbulence intensity appears at the interaction between the flow separation of the bottom of a model and the free surface. This study provides statistic flow information such as turbulence intensity, Reynolds stress and turbulence kinetic energy. Model C type (Raked transom) has low Reynolds stress and turbulence kinetic energy.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.6 s.144
• /
• pp.583-592
• /
• 2005

The turbulent free surface flow around a self-propelled KRISO 138K LNG Carrier is numerically simulated using the finite volume based multi-block RANS code, WAVIS developed at HRISO. The realizable k-$\varepsilon$ turbulence model with a wail function is employed for the turbulence closure. The free surface is captured with the Level-Set method and body forces are used to model the effects of a propeller without resolving the detail blade flow. In order to obtain an accurate free surface solution and stable convergence, the computations are executed with a proper fine grid refinement around the free surface and with an adoption of implicit discretization scheme for the Level-Set formulation. The computed velocity vectors at the several stations and wave patterns show a good agreement with the experimental results measured at the KRISO towing tank.

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

저수지나 하천 사면에서 발생하는 산사태와 토석류는 저수지와 하천 수체에 충격을 가한다. 이로 인해 발생하는 수면 충격파는 전파되어 반대편 제방으로 파의 처오름 또는 댐 제체위로의 물넘이로 큰 피해를 줄 수 있다. 최근 외국에서는 2차원 충격파 생성 및 전파의 기본 과정을 구명하기 위한 실험적 연구가 이뤄지고 있으며, 이들 연구들은 충격파의 발생과 전파, 사면활동 물질과 수체의 상호작용 그리고 자유 수면과 유속분표의 발달에 대한 자세한 관측 자료를 제시하고 있다. 아울러 충격파에 영향을 주는 지배 매개변수를 제시하고 있다. 하지만, 이러한 실험적 연구의 최근 진보에도 불구하고, 이들 지배 매개변수를 고려한 충격파 지배공식들은 대상 지역의 복잡한 바닥 지형이나, 평면적 지형 변화를 단순한 추정치로만 고려하게 된다. 따라서 복잡한 지형조건에서 토석류와 수체의 상호작용과 수면 충격파의 전파를 합리적으로 해석하는 데는 한계가 있다. 이 경우 수치모델링 기법을 대안으로 적용할 수 있으나, 수치모델링은 수면에서 충격파의 전파와 수중에서 토석류의 전파를 동시에 모의해야 하고, 뉴턴 유체와 비뉴턴 유체의 특성을 동시에 고려해야하므로 수치해석 연구자들에게는 하나의 큰 도전사항이다. 이 연구는 경계면 포착기법을 이용한 계산유체동력학 기법을 이용하여 사면활동과 이로 인한 정지 수역에서의 충격파의 발생 및 전파를 재현하기 위한 수치 모델링 기법을 개발하는 것이 목적이다. 사면활동과 수면의 경계면을 포착하고 위치를 정립하기 위해서 VOF (volume of fluid) 경계면 재구축 기법을 이용한다. 지배 방정식은 비압축성(incompressible) 질량 보존방정식과 나비어-스톡스(Navier-Stokes) 방정식이며, 서로 다른 유체의 상(phase)애 대한 체적분할이송방정식을 이용한다. 큰와 모의 계열의 난류 모델링 기법을 적용하여 충격파의 전파와 붕괴에 대한 난류의 영향을 고려하였다. 토석류는 비뉴턴 흐름저항 관계식을 적용하여 그 흐름특성을 재현하였다. 이들 지배방정식은 2차 정확도의 유한체적법(finite volume method)을 이용하여 해석한다. 외국의 연구자들이 관측하여 제시한 길이 11 m 그리고 폭 0.5 m의 수로에서 발생한 충격파를 수치적으로 재현하여 개발된 모형의 실제 문제에 대한 적용성을 보여준다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.30 no.2B
• /
• pp.179-189
• /
• 2010

This paper presents a numerical investigation of vorticity generation in fully vegetated open-channel flows. The Reynolds stress model is used for the turbulence closure. Open-channel flows with rough bed-smooth sidewalls and smooth bed-rough sidewalls are simulated. The computed vectors show that in channel flows with rough bed and rough sidewalls, the free-surface secondary currents become relatively smaller and larger, respectively, compared with that of plain channel flows. Also, open-channel flows over vegetation are simulated. The computed bottom vortex occupies the entire water depth, while the free-surface vortex is reduced. The contours of turbulent anisotropy and Reynolds stress are presented with different density of vegetation. The budget analysis of vorticity equation is carried out to investigate the generation mechanism of secondary currents. The results of the budget analysis show that in plain open-channel flow, the production by anisotropy is important in the vicinity of the wall and free-surface boundaries, and the production by Reynolds stress is important in the region away from the boundaries. However, this rule is not effective in vegetated channel flows. Also, in plain channel flows, the vorticity is generated mainly in the vicinity of the free-surface and the bottom, while in vegetated channel flows, the regions of the bottom and vegetation height are important to generate the vorticity.