• 한국수자원학회논문집
• /
• 제30권3호
• /
• pp.201-210
• /
• 1997

$textsc{k}$-$\varepsilon$난류모델을 이용한 개수로 급축소부 흐름의 해석을 위한 수치모형을 개발하였고, 그 수치실험결과는 실험결과와 잘 일치하였다. 이는 난류모델에 의한 적절한 난류 와점성계수의 산정이라 생각된다. 유함수 및 유속분포의 분석을 통한 축소부 통과후 급변화 흐름의 수리특성을 분석하였고, 또한 난류 와점성계수의 분포를 분석하였다. 따라서 천수방정식의 점성항에 경험적인 유효점성계수의 도입보다는 흐름의 양상과 장소에 따라서 변화하는 난류와점성계수의 산정에 의한 급변화흐름의 해석이 필요하다.

• 한국전산유체공학회지
• /
• 제20권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.

• 한국안전학회지
• /
• 제32권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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 춘계학술대회
• /
• 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년 춘계학술대회논문집
• /
• 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).

• 해양환경안전학회지
• /
• 제20권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.

• 한국해안해양공학회지
• /
• 제6권3호
• /
• pp.205-215
• /
• 1994

난류모형을 도입하여 풍성류에 대한 3차원 수치모형을 수립하고, 풍성류 계산에 있어서 연직 와점성계수 산정방법들에 대한 비교평가를 수행하였다. 검토된 방법은 함수형(0-방정식), 1-방정식, 2개의 2-방정식 난류모형($ extsc{k}$-$\varepsilon$과 $textsc{k}$-ι)이며, 난류모형은 난류확산 특성길이의 연직분포에 대한 비교를 통하여 각 방법의 특징이 검토되었다. 양단이 막힌 수로와 순환수로에서 유속의 연직분포 계산결과를 수리실험자료와 비교검증하였다. 분석결과에 의하면 1-방정식 난류모형에 적합한 난류확산 특성길이 산정식은 포물선형이었으며, 2-방정식 난류모형($textsc{k}$-ι 모형과 $textsc{k}$-$\varepsilon$모형)은 대체로 실험치와 일치하는 경향을 보였다. 바람에 의한 영향이 전수심에 미치지 않는 경우에 가정된 연직 와점성계수의 분포는 바람이 영향을 미치는 수심까지만 적용되며, 영향이 없는 수심에서는 연직확산이 거의 일어나지 않아 적정성계수의 크기가 0에 가까웠다.

• 해양환경안전학회지
• /
• 제25권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년도 유체기계 연구개발 발표회 논문집
• /
• 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.

• 한국전산유체공학회지
• /
• 제15권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).