• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.566-570
• /
• 2012

본 연구에서는 비선형 k-$\varepsilon$ 모형을 이용하여 직사각형 개수로에서 평균흐름과 난류구조를 모의하였다. 표준 k-$\varepsilon$ 난류모형은 난류의 등방성을 가정하여 국부적 평형상태에서 계산하기 때문에 유선에 따른 레이놀즈 응력의 변형이 큰 경우나 이방성이 강한 경우 이를 계산하지 못한다. 이를 보완하기 위하여 제시된 것이 비선형 k-$\varepsilon$ 난류모형이다. 본 연구에서는 표준 k-$\varepsilon$ 모형과 비선형 k-$\varepsilon$ 모형에 의한 모의결과를 비교하였다. 난류모형을 검증하기 위하여 직사각형 개수로에 흐름을 완전 발달된 등류로 가정하여 해석하였다. 지배방정식을 해석하기 위해 Patankar와 Spalding (1972)이 제시한 SIMPLER 알고리즘을 사용하였고 유한체적법을 이용하여 이산화하고 엇갈린 격자체계를 사용하여 계산에서 발생하는 과도한 진동을 줄였다. 또한 차분기법은 Patankar (1980)가 제시한 Power-law 기법을 채택하였으며 경계조건으로 2층 벽법칙 모형과 Hossain과 Rodi (1993)의 모형을 이용하였다. 두 모형의 적용성을 검증하기 위하여 실측자료를 이용하여 비교하였고 그 결과 비선형 k-$\varepsilon$ 모형이 표준 k-$\varepsilon$ 모형에 비해 좀 더 실측지에 가깝게 모의하는 것을 볼 수 있었다.

• Journal of Korea Water Resources Association
• /
• v.32 no.3
• /
• pp.225-235
• /
• 1999

For a comparative evaluation of three turbulence models in the analyses of thermal discharge behavior into a crossflow, a 2-dimemsional near-field numerical model is developed. The selected models are k-$\varepsilon$ and k-ι turbulence models as a 2-equation turbulence model and a 4-equation turbulence model in which the transport equations for mean of the temperature fluctuation squared and its dissipation rate for the consideration of buoyancy production and turbulent heat flux terms are added to a k-$\varepsilon$ turbulence model. The developed models are applied to a steady flow in an open channel with simple geometry and the numerical results agree with the existing experimental data. Numerical results of buoyancy induced gravitational lateral spreading by 4-equation turbulence model agree with the experimental data better than those of 2-quation turbulence models. The flow patterns by 4 and 2-equation turbulence models are similar.

• Journal of Korea Water Resources Association
• /
• v.44 no.3
• /
• pp.189-198
• /
• 2011

The aim of this study is to compare k-$\varepsilon$ and k-$\omega$ closures under the condition of oscillatory layer flow at high Reynolds number. A one dimensional vertical model incorporated with flow momentum equations and turbulence models (k-$\varepsilon$ and k-$\omega$) is applied to the laboratory measurements in the turbulent oscillatory boundary layer. The numerical simulation reveals that both turbulence models calculate similar velocity profiles and turbulent kinetic energy (TKE). In addition, both deliver high accuracy under the condition of negligible spanwise pressure gradient. Therefore, it is recommended in this study to use k-$\varepsilon$ closure, of which numerical coefficients have been calibrated from many studies, for the cases of straight channel, estuary, and coastal environment where the spanwise pressure gradient is not significant.

• Journal of Korea Water Resources Association
• /
• v.50 no.9
• /
• pp.637-646
• /
• 2017

This study presents a numerical model for simulating dense interflows. The governing equations are provided and the finite difference method is used with the $k-{\varepsilon}$ turbulence model. The model is used to simulate a dense interflow established in a deep ambient water, resulting velocity and excess density profiles. It is observed that velocity decreases in the longitudinal direction due to water entrainment in the vicinity of the outlet and rarely changes for increased Richardson number. Similarity collapses of velocity and excess density are obtained, but those of turbulent kinetic energy and dissipation rate are not. A shape factor for the dense interflow is obtained from the simulated profiles. The value of this shape factor can be used in the layer-averaged modeling of dense interflows. In addition, a buoyancy-related parameter ($c_{3{\varepsilon}}$) for the $k-{\varepsilon}$ model and the volume expansion coefficient (${\beta}_0$) are obtained from the simulated results.

• Journal of Navigation and Port Research
• /
• v.30 no.5 s.111
• /
• pp.369-374
• /
• 2006

The flow analysis was made by applying the turbulent models in the complicated fuel chamber of engine. The $k-\varepsilon,\;k-\omega$, Spalart-Allmaras and reynolds stress models are used in which the hybrid grid is applied for the simulation. The velocity vector, the pressure contour, the change of residual along the iteration number, and the dynamic head are simulated for the comparison of four example cases. Computational results are compared with others. For the code's validation, 2-D bodies were simulated in advance by predicting the drag coefficients.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2008.05a
• /
• pp.159-163
• /
• 2008

