• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.8
• /
• pp.1088-1094
• /
• 2002

A channel flow with a high Reynolds number but coarse grids is numerically studied to investigate the prediction possibility of its turbulence which is three-dimensional and time-dependent. In the present paper, a Reynolds-Averaged Navier-Stokes (RANS) model, a Large Eddy Simulation (LES) and a Navier-Stokes equation with no model are tested with a new approach of hybrid RANS/LES, which reduces to RANS model in the boundary layers and at separation, and to Smagorinsky-like LES downstream of separation, and then compared with each other. It is found that the simulations of hybrid RANS/LES method sustain turbulence like those of LES and with no model, and the results are stable and fairly accurate. This indicates strongly that gradual improvements could lead to a simple, stable, and accurate approach to predict turbulence phenomena of wall-bounded flow.

• Journal of Korea Water Resources Association
• /
• v.52 no.10
• /
• pp.697-706
• /
• 2019

A multiphase flow modeling approach equipped with a hybrid turbulence modeling method is applied to compute the gravity currents in a rectangular channel. The present multiphase solver considers the dense fluid, the less-dense ambient fluid and the air above free surface as three phases with separate flow equations for each phase. The turbulent effect is simulated by the IDDES (improved delayed detach eddy simulation), a hybrid RANS/LES, approach which resolves the turbulent flow away from the wall in the LES mode and models the near wall flow in RANS mode on moderately fine computational meshes. The numerical results show that the present model can successfully reproduce the gravity currents in terms of the propagation speed of the current heads and the emergence of large-scale Kelvin-Helmholtz type interfacial billows and their three dimensional break down into smaller turbulent structures, even on the relatively coarse mesh for wall-modeled RANS computation with low-Reynolds number turbulence model. The present solutions reveal that the modeling approach can capture the large-scale three dimensional behaviors of gravity current head accompanied by the lobe-and-cleft instability at affordable computational resources, which is comparable to the LES results obtained on much fine meshes. It demonstrates that the multiphase modeling method using the hybrid turbulence model can be a promising engineering solver for predicting the physical behaviors of gravity currents in natural environmental configurations.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2019.05a
• /
• pp.105-105
• /
• 2019

보나 여수로와 같은 수공구조물의 주변에서 발생하는 흐름 거동은 구조물 모서리에서 발생하는 흐름분리(flow separation)와 이에 따른 전단층(shear layer)과 재순환(recirculation) 흐름 영역의 발달 그리고 분리된 흐름의 재부착(reattachment)이 특징이다. 공학적으로 난류의 해석에 있어서 이러한 흐름 거동들을 정확하게 예측하는 것은 수공구조물 설계에 있어서 중요하다. 이 연구에서는 흐름 분리와 재순환 영역의 발달 그리고 흐름 재부착을 포함하는 후방계단(backward-facing step) 흐름을 155,000의 레이놀즈수 조건에서 하이브리드 RANS/LES 모델을 적용하여 해석결과를 평가한다. 하이브리드 모델로는 벽에 인접한 격자의 해상도에 상대적으로 민감하지 않은 SST(shear-stress transport) 난류 모델을 이용하는 DES(detached-eddy simulation) 기법을 적용하였다. 계단 높이가 h인 계산영역은 흐름방향 길이가 34h, 높이는 계단 상류와 하류에서 각각 1h와 2h 그리고 폭은 $2{\pi}$이다. 계단은 상류단으로부터 10h 하류부 지점에 위치한다. 경계조건으로 상부와 하부 벽면에 대해서는 비활조건을 적용한다. 상류부 수로에서 완전 발달한 흐름을 재현하기 위해서 유입경계조건은 유입부 하류 $2{\pi}h$ 지점에서 계산된 유속과 난류량을 매핑(mapping)기법을 이용하여 반복적으로 적용한다. 총 3.1백만개와 7.3백만개의 셀로 계산영역을 구현한 두 개의 계산격자 그리고 약 3.1백만개의 셀을 이용했지만 벽면 근처에서의 격자 구성을 다른 방식으로 설정한 두 가지 격자를 이용하여 격자 해상도가 DES 수치해석 결과에 미치는 영향을 분석하였다. 수치해석결과는 본 연구에서 상류단 조건으로 적용한 매핑기법이 대상 수로에서 완전 발달한 흐름을 잘 재현함을 보여주며, 합리적인 DES 해석 결과를 얻기 위해서는 벽에 수직한 방향으로 적절한 격자의 해상도와 분포가 필요함을 보여준다.

• Journal of the Society of Naval Architects of Korea
• /
• v.59 no.3
• /
• pp.141-148
• /
• 2022

