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

• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 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.

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

• 한국전산유체공학회지
• /
• 제13권3호
• /
• pp.1-7
• /
• 2008

A supersonic base flow is computed to investigate the effect of the eddy viscosity coefficient to the linear ${\kappa}-{\varepsilon}$ turbulence models. Slight modifications to the eddy viscosity coefficient, which are based on the realizability condition, are given to the Launder-Sharma turbulence model so that present models satisfy the realizability condition. Numerical results for supersonic base flow show that turbulence models with the weaky-nonlinear eddy viscosity coefficient can lead to reasonable enhancements in the prediction of the velocity and turbulent kinetic energy profiles.

• 에너지공학
• /
• 제8권2호
• /
• pp.293-297
• /
• 1999

Direct comparisons are made between curvature-corrected eddy viscosity models and the present experimental data. The results show that the curvature effects can be quantified through a curvature parameter R$\sub$c/ or S$\sub$c/ and a non-equilibrium value of p/$\varepsilon$. The data reveal a significant dependence of the eddy viscosity on the curvature and strain history for a fluid in a stabilizing curvature field, S$\sub$c/>1.0. Especially, experimental result shows that the eddy viscosity coefficient ratio at S$\sub$c/=3 changes from 10 to -10 although shear rate preserved constant. It is therefore suggested that proper curvature modifications, particularly the strain history effect, must be introduced into current eddy viscosity models for their application to turbulent flows subjected to curvature straining field for a non-negligible period of time.

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

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

• 대한기계학회논문집B
• /
• 제27권7호
• /
• pp.853-858
• /
• 2003

The suitable estimation of the filter width in the dynamic eddy viscosity model were investigated in high Reynolds number channel flow. In this study, the improvement on matters by optimizing the test filter shape was attempted through the numerical experiment. The way that select optimum test filter width is recommended. Some test filters, one is based on a discrete representation of the top-hat filter and another are based on a high-order filtering operation, are evaluated in simulations of the turbulent channel flow at Reynolds number 1020, based on friction velocity and channel half width. It appears that the estimation of test filter width practically can decrease the dissipative nature of dynamic eddy viscosity model with explicit test filter. It shows that the value of the filter width ratio used in the dynamic procedure must match the properties of the test filter actually used in the calculation.

• 대한기계학회논문집B
• /
• 제33권12호
• /
• pp.1014-1021
• /
• 2009

A new eddy viscosity equation was formulated from assumption of turbulence length scale equation and specific dissipation ratio equation. Then, a set of turbulence model equations for the turbulent kinetic energy ${\kappa}$, the viscosity ${\nu}_t$, and the intermittency factor ${\gamma}$ is proposed by considering the entrainment effect. Closure coefficients are determined by experimental data and resorting to numerical optimization. Present model has been applied to compute four representative cases of free shear flows and successfully compared with experimental data. In particular, the spreading rate, the centreline mean velocity and the profiles of intermittency are calculated with improved accuracy. Also, the proposed ${\nu}_t-{\kappa}-{\gamma}$ model was applied to channel flow by considering the wall effect and the results show good agreements with the Direct Numerical Simulation data.

• 대한기계학회논문집B
• /
• 제25권3호
• /
• pp.305-314
• /
• 2001

A modified 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 universal 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 mean 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.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 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).