• 대한기계학회논문집B
• /
• 제27권12호
• /
• pp.1681-1690
• /
• 2003

The effect of roughness is a change in the velocity and turbulence distributions near the surface. Turbulence models with surface roughness effect are applied to the fully developed flow in a two-dimensional, rough wall channel. Modified wall function model, low-Reynolds number k-$\varepsilon$ model, and k-$\omega$ model are selected for comparison. In order to make a fair comparison, the calculation results are compared with the experimental data. The modified wall function model and the low-Reynolds number k-$\varepsilon$ model require further refinement, while the k-$\omega$ model of Wilcox performs remarkably well over a wide range of roughness values.

• Wind and Structures
• /
• 제17권4호
• /
• pp.361-377
• /
• 2013

Wall functions have been widely used in computational fluid dynamics (CFD) simulations and can save significant computational costs compared to other near-wall flow treatment strategies. However, most of the existing wall functions were based on the asymptotic characteristics of near-wall flow quantities, which are inapplicable in complex and non-equilibrium flows. A modified wall function is thus derived in this study based on flow over a plate at zero-pressure gradient, instead of on the basis of asymptotic formulations. Turbulent kinetic energy generation ($G_P$), dissipation rate (${\varepsilon}$) and shear stress (${\tau}_{\omega}$) are composed together as the near-wall expressions. Performances of the modified wall function combined with the nonlinear realizable k-${\varepsilon}$ turbulence model are investigated in homogeneous equilibrium atmosphere boundary layer (ABL) and flow around a 6 m cube. The computational results and associated comparisons to available full-scale measurements show a clear improvement over the standard wall function, especially in reproducing the boundary layer flow. It is demonstrated through the two case studies that the modified wall function is indeed adaptive and can yield accurate prediction results, in spite of its simplicity.

• 한국전산유체공학회지
• /
• 제21권4호
• /
• pp.78-84
• /
• 2016

In this paper, turbulence models for considering roughness in the open source code(OpenFOAM) was investigated. Wall function in the RANS(Reynolds-averaged Navier - Stokes) turbulence model was modified considering roughness on the flat plate by using roughness function. Correlation between the first layer height in the CFD model and roughness height of the plate was observed, and the most proper roughness function, and the first layer height from the plate wall in the CFD analysis was suggested in this paper.

• Journal of Advanced Marine Engineering and Technology
• /
• 제23권5호
• /
• pp.662-670
• /
• 1999

A numerical prediction was performed to clarify the air motion in the cylinder of an axisymmet-ric four-stroke reciprocating engine at its intake and compression stage. A scheme of finite volume method is used for the calculation. Modified $k-{\varepsilon}$ turbulence model is adopted and wall function is applied to the grids near the wall. The predicted mean velocity and rms velocity profiles showed a reasonable agreement with an available experimental data at its intake and compression stage. The predicted in-cylinder flow fields show that a strong turbulent twin vortex structure is pro-duced during induction but it commences to decay rapidly around inlet valve closure. The mean velocity continues to fall to a low level during compression but the turbulence intensity attains an approximate constant level.

• 대한기계학회논문집B
• /
• 제20권1호
• /
• pp.304-320
• /
• 1996

A multiple production .epsilon. equation model was developed in the low Reynolds number $\kappa$-$\varepsilon$ model with the aids of DNS data. We derived the model theoretically and avoided the use of empirical correlations as much as possible in order for the model to have generality in the prediction of complex turbulent flow. Unavoidable model constants were, however, optimized with the aids of DNS data. All the production and dissipation models in the $\varepsilon$ equation were modified with damping functions to satisfy the wall limiting behavior. A new $f_{\mu}$ function, turbulent diffusion and pressure diffusion model for the k and .epsilon. equations were also proposed to satisfy the wall limiting behavior. By, computational investigation on the plane channel flows, we found that the multiple production model for .epsilon. equation could improve the near wall turbulence behavior compared with the standard production model without the complicated empirical modification. Satisfication of the wall limiting conditions for each turbulence model term was found to be most important for the accurate prediction of near wall turbulence behaviors.

• 한국대기환경학회지
• /
• 제13권1호
• /
• pp.65-77
• /
• 1997

