• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2003년도 추계 학술대회논문집
• /
• pp.158-163
• /
• 2003

This paper assesses the two-equation turbulence models available in a commercial code, FLUENT, for heat transfer in a turbulent heated pipe flow. In case of flow under $Re_D=10,000$, Standard $\kappa-\epsilon$ and Realizable $\kappa-\epsilon$ models overpredict the Nusselt number about $20\%$ compared with the experimental correlation, and RNG $\kappa-\epsilon$ model overpredicts about $30\%$ when the two-layer zonal method is employed. When wall function method is adopted, all $\kappa-\epsilon$ models show better predictions. Standard $\kappa-\omega$ and SST $\kappa-\omega$ models have the dependency on the first grid point ($0.3). As Reynolds number becomes high, the predictions of all $\kappa-\epsilon$ and $\kappa-\omega$ models are in a good agreement with the experimental correlation.

• 한국대기환경학회지
• /
• 제21권5호
• /
• pp.537-545
• /
• 2005

A new dispersion model for dense gas is constructed in the Lagrangian framework. Prediction of concentration by the proposed model is compared with measure data obtained in the experiment conducted in Thorney Island in 1984. Two major effects of dense gas dispersion, gravity slumping and stratification effect, are successfully incorporated into LDM (Lagrangian dense gas model). Entrainment effect is naturally modelled by introducing stochastic dispersion model with the effect of turbulence suppression by stratification. Not only various releasing conditions but also complex terrain can be extended to, although proposed model is appropriate for flat terrain.

• 한국환경과학회지
• /
• 제23권7호
• /
• pp.1223-1232
• /
• 2014

Density currents have been easily observed in environmental flows, for instance turbidity currents and pollutant plumes in the oceans and rivers. In this study, we explored the propagation dynamics of density currents using the FLOW-3D computational fluid dynamics code. The renormalization group (RNG) $k-{\varepsilon}$ scheme, a turbulence numerical technique, is employed in a Reynold-averaged Navier-Stokes framework (RANS). The numerical simulations focused on two different types of intrusive density flows: (1) propagating into a two-layer ambient fluid; (2) propagating into a linearly stratified fluid. In the study of intrusive density flows into a two-layer ambient fluid, intrusive speeds were compared with laboratory experiments and analytical solutions. The numerical model shows good quantitative agreement for predicting propagation speed of the density currents. We also numerically reproduced the effect of the ratio of current depth to the overall depth of fluid. The numerical model provided excellent agreement with the analytical values. It was also clearly demonstrated that RNG $k-{\varepsilon}$ scheme within RANS framework is able to accurately simulate the dynamics of density currents. Simulations intruding into a continuously stratified fluid with the various buoyancy frequencies are carried out. These simulations demonstrate that three different propagation patterns can be developed according to the value of $h_n/H$ : (1) underflows developed with $h_n/H=0$ ; (2) overflows developed when $h_n/H=1$ ; (3) intrusive interflow occurred with the condition of 0 < $h_n/H$ < 1.

• 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.

• 한국해안·해양공학회논문집
• /
• 제32권1호
• /
• pp.55-68
• /
• 2020

너울이 우세한 온화한 해양환경에서 출현하는 beach cusp에서의 경계층 streaming 수리특성을 살펴보기 위해 edge waves의 천수 과정을 수치 모의하였다. Beach cusp을 유지하는 것으로 알려진 synchronous edge waves는 같은 주기와 파고를 지니는 두 개의 Cnoidal wave가 전면해역에서 비스듬히 조우 되도록 조파하여 재현하였다. Beach cusp의 진폭 A_B과 파장 L_B은 맹방 해변에서 수행된 관측결과를 토대로 각각 1.25 m, 18 m로 선정하였다. 모의결과 천수 각 단계에서 예외 없이 경계층 streaming을 관측할 수 있었으며 최대 경계층 streaming은 사주 정점에서 발생하였다. 주기가 가장 짧은 RUN 1의 경우 그 세기는 약 0.32 m/s 내외에 분포하며 이러한 수치는 free stream 유속 u_∞ 진폭의 두 배에 달하는 것으로 wave Reynolds 응력에 기반한 Longuet-Higgins(1957)의 해석 해와는 상당한 차이를 보였다. 수치 모의과정에서 온화한 해양환경에서 해빈이 복원되는 과정을 특정할 수 있었으며 이 과정을 정리하면 다음과 같다: 너울로 구성된 파랑 무리에서 성분 파랑 간의 공진성 상호작용으로 생성된 외 중력파가 쇄파선 인근에 도달하는 경우 중력으로 인한 가속이 더해진 Phase II 파랑 궤도 운동으로 수면 가까이 상승한 많은 모래가 쇄파 시 발생하는 파 마루로부터 시작된 up-rush에 의해 전 빈 정점 가까이 이동하며 이 과정에서 발생하는 침투로 인해 퇴적되는 것으로 모의 되었다.

