• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.9
• /
• pp.1290-1298
• /
• 2003

The stirred tank reactor is one of the most commonly used devices in industry for achieving mixing and reaction. Here we report on results obtained from the large eddy simulations of flow inside the tank performed using a spectral multi-domain technique. The computations were driven by specifying the impeller-induced flow at the blade tip radius. Stereoscopic PlY measurements (Hill et al. $^{(1)}$) along with the theoretical model of the impeller-induced flow (Yoon et al. $^{(2)}$) were used in defining the impeller-induced flow as superposition of circumferential, jet and tip vortex pair components. Large eddy simulation of flow in a stirred tank was carried out for the three different Reynolds numbers of 4000, 16000 and 64000. The effect of different Reynolds numbers is well observed in both instantaneous and time averaged flow fields. The instantaneous and mean vortex structures are identified by plotting an isosurfaces of swirling strength for all Reynolds numbers. The Reynolds number dependency of the non-dimensional eddy viscosity, resolved scale and subgrid scale dissipations is clearly shown in this study.

• Nuclear Engineering and Technology
• /
• v.53 no.11
• /
• pp.3625-3634
• /
• 2021

The crossflow is the key phenomenon in turbulent flow inside rod bundles. In order to establish confidence on application of computational fluid dynamics (CFD) to simulate the crossflow in rod bundles, three Reynolds-Averaged Navier Stokes (RANS) models i.e. the realizable k-ε model, the k-ω SST model and the Reynolds stress model (RSM), and the Large Eddy simulations (LES) with the Wall-Adapting Local Eddy-viscosity (WALE) model are validated based on the Particle Image Velocimetry (PIV) flow measurement experiment in a 5 × 5 rod bundle. In order to investigate effects of periodic boundary condition in the gap, the numerical results obtained with four inner subchannels are compared with that obtained with the whole 5 × 5 rod bundle. The results show that periodic boundaries in the gaps produce strong errors far downstream of the spacer grid, and therefore the full 5 × 5 rod bundle should be simulated. Furthermore, it can be concluded, that the realizable k-ε model can only provide reasonable results very close to the spacer grid, while the other investigated models are in good agreement with the experimental data in the whole downstream flow in the rod bundle. The LES approach shows superiority to the RANS models.

• Journal of Mechanical Science and Technology
• /
• v.17 no.3
• /
• pp.440-448
• /
• 2003

In the present paper, numerical simulations of buoyant plume dispersion in a neutral and stable atmospheric boundary layer have been carride out. A Lagrangian Stochastic Model (LSM) with a Non-Linear Eddy Viscosity Model (NLEVM) for turbulence is used to generate a Reynolds stress field as an input condition of dispersion simulation. A modified plume-rise equation is included in dispersion simulation in order to consider momentum effect in an initial stage of plume rise resulting in an improved prediction by comparing with the experimental data. The LSM is validated by comparing with the prediction of an Eulerian Dispersion Model (EDM) and by the measured results of vertical profiles of mean concentration in the downstream of an elevated source in an atmospheric boundary layer. The LSM predicts accurate results especially in the vicinity of the source where the EDM underestimates the peak concentration by 40% due to inherent limitations of gradient diffusion theory. As a verification study, the LSM simulation of buoyant plume dispersions under a neutral and stable atmospheric condition is compared with a wind-tunnel experiment, which shows good qualitative agreements.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.17 no.6
• /
• pp.861-867
• /
• 2014

In the present study, the flow around a circular cylinder at $Re=3.6{\time}10^6$ is numerically simulated using URANS approach. The objective of this study is to evaluate the turbulence models(Realizable k-${\varepsilon}$, RNG k-${\varepsilon}$) through the prediction of the unsteady flow characteristics around the cylinder. The time-averaged drag coefficients and vortex shedding phenomenon in the wake region are compared to available experimental data and other numerical results. The simulation with Realizable k-${\varepsilon}$ model is found to be more dissipative due to large eddy viscosity predicted in the wake region while the simulation with RNG k-${\varepsilon}$ model predicts a complex vortex shedding phenomenon with more coherent structures realistically.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.3
• /
• pp.256-267
• /
• 1996

In estuarine 3-dimensional numerical modeling. it is very important to calculate vertical eddy viscosity accurately. Various turbulence models employing eddy viscosity concept were applied to the steady flow in an open-channel and the tidal flow in long tidal channel and compared. The evaluations include the verification tests against experimental data sets for steady and tidal flows. The simulation results have shown that the compared models are in good agreements with experimental data of steady flow while only $textsc{k}$-$\varepsilon$ model, $textsc{k}$-ι model, and 1-equation model with well-defined mixing length profile give good agreements with experimental data of tidal flow.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.22 no.5
• /
• pp.685-697
• /
• 1998

