• Journal of Korea Water Resources Association
• /
• v.30 no.3
• /
• pp.201-210
• /
• 1997

To analyze the flow through open-channel constrictions using $\kappa$-$\varepsilon$ turbulence mode, a numerical model is developed. The simulated results agree well with existing experimental data which attributes to the adequate input of turbulent eddy-viscosity by turbulence model. A stream function and velocity distributions enable the analysis of flow characteristics at the downstream of constriction. Turbulent eddy viscosities over channel are spatially varied with stream pattern. For the evaluation of rapidly varied flow, the eddy-viscosity input by turbulence model is required instead of the empirical effective viscosity to solve a shallow water equation.

• Journal of computational fluids engineering
• /
• v.20 no.3
• /
• pp.94-103
• /
• 2015

The turbulent flow characteristics in the channel flow are investigated using large eddy simulation(LES) of FDS code, built in NIST(USA), in which the near-wall flow is solved by Werner-Wengle wall function. The periodic flow condition is applied in streamwise direction to get the fully developed turbulent flow and symmetric condition is applied in lateral direction. The height of the channel is H=1m, and the length of the channel is 6H, and the lateral length is H. The total grid is $32{\times}32{\times}32$ and $y^+$ is kept above 11 to fulfill the near-wall flow requirement. The Smagorinsky model is used to solve the sub-grid scale stress. Smagorinsky constant $C_s$ is 0.2(default in FDS). Three cases of Reynolds number(10,700, 26,000, 49,000.), based on the channel height, are analyzed. The simulated results are compared with direct numerical simulation(DNS) and particle image velocimetry(PIV) experimental data. The linear low-Re eddy viscosity model of Launder & Sharma and non-linear low-Re eddy viscosity model of Abe-Jang-Leschziner are utilized to compare the results with LES of FDS. Reynolds normal stresses, Reynolds shear stresses, turbulent kinetic energys and mean velocity flows are well compared with DNS and PIV data.

• Journal of the Korean Society of Safety
• /
• v.32 no.2
• /
• pp.112-116
• /
• 2017

Present study investigated the influence of time step size, turbulent eddy viscosity, and the number of layer on rapid and unsteady propagation of dam break flow. When the time step size had a value such that it resulted in Cr of 0.89, a significant numerical oscillation was observed in the vicinity of the wave front. Higher turbulent viscosity ensured smooth and mild slope of velocity and water stage compared with the flow behavior by no viscosity. The vertical velocity at the lower layer positioned near the bottom showed lower velocity compared with other layers.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1806-1813
• /
• 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 PIV 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 nondimeansional eddy viscosity, resolve scale and subgrid scale dissipations is clearly shown in this study.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.407-413
• /
• 2011

VIC (Vortex-In-Cell) method for viscous incompressible flow is presented to simulate the wake behind a modified NACA16 foil. With uniform rectangular grid, the velocity in field is calculated using streamfunction from vorticity field by solving the Poisson equation in which FFT(Fast Fourier Transform) is combined with 2nd order finite difference scheme. Here, LES(Large Eddy Simulation) with Smagorinsky model is applied for turbulence calculation. Effective viscosity is formulated using magnitude of strain tensor(or vorticity). Then the turbulent diffusion as well as viscous diffusion becomes particle strength exchange(PSE) with averaged eddy viscosity. The well-established panel method is combined to obtain the irrotational velocity and to apply the no-penetration boundary condition on the body panel. And wall diffusion is used for no-slip condition numerical results of turbulent stresses are compared with experimental results (Bourgoyne, 2003). Before comparing process, LES(Large Eddy Simulation) SGS(Subgrid scale) stress is transformed Reynolds averaged stress (Winckelmans, 2001).

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.20 no.3
• /
• pp.323-333
• /
• 2014

In this paper, the wakes behind a square cylinder were simulated using two kinds of different turbulence models for the eddy viscosity concept such as the zero- and the one-equation model in which the former is the mixing length model and the latter is the k-equation model. For comparison between numerical and analytical solutions, we employed three skill assessments: the correlation coefficient(r) for the similarity of the wake shape, the error of maximum velocity difference(EMVD) for the accuracy of wake velocity and the ratio of drag coefficient(RDC) for the pressure distribution around the structure. On the basis of the numerical results, the feasibility of each model for wake simulation was discussed and a suitable value for the empirical constant was suggested in these turbulence models. The zero-equation model, known as the simplest turbulence model, overestimated the EMVD and its absolute mean error(AME) for r, EMVD and RDC was ranging from 20.3 % to 56.3 % for all test. But the AME by the one-equation model was ranging from 3.4 % to 19.9 %. The predicted values of the one-equation model substantially agreed with the analytical solutions at the empirical mixing length scale $L=0.6b_{1/2}$ with the AME of 3.4 %. Therefore it was concluded that the one-equation model was suitable for the wake simulation behind a square cylinder when the empirical constant for eddy viscosity would be properly chosen.

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