• Wind and Structures
• /
• 제34권5호
• /
• pp.421-435
• /
• 2022

To investigate the non-Gaussian properties of fluctuating wind pressures and the error margin of extreme wind loads on a long-span curved roof with matching and mismatching ratios of turbulence integral scales to depth (L_u^x/D), a series of synchronized pressure tests on the rigid model of the complex curved roof were conducted. The regions of Gaussian distribution and non-Gaussian distribution were identified by two criteria, which were based on the cumulative probabilities of higher-order statistical moments (skewness and kurtosis coefficients, S_k and K_u) and spatial correlation of fluctuating wind pressures, respectively. Then the characteristics of fluctuating wind-loads in the non-Gaussian region were analyzed in detail in order to understand the effects of turbulence integral-scale. Results showed that the fluctuating pressures with obvious negative-skewness appear in the area near the leading edge, which is categorized as the non-Gaussian region by both two identification criteria. Comparing with those in the wind field with matching L_u^x/D, the range of non-Gaussian region almost unchanged with a smaller L_u^x/D, while the non-Gaussian features become more evident, leading to higher values of S_k, Ku and peak factor. On contrary, the values of fluctuating pressures become lower in the wind field with a smaller L_u^x/D, eventually resulting in underestimation of extreme wind loads. Hence, the matching relationship of turbulence integral scale to depth should be carefully considered as estimating the extreme wind loads of long-span roof by wind tunnel tests.

• 대한기계학회논문집
• /
• 제9권4호
• /
• pp.403-413
• /
• 1985

본 논문에서는 수직벽 하류에 형성되는 박이 전단층의 발전과 재부착 그리고 재발전 경계층에 대해 평균 속도, 벽면의 압력 분포, 난류 강도, 레이놀즈 전단 압력 및 아직 수직벽에 대해서는 보고된바 없는 난류 떨림 속도 성분들의 3승곱 통계치를 측정하여 난류 구조의 변화를 분석하고 이를 수치적 계산 모델개발의 자 료로 제공하고자 함이 이 연구의 목적이다.

• 한국전산유체공학회지
• /
• 제10권4호통권31호
• /
• pp.51-58
• /
• 2005

The unsteady supersonic flow over two- and three-dimensional cavities has been analyzed by the integration of unsteady Reynolds-Averaged Navier-Stokes(RANS) with the k-$\omega$ turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in the cavity. An explicit 4th order Runge-Kutta scheme and an upwind TVD scheme based on the flux vector split with the van Leer limiters are used for time and space discritizations, respectively. The cavity has a L/D ratio of 3 for two-dimensional case, and same L/D and W/D ratio of I for three-dimensional case. The Mach and Reynolds numbers are 1.5 and 450000 respectively. In the three-dimensional flow, the field is observed to oscillate in the 'shear layer mode' with a feedback mechanism that follows Rossiter's formula. In the two-dimensional simulation, the self-sustained oscillating flow has more violent fluctuation inside the cavity. The primary fluctuating frequencies of two- and three- dimensional flow agree very well with the 2nd mode of Rossiter's frequency. In the three-dimensional flow, the 1st mode of frequency could be seen.

• 설비공학논문집
• /
• 제17권9호
• /
• pp.789-798
• /
• 2005

This paper is a numerical study concerning how the interactions between a pair of the vortices effect flow field and heat transfer. The flow field (common flow down) behind a vortex generator is modeled by the information that is available from studies on a half-delta winglet. Also, the energy equation and the Reynolds-averaged Wavier-Stokes equation for three-dimensional turbulent flows, together with a two-layer turbulence model to resolve the near-wall flow, are solved by the method of AF-ADI. The present results predict that the boundary layer is thinned in the regions where the secondary flow is directed toward the wall and thickened where it Is directed away from the wall. Although some discrepancies are observed near the center of the vortex core, the overall performance of the computational model is found to be satisfactory.

• 한국전산유체공학회지
• /
• 제4권2호
• /
• pp.17-30
• /
• 1999

The unsteady transitional boundary layer due to rotor-stator interaction was studied at two operation points, the design and one off design points. The off design point leads to lower blade loading and lower Reynolds number. A Navier-Stokes code developed in the previous study was parallelized to expedite computations. A low Reynolds number turbulence model was used to close the momentum equations. All computations show good agreement with experimental data. The wake induced transitional strip on the suction side of the stator is clearly captured at design point operation. There is no noticeable change in shape and phase angle of the wake induced strip even in the laminar sublayer. The wake induced transitional strip at off design point shows more complex structure. The wake induced transitional strip is observed only in the turbulent layer, and becomes obscure in the laminar sublayer and buffer layer. This behavior is probably consequent upon that the transition is governed by both wake induced strip and natural transition mechanism by Tollmien-Schlichting wave.