A hybrid turbulence model has developed by combining a sub-grid scale model using dynamic k equation in LES with k-𝜔 SST model of RANS equation. To ascertain potential applicability of the hybrid turbulence model, fully developed turbulent channel flows at Re_𝜏=180 have been simulated of which computational domain has a top wall with coarse cells and a bottom wall with fine cells. The streamwise mean velocity and turbulent intensity profiles showed a good agreement with DNS data when using the hybrid model rather than using a single model in k-𝜔 SST or dynamic k equation models. Computational simulations of turbulent flows around KVLCC2 with a pre-swirl duct have been mainly performed using the hybrid turbulence model. Compared to the results obtained from RANS simulation with k-𝜔 SST model as well as LES with dynamic k equation SGS model, turbulent wakes of the duct in the present simulation using the hybrid turbulence model were very similar to that of LES. Also, the resistances acting on hull, rudder and duct in hybrid turbulence model were similar to those in RANS simulation whereas the viscous forces acting on the hull in LES had a significant error due to coarse cells inappropriate to the sub-grid scale model.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.244-244
• /
• 2021

난류흐름 거동은 지형이나 수공구조물과 같은 고체 경계면의 변화에 민감하게 반응하며 특징 또한 다양하다. 보나 여수로 등과 같은 단차 구조물을 통과하는 흐름은 구조물의 모서리 같은 흐름 경계면이 급변하는 지점에서는 흐름분리(flow separation)가 발생하는 것이 특징이다. 이러한 흐름분리로 인해 전단층이 발생하며 흐름 재순환(recirculation)이 구조물 하류부에 형성된다. 이 연구에서는 낙차공 형식의 단차 구조물 하류부에서의 흐름 거동을 이해하기 위해 CFD모델링을 통하여 계산된 3차원 유동장을 분석한다. 난류 모의는 하이브리드 LES(large-eddy simulation)/RANS 계산 기법인 IDDES(improved delayed detached-eddy simulation)기법을 적용한다. IDDES의 기본 모형으로는 k-ω SST모형과 Spalart-Allmaras모형을 이용하여 두 모형의 성능을 평가한다. 자유수면의 변동은 VoF(volume of fluid)기법을 이용하여 계산하며, 각 지배방정식은 최소의 수치분산을 유지하면서 수치해의 안정성을 확보할 수 있는 2차 정확도의 유한체적법을 이용하여 이산화하였다. 수치해석 결과는 레이놀즈수 23,400과 후르드수 0.22의 조건에서 기존에 계측된 자료와 비교하여 수치모형의 정확도를 평가하고 하상 단차 하류부에서의 흐름 거동 특성을 분석한다. 계산 결과는 공학적으로 널리 사용되는 RANS 수치모의에서 볼 수 없는 전단층과 난류구조의 동적 거동 특성과 이에 따른 레이놀즈 응력분포의 특성을 설명해준다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.37 no.3
• /
• pp.575-583
• /
• 2017

Various types of flow patterns around the stepped weir and spillway, such as the skimming flow over such structures and the wave-type flow with a standing undular hydraulic jump and roller downstream of the structures, are developed in open channels. Unsteady three-dimensional numerical simulations are carried out using a hybrid RANS-LES turbulence modeling approach and the volume of fluid method for resolving free surface fluctuations to represent the turbulent flow including the skimming flow and wave-type flow over a stepped weir installed in a rectangular channel. The comparison of numerical results with an existing experimental measurement reveals that the present numerical simulations reasonably well reproduce the turbulent flow passing the stepped weir, in terms of time-averaged velocity profiles at selected locations downstream of the weir, flow topology characterized by the wave-type and skimming flows, the maximum height and length of the standing wave and the length of reattachment of recirculating zone. The numerical result further elucidates the distinct flow behaviors of the wave-type and skimming flow by presenting instantaneous intense variations of free surface and velocity vectors, the distributions of Reynolds shear stress and turbulent kinetic energy and three-dimensional complex features of coherent structures and total pressure distribution.

• Journal of Korea Water Resources Association
• /
• v.52 no.10
• /
• pp.743-752
• /
• 2019

The flow in the meandering channel is characterized by the spiral motion of secondary currents that typically cause the erosion along the outer bank. Hydraulic structures, such as spur dike and groyne, are commonly installed on the channel bottom near the outer bank to mitigate the strength of secondary currents. This study is to investigate the effects of submerged vanes installed in a $90^{\circ}$ meandering channel on the development of secondary currents through three-dimensional numerical modeling using the hybrid RANS/LES method for turbulence and the volume of fluid method, based on OpenFOAM open source toolbox, for capturing the free surface at the Froude number of 0.43. We employ the second-order-accurate finite volume methods in the space and time for the numerical modeling and compare numerical results with experimental measurements for evaluating the numerical predictions. Numerical results show that the present simulations well reproduce the experimental measurements, in terms of the time-averaged streamwise velocity and secondary velocity vector fields in the bend with submerged vanes. The computed flow fields reveal that the streamwise velocity near the bed along the outer bank at the end section of bend dramatically decrease by one third of mean velocity after the installation of vanes, which support that submerged vanes mitigate the strength of primary secondary flow and are helpful for the channel stability along the outer bank. The flow between the top of vanes and the free surface accelerates and the maximum velocity of free surface flow near the flow impingement along the outer bank increases about 20% due to the installation of submerged vanes. Numerical solutions show the formations of the horseshoe vortices at the front of vanes and the lee wakes behind the vanes, which are responsible for strong local scour around vanes. Additional study on the shapes and arrangement of vanes is required for mitigate the local scour.