A 3-dimensional numerical model of wind-induced flows has been established. and comparative evaluation of determination methods of vertical eddy viscosity has been performed. The model uses turbulence models to calculate vertical eddy viscosity. The examined methods arp 0-equation model of functional form, 1-equation model of turbulence kinetic energy, and two 2-equation models ($textsc{k}$-$\varepsilon$ and $textsc{k}$-ι models). The evaluation includes the verification tests against experimental data for wind-driven current On a closed one-dimensional channel and a recirculating one-dimensional channel. Comparative study of turbulence models has shown that the proper distribution of turbulence scale is parabolic and the eddy viscosity is depending strongly on mixing depth due to wind.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.25 no.3
• /
• pp.344-353
• /
• 2019

The wakes behind a square cylinder were simulated using two-equation turbulence models, $k-{\varepsilon}$ and RNG $k-{\varepsilon}$ models. For comparisons between the model predictions and analytical solutions, we employed three skill assessments:, the correlation coefficient for the similarity of the wake shape, the error of maximum velocity difference (EMVD) of the accuracy of wake velocity, and the ratio of drag coefficient (RDC) for the flow patterns as in the authors' previous study. On the basis of the calculated results, we discussed the feasibility of each model for wake simulation and suggested a suitable value for an eddy viscosity related constant in each turbulence model. The $k-{\varepsilon}$ model underestimated the drag coefficient by over 40 %, and its performance was worse than that in the previous study with one-equation and mixing length models, resulting from the empirical constants in the ${\varepsilon}-equation$. In the RNG $k-{\varepsilon}$ model experiments, when an eddy viscosity related constant was six times higher than the suggested value, the model results were yielded good predictions compared with the analytical solutions. Then, the values of EMVD and RDC were 3.8 % and 3.2 %, respectively. The results of the turbulence model simulations indicated that the RNG $k-{\varepsilon}$ model results successfully represented wakes behind the square cylinder, and the mean error for all skill assessments was less than 4 %.

• 유체기계공업학회:학술대회논문집
• /
• 2000.12a
• /
• pp.337-344
• /
• 2000

A modified k-$\epsilon$ model is proposed for calculation of transitional boundary layer flows. In order to develop the eddy viscosity model for the problem, the flow is divided into three regions; namely, pre-transition region, transition region and fully turbulent region. The pre-transition eddy-viscosity is formulated by extending the mixing Length concept. In the transition region, the eddy-viscosity model employs two length scales, i.e., pre-transition length scale and turbulent length scale pertaining to the regions upstream and the downstream, respectively, and a university model of stream-wise intermittency variation is used as a function bridging the pre-transition region and the fully turbulent region. The proposed model is applied to calculate three benchmark cases of the transitional boundary layer flows with different free-stream turbulent intensity ( $1\%{\~}6\%$ ) under zero-pressure gradient. It was found that the profiles of mom velocity and turbulent intensity, local maximum of velocity fluctuations, their locations as well as the stream-wise variation of integral properties such as skin friction, shape factor and maximum velocity fluctuations are very satisfactorily Predicted throughout the flow regions.

• Journal of computational fluids engineering
• /
• v.15 no.3
• /
• pp.39-45
• /
• 2010

The turbulent flows in a tunnel mock-up($10L{\times}0.5W{\times}0.25H$ m3 : scale reduction 1/20) with rectangular cross section were investigated. The instantaneous velocity fields of Re = 49,029, 89,571 were measured by the 2-D PIV system which is consisted of double pulsed Nd:Yag laser and the tracer particles in the straight-duct mock-up where the flows were fully developed. The mean velocity profiles were taken from the ensemble averages of 1,000 instantaneous velocity fields. Simultaneously, numerical simulations(RANS) were performed to compare with experimental data using STREAM code. Non-linear eddy viscosity model (NLEVM : Abe-Jang-Leschziner Eddy Viscosity Model) was employed to resolve the turbulent flows in the duct. The calculated mean velocity profiles were well compared with PIV results. In the log-law profiles, the experimental data were in good agreement with numerical simulations all the way to the wake region except the viscous sub-layer (near wall region).