본 연구에서는 2-방정식 난류모형을 이용하여 사면을 따라 발달하는 하층밀도류의 시간에 따른 이동현상에 대해 살펴보았다. 이를 위해 타원형의 편미분 방정식을 지배방정식으로 구성하고, 난류 완결을 위해 k-$\varepsilon$ 난류모형을 이용하였다. 개발된 모형을 이용하여 경사의 사면을 따라 진행하는 연속 유입 밀도류를 수치모의 하였다. 완전 발달된 하층밀도류의 거리에 따른 주흐름방향 유속 분포, 체적 농도를 계산하였고, 이를 기존의 실험결과와 비교하였다. 실험과 수치모의 결과가 잘 일치함을 확인하였다. 또한, 불연속 유입 밀도류의 시간에 따른 진행 상황을 수치모의하여 밀도류와 주변수체의 경계부에서 Kelvin-Helmholtz 불안정에 의한 와(渦)가 형성되는 것을 확인하였으며, 밀도류 선단부의 진행 속도와 주변수체의 유입에 대해 고찰하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.545-549
• /
• 2009

본 연구에서는 주흐름방향으로 식생 영역이 교차적으로 존재하는 개수로 흐름에 대한 3차원 수치 모의를 수행하였다. 지배방정식의 난류 폐합을 위해 $k-{\varepsilon}$ 난류모형을 이용하였다. 먼저, 하상의 일부만 식재된 부분 식생 수로를 수치모의 하고 기존의 실험 결과와 비교하였다. 그 결과 본 모형이 평균유속 분포를 매우 잘 예측하는 것으로 나타났으나, 레이놀즈응력 분포는 실험 결과에 비해 비식생영역에서는 다소 과소 산정하고 식생영역에서는 과대 산정하는 것으로 나타났다. 이는 본 모형이 등방성 모형이기 때문에 식생 경계부에서 발생되는 난류의 비등방성 효과를 정확히 예측 할 수 없기 때문인 것으로 판단된다. 또한 주흐름방향으로 식생 영역이 교차적으로 존재하는 대응 식생 수로를 수치모의하고, 계산 결과를 기존의 실험 결과와 비교하였다. 그 결과 본 모형이 대응 식생 수로에서의 유속 분포를 매우 잘 예측하는 것으로 나타났다. 또한 식생 밀도가 증가함에 따라 식생이 흐름 방향을 변화시켜 점차 만곡수로와 유사한 형태의 흐름이 형성되는 것으로 나타났다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.5 no.4
• /
• pp.357-368
• /
• 1993

With the assumption of the diffusion dominated flow, a numerical model has been developed for undertow and turbulence structure under the breaking wave by using the $textsc{k}$-$\varepsilon$ turbulence model. Undertow is a strong mean current which moves seqwards below the level of wave trough in the surf zone. The turbulence, generated by wave breaking in the roller, spreads and dissipates downwards. The governing equations are composed of the equation of motion with the period-averaged shear stress due to waves; $textsc{k}$- and $\varepsilon$-equations with the turbulence energy Production due to wave breaking. They are discretised by the three-level fully implicit scheme, which can be solved by using Thomas algorithm. The model gives good agreements with measurements except for the station that is closest to the breaking point.

• Journal of Korea Water Resources Association
• /
• v.50 no.2
• /
• pp.89-97
• /
• 2017

Turbulent flow structure in the high amplitude meandering channel is complex due to secondary recirculation with helicoidal motions and shear layers formed by flow separation from the curved sidewall. In this work, the secondary flow and the superelevation of the water surface produced in the high-amplitude Kinoshita channel are reproduced by the unsteady Reynolds-averaged Navier-Stokes (RANS) computations using the VOF technique for resolving the variation of water surface elevation and three statistical turbulence models ($k-{\varepsilon}$, RNG $k-{\varepsilon}$, $k-{\omega}$ SST). The numerical results computed by a second-order accurate finite volume method are compared with an existing experimental measurement. Among applied turbulence models, $k-{\omega}$ SST model relatively well predicts overall distribution of the secondary recirculation in the Kinoshita channel, while all three models yield similar prediction of water superelevation transverse slope. The secondary recirculation driven by the radial acceleration in the upstream bend affects the flow structure in the downstream bend, which yields a pair of counter-rotating vortices at the bend apex. This complex flow pattern is reasonably well reproduced by the $k-{\omega}$ SST model. Both $k-{\varepsilon}$ based models fail to predict the clockwise-rotating vortex between a pair of counter-rotating vortices which was observed in the experiment. Regardless of applied turbulence models, the present computations using the VOF method appear to well reproduce the superelevation of water surface through the meandering channel.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.3
• /
• pp.813-820
• /
• 2014

This study presents a numerical modeling of mean flow and turbulence structures of partly-vegetated open-channel flows. For this, Reynolds-averaged Navier-Stokes equations with vegetation drag terms are solved numerically using the non-linear k-${\varepsilon}$ model. The numerical model is applied to laboratory experiments of Nezu and Onitsuka (2001), and simulated results are compared with data from measurement and computations by Kang and Choi's (2006) Reynolds stress model. The simulation results indicate that the proposed numerical model simulates the mean flow well. Twin vortices are found to be generated at the interface between vegetated and non-vegetated zones, where turbulence intensity and Reynolds stress show their maximums. The model simulates the pattern of the Reynolds stress well but under-predicts the intensity of Reynolds stress slightly.