A theoretical and numerical investigation on the boundary-layer flow over a two- or three-dimensional hill is presented. The numerical model is based on the finite volume method with boundary-fitted coordinates. The k-$\varepsilon$ turbulence model with modified wall function and the low-Reynolds-number model are employed. The hypothesis of Reynolds number independency for the atmospheric boundary-layer flow over aerodynamically rough terrain is confirmed by the numerical simulation. Comparisons of the mean velocity profiles and surface pressure distributions between the numerical predictions and the wind-tunnel experiments on the flow over a hill show good agreement. The linear theory provides generally good prediction of speed-up characteristics for the gentle-sloped hills. The flow separation occurs in the hill slope of 0.5 and the measured reattachment points are compared with the numerical prediction. It is found that the k- $\varepsilon$ turbulence model is reasonably accurate in predicting the attached flow, while the low- Reynolds-number model is more suitable to simulate the separated flows.ows.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 춘계학술대회논문집A
• /
• pp.900-907
• /
• 2000

The three-dimensional turbulent flow in curved pipes susceptible to flow-accelerated corrosion has been analyzed numerically to predict the pressure and shear stress distributions on the inner surface of the pipes. The analysis employs the body-fitted non-orthogonal curvilinear coordinate system and a standard $ {\kappa}-{\varepsilon}$ turbulence model with wall function method. The finite volume method is used to discretize the governing equations. The convection term is approximated by a high-resolution and bounded discretization scheme. The cell-centered, non-staggered grid arrangement is adopted and the resulting checkerboard pressure oscillation is prevented by the application of a modified version of momentum interpolation scheme. The SIMPLE algorithm is employed for the pressure and velocity coupling. The numerical calculations have been performed for two curved pipes with different bend angles and curvature radii, and discussions have been made on the distributions of the primary and secondary flow velocities, pressure and shear stress on the inner surface of the pipe to examine applicability of the present analysis method. As the result it is seen that the method is effective to predict the susceptible systems or their local areas where the fluid velocity or local turbulence is so high that the structural integrity can be threatened by wall thinning degradation due to flow-accelerated corrosion.

• 대한조선학회논문집
• /
• 제42권6호
• /
• pp.593-600
• /
• 2005

The finite volume based multi-block RANS code, WAVIS developed at KRISO, is used to simulate the turbulent flow field around the KRISO container ship (KCS) and the modified KRISO tanker (KVLCC2M). The realizable k-$\varepsilon$ turbulence model with a wall function is employed for the turbulence closure. The free surface flow with and without propeller is mainly investigated for the KCS and the double model flow is concerned for the KVLCC2M which is obliquely towed in still water. The computed results are compared with the experimental data provided by CFD Tokyo Workshop 2005 in terms of wave profiles, hull surface pressure and wake distribution with and without propeller for the HCS and wake distribution and hydrodynamic forces and moments with various drift angles for the KVLCC2M.

• 대한조선학회논문집
• /
• 제31권4호
• /
• pp.58-65
• /
• 1994

계산격자크기의 변화에 따른 해의 특성을 이해하는 것은 중요하며, 특히 벽함수를 사용하지 않는 난류모형을 사용할 때 더욱 그러하다. 본 논문에서는 배주위의 3차원 난류유동장에 대한 수치해의 격자의존성에 대한 수치계산적인 조사를 수행한다. 본 연구에서는 수정된 sub-grid-scale 난류모형과 함께 유한체적법을 사용하며 복잡한 배의 기하학적 형상에 적합한 비직교의 곡선좌표계를 수치적으로 만들어 사용한다. 그리고 수학선형인 Wigley 선형과 Series 60($C_B=0.8$) 선형에 대하여 수치계산적인 연구를 수행하고, 수치해의 격자의존성을 보이기 위하여 여러가지 격자크기에 대한 계산결과들을 서로 비교하였으며 실험결과와도 비교해 보인다.

• 설비공학논문집
• /
• 제3권1호
• /
• pp.11-25
• /
• 1991

The steady, incompressible developing 3-dimensional turblent flow in a square sectioned curved duct has been investigated by using partially-parabolic equation and Finite Analytic Method. The calculation of turbulent flow field is performed using 2-equation K-$\epsilon$ turbulence model, modified wall function, simpler algorithm and numerically generated body fitted coordinates. Iso-mean velocity contours at the various sections are compared with the existing experimental data and elliptic solutions by other authors. In the region of $0^{\circ}<{\theta}<71^{\circ}$, present results agree with the experimental data much better than the elliptic solution for the similar number of grid points. Furthermore, for the same tolerance, the present solution converges four times faster than the elliptic solution.