• 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 Korea Water Resources Association
• /
• v.53 no.1
• /
• pp.1-13
• /
• 2020

The most common approach for computing engineering flow problems at high Reynolds number is still the Reynolds-averaged Navier-Stokes (RANS) computations based on turbulence models with wall functions. The recently developed generalized wall functions blending between the wall-limiting viscous and the outer logarithmic relations ensure a smooth transition of flow quantities across two regions. The performances and convergence properties of widely used turbulence models with wall functions that are applicable for turbulence kinetic energy (TKE), turbulent and specific dissipation rates, and eddy viscosity are presented through a series of near wall flow simulations. The present results show that RNG k-𝜖 model should be carefully applied with small tolerance to get the stable solution when the first grid lies in the buffer layer. The standard k-𝜖 and RNG k-𝜖 models are not sensitive to the selection of wall functions for both TKE and eddy viscosity, while the k-ω SST model should be applied together with kL-wall function for TKE and nutUB-wall functions for eddy viscosity to ensure accurate and stable boundary conditions. The applications to a backward-facing step flow at Re=155,000 reveal that the reattachment length is reasonably well predicted on appropriately refined mesh by all turbulence models, except the standard k-𝜖 model which about 13% underestimates the reattachment length regardless of the grid refinement.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.19 no.2
• /
• pp.162-169
• /
• 2007

This study is focused on the comparison of a 3 dimensional numerical and hydraulic model experiment for the flow phenomenon when a lock gate is opened. The lock gate is designed to discharge the flood flow rate at $218m^3/s$ of Solicheon at the Kun Jang national industry complex. The three dimensional ${\kappa}-{\epsilon}$ turbulent model of ANSYS CFX-10 of the computational fluid dynamics(CFD) program was used. The characteristics of CFX-10 are able to be simulated effectively for turbulent flow, especially the flow separation of the boundary layer of the two phase interface of air and water. The velocity and the flow pattern of the numerical model was showed to be similar to the results of the hydraulic model experiment.

• Journal of Ocean Engineering and Technology
• /
• v.10 no.3
• /
• pp.44-49
• /
• 1996

후방계단(backward facting step) 주위의 난류 유동 특성을 수치 해석을 통해 파악하고자 하였다. 지배방정식은 2차 정도의 유한 차분 기법으로 이산화하였으며 비교차격자계를 사용하여 양해법으로 계산하였다. 난류 모형으로는 이층 모형(two-layer)을 사용하였고 압력 Poisson 방정식을 이용하여 압력과 속도를 연성 시켰다. Re=44,000인 경우에 대해 계산 결과로 부터 후방 계단 뒤의 속도 벡터, 유선, 압력 및 속도 분포, 재부착 길이(reattachment length)등의 실험치와 비교하였다. 본 계산에 사용한 수치 해석 기법은 박리등이 포함된 복잡한 난류 유동 현상을 잘 재현할 수 있음을 확인할 수 있었다.

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

The non-linear k-.epsilon. model developed by Speziale was employed for the prediction of developing turbulent flow in a square duct. The numerical procedure incorporated a finite volume method using a strong conservation form of the partially-parabolized Navier-Stokes equation. Results of the calculation were compared with available experimental data on the mean velocity field and turbulent kinetic energy, and was found to be in favorable agreement.

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

본 연구에서는 레이놀즈응력모형(RSM: Reynolds Stress Model)을 이용하여 부분 식생된 개수로 흐름을 수치모의 하였다. 부분 식생된 개수로 흐름에서의 평균유속 및 난류구조를 수치모의 하고 기존의 실험결과와 비교하였다. 그 결과 개발된 모형이 식생된 개수로 흐름을 매우 잘 예측하는 것으로 나타났다. 특히, 이차흐름 벡터도를 수치모의 한 결과 식생구간과 비식생 구간에서 방향이 서로 다른 새로운 이차흐름 구조가 형성되는 것으로 나타났다. 또한 주흐름방향으로의 최대유속이 비식생 영역의 수면 아래에서 발생되고, 식생 및 비식생 영역의 경계면에서 난류량이 최대 값을 갖는 것을 확인하였다.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.1
• /
• pp.15-23
• /
• 1999

