• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.1
• /
• pp.304-320
• /
• 1996

A multiple production .epsilon. equation model was developed in the low Reynolds number $\kappa$-$\varepsilon$ model with the aids of DNS data. We derived the model theoretically and avoided the use of empirical correlations as much as possible in order for the model to have generality in the prediction of complex turbulent flow. Unavoidable model constants were, however, optimized with the aids of DNS data. All the production and dissipation models in the $\varepsilon$ equation were modified with damping functions to satisfy the wall limiting behavior. A new $f_{\mu}$ function, turbulent diffusion and pressure diffusion model for the k and .epsilon. equations were also proposed to satisfy the wall limiting behavior. By, computational investigation on the plane channel flows, we found that the multiple production model for .epsilon. equation could improve the near wall turbulence behavior compared with the standard production model without the complicated empirical modification. Satisfication of the wall limiting conditions for each turbulence model term was found to be most important for the accurate prediction of near wall turbulence behaviors.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3B
• /
• pp.247-257
• /
• 2009

In the present paper, turbulent open-channel flows with alternate vegetated zones are numerically simulated using threedimensional model. The Reynolds-averaged Navier-Stokes Equations are solved with the ${\kappa}-{\varepsilon}$ model. The CFD code developed by Olsen(2004) is used for the present study. For model validation, the partly vegetated channel flows are simulated, and the computed depth-averaged mean velocity and Reynolds stress are compared with measured data in the literature. Comparisons reveal that the present model successfully predicts the mean flow and turbulent structures in vegetated open-channel. However, it is found that the ${\kappa}-{\varepsilon}$ model cannot accurately predict the momentum transfer at the interface between the vegetated zone and the non-vegetated zone. It is because the ${\kappa}-{\varepsilon}$ model is the isotropic turbulence model. Next, the open channel flows with alternate vegetated zones are simulated. The computed mean velocities are compared well with the previously reported measured data. Good agreement between the simulated results and the experimental data was found. Also, the turbulent flows are computed for different densities of vegetation. It is found that the vegetation curves the flow and the meandering flow pattern becomes more obvious with increasing vegetation density. When the vegetation density is 9.97%, the recirculation flows occur at the locations opposite to the vegetation zones. The impacts of vegetation on the flow velocity and the water surface elevation are also investigated.

• Water for future
• /
• v.29 no.2
• /
• pp.221-233
• /
• 1996

To investigate the hydrodynamic characteristics of a two-dimensional buoyant surface jet, the most important factors of the numerical analysis are the evaluation of the free surface and the turbulence transportation under the stratification. In present study, a numeriacal simulation model used with the semi-implicit method for pressure-linked equations (SIMPLE), the non-hydrostatic approximation and the algebraic stress model (ASM) is applied to investigate the vertical structure of internal flow hydrodynamically. The ASM enables to take account of anisotropy of turbulence, the damping effects of the density interface, and the free surface on the turbulence structure accurately. The ASM tested produces better agreement than the $\kappa-\varepsilon$ model with measurements by Nakatsuji (1984) on the flow development and turbulence structure. Applicability of the ASM to a two-dimensional buoyant surface jet is examined through comparison with experimental data.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3B
• /
• pp.231-237
• /
• 2009

This study presents a numerical model to simulate density currents developing two dimensionally. The ${\kappa}-{\varepsilon}$ model is used for the turbulence closure. Elliptic flow equations are solved by the finite volume method. In order to investigate the applicability of the numerical model, discontinuous density currents are simulated numerically. The vortices due to the instability at the interface are simulated, showing a good agreement with the experimental visualizations in the literature. It is also investigated that the transition from slumping phase to inertial phase occurs when a bore generated at the end wall overtakes the front. However, the propagation of the density current is retarded compared with the experimental results. Two-dimensional modeling seems to have an effect on underestimating the front velocity of the density current.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.5
• /
• pp.1052-1060
• /
• 1989

