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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.8
• /
• pp.1046-1055
• /
• 1997

An anisotropic k-.epsilon. turbulence model for predicting the rotating flows is proposed with the simple inclusion of a new parameter dealing with the extra straining effects in the .epsilon.-equation. This model is employed to compute the effects of Coriolis forces on fully-developed flow in a rotating channel. The predicted results indicate that the present model captures fairly well the striking rotational-induced effects on the Reynolds stresses and the mean flow distributions, including the argumentation of turbulent transport on the unstable side (pressure surface) of the channel and its damping on the stable side (suction surface).

• Journal of Korea Water Resources Association
• /
• v.33 no.5
• /
• pp.581-592
• /
• 2000

Turbulence structure and suspended sediment transport capacity in vegetated open-channel flows are investigated numerically in the present paper. The $\textsc{k}-\;\varepsilon$ model is employed for the turbulence closure. Mean velocity and turbulence characteristics including turbulence intensity, Reynolds stress, and production and dissipation of turbulence kinetic energy are evaluated and compared with measurement data available in the literature. The numerical results show that mean velocity is diminished due to the drag provided by vegetation, which results in the reduction of turbulence intensity and Reynolds stress. For submerged vegetation, the shear at the top of vegetation dominates turbulence production, and the turbulence production within vegetation is characterized by wakes. For emergent condition, it is observed that the turbulence generation is dominated by wakes within vegetation. In general, simulated profiles compares favorably to measured data. Computed values of eddy viscosity are used to solve the conservation equation for suspended sediment, yielding sediment concentration more uniform over the depth compared with the one in the plain channel. The simulation reveals that the suspended load decreases as the vegetation density increases and the suspended load increases as the particle diameter decreases for the same vegetation density.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.34 no.E
• /
• pp.30-44
• /
• 1992

A modification of the TEACH-like computer program based on the k-$\varepsilon$ turbulence transport was applied for predicting air mixing patterns and temperature distributions in a rectangular, slot-ventilated enclosure having obstructions ; a rectangular obstacle with heat flux, solid walls separates the passage and the pig pens, and purlins beneath the ceiling. Air flow patterns were calculated for the cases with and without the purlin, extending 300mm beneath the ceiling. Comparisons of prediction data of Randall & Battams(1976) showed air flow pattern predicted well for the case without the purlin. Heat was accumulated at the corner of the left side of the solid wall and the right-upper region of the simulated pigs. However the air distribution pattern was completely different from data for the case with the purlin. The deviation from the observation may be attributed to the difference of the geometric configuation. Exploring the cause of the deviation should be conducted in a further study. Temperature stratification was also observed due to incomplete mixing. The obstruction in the route of the inlet air jet at inlet should be avoided since most of kinetic energy dissipates at the abstacle duet to impingement.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.10a
• /
• pp.126-131
• /
• 2004

A numerical code for cavitation is developed based on pressure-based algorithm. The k-\varepsilon$ model (with wall function) is used for turbulence, and volume transport equation is used for cavitation model. The compressibility is not considered for the flow field is low speed.

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

• Water Engineering Research
• /
• v.3 no.2
• /
• pp.75-84
• /
• 2002

An Euler-Lagrange two-phase flow model is presented fur simulation sheet-flow transport under wave and current. The flow is computed by solving the Reynolds Averaged Navier-Stokes equation in conjunction with the k-$\varepsilon$ turbulence model for turbulence closure. The sediment transport is introduced as a motion of granular media under the action of unsteady flow from the Lagragian point of view. In other word, motion of every single particle is numerically traced with Movable Bed Simulator (MBS) code based on the Distinct Element Method (DEM), in which the frequent interparticle collision of the moving particles during the sheet-flow transport is sophisticatedly taken into account. The particle diameter effect on time-dependent developing process of sheet-flow transport is investigated, by using three different diameter sizes of sediment. The influence of an imposed current on oscillatory sheet-flow transport is also investigated. It is concluded that the sediment transport rate increases due to the relaxation process related to the time-lag between flow velocity and sediment motion.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.3 no.2
• /
• pp.81-91
• /
• 1991

A k-$\varepsilon$ equation model was established to investigate the behaviours of two-dimensional surface buoyant jets. Its computational results were compared with experimental data on the mean flow and the turbulent transport. The model was proved to predict the flow characteristics reasonably. The influence of the values of k and $\varepsilon$ given in the inlet on the evaluation of surface buoyant jets was examined to determine them quantitatively. Computations for several values of buoyancy production coefficient $C\varepsilon$$_3$ in the $\varepsilon$ equation, which has been neglected by many researchers. were carried out to evaluate its effect on the flow development. Computational results of the two-dimensional surface buoyant jets were presented and briefly discussed.

• Wind and Structures
• /
• v.17 no.1
• /
• pp.87-105
• /
• 2013

Modeling an equilibrium atmospheric boundary layer (ABL) in an empty computational domain has routinely been performed with the k-${\varepsilon}$ turbulence model. However, the research objects of structural wind engineering are bluff bodies, and the SST k-${\omega}$ turbulence model is more widely used in the numerical simulation of flow around bluff bodies than the k-${\varepsilon}$ turbulence model. Therefore, to simulate an equilibrium ABL based on the SST k-${\omega}$ turbulence model, the inlet profiles of the mean wind speed U, turbulence kinetic energy k, and specific dissipation rate ${\omega}$ are proposed, and the source terms for the U, k and ${\omega}$ are derived by satisfying their corresponding transport equations. Based on the proposed inlet profiles, numerical comparative studies with and without considering the source terms are carried out in an empty computational domain, and an actual numerical simulation with a trapezoidal hill is further conducted. It shows that when the source terms are considered, the profiles of U, k and ${\omega}$ are all maintained well along the empty computational domain and the accuracy of the actual numerical simulation is greatly improved. The present study could provide a new methodology for modeling the equilibrium ABL problem and for further CFD simulations with practical value.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.10
• /
• pp.1940-1954
• /
• 1992

Fine grid computations were attempted to analyze the turbulent flows in the near wall low Reynolds number region and the numerical analyses were incorporated by a finite-volume discretization with full find grid system and low Reynolds number k-.epsilon. model was employed in this region. For the improvement of low Reynolds number k-.epsilon. model, modification coefficient of eddy viscosity $f_{\mu}$ was derived as a function of turbulent Reynolds number $R_{+}$ and nondimensional length $y^{+}$ from the concept of two length scales of dissipation rate of turbulent kinetic energy. The modification coefficient $f_{\epsilon}$ in .epsilon. transport equation was also derived theoretically. In the turbulent kinetic energy equation, pressure diffusion term was added in order to consider low Reynolds number region effect. The main characteristics of this low Reynolds number k-.epsilon. model were founded as : (1) In high Reynolds number region, the present model has limiting behavior which approaches to the high Reynolds number model. (2) Present low Reynolds number k-.epsilon. model dose not need additional empirical constants for the transport equations of turbulent kinetic energy and dissipation of turbulent kinetic energy in order to consider wall effect. Present low Reynolds number turbulence model was tested in the pipe flow and obtained improved results in velocity profiles and Reynolds stress distributions compared with those from other k-.epsilon. models.s.s.