K-$\omega$ model is used for simulation of flowfield around the cylinderical afterbody. In addition to two-equation turbulence model, modification terms for the compressibility effects are applied to the simulation. Although the estimations of the skin friction and the surface pressure distribution at hypersonic ramp flowfield were satisfactory, the result of the simulation with the modifications for this flowfield is worse than that of the original K-$\omega$ model. The compressiblility modification terms do negative effects on the estimation. The basic research on the turbulence model for the compressible flowfield has to be further conducted.

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

하천에서 유사이동은 하천환경과 하천형상을 결정하는 주요 요소이므로 이를 해석하는 것은 매우 중요하다. 그러나 유사이동은 일반적으로 이상흐름 (two-phase flow)이며 난류를 동반하기에 이를 해석하기에는 쉽지 않다. 이상흐름을 해석하는 방법으로는 유사를 연속상인 유사구름(sediment cloud)으로 표현하여 해석하는 Euler-Euler 모형이 있으며 입자를 직접 추적하여 해석하는 Euler-Lagrange 모형이 있다. 본 연구에서는 유사이동 해석을 위하여 Euler-Lagrange 모형을 사용하였으며 흐름의 진동성분을 고려하기 위하여 EIM (Eddy Interaction Model)을 사용하였다. 유체의 유속은 Dou (1987)가 제시한 경험식을 사용하였고 난류운동에너지와 소산률은 Nezu and Nakagawa (1993)가 제시한 식을 사용하였다. EIM에서 입자에 발생하는 와의 영향시간(eddy interaction time)을 계산하기 위해 Gosman and Ioannides (1983)가 제시한 eddy lifetime과 eddy crossing time을 사용하였다. 유사입자는 입자의 운동량방정식을 풀어 그 거동을 추적하였으며 일정 시간 후 입자의 수를 이용하여 농도를 계산하였다. 유체에 발생하는 유속의 진동성분에 의해 입자가 부상하고 중력에 의해 흐름에 따른 일정한 농도분포 형태를 가지는 것을 확인하였다. 유사의 입자크기와 흐름에 따른 농도분포를 계산하였으며, 이를 측정치와 비교하여 EIM의 적용성을 확인하였다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.3
• /
• pp.234-244
• /
• 1994

Longshore currents driven by monochromatic waves have been described using 2-equation $textsc{k}$-ι turbulence transport model. When using $textsc{k}$-ι closure both profiles of eddy viscosity and current velocity are found to be satisfactory. Several terms of ι equation are related to various variables concerned with turbulence mechanism. New form of turbulence frequency used in ι equation is suggested in the present approach, and non-dimensional parameters are evaluated by comparing the computational results with the laboratory measurements. Various values of a large range are applied to the non-dimensional parameters for the sensitivity test and in order to improve the predictability common values of constants are deduced, which produce similarly good computational results for the well-controlled laboratory measurements.

• Water for future
• /
• v.29 no.4
• /
• pp.187-198
• /
• 1996

A numerical model for unsteady dispersion of horizontal line source in turbulent shear flow is developed. A fractional step finite difference method is used which splits the unsteady two-dimensional advective diffusion equation into the longitudinal advection and the vertical diffusion equations, and solves them alternately for half time intervals by the Holly-Preissmann scheme and the Crank-Nicholson scheme, respectively. The developed numerical model is verified using a semi-analytic solution for steady dispersion in turbulent shear flow. Dispersion of an instantaneous plane source in turbulent shear flow is analyzed using the model. The degree of mixing at the same dimensionless time is almost the same regardless of the friction factor, and the travel distance required to reach a certain degree of mixing is inversely proportional to the square root of the friction factor.