Series of recent k-.epsilon. model modification have been carried out with the aid of DNS data to include the effect of near wall. Though these methods opened new way of turbulence modelings, newly developed turbulence models of its kind had yet shortcomings in prediction for the turbulent flows with various Reynolds numbers and various geometric conditions. As a remedy for these shortcomings, a new k-.epsilon. model proposed here by improving the dissipation rate equation and the damping function for eddy viscosity model. The new dissipation rate equation was modeled based on the energy spectrum and magnitude analysis. The damping function for eddy viscosity was also formulated on the ground of distribution of dissipation rate length scales near a wall and the DNS data. The new k-.epsilon. model was applied to the fully developed turbulent flows in a channel and a pipe with a wide range of Reynolds numbers. Prediction results showed that the present model represents properly the turbulence properties in all turbulent regions over a wide range of Reynolds numbers.

• Journal of the Korean Society of Safety
• /
• v.28 no.5
• /
• pp.54-60
• /
• 2013

In this study, a finite element model based on the SU/PG scheme was developed to solve shallow-water equations and the influences of parameters and internal boundary conditions on depth-averaged flow behavior were investigated. To analyze the effect of roughness coefficient and eddy viscosity on flow characteristics, the developed model was applied to rectangular meandering channel with two bends, and transverse velocities and water depth distributions were examined. As the roughness coefficient adjacent to wall increased, the velocities near the wall decreased, and the reduced velocities were compensated by the expanding mid-channel velocities. In addition, the flow characteristics around a circular cylinder were analyzed by varying the internal boundary conditions as free slip and no slip. The assignment of slip condition changed the velocity distribution on the cylinder surface and reduced the magnitude of the shear stress up to one third.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.267-277
• /
• 1993

A finite element method(FEM) is presented and applied to the two-dimensional bottom turbulent boundary layer. The time-dependent incompressible motion of a viscous fluid is formulated by using the well-known Navier-Stokes equations and vorticity equation in terms of the velocity and pressure fields. The general numerical formulation is based on Galerkin method and solved by introducing the mixing length theory of Prandtl for eddy kinematic viscosity of a turbulent flow field. Numerical computations of the transport of sediment on an arbitrary sea-bed due to wave motion in the turbulent boundary layer are carried out. The results obtained by the FEM made clear the difference in characteristic features between the boundary layer due to oscillatory flow and the boundary layer due to wave motion.

• Journal of computational fluids engineering
• /
• v.21 no.3
• /
• pp.24-30
• /
• 2016

In the water purification plant, the mixture of water and chemical from the mixing basin enters the flocculation basin. The rotating flocculators are generally used for the efficient flocculation of dregs. In this paper, the performance of flocculators of a vertical paddle type, widely used in the typical flocculation basins, and a hydro-foil type, recently disseminated in the field, are compared with each other by use of the numerical method. Also the characteristics and the efficiency are analyzed with CFD techniques. The strain rate and the eddy viscosity are compared for two types to predict the mixing efficiency, and the maximum speed and its location are pursued from the computed data. The hydrofoil type shows that the eddy viscosity is enhanced 1.66 to 3.03 times larger than that of vertical paddle type, and also produced 1.87 to 1.95 times larger flocs for each stage. However, the rapid rotation of hydrofoil may chop the floc to small size due to the higher turbulence intensity. From the result of computation, the strong and weak points of each type have been analyzed for the decision making.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2011.05a
• /
• pp.132-136
• /
• 2011

수제는 유수에 의한 하안 침식을 방지하기 위한 목적의 수공구조물이다. 최근에는 하천복원 및 생태계복원에 대한 관심이 증가되면서, 하안 침식 방지뿐만 아니라 수제에 의한 흐름분리 및 재순환영역의 생태 서식처 제공에 대한 역할이 부각되고 있다. 따라서 수제 설계 시, 수제 설치에 따른 수리특성 변화를 미리 파악하는 것이 필요하지만, 수제 간격, 설치 방향, 수제 높이, 수제 길이, 하도 특성 및 접근 유수 특성 등 다양한 인자에 의해 수리특성 변화를 예측하기란 쉽지 않으며, 비용과 시 공간 측면에서 물리 모형을 이용한 접근법도 용이하지 못하다. 따라서 본 연구에서는 2차원 수치모형을 이용하여 수제 설치에 따른 흐름특성 변화를 모의하고, 모형의 적용성을 판별하고자 하였다. 지배 방정식은 기본적으로 천수방정식을 적용하게 되는데, 수제 주변 흐름 특성은 난류에 의한 횡방향 확산이 중요한 인자로작용하게 된다. 상수 와점성(constant eddy viscosity) 모형, 포물선형 와점성 (parabolic eddy viscosity) 모형을통해 개발 모형의예측 결과와 실험수로 결과를 비교하였다. 수제하류의 재순환영역의 길이가 실험수로의 결과에 비해서 5~6 배까지 확대된 결과를 나타내고 있어, 만족할 만한 결과를 현재까지 얻지 못하였다. 횡방향 확산이 중요한 만큼 k-$\varepsilon$ 모형 등 다양한 난류 모형의 적용이 고려되어질 필요가 있으며, 또한 1차원 요소가 아닌 2차원 요소의 적용을 고려할 필요도 있다.