Using a vortex stretching invariant term, the two-layer k-.epsilon. model has been modified to account for the extra staining of turbulence due to the mean-flow convergence and divergence. The calculations of turbulent boundary layers in a plane of symmetry are compared for experimental cases which are an axisymmetric body at an incidence of 15.deg.. The comparisons between the calculations and experimental data show that additional modifications to the dissipation rate equation have brought the significant improvement to the prediction of plane of symmetry boundary layers in the strong mean-flow convergence and divergence.

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

바닥이 언덕과 저면으로 이루어진 연속적인 횡방향 바닥형상을 갖는 개수로 흐름의 부유사 농도분포를 수치모의 하였다. 비선형 ${\kappa}-{\varepsilon}$ 모형을 이용하여 흐름의 지배방정식인 곡선 직교좌표계에 대한 Navier-Stokes 방정식을 해석하였으며, 와점성계수 개념을 이용하여 부유사 수송 방정식을 해석하였다. 기존의 실험결과와 비교하여 모형이 격자형 이차흐름을 비교적 잘 예측하는 것을 확인하였으며, 동일한 흐름에 대하여 부유사 농도 분포를 계산하였다. 부유사량의 계산 결과 언덕 위의 부유사 농도분포가 저면 위의 농도분포보다 크게 나타났다.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.19-26
• /
• 1991

The complete, fully-elliptic Reynolds-averaged Navier-Stokes equations have been solved using a two-layer model, in the $\kappa-\varepsilon$ turbulence model, for the axisymmetric body. Numerically generated boundary-fitted coordinate system and the finite analytic methods are used to solve the governing equations. Calculations are started after the middle body with given inlet conditions. The velocities and the turbulent quantities at the inlet section are specified by solving the boundary layer equations or by standard flat-plate boundary profiles. The effects of the inlet conditions on the solution are investigated.

• Journal of Korea Water Resources Association
• /
• v.37 no.6
• /
• pp.461-466
• /
• 2004

This study presents a combined experimental and numerical effort to investigate wave breaking of sinusoidal waves in a surf zone. Numerical predictions are verified by comparing to laboratory measurements. The model solves the Reynolds equations and$textsc{k}$-$\varepsilon$ models for the turbulence analysis. To track the free surface displacement, the volume of fluid method is employed. As the height of incident wave increases, the wave breaking occurs at a closer point of the slope in the numerical model and laboratory experiments with the same depth and period. When a wave breaking occurs, the ratio of wave height becomes larger, with the same wave height and depth, as the period increases.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.6 no.4
• /
• pp.325-336
• /
• 1994

A numerical simulation of velocity and temperature fields and Nusselt number distributions is performed by using the algebraic stress model (ASM) for the velocity profiles and low Reynolds number ${\kappa}-{\varepsilon}$ model and the algebraic heat flux model(AHFM) for turbulent heat transfer in a $180^{\circ}$ bend with a constant wall heat flux. In the low Reynolds number ${\kappa}-{\varepsilon}$ model, turbulent Prandtl number is modified by considering the streamline curvature effect and the non-equilibrium effect between turbulent kinetic energy production and dissipation rate. Every heat flux term presented in the transport equation of turbulent heat flux is reduced to algebraic expressions in a way similar to algebraic stress model. Also. in the wall region, low Reynods number algebraic heat flux model(AHFM) is applied.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.10 no.4
• /
• pp.151-158
• /
• 1990

Various factors may contribute on the mixing processes in the surf zone formed by irregular waves. The turbulence motion driven by wave breaking may be one of the major causes, the effect due to spatial variation on current velocity be a secondary one, and the additional process may result from the irregular superposition of radiation stresses or wave breaking dissipation incurred by random breaking waves in a broadened surf zone. In the present study a numerical model of spectral waves and induced currents was developed using a superposition technique with ${\kappa}-{\varepsilon}$ closure for mixing process and applied to a field situation of longshore current generated by spectral waves on a uniform beach. It was found from the application that the surf-zone mixing processes formed by irregular waves can be well described by using ${\kappa}-{\varepsilon}$ equations if the source of ${\kappa}$ is properly represented. The nonlinear energy transfer was also found to have some influence on the velocity profile of longshore current particularly in very shallow water region near coast.